Решение краевых задач для обыкновенных дифференциальных уравнений и систем. Основные положения - umotnas.ru o_O
Главная
Поиск по ключевым словам:
страница 1страница 2
Похожие работы
Название работы Кол-во страниц Размер
Решение задачи Коши для обыкновенных дифференциальных уравнений 1 33.74kb.
Iv. Численные методы решения обыкновенных дифференциальных уравнений. 1 78.46kb.
Сингулярная краевая задача типа николетти с кусочно-непрерывными... 1 105.68kb.
1. Качественный анализ краевых задач и экстремальных задач для уравнений... 1 130.19kb.
Программа курса "Дифференциальные уравнения" 1 108.99kb.
4 Вычислительные методы решения дифференциальных уравнений 1 110.22kb.
Не разрешенные относительно производной 1 85.76kb.
Отчет по лр№1: «Решение систем линейных алгебраических уравнений... 1 64.28kb.
Задача Дирихле: а внутренняя задача 1 39.4kb.
Корректность краевых задач для сингулярных параболических уравнений... 1 150.81kb.
Определение функций источника систем уравнений составного типа для... 1 48.02kb.
Методические указания к выполнению контрольной работы Направление... 1 286.09kb.
Викторина для любознательных: «Занимательная биология» 1 9.92kb.

Решение краевых задач для обыкновенных дифференциальных уравнений и систем. Основные - страница №1/2


  1. Оптимизация. Критерии оптимальности; целевая функция. Классификация. Модели оптимизационных задач. Примеры.

  2. Решение краевых задач для обыкновенных дифференциальных уравнений и систем. Основные положения.

  3. Метод стрельбы. Схемы дихотомии и секущих.

  4. Метод стрельбы для линейной краевой задачи.

  5. Метод конечных разностей, или метод сеток.

  6. Полуаналитические методы решения краевой задачи. Метод коллокации.

  7. Полуаналитические методы решения краевой задачи. Метод Галеркина.

  8. Разностные схемы для уравнений в частных производных. Основные понятия.

  9. Сходимость, аппроксимация и устойчивость разностных схем

  10. Разностные схемы для уравнений параболического типа. Решение задачи Коши.

  11. Устойчивость двухслойных разностных схем для уравнений параболического типа.

  12. Построение разностной аппроксимации для уравнения Пуассона.

  13. Различные краевые задачи и аппроксимация граничных условий.

  14. Построение разностной схемы в случае задачи Дирихле для уравнения Пуассона.

  15. Метод матричной прогонки. Итерационный метод решения разностной схемы для задачи Дирихле.

  16. Разностные схемы для уравнений гиперболического типа. Решение задачи Коши.

  17. Разностные схемы для уравнений гиперболического типа. Решение смешанной задачи.

  18. Общие понятия метода конечных элементов.

  19. Дискретизация области и нумерация узлов.

  20. Линейные интерполяционные полиномы.

  21. Одномерный симплекс-элемент.

  22. Двумерный симплекс-элемент.

  23. Местная система координат.

  24. Двумерные L-координаты.

  25. Объединение элементов в ансамбль.

  26. Вывод уравнений для элементов с помощью метода Галеркина.

  27. Пример расчета одномерного температурного поля в однородном стержне.




  1. Оптимизация. Критерии оптимальности; целевая функция. Классификация. Модели оптимизационных задач. Примеры.


Оптимизация - целенаправленная деятельность, заключающаяся в получении наилучших результатов при соответствующих условиях.

Критерий оптимальности — характерный показатель решения задачи, по значению которого оценивается оптимальность найденного решения, то есть максимальное удовлетворение поставленным требованиям. В одной задаче может быть установлено несколько критериев оптимальности.

Функция, связывающая цель (оптимизируемую переменную) с управляемыми переменными в задаче оптимизации.

В широком смысле целевая функция есть математическое выражение некоторого критерия качества одного объекта (решения, процесса и т.д.) в сравнении с другим. 

Методы оптимизации классифицируют в соответствии с задачами оптимизации:



  • Локальные методы: сходятся к какому-нибудь локальному экстремуму целевой функции. В случае унимодальной целевой функции, этот экстремум единственен, и будет глобальным максимумом/минимумом.

  • Глобальные методы: имеют дело с многоэкстремальными целевыми функциями. При глобальном поиске основной задачей является выявление тенденций глобального поведения целевой функции.

