Численные методы, ч.1 (стр. 19 )

Пятидиагональная матрица имеет следующее строение:

Вектор правой части СЛАУ задается в данной задаче по закону:

,

где - число, определяющее в индексах вектора решения начало последнего слоя внутренних узлов по оси , на которых учитывается заданное граничное условие. Заданная функция из уравнения, в данной задаче являющаяся константой, , входит в правую часть СЛАУ в виде слагаемого . Значения вектора правой части в данной задаче:

На рис.3 показано распределение функции решения аналогичной краевой задачи в двумерной области при порядке СЛАУ . Решение получено комбинированным методом Зейделя-ОСП при оптимальном параметре за итераций с относительной точностью решения в . Обычный метод Зейделя сходится здесь лишь за итераций и сопоставим по времени решения с прямым методом. Еще большее число требуемых итераций показывают в данной задаче метод ОСП с матрицей (1.2) - .

Отметим, что матрица задачи при и заданном постоянна и не зависит от краевых условий и источников , которые входят в правую часть СЛАУ. Соответственно задачи с различными краевыми условиями и источниками могут решаться с тем же самым оптимальным параметром, найденным один раз для данной сетки.

рис.3

7.5.  Лабораторные задания к теме «Численное решение уравнений в частных производных»

Лабораторные работы по теме могут быть выполнены с помощью математических пакетов программ Mathcad или Matlab. В результате работы должна быть представлена искомая сеточная функция в виде матрицы значений в узлах сетки либо в виде послойного по времени распределения значений сеточных векторов. Следует также привести графическое представление результатов.

7.5.1. Гиперболические уравнения

Варианты заданий для одномерного волнового уравнения с граничными и начальными условиями (см. 7.2).

№ п/п									Метод
1.1	1	1	0	0			2	0.1, 0.05	1
1.2	1	2	0	0			1	0.1, 0.05	1
1.3		1	0	0			2	0.1, 0.05	1
1.4	1	2	0	0			1	0.1, 0.05	1
1.5		1	0	0			2	0.1, 0.05	1
1.6	1	2	0	0			1	0.1, 0.05	1
1.7		1	0	0			2	0.1, 0.05	1
1.8	1	2	0	0			1	0.1, 0.05	1
1.9		1	0	0			2	0.1, 0.05	1
1.10	1	2	0	0			1	0.1, 0.05	1

7.5.2. Параболические уравнения

Варианты заданий для одномерного уравнения теплопроводности с граничными и начальными условиями (см. 7.3).

№ п/п									Метод
2.1	1	0.1	0	0			1	0.01	1
2.2	1	1	0	0			1	0.01	2
2.3	1	0.1	0	0			1	0.01	1
2.4	1	2	0	0			1	0.01	2
2.5		0.1	0	0			1	0.01	1
2.6		1	0	0			1	0.01	2
2.7		0.1	0	0			1	0.01	1
2.8		2	0	0			1	0.01	2
2.9	1	0.1	0	0			1	0.01	1
2.10	1	1	0	0			1	0.01	2
2.11		0.1	0	0			1	0.01	1
2.12		2	0	0			1	0.01	2

Методы: 1- явный; 2 - неявный, сведением к СЛАУ и последующим решением стандартным методом. Сравнить со строгим решением.