Метод прогону. Розділимо область інтегрування [a,b] на досить велике число рівних частин точками xi = x0 + ih, i = 1 N

Розділимо область інтегрування [a,b] на досить велике число рівних частин точками x_i = x₀ + ih, i = 1 … N... Перетворимо систему (1), (2) до виду

y_i+1 + m_iy_i + n_iy_i-1 = h²f(x_i); i = 1 … N-1, (3)

де r_i = 1/(1 + h/2p_i); p_i = p(x_i); q_i = q(x_i);

m_i = r_i(h²q_i – 2); n_i = r_i(1 – h/2p_i); ¦ˆ_i = h²r_if(x_i). (4)

Розв’язавши рівняння (3) відносно у _i, одержуємо

y _i = c_i(d_i - y_i+1), i = N, N-1, …, 1; (5)

де c_i = 1/(m_i - n_ic_i-1), d_i = ¦ˆ_i - n_ic_i-1d_i-1, i = 1 … N... (6)

c₀ = a₁/(a₀h - a₁), d₀ = h/a_1.

По формулах (6) обчислюються коефіцієнти c_i, d_i – прямий хід прогону. По формулі (5) обчислюємо y_i, i = N, N-1, …, 1 – зворотний хід прогону. При цьому з (3) і (4) знаходимо

y₊₁ = (Bh + b₁c_N-1d)/(b₀h + b₁(c+1)).

Проекційні методи (на прикладі методу Гальоркіна).

Сутність проекційних методів обчислювальної математики полягає в представленні розв’язку задачі множиною проекцій (відліків) у визначеній системі координатних функцій.

У традиційному методі, запропонованому Б. П. Гальоркіним і розвинутому у роботах М. В. Келдиша, наближений розв’язок у (х) відшукується у вигляді

y_n = a₁u_1(x) + a₂u₂ + _... + a_nu_n(x), (7)

де j_1(x), j_2,,,,, j_n(x), – система базисних функцій, що задовольняють вихідним граничним умовам; a_1, a_2,..., a_n – невідомі постійні коефіцієнти.

При підстановці (7) у (1) одержуємо функцію

R(x; a₁; a₂;...; a_n) = L[y_n(x)] – f(x),

яку називають нев'язкою розв’язку. Очевидно, для точного розв’язку задачі

R(x; a₁; a₂;...; a_n) = 0.

Коефіцієнти а_i визначимо з умови ортогональності нев'язки першим n функціям деякої системи функцій {j _i, i = 0... n }:

(8)

Такий метод розв’язку задачі називається методом моментів.

Якщо при цьому j_i=u_i, то виходить метод Гальоркіна. Система рівнянь (8) являє собою систему лінійних рівнянь щодо вектора з матрицею

.

Розв’язок цієї системи є каркасом наближеного розв’язку крайової задачі.

Система (8) може бути представлена у вигляді

. (9)

Коефіцієнти c_ki, d_i, обчислюються по формулах

Для рівняння (1) L(y) = y¢¢ - p(x)y¢ - q(x)y.

Тому

c_ik = c_ki

Функцію u₀(x) можна вибирати довільно, але так щоб u₀(a)=A, u₀(b)=B. Нехай u₀(x) = a + bx. Тоді з (2) одержуємо:

Інші функції u_i(x) можна обчислити по кожному з правил:

u_i(x) = (x – a)ⁱ(x – b);

u_i(x) = (x – a)(x – b)ⁱ.

ОБЧИСЛЮВАЛЬНІ СХЕМИ

Дата добавления: 2014-01-07