Модели, способные воспроизводить только качественную сторону изучаемого процесса, называются качественными. Более сложными моделями являются количественные, или математические, модели, с помощью которых можно предсказывать конкретные численные значения характеристик явления. Так, например, модель погоды, основанная на народных приметах, будет качественной. Модели же, используемые в гидрометцентрах и позволяющие предсказывать численные значения характеристик погоды (температура, скорость ветра и т.п.), относятся к количественным моделям. Математическая модель – это совокупность уравнений (алгебраических, дифференциальных, интегральных), описывающих процессы, происходящие в моделируемом явлении. Каждое такое уравнение само по себе есть модель процесса более низкого уровня. Модели же сложных явлений часто строятся, как из кирпичиков, из моделей более простых. Результатом математического моделирования являются формулы (если удается получить аналитическое решение уравнений модели), которые позволяют рассчитать характеристики моделируемого процесса или таблицы значений этих характеристик (если аналитическое решение получить не удается, либо оно не удобно). Чаще всего встречается именно последний случай, когда решение можно получить лишь в виде таблицы с помощью каких-либо вычислительных процедур. Такого рода модели носят специальное название  численных моделей. Специфика этих моделей состоит в том, что они позволяют получить принципиально приближенное решение, что связано с дискретностью вычислительных устройств, ограниченным объемом памяти, конечным быстродействием и т.п.




  1. Решение краевых задач для обыкновенных дифференциальных уравнений и систем. Основные положения.


Краевая задача – это задача отыскания частного решения системы

, (2.1)

на отрезке , причем дополнительные условия налагаются на значения функций более чем в одной точке этого отрезка.

Заметим, что дифференциальное уравнение порядка р

(2.3)

где − производная порядка

может быть сведено к системе дифференциальных уравнений вида (2.1) заменой переменных

(2.4)

Действительно, по замене (2.4)




и уравнение (2.3) сведется к следующей системе вида (2.1):


Здесь последнее уравнение получено подстановкой (2.4) в (2.3).

Заметим, что общая краевая задача (2.1) может:


  1. не иметь решений;

  2. иметь единственное решение;

  3. иметь несколько и даже бесконечно много решений.

Рассмотрим более подробно важный частный случай, когда дифференциальное уравнение и краевые условия линейны. Такая краевая задача называется линейной краевой задачей.

Линейное дифференциальное уравнение n-го порядка сокращенно можно записать в виде



(2.6)

Где

причем обычно предполагается, что − известные непрерывные функции на данном отрезке .

Для простоты будем предполагать, что в краевые условия входят две абсциссы и − концы отрезка . Такие краевые условия называются двухточечными. Краевые условия называются линейными, если они имеют вид

(2.7)

Где и − заданные постоянные, причем


Например, краевые условия, приведенные в предыдущих примерах, линейны.



Линейными краевыми условиями являются также условия периодичности, которые в случае дифференциального уравнения второго порядка имеют вид




Линейная краевая задача называется однородной, если:

  • во-первых, то есть дифференциальное уравнение (2.6) однородно, и,

  • во-вторых, то есть имеют место однородные краевые условия.

В противном случае краевая задача (2.6)-(2.7) называется неоднородной.


  1. Метод стрельбы. Схемы дихотомии и секущих.


Это численный метод, заключающийся в сведении краевой задачи к решению последовательности задач Коши для той же системы дифференциальных уравнений. Рассмотрим его на примере простейшей задачи для системы двух уравнений первого порядка с краевыми условиями достаточно общего вида

(2.15а)

; . (2.15б)

Выберем произвольно значение , рассмотрим левое краевое условие как алгебраическое уравнение




и найдем из него . Возьмем значения в качестве начальных условий задачи Коши для системы (2.15а) и проинтегрируем эту задачу Коши любым численным методом. При этом получим решение, зависящее от , как от параметра.

Значение выбрано так, что найденное решение удовлетворяет левому краевому условию (2.15б). Однако правому краевому условию это решение, скорее всего, не удовлетворяет: при его подстановке левая часть краевого условия в точке b, рассматриваемая как функция параметра

, (2.16)

не обратится в нуль.

Необходимо каким-либо способом менять параметр , пока не подберем такое значение, для которого с требуемой точностью. Таким образом, решение краевой задачи (2.15) сводится к нахождению корня уравнения

. (2.17)

Простейшим методом его решения является метод дихотомии (деления отрезка пополам).

Делают пробные "выстрелы" – расчеты с наудачу выбранными значениями до тех пор, пока среди величин не окажется разных по знаку. Пара таких значений образует "вилку". Деля ее последовательно пополам до получения нужной точности, производим "пристрелку" параметра . Благодаря этому процессу весь метод получил название стрельбы.

Однако нахождение каждого нового значения функции требует численного интегрирования системы (2.15а), то есть достаточно трудоемко.

Поэтому корень уравнения (2.17) желательно находить более быстрым методом.

Попробуем сделать это методом Ньютона:

. (2.18)

Однако вычисление производной затруднительно и лучше ее заменить разностным отношением

. (2.19)

Подставляя (2.19) в (2.18), получим итерационную формулу метода секущих:

