- •Численные методы II
- •Содержание
- •1. Решение линейной краевой задачи для обыкновенных дифференциальных уравнений высших порядков и их систем Справочная информация
- •Метод стрельбы
- •Программное обеспечение
- •Пример решения на пэвм
- •Контрольные задания
- •2. Алгебраическая задача на собственные значения Справочная информация
- •Программное обеспечение
- •Пример решения на пэвм
- •Собственные значения - первые 4
- •3. Решение систем нелинейных алгебраических уравнений Справочная информация
- •Программное обеспечение
- •Пример решения на пэвм
- •Контрольные задания
- •4. Применение метода конечных элементов для решения эллиптического уравнения Справочная информация
- •Вариационная постановка краевой задачи
- •Программное обеспечение
- •Пример решения на пэвм
- •Контрольные задания
- •5. Применение метода конечных элементов для решения параболического уравнения Справочная информация
- •Вариационная постановка начально-краевой задачи
- •Программное обеспечение
- •Пример решения на пэвм
- •Контрольные задания
- •6. Применение метода конечных элементов для решения гиперболического уравнения Справочная информация
- •Вариационная постановка начально-краевой задачи
- •Программное обеспечение
- •Пример решения на пэвм
- •Контрольные задания
- •7. Решение экстремальных задач Справочная информация
- •Программное обеспечение
- •Пример решения на пэвм
- •Контрольные задания
- •Список литературы
Программное обеспечение
Для решения начально-краевой задачи для параболического уравнения в Matlab’е используется описанная ранее среда PDE Tool. Процесс построения решения начально-краевой задачи для параболического уравнения с её помощью показан ниже на примере решения конкретной задачи.
Пример решения на пэвм
Условия задачи. Найти решение начально-краевой задачи для параболического уравнения
,
,
где область S поиска решения показана на рис.2. На границе Г области S заданы следующие граничные условия
,
,
а
Рис.2.
.
Решение. Процесс построения решения начально-краевой задачи для параболического уравнения в среде PDE Tool аналогичен описанному ранее для краевой задачи с эллиптическим уравнением.
Сначала надо задать систему глобальных координат и нарисовать область S поиска решения, как это показано на рис.3.
Рис.3.
Рисование области S, осуществляемое с помощью мыши, характеризуется существенными погрешностями. Однако их можно избежать. Для этого достаточно сделать двойной щелчок по нарисованному объекту. Это выведет на экран окно его свойств (см. рис.4). Так для эллипса можно точно задать координаты его центра (X-center, Y-center) и радиусы полуосей (A-semiaxes, B-semiaxes).
Рис.4.
После этого надо привести параболическое уравнение и его граничные условия к используемому в среде PDE Tool виду
,
,
и задать в пункте меню PDE Specification тип уравнения – parabolic и коэффициенты c=1, a=0, f=0, d=10, а в пункте Specify Boundary Conditions – условия для границы Г1 – тип Neumann, g=1, q=0.5 и условия для границ Г2, Г3 – тип Dirichlet, h=1, r=2+10*t.*exp(-t).
Следующим шагом является формирование конечно-элементной сетки в области поиска решения c помощью меню Mesh. После подбора степени густоты сетки в рабочем окне получается изображение, показанное на рис.5.
Последним шагом в этой процедуре является построение конечно-элементного решения с помощью пункта меню Solve, но с одним отличием. Оно заключается в том, что до решения задачи Коши, к которой сводится начально-краевая задача для параболического уравнения, необходимо установить диапазон изменения времени t и задать
Рис.5.
н
Рис.6.
Для установки необходимой формы представления результатов расчёта, отличной от той, которая задана по умолчанию, следует воспользоваться пунктом меню Plot. Если есть необходимость увидеть процесс видоизменения решения задачи во времени, то в окне Plot Solution надо активизировать строки Height(3-D plot), Animation и Show mesh. В этом случае в графическом окне можно получить мультипликацию трёхмерного изображения решения u(x, y, t), некоторые кадры которой для разных моментов времени показаны на рис.7.
Рис.7.