- •РАЗДЕЛ I. МЕТОДЫ ОПТИМИЗАЦИИ
- •1. ПОСТАНОВКА ЗАДАЧ ОПТИМИЗАЦИИ И ИХ КЛАССИФИКАЦИЯ.
- •1.1. Постановка задачи
- •1.2. Методы нахождения минимума функции одной переменной
- •Метод деления отрезка пополам
- •Метод золотого сечения
- •Метод последовательного перебора
- •Метод квадратичной параболы
- •Метод кубической параболы
- •1.3. Рельеф функции
- •1.4. Представление функции в окрестности минимума
- •1.5. Классификация методов оптимизации
- •2.МЕТОДЫ НУЛЕВОГО ПОРЯДКА
- •2.1. Описание общего алгоритма методов покоординатного спуска.
- •2.2 Метод ГАУССА-ЗЕЙДЕЛЯ
- •2.3. Метод ПАУЭЛЛА
- •2.4. Метод ДСК
- •2.5. Метод РОЗЕНБРОКА
- •2.6. Метод покоординатного последовательного перебора
- •2.7. Метод ХУКА-ДЖИВСА
- •2.8. Метод НЕЛДЕРА-МИДА
- •3. МЕТОДЫ ПЕРВОГО ПОРЯДКА И МЕТОДЫ С ПЕРЕМЕННОЙ МЕТРИКОЙ
- •3.1. Метод спуска по градиенту
- •3.3. Метод Давидона-Флетчера-Пауэлла
- •3.4. Метод проективного градиента
- •3.5. Метод Мак-Кормика 1
- •3.6. Метод Мак-Кормика 2
- •3.7. Метод Гольдфарба
- •3.8. Метод Гринстадта
- •4. МЕТОДЫ УСЛОВНОЙ ОПТИМИЗАЦИИ ФУНКЦИИ N ПЕРЕМЕННЫХ.
- •4.1. Метод штрафных функций
- •РАЗДЕЛ II. ЧИСЛЕННЫЕ МЕТОДЫ УРАВНЕНИЙ МАТЕМАТИЧЕСКОЙ ФИЗИКИ
- •1. МАТЕМАТИЧЕСКИЕ МОДЕЛИ И ЧИСЛЕННЫЕ МЕТОДЫ
- •1.1. Как исследуются физические явления и решаются задачи
- •1.2. Как оценивается погрешность вычислений?
- •1.3. Откуда возникают погрешности расчетов?
- •1.4. Итерационные методы решения задач
- •2. КЛАССИФИКИЦИЯ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ
- •2.1. Определение дифференциальных уравнений
- •2.2. Постановка задач для обыкновенных ДУ
- •2.4. Как получают дифференциальные уравнения
- •2.5. Подобие физических явлений, безразмерные переменные
- •3. ТЕОРЕТИЧЕСКИЕ ОСНОВЫ МЕТОДА СЕТОК
- •3.1. Теоретические основы метода сеток
- •3.2 Решение систем обыкновенных дифференциальных уравнений
- •4. РЕШЕНИЕ КРАЕВОЙ ЗАДАЧИ ДЛЯ ДИФФЕРЕНЦАЛЬНЫХ УРАВНЕНИЙ 2-ГО ПОРЯДКА МЕТОДОМ ПРОГОНКИ
- •5. РЕШЕНИЕ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ В ЧАСТНЫХ ПРОИЗВОДНЫХ
- •5.1. Одномерное нестационарное уравнение теплопроводности
- •5.2. Задача Дирихле для двухмерного уравнения Пуассона
2. КЛАССИФИКИЦИЯ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ
2.1. Определение дифференциальных уравнений
Любой физический закон наиболее точно может быть выражен только в виде некоторого (часто дифференциального) уравнения и его решение при заданных исходных данных позволяет предсказывать всевозможные физические явления и результаты эксперимента. Дифференциальными уравнениями можно описать широкий круг явлений, таких, как поведение планет и галактик, движение воздушных потоков в атмосфере, течения в океанах, распространение радиоволн, протекание химических реакций, процессы выращивания кристаллов, поведение электронов в атомной решетке полупроводника.
В общем случае дифференциальное уравнение (ДУ) – это уравнение, связывающее искомую функцию u(x, y, z,t) , ее частные производные и независи-
мые переменные x, y, z,t : |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Φ x, y, z,t,u(x, y, z,t), |
∂u |
, |
∂u |
, |
∂u |
, |
∂u , |
∂2u2 |
, |
∂2u2 |
, |
∂2u2 |
, |
∂2u2 |
= f (x, y, z,t) . |
(2.1) |
|
∂x |
|
∂y |
|
∂z |
|
∂t |
∂x |
|
∂y |
|
∂z |
|
∂t |
|
|
Максимальный порядок входящих в это уравнение производных называется порядком ДУ.
Описание конкретной прикладной задачи может приводить к одному ДУ
порядка m или к эквивалентной системе ДУ, суммарный порядок которой равен m. Установлены теоремы о том, что систему ДУ порядка m можно привести к эквивалентной системе из m ДУ первого порядка, причем если m - четное, то и к системе из m ДУ второго порядка. В настоящее время теория систем ДУ первого и второго порядка детально разработана. Для решения таких ДУ имеются хорошо апробированные как приближенные, так и численные методы. Это об-
стоятельство позволяет при создании прикладных пакетов программ сосредо-
точить усилия на разработке программ решения систем ДУ первого и второго
порядка.
По количеству независимых переменных ДУ разделяют на два класса: обыкновенные ДУ (ОДУ), которые описывают процессы, зависящие от одной
независимой переменной u(x) , и ДУ в частных производных (ДУЧП), которые
описывают поля физических величин, зависящих более чем от одной независимой переменной.
2.2. Постановка задач для обыкновенных ДУ
Систему ОДУ можно привести к системе, разрешенной относительно пер-
вой или второй производной. Для них имеются хорошо разработанные методы
и стандартные программы.
Система ОДУ первого порядка:
30
du1 |
= |
f (x,u ,u ,...,u |
|
); |
|
|
|
|
|
|
|
|
|
|
|
|||||||||
dx |
|
1 |
|
1 |
2 |
|
|
|
m |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
.......................................... |
|
|
|
|
|
|
|
|
|
|
|
|
|
или коротко du = fr(x,ur). |
|
(2.2) |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dx |
|
|
|
|
|
du |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
= fm (x,u1,u2 ,...,um ). |
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
m |
|
|
|
|
|
|
|
|
|
|
||||||||||||
dx |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Система ОДУ второго порядка: |
|
|
|
|
|
|
|
|
|
|
||||||||||||||
∂ |
(g |
∂u1 ) + q |
∂u1 |
+ p u = f (x,u |
|
,...,u |
|
, ∂u2 |
,..., ∂um ); |
|
|
|||||||||||||
|
∂x |
|
m |
|
|
|||||||||||||||||||
|
1 |
∂x |
|
1 ∂x |
|
|
|
|
1 1 |
|
2 |
|
∂x |
|
∂x |
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(2.3) |
..................................................................................... |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂ |
|
|
∂u |
|
|
|
∂u |
|
|
|
|
|
|
|
|
|
|
∂u |
|
∂u |
|
|
|
|
|
(gm |
|
|
) + qm ∂x |
|
+ pmum = f (x,u1,...,um−1 |
|
|
|
|
). |
||||||||||||
|
∂x |
∂x |
|
, ∂x |
,..., ∂x |
|||||||||||||||||||
|
|
|
|
|
m |
|
|
|
|
m |
|
|
|
|
|
|
|
|
1 |
|
|
m−1 |
|
|
Требуется найти ur ={u (x),...,u |
m |
(x)} |
для значений х из заданного интервала |
|||||||||||||||||||||
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a ≤ x ≤ b .
Известно, что системы (2.2) и (2.3) имеют бесконечное множество решенийr , семейство которых зависит от набора произвольных параметров
C = (c1,...,cm ) , их количество равно порядку системы m для 2.2) и 2m для (2.3).
Это семейство решений описывает многообразие реализаций физического процесса, математической моделью которого являются системы ОДУ (2.2), (2.3). Для выделения одной искомой реализации среди этого многообразия необходимо наложить дополнительные условия на функцию ur(x) , количество
этих условий равно порядку системы ОДУ.
В зависимости от способа постановки дополнительных условий можно выделить два основных типа задач для ОДУ.
Задача Коши (задача с начальным условием), когда все условия заданы в начале отрезка интегрирования [a,b] (при x=a). Эта задача чаще всего ставится
для системы (3.2) в виде |
|
|
или коротко ur(a) = ur0 . |
|
|
||||
u (a) = u0 |
; ... u |
m |
(a) = u0 |
|
(2.4) |
||||
1 |
1 |
|
m |
|
|
|
|
|
|
Краевая задача, когда условия заданы на обоих концах отрезка [a,b]. Эта |
|||||||||
задача обычно ставится для ДУ второго порядка (3.3) в общем виде: |
|
||||||||
αa dui (a) |
+ βau (a) = γ a |
; αb dui (b) |
+ βbu (b) = γ b ; i =1..m. |
(2.5) |
|||||
i |
dx |
i i |
i |
i |
dx |
i i |
i |
|
|
|
|
|
|
|
|
|
|
||
Имеется ряд установленных теорем о том, что при определенных условиях как задача Коши (2.2), (2.4), так и краевая задача (2.3), (2.5) имеют единственное решение. Прежде чем приступать к решению ОДУ, необходимо удостовериться в выполнении этих условий. Обычно, если физическая задача поставлена правильно, то и математическая ее постановка в «разумном» диапазоне параметров будет иметь решение.
2.3. Постановка задач для ДУ в частных производных (ДУЧП)
Наш мир устроен так удивительно, что большинство окружающих нас физических полей описывается ДУЧП второго порядка или, если имеется не-
31
сколько взаимосвязанных полей, - системами таких ДУ. Более того, если процессы не слишком интенсивны, то при их описании можно ограничиться линейными дифференциальными уравнениями второго порядка. Самый общий вид такого уравнения:
∑aij |
∂2u |
|
+∑bi |
∂u |
+ cu = f ; x1 = x, x2 = y, x3 = z, x4 = t. |
(2.6) |
∂x ∂x |
|
∂x |
||||
ij |
j |
i |
|
|
||
|
i |
|
i |
|
|
Оказывается, что свойства решений (2.6) существенно зависят от коэффициентов, стоящих при старших производных. Причем с помощью соответствующей замены независимых переменных уравнение (2.6) может быть приведено к одному из следующих трех типов.
1. Параболические
∂u = ∂2u + ∂2u + ∂2u + f . ∂t ∂x2 ∂y2 ∂z2
Более общий случай
a |
∂u |
= |
∂ |
g |
∂u + |
∂ |
g |
∂u + |
∂ |
g |
∂u + f . |
(2.7) |
∂t |
|
|
|
|||||||||
|
|
∂x |
∂x |
∂y |
∂y |
∂z |
∂z |
|
||||
К этому типу сводятся уравнения для физических полей, описывающих диссипативные процессы: теплопроводность – поле температуры Т (x, y, z, t); диффузия – поле концентрации С (x, y, z, t) (например, сахара в растворе). Уравнение (3.7) получило название уравнение теплопроводности.
Процесс распространения возмущения, возникающего в таком поле, имеет затухающий характер и для одномерного случая этот процесс иллюстрирует рис. 2.1:
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
u |
|
|
|
|
t=0 |
|
|
|
|
|
∂u |
|
|
∂2u |
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
= |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
∂t |
∂x2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t=t1>0 |
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t=t2>t1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
0 |
z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
2. Гиперболические |
|
|
|
|
|
Рис. 3.1 |
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
∂2u |
= |
∂2u |
+ |
∂2u |
+ |
∂2u |
+ f . |
|
|
|
(2.8) |
||||||||||
|
|
∂t2 |
|
∂x2 |
∂y2 |
∂z2 |
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
Более общий случай |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
a |
∂2u |
|
∂ |
|
|
|
∂u |
|
+ |
∂ |
|
∂u |
∂ |
∂u |
|
|
|||||||
∂t |
2 = |
∂x |
g |
∂x |
|
|
|
g |
|
+ |
g |
+ f . |
|
|
|||||||||
|
|
|
|
|
|
|
∂y |
∂y |
∂z |
∂z |
|
|
|||||||||||
К этому типу сводятся уравнения для физических полей, описывающих волно- |
|||||||||||||||||||||||
вые процессы – |
электромагнитные волны E, Hr (x, y, z,t) , |
звук, поле вероятно- |
|||||||||||||||||||||
стей в |
квантовой |
|
физике. |
Уравнение |
(2.8) получило |
название |
уравнение |
||||||||||||||||
Д’Аламбера, или неоднородное волновое уравнение. |
|
|
|||||||||||||||||||||
32
|
Процесс распространения возмущения, возникающего в таком поле, имеет |
|||||||||||||||||
вид бегущих в разных направлениях возмущений, причем в одномерном случае, |
||||||||||||||||||
представленном на рис. 2.2, возмущения не затухают. |
|
|||||||||||||||||
|
∂2u |
= |
∂2u |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
∂t2 |
∂z2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
u |
|
|
|
|
|
|
|
|
t=0 |
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
t=t2 |
t=t1 |
t=t1 |
t=t2 |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рис. 3.2 |
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
3. Эллиптические (не изменяющиеся во времени процессы): |
|
|||||||||||||||||
|
|
|
|
|
|
∂2u |
+ |
∂2u |
+ |
∂2u |
= f . |
|
(2.9) |
|||||
|
|
|
|
|
|
∂x2 |
∂y2 |
|
∂z2 |
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
В более общем случае: |
|
|
|
|
|
|
|
|
||||||||||
|
∂ |
|
∂u |
+ |
∂ |
|
|
∂u |
|
|
+ |
∂ |
|
∂u |
|
|
||
|
|
g |
|
∂y |
g |
∂y |
|
|
g |
= f . |
|
|
||||||
|
∂x |
|
∂x |
|
|
|
|
|
|
∂z |
|
∂z |
|
|
||||
К этому типу сводятся уравнения для многих физических полей: электростати- |
||||||||||||||||||
ческие и магнитостатические поля, установившееся температурное поле и поле |
||||||||||||||||||
концентраций и др. |
|
|
|
|
|
|
|
|
|
|
||||||||
|
Уравнение (2.9) получило название уравнение Пуассона; если |
f =0 – |
||||||||||||||||
уравнение Лапласа. |
|
|
|
|
|
|
|
|
|
|||||||||
|
При постановке задач для ДУ в частных производных обычно требуется |
|||||||||||||||||
найти распределение u (x, y, z,t) , удовлетворяющее ДУ в некоторой области |
||||||||||||||||||
пространства Ω с границей Г, причем область может быть как ограниченной, |
||||||||||||||||||
так и полуограниченной. |
|
|
|
|
|
|
|
|
||||||||||
|
Оказывается, что общее решение ДУЧП в отличие от ОДУ не зависит от |
|||||||||||||||||
конечного набора параметров, а содержит произвольные дифференцируемые |
||||||||||||||||||
функции. Так, легко показать, что решением одномерного волнового уравнения |
||||||||||||||||||
∂2u |
= |
∂2u |
является произвольная, дважды дифференцируемая функция с аргу- |
|||||||||||||||
∂t2 |
∂z |
2 |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
ментом t ± z : u = f (t ± z) |
(бегущие влево и вправо волны, см. рис. 2.2). Это зна- |
|||||||||||||||||
чит, что класс функций, являющихся решением ДУЧП, гораздо шире, чем у |
||||||||||||||||||
ОДУ. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
Для выделения единственного искомого решения из класса решений необ- |
|||||||||||||||||
ходимо задать дополнительные условия. Чаще всего ими являются начальные |
||||||||||||||||||
условия (для нестационарных полей, описываемых уравнениями (2.7), (2.8)) и |
||||||||||||||||||
граничные (краевые) условия. |
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
33 |
|
|
|