. (2.20)

В методе секущих первые два расчета делают с наудачу выбранными близкими значениями и , а следующие значения параметра вычисляют по формуле (2.20).

Заметим, что этот метод быстро сходится вблизи корня. Сходимость вдали от корня зависит от того, насколько удачно выбрано нулевое приближение.




  1. Метод стрельбы для линейной краевой задачи.


Рассмотрим теперь линейную краевую задачу, решение которой методом стрельбы особенно просто:

(2.21а)

(2.21б)

Воспользуемся известным результатом из теории дифференциальных уравнений: общее решение линейной неоднородной системы равно сумме ее какого-нибудь частного решения и общего решения соответствующей однородной системы.

Найдем частное решение неоднородной системы (2.21а), положив в левом условии (2.21б) . Обозначим это частное решение через и заметим, что .

Рассмотрим теперь соответствующую однородную систему




с однородными начальными условиями

.

Вычислим решение этой задачи Коши и обозначим его через . Рассмотрим функции и . Очевидно, что в точке a эти функции удовлетворяют краевому условию:




Поэтому общее решение неоднородной задачи Коши, удовлетворяющее левому краевому условию (2.21б) дается следующим однопараметрическим семейством

(2.22)

Значение параметра C выбираем так, чтобы удовлетворить правому краевому условию (2.21б):

. (2.23)

Искомое решение краевой задачи (2.21) тогда находится по формуле (2.23).

Итак, решение линейной краевой задачи требует только двух “выстрелов” – вспомогательные задачи Коши решаются дважды.

  1. Метод конечных разностей, или метод сеток.


Рассмотрим линейную краевую задачу

(2.24)

(2.25)

,

где , , и непрерывны на [a, b].

Разобьем отрезок [a, b] на n равных частей длины, или шага

.

Точки разбиения , называются узлами, а их совокупность – сеткой на отрезке [a, b]. Значения в узлах искомой функции и ее производных обозначим соответственно через

.

Введем обозначения




Заменим производные так называемыми односторонними конечно-разностными отношениями:

(2.26)

Формулы (2.26) приближенно выражают значения производных во внутренних точках интервала [a, b].

Для граничных точек положим

. (2.27)

Используя формулы (2.26), дифференциальное уравнение (2.24) при , (i=1, 2,..., n–1) приближенно можно заменить линейной системой уравнений

(2.28)

Кроме того, в силу формул (2.27) краевые условия (2.25) дополнительно дают еще два уравнения:

. (2.29)

Таким образом, получена линейная система n+1 уравнений с n+1 неизвестными , представляющими собой значения искомой функции в узлах сетки. Система уравнений (2.28), (2.29), заменяющая приближенно дифференциальную краевую задачу (2.24), (2.25) обычно называется разностной схемой. Решить эту систему можно каким-либо общим численным методом. Однако схема (2.28), (2.29) имеет специфический вид и ее можно эффективно решить специальным методом, называемым методом прогонки. Специфичность системы заключается в том, что уравнения ее содержат три соседних неизвестных и матрица этой системы является трехдиагональной.

Преобразуем уравнения (2.28):

. (2.30)

Введя обозначения




Получим , (i=0, 1,..., n-2). (2.31)

Краевые условия по-прежнему запишем в виде

. (2.32)

Метод прогонки состоит в следующем.

Разрешим уравнение (2.31) относительно :

. (2.33)

Предположим, что с помощью полной системы (2.31) из уравнения исключен член, содержащий. Тогда уравнение (2.33) может быть записано в виде

, (2.34)

где и должны быть определены. Найдем формулы для этих коэффициентов. При i=0 из формулы (2.33) и краевых условий (2.32) следует, что




Исключая из этих двух уравнений , найдем

.

Выразим теперь отсюда :

(2.35)

Но, согласно формуле (2.34),

(2.36)

Сравнивая теперь (2.35) и (2.36), найдем, что

(2.37)

Пусть теперь i >0, то есть i=1, 2,..., n2. Выражая по формуле (2.34), получим:

.

Подставляя это в формулу (2.33), будем иметь

.

Разрешая полученное уравнение относительно, находим

, или

. (2.38)

Отсюда, сравнивая формулы (2.34) и (2.38), получаем для коэффициентов и рекуррентные формулы:

(2.39)

Так как и уже определены по формулам (2.37), то, используя формулы (2.39), можно последовательно определить коэффициенты и до и включительно. Эти вычисления называются прямым ходом метода прогонки.

Из формулы (2.33) при i=n2 и второго краевого условия (2.32) получаем






Разрешая эту систему относительно, будем иметь

. (2.40)

Теперь, используя (2.34) и первое краевое условие (2.32), мы можем последовательно найти . Это − обратный ход метода прогонки.

