Metody matematicheskoj fiziki v obrabotke signalov i izobrazhenij
Методы математической физики в обработке сигналов и изображений
Textbook, 2009
70 Pages
Dr. Eugene Postnikov (Author)
I. FLEXPDE: СРЕДА ОПИСАНИЯ И РЕШЕНИЯ ЗАДАЧ МАТЕМАТИЧЕСКОЙ ФИЗИКИ
Интерфейс системы
Структура описания задачи
Задание области решения на плоскости
Задание граничных условий
Графический вывод решения
Контроль вида графического вывода
II. ОБРАБОТКА ИЗОБРАЖЕНИЙ И СИГНАЛОВ
1. Изображение как поле яркости
2. Простая диффузия: сглаживание деталей и подавление шума
3. Метод разделения переменных, собственные числа и собственные функции
4. Градиент, поток и плотность потока его источников
5. Выделение контуров и анизотропная диффузия
6. Кратномасштабные методы
7. Глобальный и локальный спектральный анализ
ВВЕДЕНИЕ
Одной из важных сторон работы с сигналами, в частности, с изображениями, является предварительная обработка, имеющая своей целью улучшение их качества за счет подавления шума, подчеркивания контуров и границ или же, наоборот, сглаживания контрастных деталей. Помимо улучшения визуального восприятия, подобные процедуры важны также для задач сжатия исходных данных путем выделения характерных особенностей образа.
Разработка соответствующей математической основы существенным образом базируется на теории дифференциальных уравнений математической физики. При таком подходе яркость изображения (или же интенсивность его цветовых составляющих) рассматривается как скалярная функция координат экрана, подчиняющаяся тем же закономерностям, что и, например, поля температуры или концентрации диффундирующих частиц. Основанный на такой математической аналогии с уравнением простой диффузии кратномасштабный анализ изображений был обоснован в конце 1970-х – начале 1980-х годов Марром, Хильдретом, Кёндеринком, Виткином (Marr, Hildreth, Koenderink, Witkin). Следующим заметным шагом стала концепция анизотропной диффузии, предложенная в 1990 году Перона и Маликом (Perona, Malik) [1] и позволяющая эффективно учесть при обработке локальные особенности изображения. Эта работа положила начало целому направлению в математической теории обработки изображений – применению нелинейных дифференциальных уравнений в частных производных и разработке прикладных средств на их основе [2] .
Еще одно важное направление, возникшее в середине 1980-х годов и активно развивающееся в настоящее время – вейвет-анализ, имеющее важные приложения как для локального анализа сигналов и изображений, как и для их обработки и сжатия. Функции-вейвлеты представляют весьма широкое семейство, различные представители которого адаптированы к решению разнообразных задач. Их представители из действительнозначного гауссова семейства позволили Малла (Mallat) развить кратномасштабные методы Марра–Хильдрета–Кёндеринка–Виткина, а использование комплексного вейвлета Морле стало мощным обобщением спектрального анализа сигналов.
Основная цель данного учебного пособия – дать обзор основных методов обработки сигналов и изображений, базирующихся на решении дифференциальных уравнений в частных производных. Изложение материала рассчитано на то, что читатель владеет основами математического анализа (дифференцирование, интегрирование, основы векторной алгебры), все же вопросы, относящиеся к методам математической физики, вводятся и поясняются по ходу изложения в объеме, необходимом для понимания сущности изучаемых процессов. Естественно, это вызывает определенное упрощение математической стороны проблемы, основное внимание уделяется пояснению на примерах, а не строгим доказательствам. Читатели, заинтересованные в более глубоком знании чисто-математической стороны вопроса могут обратиться к учебникам математической физики [3] . Метод вычисления комплексного интегрального вейвлет-преобразования с вейвлетом Морле описан на основе статьи автора [4] .
Так как пособие ориентировано, прежде всего, на пояснение сущности методов, основанных на решении дифференциальных уравнений в частных производных, то перед читателем прежде всего ставится задача правильно формулировать соответствующую математическую проблему. Нахождение ее же решение предполагается численно, с использованием система решения задач математической физики FlexPDE, разработанной фирмой PDE solutions [5] . Она представляет собой совокупность среды для описания задач математической физики, решателя, базирующегося на методе конечных элементов (Галеркина), и средств графического вывода решения.
Применяемый во FlexPDE язык описания задачи является полностью декларативным и весьма близок к естественному языку: пользователь должен только описать конфигурацию области решения и постановку граничной задачи, не производя программирования самого процесса решения. Написание операторов, используемых в дифференциальных уравнениях, практически совпадает со стандартным математическим. FlеxPDE может решать системы дифференциальных уравнений первого и второго порядка на прямой, плоскости (в декартовых и полярных координатах) и в пространстве (в декартовых координатах). Уравнения могут быть как линейными, так и нелинейными, численное решение осуществляется модифицированным методом Ньютона-Рафсона.
Студенческая версия FlexFDE (с ограничением узлов сетки до 100 узлов в одномерных, 800 – в двумерных и 1600 – в трехмерных задачах) может быть бесплатно скачана с сайта фирмы-разработчика. Там же имеется полный комплект документации и электронные варианты книг, посвященных применению системы для решения задач математической физики.
Первую часть данного пособия составляет обзор среды описания задач, а также принципов программирования и основных функций FlexPDE . Каждый параграф второй части состоит из изложения теоретического материала, примера, в котором приведен подробный разбор написания программы, реализующий описанный алгоритм, упражнения с графическими ответами и вариантов контрольных заданий, которые могут использоваться в учебных группах.
Курс построен на основе лекций и практических занятий по предмету «Методы математической физики» для специальности «Информатика» в Курском государственном университете.
I. FLEXPDE: СРЕДА ОПИСАНИЯ И РЕШЕНИЯ ЗАДАЧ МАТЕМАТИЧЕСКОЙ ФИЗИКИ
ИНТЕРФЕЙС СИСТЕМЫ
illustration not visible in this excerpt
При запуске FlexPDE открывается редактор среды:
illustration not visible in this excerpt
Области редактора
illustration not visible in this excerpt
illustration not visible in this excerpt
Интерактивное управление видом графиков решения
Двойной клик левой клавишей мыши по любому из окон графиков приводит к его увеличению до размеров полного экрана. Повторный двойной клик возвращает многооконный режим.
Однократный клик правой клавишей мыши вызывает появление меню, в котором управлению видом графика отвечают пункты: Maximize, Restore (описанное выше увеличение – уменьшение окна), Rotate * .
Выполнение команды Rotate вызывает появление куба с координатными осями и пирамидой, указывающей направление вверх (см. рис) на месте области контроля сетки разбиения ** . Смещением движков можно вращать график (направление координатных осей можно контролировать по кубу). Движок Az определяет вращение вокруг вертикальной оси (проходящей через вершину пирамиды перпендикулярно ее основанию); движок Po осуществляет поворот вокруг горизонтальной оси (перпендикулярной оси Az). Переключатели режимов отображения вращения: при включенном Live Redraw график вращается одновременно с кубом, при Static Redrow – вращается только куб, а для поворота графика требуется нажать Redrow.
Интерактивное сохранение графика
В меню, вызываемом однократным кликом правой кнопкой мыши по графику, пункт Export позволяет сохранить содержимое окна в файл в различных графических форматах и используя различную толщину линий (см. рис.)
illustration not visible in this excerpt
Пункт Print отправляет на печать изображение, отображаемое в данном графическом окне.
СТРУКТУРА ОПИСАНИЯ ЗАДАЧИ
Описание задачи для FlexPDE состоит из фиксированной последовательности разделов, снабженных заголовками. После заголовка указываются операторы, задающие соответствующую часть задачи. В том случае, когда какой-либо раздел в данной задаче не используется, соответствующий заголовок не пишется. Текст может быть снабжен комментариями в фигурных скобках (<>), с употреблением только латинских букв. Кроме того, восклицательный знак (!) обращает в комментарий стоящий после него текст строки.
TITLE ‘Имя’– имя данной задачи (латиницей), указывается в апострофах.
В этом разделе указывается значение контролирующих параметров.
В случае задачи на собственные значения тут помещается обязательный оператор: modes=#. Здесь # – рассчитываемое число мод; одновременно инициализируется системная переменная LAMBDA, которой будут присваиваться собственные числа задачи.
Выбор системы координат: декартовой прямоугольной (по умолчанию) или цилиндрической. Команды: YCYLINDER – цилиндрическая система координат, в которой ось Oz направлена вертикально, а радиальное направление – горизонтально; XCYLINDER – система координат, к которой ось Oz горизонтальна, а Or – вертикальна. В декартовых координатах командой CARTESIAN(‘name1’, ‘name2’) можно изменить имена координат с x и y на name1 и name2 (в этом случае эти имена должны использоваться и в дифференциальных операторах, например dname1() вместо dx()).
Задание имен неизвестных функций, которые ищутся в качестве решения дифференциальных уравнений. Именем переменной может быть одна или несколько латинских букв, а также комбинация букв, цифр и подчеркиваний (имя должно начинаться с буквы). Заглавные и строчные буквы не различаются.
В случае сложных расчетов следует указывать в скобках после имени переменной задается параметр threshold, который определяет уровень точности определения переменной (ее значение, меньше которого не следует проводить уточнение деталей решения).
Например, u(threshold=1e-3), v(threshold=0.01).
Инициализация пользовательских констант (присвоение им значений) и функций (задание функции ее выражением через координаты, время или другие переменные)
Задание значения константы:
Задание функции: =. Выражение связывает координаты, время, векторы, заданные ранее константы и функции при помощи арифметических, векторных, дифференциальных операторов.
Например: f=x^2, g=sin(x)+exp(y), f1=1/r, v=-k*grad(phi) – в последнем случае предварительно должны быть заданы k (например, k=5) и phi (например, phi=x^2+y^2).
При использовании функции в дальнейших выражениях, достаточно указывать ее имя, аргумент в скобках не указывается.
Запись системы дифференциальных уравнений, которым удовлетворяют неизвестные функции, указанные в разделе VARIABLES. Число уравнений должно равняться числу неизвестных, перед каждым уравнением указывается имя одной входящей в него переменной c двоеточием, например, u: dt(u)=dxx(u).
Значения неизвестных функций в начальный момент времени или функция, служащая начальным приближением, для нелинейных задач.
Описание геометрии области, в которой будет проводиться решение, а также указание граничных условий для дифференциальных уравнений.
Указание границ изменения времени для нестационарных задач:
FROM time1 TO time 2, где time1 – начальное время, time 2 – конечное время. По умолчанию, начальный временной шаг равен
(time2– time1)/104. Это значение можно изменить на другое, командой в форме: FROM time1 BY timestep TO time 2.
Команды графического вывода текущего состояния решения.
Команды вывода графиков полученного решения, а также сохранения численных результатов в файл.
Графический вывод (и сохранение в файл) значений искомой функции в заданных точках, с течением времени в нестационарных задачах.
END – команда указывает на завершение текста программы.
ЗАДАНИЕ МАТЕМАТИЧЕСКИХ ВЕЛИЧИН: ПЕРЕМЕННЫХ, ФУНКЦИЙ, ОПЕРАТОРОВ
Зарезервированные системные переменные
(не требуют декларации, не могут использоваться как функции)
illustration not visible in this excerpt
Зарезервированные системные функции
illustration not visible in this excerpt
illustration not visible in this excerpt
illustration not visible in this excerpt
illustration not visible in this excerpt
Операторы интегралов и суммы
illustration not visible in this excerpt
Неаналитические функции
illustration not visible in this excerpt
ЗАДАНИЕ ОБЛАСТИ РЕШЕНИЯ НА ПЛОСКОСТИ
Двумерные задачи
Область, к которой осуществляется решение, задается прочерчиванием ее контура и может содержать несколько подобластей. Описание контура каждой области (подобласти) должно начинаться с команды REGION. При этом области (подобласти) может быть присвоен номер: REGION #, где через # обозначено число – номер (области) подобласти, или имя: REGION ‘name’. Контур должен быть замкнутым и рисуется последовательным соединением вершин, идущим от начальной точки START (x0,y0) (между START и значением координат может быть задано имя контура в кавычках: “con”) линиями или дугами к дальнейшим точкам с координатами (x1, y1), (x2, y2) и т.д. до точки (xn, yn), которая соединяется с начальной, что определяется указанием команды CLOSE вместо координат. Имена областей (контуров) могут быть использованы в интегральных оперторах. Чтобы подобласть была вырезана из области, необходимо начать ее описание не с команды REGION, а с команды EXCLUDE. Вырезаемая область описывается после той, из которой она вырезается.
Построение отрезка
LINE TO (xk,yk) проводит отрезок от текущей точки до точки с координатами (xk,yk).
Построение дуги окружности
ARC (radius=R) TO (xk,yk) проводит дугу окружности радиусом R от текущей точки до точки с координатами (xk,yk).
ARC (center=xc,yc) TO (xk,yk) проводит дугу окружности c центром в точке (xc,yc) от текущей точки до точки с координатами (xk,yk).
ARC (center=xc,yc) angle=A проводит дугу окружности c центром в точке (xc,yc) от текущей точки на угол А, выраженный в градусах. Положительное значение А соответствует повороту против часовой стрелки.
ARC (center=xc,yc) angle=360 строит окружность с центром в точке (xc,yc).
Построение линии, заданной функцией
Контур области в виде линии, заданной функцией Abbildung in dieser Leseprobe nicht enthalten строится из линий, последовательно соединяющих точки с координатами Abbildung in dieser Leseprobe nicht enthalten:
REPEAT xi=x0 BY step TO xn
Здесь x0 – координата х начальной точки, step – шаг по координате x, xe – координата конечной точки. Функция f(xi) может быть записана непосредственно в виде выражения, содержащего арифметические действия и зарезервированные функции FlexPDE, либо может быть предварительно описана в разделе DEFINITIONS в формате f(имя аргумент)=действия над аргументом (например, f(i)=i^2). Имя аргумента не должно совпадать с зарезервированными переменными.
Построение границы по массиву координат ее вершин
Часть контура области, состоящая из линий, последовательно соединяющих точки с координатами (xc[i], yc[i]), где xc[i] и yc[i] – числа, содержащиеся в массивах xc и yc под номером i, организуется в виде цикла:
illustration not visible in this excerpt
Одномерные задачи
В разделе SELECT указывается соотвтствующая спецификация, например, coordinate=cartesian1 для обыкновенных декартовых координат.
В разделе BOUNDARIES линии задаются аналогично двумерному случаю, за исключением того, что не производится замыкания (отсутствует команда to close) и координата указывается одна. Задание отрезка [x0, x1]:
START (x0) LINE TO (x1)
ЗАДАНИЕ ГРАНИЧНЫХ УСЛОВИЙ
Граничные условия указываются для всех искомых функций задачи при вычерчивании области в разделе BOUNDARIES. Условие действует на всех следующих за ним линиях, до тех пор, пока не будет указано иное, либо не встретится оператор nobc(имя_функции), который прекращает действие граничного условия для указанной функции.
illustration not visible in this excerpt
Условия в точках и на незамкнутых линиях
Чтобы выделить точку с координатами (xp,yp), надо после описания области решения указать оператор fixed point(xp,yp). После него можно задать значение функции в этой точке point value(f)=g или значение величины плотности источника поля point load(f)=g (положительное значение g соответствует поглощению поля в точке). Эти же условия употребляются при задании граничных условий для одномерных задач (без оператора fixed point на границах интервала), point load играет для одномерных задач роль natural для двумерных.
Незамкнутые линии строятся в разделе BOUNDARIES после описания области решения. Заданию кривой предшествует заголовок FEATURE, далее – описанное выше построение линии, начинающейся с оператора START операторами LINE и ARC; построение не должно заканчиваться оператором FINISH. Граничные условия указываются так же, как и на границе области расчета. Если линии присвоено имя, то оно может быть использовано для вычисления интеграла по контуру, совпадающему с этой линией.
ГРАФИЧЕСКИЙ ВЫВОД РЕШЕНИЯ
Графики решения строятся в разделах MONITORS и PLOTS. Для стационарных задач в разделе MONITORS отображается процесс приближения к решению, в PLOTS – окончательный результат (который записывается в файл с расширением .pg4 и может быть просмотрен после окончания решения по команде верхнего меню P lots). В нестационарных задачах оба раздела работают одинаково. По умолчанию двумерные графики строятся на всей области решения. Для построения графика в выделенной подобласти с именем ‘name’ следует указать после оператора построения графика on ‘name’. Аналогично строится график типа elevation вдоль именованной линии.
Для нестационарных задач необходимо после заголовка использовать один из операторов, указывающих, когда следует перерисовывать график: for cycle=# (на каждом #-ом цикле пересчета); for time=time1,time2 (в моменты времени time1, time2 и т.д.); for time time1 by time_step to time2 (от момента времени time1 до time2 с шагом time_step).
illustration not visible in this excerpt
Каждый двумерный график строится в отдельном окне; вдоль линии можно одновременно построить несколько графиков, используя перечень аргументов: elevation (f1,f2,f3) from (x1,y1) to (x2,y2).
[1] Perona P., Malik J. Scale-Space and Edge Detection Using Anisotropic Diffusion // IEEE Transactions on Pattern Analysis and Machine Inteligence.– 1990.– V. 12.– No. 7.– P. 629–639.
[2] См. широкий обзор различных методов в Special Issue on Partial Differential Equations and Geometry-Driven Diffusion in Image Processing // IEEE Transactions on Image Processing.– 1998.– V. 7.– No. 3.
[3] Например, можно рекомендовать такие учебники, сочетающие математическую строгость с ясностью изложения, как
Арсенин В.Я. Математическая физика. Основные уравнения и специальные функции.- М.: Наука, 1966 (переиздание: Арсенин В.Я. Методы математической физики и специальные функции.- М.: Наука, 1974.); Кошляков Н.С., Глинер Э.Б., Смирнов М.М. Уравнения в частных производных математической физики.- М.: Высшая школа.- 1970; Мартинсон Л.К., Малов Ю.И. Дифференциальные уравнения математической физики.– М.: Изд-во МГТУ им. Н.Э. Баумана, 2002.
[4] Постников Е.Б. Вычисление непрерывного вейвлет–преобразования как решение задачи Коши для системы дифференциальных уравнений в частных производных // Журнал вычислительной математики и математической физики.– 2006.– Т. 46.– №1.– С. 77–82.
* Команда применима только для трехмерных графиков (поверхностей).
** Редактор написан для разрешения экрана 1024´768 точек и не масштабируется (за исключением окна графического вывода), поэтому при меньшем разрешении экрана (например, 800´600) движки не видны.
Основы работы в конечно-элементном пакете FlexPDE. Часть I
Пособие посвящено основам работы в пакете FlexPDE, который предназначен для построения численных решений дифференциальных уравнений в частных производных с помощью метода конечных элементов. Описаны основные разделы типовой программы. Для демонстрации и эффективного освоения базового функционала пакета представлены примеры скриптов и задания для самостоятельной работы. Учебное пособие предназначено для широко круга читателей, в том числе для студентов бакалавриата, обучающихся по направлениям подготовки 01.03.02 прикладная математика или 01.03.03 механика и математическое моделирование.
Публикуется в авторской редакции.
Год выпуска: 2022
Авторы: В. В. Дударев, О. Г. Пустовалова
ISBN: 978-5-9275-4121-8
Похожие записи:
- Математика и анализ данных с поддержкой MS Excel и языка R
- Российская история рекламы
- Итоговая государственная аттестация. Методические указания к выполнению выпускной квалификационной работы магистра
- Психология социальных коммуникаций и рекламы. Часть 1. Психологические основы рекламы и массовых коммуникаций
Добавить комментарий Отменить ответ
«Словарь терминов по изобразительному искусству. Живопись. Графика. Скульптура» включает более четырехсот искусствоведческих терминов и
Для Джин Хэйнс краски – это сокровища, а кисть – волшебная палочка, с помощью
«Из всех богов лишь Смерть к дарам бесчувственна…» однако некоторые люди умеют ей угодить.
Новое издание известного учебника по органической химии отражает последние достижения в теории и изучении
Основы работы в конечно-элементном пакете FlexPDE. Часть I
Также данная книга доступна ещё в библиотеке. Запишись сразу в несколько библиотек и получай книги намного быстрее.
Как читать книгу после покупки
- Чтение только в Литрес «Читай!»
Посоветуйте книгу друзьям! Друзьям – скидка 10%, вам – рубли
По вашей ссылке друзья получат скидку 10% на эту книгу, а вы будете получать 10% от стоимости их покупок на свой счет ЛитРес. Подробнее
Стоимость книги: 198 ₽
Ваш доход с одной покупки друга: 19,80 ₽
Чтобы посоветовать книгу друзьям, необходимо войти или зарегистрироваться Войти
- Объем: 133 стр.
- Жанр:п рограммы, у чебники и пособия для вузов
- Теги:б акалавриат, д ифференциальные уравнения, м атематические модели, м атематическое и компьютерное моделирование, п рикладная математика, п рикладная механикаРедактировать
Эта и ещё 2 книги за 399 ₽
По абонементу вы каждый месяц можете взять из каталога одну книгу до 700 ₽ и две книги из специальной подборки. Узнать больше
Оплачивая абонемент, я принимаю условия оплаты и её автоматического продления, указанные в оферте
Оплатить Отмена
Описание книги
Пособие посвящено основам работы в пакете FlexPDE, который предназначен для построения численных решений дифференциальных уравнений в частных производных с помощью метода конечных элементов. Описаны основные разделы типовой программы. Для демонстрации и эффективного освоения базового функционала пакета представлены примеры скриптов и задания для самостоятельной работы. Учебное пособие предназначено для широко круга читателей, в том числе для студентов бакалавриата, обучающихся по направлениям подготовки 01.03.02 прикладная математика или 01.03.03 механика и математическое моделирование. Публикуется в авторской редакции.
Подробная информация
Возрастное ограничение: 0+ Дата выхода на ЛитРес: 12 января 2023 Дата написания: 2022 Объем: 133 стр.
ISBN: 978-5-9275-4121-8 Общий размер: 4 MB Общее кол-во страниц: 133 Размер страницы: Правообладатель: Южный Федеральный Университет
«Основы работы в конечно-элементном пакете FlexPDE. Часть I» — читать онлайн бесплатно фрагмент книги. Оставляйте комментарии и отзывы, голосуйте за понравившиеся.
Моделирование электромагнитных полей c помощью программы FlexPDE
FlexPDE – программа, предназначенная для построения сценарных моделей решения дифференциальных уравнений методом конечных элементов, и специализирующаяся на моделировании векторных полей, включая электромагнитные. FlexPDE по сценарию, написанному пользователем, производит операции, необходимые для того, чтобы преобразовать описание системы дифференциальных уравнений в частных производных в модель для расчета методом конечных элементов, найти решение для этой системы и представить результаты в графической форме. Модели могут быть одномерными, двухмерными и трехмерными. Студенческая версия свободно скачивается с официального сайта программы http://www.pdesolutions.com , обладает таким же функционалом, как и промышленная версия, но имеет ограничения: пять дифференциальных уравнений и 100, 800 или 1600 узлов сетки для одномерных, двухмерных и трехмерных задач соответственно. FlexPDE не ограничивает пользователя заранее заданным списком прикладных задач или видов уравнений. Выбор вида дифференциальных уравнений в частных производных полностью зависит от пользователя. Язык сценария позволяет пользователю описывать в естественном формате математический аппарат его системы дифференциальных уравнений в частных производных и структуру области решений в целом. Например, в сценарии имеется раздел EQUATIONS (УРАВНЕНИЯ), в котором уравнение Лапласа можно представить как Div(grad(u)) = 0 Аналогично, в сценарии имеется раздел BOUNDARIES (ГРАНИЧНЫЕ УСЛОВИЯ), в котором геометрические границы двумерной области решений описываются по кривой периметра, Start(x1,y1) line to (x2,y1) to (x2,y2) to (x1,y2) to finish. Следующее описание дается для версии FlexPDE 6. Составление сценария расчета Сценарий описания задачи представляет собой текстовый файл. Содержание этого файла представляет собой ряд разделов, каждый из которых идентифицируется при помощи заголовка. Чаще всего используются следующие разделы: TITLE – заголовок программы, COORDINATES, VARIABLES, SELECT, DEFINITIONS, EQUATIONS, BOUNDARIES, PLOTS, SUMMARY, END – обозначает конец программы. В любом месте скрипта в фигурных скобках <> и после восклицательного знака ! можно оставлять комментарии. Более подробно о некоторых секциях можно прочитать ниже. С полной спецификацией используемых команд можно ознакомиться в руководстве пользователя. COORDINATES Описывает систему координат (декартову, цилиндрическую или сферическую). Двухмерная и трехмерная декартовы системы координат описываются командами: cartesian1, cartesian2, cartesian3. VARIABLES Данный раздел используется задания переменных задачи. Для моделирования электрического поля используется электрический потенциал. Для моделирования магнитного поля используется магнитный векторный потенциал. SELECT Данный раздел предназначен для задания параметров сетки, точности решения задачи и т.п. DEFINITIONS 10
Это раздел задания вспомогательных переменных и определений задачи. В роли определений выступает задание параметров материалов, констант, функций и т.д. При наборе определений следует учесть, что некоторые имена зарезервированы (например, нельзя использовать букву t, поскольку она зарезервирована как переменная времени). Все выражения, используемые в данной секции, являются вспомогательными к основным уравнениям, которые описываются в следующей секции. Примеры выражений приведены в таблицах 1-3. Таблица 1
| Определение | Пример |
| Eps_metal=1e6 | |
| Относительная диэлектрическая | Eps_diel = 88 |
| проницаемость ε r ε 0 материалов и сред | Eps_air=1 |
| Eps0 = 8.854e-12 | |
| Относительная магнитная | mu = 1 |
| проницаемость μμ 0 | |
| Параметры расчетной области, | r_provod = 0.02 |
| R = 0,10 | |
| геометрические размеры модели | |
| R0 = 10*R | |
| Данные о заданных токах, зарядах, | U_provod = 500 |
| напряжений и т.п. | I = 1 |
| S=800e-6 | |
| J = I/S |
Таблица 2. Вспомогательные формулы расчета физических величин электрического поля
| Определение | Пример | |
| Ex=-dx(U) | ||
| Напряженность электрического поля | Ey=-dy(U) | |
| Eabs=sqrt(Ex^2+Ey^2) | ||
| E = -grad(U) | ||
| Энергия электрического поля | W = integral(0.5*eps*eps0*E^2) | |
| Электрический заряд. | ||
| Следует учесть, что из теоремы | Q = eps*eps0*sintegral(normal(grad(U)), | |
| ‘ребро или область1’, ‘ребро или | ||
| Остроградского-Гаусса следует, что | ||
| область2′,…) | ||
| Электрическая емкость. | ||
| Следует учесть, что конденсатор | C = q/(U1-U0) | |
| накапливает энергию | C = 2*W/(U1-U0)^2 | |
| W = CU 2 /2= qU /2= q 2 /(2 C ) | ||
Таблица 3. Вспомогательные формулы расчета физических величин магнитного поля
| Определение | Пример |
| B = curl(Ax,Ay,0) | |
| Индукция магнитного поля | Bxx= -dz(Ay) |
| Byy= dz(Ax) | |
| Bzz= dx(Ay)-dy(Ax) | |
| Энергия магнитного поля | W = integral(0.5*mu*m0*B^2) |
| Потокосцепление, магнитный поток | F = sintegral(B, ‘имя ребра или области |
| 1′, ‘имя ребра или области 2’,…) | |
| 11 |
| Индуктивность. | |
| Следует учесть, что катушка | L = F/I |
| индуктивности накапливает энергию | L = 2*W/I^2 |
| W = LI 2 /2 |
Кроме того, в разделе определений можно вычислить физические величины по заранее известным формулам для дальнейшего сравнения результата численного решения с теорией. Например, известно, что электрическая емкость плоского конденсатора вычисляется по формуле С = ε r ε 0 S / d , где ε r – относительная диэлектрическая проницаемость, S – площадь пластин, d – расстояние между ними. После построения модели плоского конденсатора и численного расчета электростатического поля программа сравнит емкость идеального и реального конденсатора. EQUATIONS В данном разделе записываются основным уравнения в частных производных для каждой переменной. Например, для электростатического поля это уравнение Лапласа, оно будет выглядеть: div( eps * grad( v)) = 0. Для магнитного поля также можно записать уравнение Лапласа: Ax: div(grad(Ax)/mu) = 0 Ay: div(grad(Ay)/mu) = 0 Здесь запись Ax и Ay требуется для указания соответствия уравнения переменной. BOUNDARIES Раздел предназначен для задания геометрии и граничных условий. Каждый блок (область) задачи начинается с команды REGION N ‘Имя области (латиницей)’, где N – порядковый номер области. Далее регион описывается точкой, отрезками от точки до точки, дугами и другими более сложными объектами (таблица 4). Таблица 4. Геометрические объекты
| Объект | Пример |
| Точка | NODE POINT (x, y) |
| Прямоугольник, полигон | LINE TO (x1,y1) TO (x2,y2) TO (x3,y3) TO CLOSE |
| Дуга, окружность, эллипс | ARC TO (x1,y1) to (x2,y2) |
| ARC ( RADIUS = R ) to (x,y) | |
| ARC ( CENTER = x1,y1 ) to (x2,y2) | |
| ARC ( CENTER = x1,y1 ) ANGLE=angle_degrees |
Например, следующий сценарий описывает модель, представленную на рис. 1: DEFINITIONS eps = 1 < значение по умолчанию>BOUNDARIES REGION 1 ‘domain’ < расчетная (основная) область задачи >START(0,0) LINE TO (10,0) TO (10,10) TO (0,10) TO CLOSE REGION 2 eps = 2 < иное значение для двух областей ниже >START ‘box1’ (1,1) LINE TO (2,1) TO (2,2) TO (1,2) TO CLOSE START ‘box3’ (5,5) LINE TO (6,5) TO (6,6) TO (5,6) TO CLOSE REGION 3 < используется по умолчанию >START ‘box2’ (3,3) LINE TO (4,3) TO (4,4) TO (3,4) TO CLOSE 12
Рис. 1. Геометрическая модель Проверить геометрию модели можно с помощью кнопки 
Domain Review Граничные условия задаются для контуров (ребер) и точек геометрической модели. Команды для указания граничных условий представлены в таблице 5.
| Таблица 5. | Граничные условия | |
| Граничное условие | Пример | |
| Условие Дирихле (задание потенциала) | VALUE (U) = 220 | |
| Условие для потока переменной сквозь | NATURAL (U) = 0 1 | |
| границу, в частности условие Неймана | LOAD (U) = expression | |
| (задание производной n ). | ||
| Скачкообразное | изменение переменной, | CONTACT (U) = EDS |
| например задание ЭДС | ||
Граничные условия указываются перед командами LINE и ARC: BOUNDARIES REGION 1 START(0,0) VALUE(u) = 0 LINE TO (1,0) < значение на нижней границе >NATURAL(u)=0 LINE TO (1,1) < значение на правой границе >VALUE(u)=0 LINE TO (0,1) < значение на верхней границе >NATURAL(u)=0 LINE TO CLOSE < значение на левой границе >REGION 2 NODE POINT (0.5, 0.5) POINT VALUE(u)=250 С команды FEATURE (аналогично команде REGION) начинается описание набора линий и дуг, не являющихся частью модели и предназначенных для последующих расчетов. Формат команды следующий: FEATURE START(0,0) LINE TO(10,10) TIME Если переменные зависят от времени, то указывается раздел TIME, в котором указываются период времени, в течение которого проводится расчет. Например, следующим образом: TIME 0 TO 1 1 Однородное условие Неймана устанавливается по умолчанию для всех контуров (ребер). 13
PLOTS Раздел предназначен для вывода графических результатов. Возможные варианты графики представлены в таблице 6. Таблица 6. Графики
| График | Пример |
| Эквипотенциальные линии | CONTOUR(U) |
| REPORT(C) as «Емкость (Ф)» 1 | |
| CONTOUR(U) painted | |
| CONTOUR(Ax, Ay) painted | |
| CONTOUR (U) as ‘Потенциал около проводника’ | |
| ZOOM(x0, y0, Δx, Δy) | |
| Картина векторного поля | VECTOR(E) as ‘Напряженность электрического |
| поля’ | |
| Зависимость физической величины | ELEVATION(U) as ‘Потенциал от стены до |
| вдоль какого-либо пути | проводника’ from (x0, y0) to (x1, y1) |
| ELEVATION(U) on ‘box2’ | |
| Трехмерный график | SURFACE(U) |
Если переменные зависят от времени, то перед построением графиков сразу после команды PLOTS указываются моменты времени, для которых будут строиться графики в следующем формате: FOR t = начальный момент BY шаг TO конечный момент. Построение графиков зависимости физических величин от времени в конкретных точках осуществляется в дополнительном разделе HISTORIES командами history(u) at (x1,y1) (x2,y2). SUMMARY Раздел описывает представление результатов расчета в текстовом виде. Для этого используется команда REPORT. Например: SUMMARY REPORT(C) as «Емкость (Ф)» REPORT(C0) as «Емкость в теории (Ф)» REPORT(W) as «Энергия» Система обозначений Программа принимает стандартные команды расчета функций sqrt(), sin(), cos(), max(), min() и т.п. Обрабатываются команды для работы с векторами и комплексными числами. Поскольку программа работает с векторными полями, имеются команды для работы с векторами. Вектор может быть указан либо своей переменной (например, E), либо собран из составляющих Ex,Ey, может быть задан командой vector(Ex,Ey,Ez) (см. таблицу 7. Таблица 7. Работа с векторами
| Описание | Команда |
| Длина вектора | MAGNITUDE ( vector ) |
| MAGNITUDE ( argx, argy ) | |
| Векторное произведение векторов | CROSS (vector1, vector2) |
| Скалярное произведение векторов | DOT (vector1, vector2) |
| Нормальная и тангенциальная | NORMAL ( vector ) |
| составляющие вектора к поверхности | NORMAL ( argx, argy ) |
| границы (используется в разделе | TANGENTIAL(vector) |
| BOUNDARIES и тех местах, где имеется | TANGENTIAL ( argx, argy ) |
1 Следует отметить, что кириллица не поддерживается. Рекомендуется использовать латиницу. 14