Для БПМ13 / Численные методы лаб / 6.Решение краевой задачи

Решение краевой задачи для линейного

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

Пусть на отрезке требуется найти решение дифференциального уравнения:

, (1)

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

;	; (2)

Численное решение задачи состоит в нахождении приближенных значений искомого решения в точках . Для этого разобьем отрезок на равных частей с шагом . Полагая и вводя обозначения , , для внутренних точек отрезка , вместо дифференциального уравнения (1)–(2) получаем систему конечноразностных уравнений:

После соответствующих преобразований будем иметь

, , (3)

где

.

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

Решая уравнение (3) относительно , будем иметь

.

Предположим, что из этого уравнения исключена неизвестная . Тогда это уравнение примет вид

, (4)

где – некоторые коэффициенты.

Отсюда . Подставляя это выражение в (3), получим и, следовательно,

. (5)

Сравнивая формулы (4) и (5), получим для определения рекуррентные формулы:

.

Определим :

.

Из формулы (4) при имеем

. (6)

Поэтому

, . (7)

На основании формул (6) и (7) последовательно определяются коэффициенты до включительно (прямой ход). Обратный ход начинается с определения . Решая систему

,

получим

и по формуле (4) последовательно находим .

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

Отсюда .

Пример. Методом прогонки решить краевую задачу:

.

Решение. Пусть .

;

;

; ;

.

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

Таблица

	0	1	2	3	4	5
	0	-0,498	-0,662	-0,878	-0,890	-0,900
		0,001	0,002	0,004	0,008	0,012
	0	0,1	0,2	0,3	0,4	0,5
	0	-0,025	-0,049	-0,072	-0,078	-0,081
	0	-0,015	-0,029	-0,041	-0,050	-0,057

	6	7	8	9	10
	-0,908	-0,915	-0,921	-0,926
	0,16	0,022	0,028	0,035
	0,6	0,7	0,8	0,9	1
	-0,078	-0,070	-0,055	-0,032	0
	-0,058	-0,054	-0,044	-0,026	0

Задание. Решить краевую задачу методом прогонки.

№	Дифференциальное уравнение	Краевые условия
1				10
2				20
3				30
4				40
5				50
6				10
7				20
8				30
9				40
10				50