Итак, получаем следующую цепочку:

(2.41)

Для простейших краевых условий формулы для и упрощаются. Полагая в этом случае


из формул (2.37), (2.40), (2.41) будем иметь




Рассмотренный нами подход сводит линейную краевую задачу к системе линейных алгебраических уравнений

Более точную разностную схему можно получить, если при переходе от линейной краевой задачи к конечно-разностным уравнениям воспользоваться центральными формулами для производных:

, (2.42)

, (2.43)

i=1, 2,..., n.

то есть формула (2.42) имеет второй порядок точности относительно шага сетки h. Подставляя выражения (2.42), (2.43) в задачу (2.24), (2.25) и выполняя некоторые преобразования, получим следующую систему:



(2.44)

где .

Система (2.44) снова трехдиагональная и ее решение также можно получить методом прогонки. Его алгоритм здесь будет выглядеть так. Сначала находят коэффициенты

(2.45)

Затем определяют коэффициенты по следующим рекуррентным формулам:

(2.46)

Обратный ход начинается с нахождения :

(2.47)

После этого находим по формулам:

, (2.48)

. (2.49)

  1. Полуаналитические методы решения краевой задачи. Метод коллокации.


Пусть необходимо определить функцию, удовлетворяющую линейному дифференциальному уравнению

(2.50)

и линейными краевыми условиями

, (2.51)

причем

Выберем некоторую совокупность линейно независимых функций

(2.52)

которую назовем системой базисных функций.

Пусть функция удовлетворяет неоднородным краевым условиям

(2.53)

а остальные функции удовлетворяют соответствующим однородным краевым условиям:

. (2.54)

Если краевые условия (2.51) однородны (A=B=0), то можно положить и рассматривать лишь систему функций .

Будем искать приближенное решение краевой задачи (2.50), (2.51) в виде линейной комбинации базисных функций

. (2.55)

Тогда функция y удовлетворяет краевым условиям (2.51). В самом деле, в силу линейности краевых условий имеем




и аналогично

Составим функцию . Подставляя сюда вместо y выражение (2.55), будем иметь

.(2.56)

Если при некотором выборе коэффициентов ci выполнено равенство

при

то функция y является точным решением краевой задачи (2.50), (2.51). Однако подобрать так удачно функции и коэффициенты ci в общем случае не удается. Поэтому ограничиваются тем, что требуют, чтобы функция обращалась в нуль в заданной системе точек из интервала [a, b], которые называются точками коллокации. Сама функция R называется невязкой уравнения (2.50). Очевидно, что в точках коллокации дифференциальное уравнение (2.50) будет удовлетворено точно, и невязка в этих точках равна нулю.

Итак, метод коллокации приводит к системе линейных уравнений

. (2.57)

Из системы (2.57) в случае ее совместности можно определить коэффициенты , после чего приближенное решение краевой задачи дается формулой (2.55).


  1. Полуаналитические методы решения краевой задачи. Метод Галеркина.


Пусть дано дифференциальное уравнение с линейными краевыми условиями

, (2.62)

(2.63)

Будем искать приближенное решение этой краевой задачи в виде суммы

(2.64)

где – некоторая непрерывная функция, удовлетворяющая неоднородным краевым условиям (2.63), а – какая-то система линейно независимых функций, удовлетворяющих однородным краевым условиям

(2.65)

и, кроме того функции при образуют в классе функций c2[a, b], удовлетворяющих условиям (2.65), полную систему.

Заметим, что свойство полноты понимается следующим образом.

Обозначим через G класс функций y(x), принадлежащих c2[a, b] (то есть дважды непрерывно дифференцируемых на [a, b]) и удовлетворяющих граничным условиям (2.65). Говорят, что система функций полна в классе G, если для любого и любой функции можно указать такое n и такие параметры , что имеет место неравенство




где

Это означает, что для любой допустимой функции найдется такая функция , которая на [a, b] будет сколь угодно точно приближать функцию y(x) вместе с ее производными и .

Возвращаясь теперь к задаче (2.62), (2.63), видим, что если бы мы нашли такую функцию y(x), удовлетворяющую условиям (2.63), и чтобы было ортогонально при любых , то это означало бы, что , и задача (2.62), (2.63) была бы решена. Если же ортогональность есть только при , то в разложении по системе входят и более старшие коэффициенты, то есть

Метод Галеркина состоит в том, что решение задачи (2.62), (2.63) ищется в виде (2.64), причем требуют ортогональности к функциям полной системы для , то есть

(2.68)

где

Это дает алгебраическую систему уравнений для определения коэффициентов ak. Найдя из нее коэффициенты, получим приближенное решение. В методе Галеркина функция должна удовлетворять краевым условиям (2.63). Поэтому можно выбрать в виде

,

и коэффициенты найти как решение системы уравнений




Таким же образом отыскиваются функции . Выберем, например, полную систему в виде многочленов последовательных степеней:

.

Коэффициенты найдем из однородных краевых условий (2.65)

(2.65а)

при всех .

Так, для и условия (2.65а) принимают вид:




В этой системе из двух уравнений три неизвестных: и . Одну из них можно выбрать произвольно, положив, например, . Аналогично отыскивают коэффициенты для .

Для простых условий вида то есть функции можно вычислять по правилу




Или

Отметим, что при нелинейном краевом условии вида, например, линейная комбинация (2.64) с произвольными коэффициентами ak уже не будет удовлетворять этому краевому условию. Поэтому метод Галеркина применим только к задачам с линейными краевыми условиями, хотя допустим и нелинейный оператор L.


  1. Разностные схемы для уравнений в частных производных. Основные понятия.

Пусть D − некоторая область изменения независимых переменных x, y, ограниченная контуром Г. Говорят, что в области D задано линейное дифференциальное уравнение второго порядка для функции , если для любой точки из области D имеет место соотношение

(3.1)

где , ... − коэффициенты, − свободный член уравнения. Эти функции известны и их обычно считают определенными в замкнутой области . Обозначим Уравнение называется эллиптическим, параболическим или гиперболическим в D, если соответственно выполняются условия для всех D.

В зависимости от типа дифференциального уравнения по разному ставятся граничные и начальные условия, связанные с этим уравнением. Далее мы будем рассматривать частные случаи уравнения (3.1):

  1. уравнение Пуассона (эллиптическое уравнение)

;

  1. уравнение теплопроводности

(параболическое уравнение)

;

  1. волновое уравнение (гиперболическое уравнение)

.


  1. Сходимость, аппроксимация и устойчивость разностных схем


Пусть u есть решение дифференциального уравнения

(3.2)

заданного в области D. Рассмотрим некоторое множество Dh состоящее из изолированных точек , принадлежащих замкнутой области . Число точек в Dh будем характеризовать величиной h; чем меньше h, тем большим будет число точек в Dh. Множество Dh называется сеткой, а точки Dhузлами сетки. Функция, определенная в узлах, называется сеточной функцией.

Обозначим через U пространство непрерывных в D функций . Через Uh обозначим пространство, образованное совокупностью сеточных функций , определенных на Dh. В методе сеток осуществляется замена пространства U на пространство Uh.

Пусть − точное решение уравнения (3.2) и принадлежит U. Поставим задачу отыскания значений . Эти значения в совокупности образуют таблицу, в которой число значений равно числу точек в Dh. Точно поставленную задачу удается решить редко. Как правило, можно вычислить некоторые сеточные значения , относительно которых можно думать, что




Величины называются приближенными сеточными значениями решения . Для их вычисления строят систему численных уравнений, которую мы будем записывать в виде

(3.3)

где Lh есть разностный оператор, соответствующий оператору L, Если то образуется по F аналогично тому, как Uh образовывалось по U. Формулу (3) будем называть разностной схемой.

Пусть в линейных пространствах Uh и Fh введены соответственно нормы , которые являются сеточными аналогами норм в исходных пространствах.

Будем говорить, что разностная схема (3.3) является сходящейся, если при выполняется условие

.

Если выполняется условие

,

где c − постоянная, не зависящая от h и s>0, то говорят, что имеет место сходимость со скоростью порядка s относительно шага h.

Говорят, что разностная схема (3.3) аппроксимирует задачу (3.2) на решении , если




Величина называется погрешностью аппроксимации или невязкой разностной схемы. Если , где M − константа, не зависящая от h и , то говорят, что разностная схема (3.3) аппроксимирует задачу (3.2) на решении с погрешностью порядка относительно шага h.

Разностная схема (3.3) называется устойчивой, если существует такое , что для всех и любых выполняются условия:

1) разностная схема (3.3) имеет единственное решение;

2) где M − постоянная, не зависящая от и.

Иначе говоря, разностная схема является устойчивой, если ее решение непрерывно зависит от входных данных. Устойчивость характеризует чувствительность схемы к различного рода погрешностям, она является внутренним свойством разностной задачи, и это свойство не связывается непосредственно с исходной дифференциальной задачей, в отличие от сходимости и аппроксимации. Между понятиями сходимости, аппроксимации и устойчивости существует связь. Она состоит в том, что из аппроксимации и устойчивости следует сходимость, что отражается в следующей теореме.

Обычно применение метода сеток заключается в следующем:

  1. Вначале указывается правило выбора сетки, то есть указывается метод замены области Д и контура Г некоторой сеточной областью. Чаще всего сетка выбирается прямоугольной и равномерной.

  2. Затем указывается и строится конкретно одна или несколько разностных схем. Проверяется условие аппроксимации и устанавливается ее порядок.

  3. Доказывается устойчивость построенных разностных схем. Это один из наиболее важных и сложных вопросов. Если разностная схема обладает аппроксимацией и устойчивостью, то о сходимости судят по доказанной теореме.

  4. Рассматривается вопрос численного решения разностных схем.

В случае линейных разностных схем это будет система линейных алгебраических уравнений. Порядок таких систем может быть очень большим.


  1. Разностные схемы для уравнений параболического типа. Решение задачи Коши.


Рассмотрим задачу Коши для уравнения теплопроводности

, , , (3.5)

с условием на прямой t=0

, . (3.6)

Требуется найти функцию , которая при и удовлетворяла бы уравнению (3.5), а при выполняла бы условие (3.6).

Будем считать, что задача (3.5), (3.6) имеет в верхней полуплоскости единственное решение , непрерывное вместе со своими производными

, i=1, 2 и , k=1, 2, 3, 4.

Запишем задачу (3.5), (3.6) в виде . Для этого достаточно положить







Будем далее считать, что t изменяется в пределах . В рассматриваемом случае

, Г − объединение прямых t=0 и t=T.

Выберем прямоугольную сетку и заменим область сеточной областью . К области отнесем совокупность узлов , где

, , ,

, , , .

Заменим задачу разностной схемой вида . Обозначим через точное значение решения задачи в узле , а через – соответствующее приближенное решение. Имеем







Для замены выражений и воспользуемся формулами численного дифференцирования. Имеем:

, (3.7)

, (3.8)

, (3.9)

(3.10)

Назовем некоторую совокупность узлов, привлекаемых для замены задачи в узле , разностной схемой , шаблоном. Наиболее употребительные шаблоны изображены на рис. 3:






Рис. 3. Явный и неявный шаблоны

Рассмотрим явный двухслойный шаблон. Для него



(3.11)

Здесь мы воспользовались формулами (3.7) и (3.10) и обозначили

.

Введем обозначение

(3.12)

Теперь на основании формул (3.11), (3.12) можно записать разностную схему для задачи :

, (3.13)

где разностный оператор определяется по правилу




Аналогично, если использовать неявный двухслойный шаблон, можно получить такую разностную схему:

, (3.14)

Где




На основании формул (3.11) и (3.13) можно записать

,

где

Аналогично, используя (3.11), (3.10), (3.14), получим

,

.

Разностная схема (3.13) позволяет по значениям решения на нулевом слое, то есть по значениям вычислить значения на первом слое . Для этого достаточно в (3.13) положить n = 0 и произвести вычисления, носящие рекурсионный характер. Потом по значениям можно аналогично при n = 1 вычислить значения и т.д. В силу этого разностную схему (3.13) называют явной.

Разностная схема (3.14) такими свойствами не обладает. Действительно, если мы в (3.14) положим n = 0, то в левой части полученной формулы будет линейная комбинация из значений , в правой части будут значения известной функции и . Для вычисления значений на первом слое в этом случае необходимо решать бесконечную систему линейных уравнений. По этой причине схему (3.14) называют неявной.


  1. Устойчивость двухслойных разностных схем для уравнений параболического типа


Определим норму в пространстве по правилу

.

Рассмотрим явную разностную схему (3.13). Выясним, при каких значениях r, возможна устойчивость этой схемы.

Для доказательства устойчивости надо показать, что разностная схема однозначно разрешима и при любых

,

имеет место оценка

,

где М – постоянная, не зависящая от и и .

Разностная схема (3.13) – явная, и поэтому ее однозначная разрешимость очевидна.

Перепишем формулу в виде

, , (3.17)

.

Пусть выполнено условие

или . (3.18)

Тогда из (3.17) получим:

,

или

. (3.19)

Неравенство (3.19) означает, что при , не превосходит , то есть не возрастает с увеличением n.

Это свойство однородной разностной схемы принято называть принципом максимума. Положим в (3.19) . Это даст

,

,

. . .

.

Заметим, что есть число, независящее от m и n. Просуммировав последние неравенства и, учитывая, что , получим

(3.20)

где обозначено




На основании (3.20) можно записать

или .

Таким образом, разностная схема (3.13) при выполнении условия (3.18), налагаемого на и h, устойчива. Условие (3.18) весьма жестко, ибо из него следует, что . (3.21)

Это приводит к тому, что если мы желаем сохранить устойчивость, то при вычислениях по схеме (3.13) шаг по времени приходится выбирать очень малым.

Обратимся теперь к разностной схеме (3.14), соответствующей шаблону, изображенному на рис. 4,



Рис. 4. Неявный двухслойный шаблон

и перепишем ее в виде

(3.22)

Посмотрим, какие надо проделать вычисления, чтобы, используя формулы (3.22), можно было вычислить, например, значения на первом временном слое со значениями на нулевом временном слое. Положив в формулах (3.22) n=0, получим:

(3.23)

Формулы (3.23) представляют собой бесконечную систему линейных уравнений относительно неизвестных .

Решение таких систем является сложной и трудоемкой задачей, поэтому разностные схемы (3.14) неудобны для задач Коши на бесконечных отрезках и применяется редко. Однако если отрезок оси x, на котором рассматривается задача Коши, конечен, то есть , а на прямых x=a и x=b дополнительно заданы некоторые ограничения на решение , то разностные схемы вида (3.14) оказываются весьма эффективными. В частности, можно показать, что такие схемы являются абсолютно устойчивыми, то есть устойчивыми при любых значениях .

Если, например, на отрезках прямых x=a и x=b, заданы условия , , то вид системы (3.23) существенно изменится:

(3.24)

Формулы (3.24) представляют собой систему M+1 алгебраических уравнений относительно . Матрица этой системы трехдиагональна и ее можно решить методом прогонки. Отсюда ясно, что реализация неявных разностных схем требует больших вычислительных затрат для вычисления решения на одном временном слое, но таких слоев может быть немного из-за того, что в этом случае отсутствуют ограничения на соотношение . Если пользоваться явной разностной схемой, то вычисление решения на следующем слое осуществляется по рекурсионному правилу и связано с минимальными вычислительными затратами, однако из-за ограничения число временных слоев в случае явных схем может быть существенно большим по сравнению с числом временных слоев для неявных схем.



Рассмотрим теперь вопрос о сходимости схемы (3.13). Эта схема аппроксимирует задачу (3.5), (3.6) с погрешностью порядка и устойчива при . Поэтому схема (3.13), по теореме об аппроксимации и устойчивости, будет сходящейся. При этом погрешность для приближенного решения будет величиной порядка .


  1. Построение разностной аппроксимации для уравнения Пуассона.


Рассмотрим в некоторой области D с границей Г уравнение Пуассона

. (3.25)

Выберем прямоугольную сетку по правилу

.

К сеточной области отнесем все узлы, принадлежащие области (рис. 5).

Рис. 5. К построению разностной аппроксимации

Возьмем пятиточечный шаблон



Пользуясь расположением точек в этом шаблоне, разобьем узлы области на две категории: внутренние и граничные. Узел будем считать внутренним, если он сам и четыре соседних точки шаблона принадлежат области (эти узлы обозначены символом ). Обозначим множество внутренних узлов через . Остальные узлы назовем граничными (помечены ) и их множество обозначим через .

Таким образом, .

Очевидно, что разбиение узлов из на внутренние и граничные зависит от выбранного шаблона.

Пусть узел . Замену дифференциального уравнения (3.25) разностным будем осуществлять только во внутренних узлах.

Имеем

. (3.26)

Воспользовавшись аппроксимацией вторых производных, получим




,

(3.27)

,

, .

Пусть и ограничены по абсолютной величине в . Тогда в формуле (3.27) при достаточно малых h и l можно пренебречь членами, содержащими в качестве множителей и , и мы получим искомое разностное уравнение

, (3.28)

где

,

, .

В силу определения невязки уравнения можно получить

, (3.29)

где – точное решение в узлах,

, . (3.30)

При сделанных предположениях относительно и , как видно из (3.30), имеет место оценка

. (3.31)

Здесь М – постоянная, не зависит от .

Оценка (3.31) означает, что разностное уравнение (3.28) аппроксимирует уравнение (3.25) на решение u(x, y) с погрешностью порядка .


  1. Различные краевые задачи и аппроксимация граничных условий.

К уравнениям эллиптического типа, в частности, к уравнению Пуассона (3.25), на границе Г области D присоединяются граничные условия трех видов:

  1. граничные условия 1-го рода:

; (3.32)

  1. граничные условия 2-го рода:

, (3.33)

– производная по внешней нормали;

  1. граничные условия 3-го рода:

, (3.34)

– известные функции.

Если требуется определить функцию , которая в области D удовлетворяет уравнению (3.25), а на границе Г – одному из краевых условий, то говорят, что поставлена граничная задача для эллиптического уравнения.

Задача (3.25), (3.32) называется задачей Дирихле,

задача (3.25), (3.33) – задачей Неймана,

а задача (3.25), (3.34) – смешанной граничной задачей.

Рассмотрим, как можно заменить граничные условия первого рода разностными условиями (рис. 6). Отметим, что граничные условия заменяются условиями на множестве граничных узлов.

Пусть – некоторый узел из ; обозначим его буквой В; – внутренний узел, ближайший к В по направлению x; обозначим его буквой А. Буквой М обозначим точку контура Г, ближайшую к В по направлению x.
Рис. 6. Замена граничных условий первого рода
разностными условиями

Координаты этих точек такие:

.

По условию (3.32) имеем .

Значит, можно положить

(3.35)

для узлов.

Найдем погрешность формулы (3.35).

Имеем




Значит,

Отсюда следует, что погрешность формулы (3.35) будет иметь первый порядок относительно h в предположении, что . Если точки М и В совпадают, то формула (3.35) будет точной.

Точность вычисления можно повысить, если воспользоваться еще значением в точке А.

Имеем:

, (3.36)

. (3.37)

Исключив из (3.36) с помощью (3.37), получим

.

Отбросив здесь величину , найдем разностное граничное условие, аппроксимирующее граничное условие (3.32) в узле с погрешностью : . (3.38)

Обратимся теперь к замене граничного условия 2-го рода разностным уравнением (рис. 7).


Рис. 7. Замена граничных условий второго рода
разностными условиями


Пусть В – граничный узел с координатами , М – ближайшая к В точка контура Г, А – внутренний узел с координатами , С – граничный узел с координатами – внешняя нормаль к Г в точке М. Обозначим угол между и осью Ox через , между и осью Oy – через . Очевидно, что.

По определению имеем



.

Предположим, что в точке В направление нормали такое же, как и в точке М. Поскольку расстояние между В и М есть величина порядка , то это предположение связано с внесением погрешности того же порядка. Значит,

Поэтому окончательно получим




Используя приблизительные сеточные значения, найдем



. (3.39)

Эта формула является разностной аппроксимацией в узле граничного условия второго рода с погрешностью .

Выражения вида (3.39) должны быть записаны для всех граничных узлов , после чего будут получены разностные граничные условия, аппроксимирующие граничные условия (3.33). Процедура замены граничных условий разностными может оказаться весьма громоздкой и сложной, особенно если контур Г имеет непростую форму. Замена граничных условий третьего рода может быть осуществлена с помощью формул вида (3.35), (3.37), (3.38).



  1. Построение разностной схемы в случае задачи Дирихле для уравнения Пуассона.

Пусть в прямоугольной области задано уравнение Пуассона , (3.40)

а на границе Г области D – условие Дирихле . (3.41)

Будем считать, что задача (3.40), (3.41) имеет единственное решение в области и это решение имеет непрерывные в D производные и .Возьмем прямоугольную сетку




Используя формулы (3.28) и (3.35), запишем разностную схему, аппроксимирующую задачу (3.40)–(3.41) с погрешностью порядка :

, (3.42)

где




Разностная схема (3.42) представляет собой систему линейных алгебраических уравнений. Число уравнений этой системы равно .

Столько же и неизвестных: при

Доказательство устойчивости схемы (3.42) сводится к выявлению двух свойств.

  1. Разностная схема ,

– произвольный элемент из , однозначно разрешима.

  1. Имеет место оценка ,

где С – постоянная, не зависящая от .Здесь нормы, как и прежде, определим по правилу.

  1. Метод матричной прогонки. Итерационный метод решения разностной схемы для задачи Дирихле.

Перепишем разностную схему (3.42) так:

(3.55)




Пусть M>N. Введем обозначение

. (3.56)

Положим в формулах (3.55) n=1, ... , N1 и, используя (3.56), запишем систему уравнений (3.55) в векторной форме:

(3.57)

где A – трехдиагональная матрица порядка N 1 вида

,

,

, .

Задача (3.57) аналогична задаче, которая решалась нами методом прогонки. Ее отличие лишь в том, что она имеет векторную форму.

Приведем здесь алгоритм решения задачи (3.57), который называется матричной прогонкой.

1. По формуле полагая , вычисляем матрицы Порядок этих матриц . После этого положим вектор , а затем по формуле

, вычислим векторы




2. Положим . По формуле

вычислим последовательно искомые значения решения задачи (3.57):

.

Итерационный метод решения разностной схемы для задачи Дирихле


Выразим из схемы (3.42) значения :

(3.58)

Значения на границах известны:

(3.59)

В равенстве (3.58) значение выражается через четыре соседних значения по пятиточечному шаблону.

В итерационном методе решения полагают значения во всех внутренних точках области равным некоторым произвольным начальным значениям.

Часто полагают.

Затем с помощью формул (3.58), (3.59) вычисляют новые значения , затем и так далее, до тех пор, пока максимальное отклонение значений сеточных функций на предыдущей и текущей итерациях не станет меньше по модулю, чем некоторая заранее заданная точность . Итак, итерации прекращают, когда выполнится условие

. (3.60)

  1. Разностные схемы для уравнений гиперболического типа. Решение задачи Коши



следующая страница >>