§5. Метод скінчених різниць для граничної задачі, для лінійного диференціального рівняння другого порядку з змінними коефіцієнтами.

Нехай задане диференціальне рівняння:

(1)

і крайові умови: (2) або

(2‘)

(2 )

Потрібно знати розв‘язок диференціального рівняння у вигляді таблиці значень на .

Розіб‘ємо на частини точками:

, ,

Другу похідну в рівнянні (1) замінимо різницевим відношенням:

Тоді отримаємо систему:

(3)

В залежності від крайових умов до системи (3) будуть додані рівняння:

(4)

(4‘)

Враховуючи формули чисельного диференціювання функцій інтерполювання многочленами Ньютона.

або

(4‘‘)

Отримаємо систему рівнянь з невідомими. Дана система буде мати єдиний розв‘язок якщо відповідна однорідна система (випадок коли ) матиме лише тривіальний розв‘язок.

Введемо позначення:

Лема1. Нехай дана довільна система з точок:

Якщо , для довільного , то найбільшим додатнім числом серед може бути або .

Нехай - найбільше додатне число, .

Тоді існують числа , . Маємо:

, тоді випливає

, тобто

, що неможливо.

Лема2. Якщо задана система точок і виконується умова , , то найменшим від‘ємним числом серед може бути лише або .

Доведення аналогічне.

Доведемо, що система (3) з крайовими умовами (4), (4‘), (4‘‘) має лише тривіальний розв‘язок у випадку коли .

1. Нехай маємо (3) і умови (4). Якщо існує ненульовий розв‘язок, то серед чисел є найбільше додатне або найменше від‘ємне, а це суперечить Лемі 1 або 2.

2. Нехай маємо (3) і умови (4‘), маємо:

Якщо існує ненульовий розв‘язок, то серед чисел є найбільше додатне або найменше від‘ємне. Знову суперечність з Лемою 1 або 2.

3. Система (3) і умови (4‘‘): і хоча б одне з них не дорівнює 0. Тоді в (3) підставимо :

Визначивши з ( ) і підставимо в ( ).

(5)

Аналогічне підставивши в систему (3) маємо:

З останньої системи можна визначити:

(6)

Розглянемо рівності (5) і (6). Нехай знову існує деякий нетривіальний розв‘язок системи, тоді найбільше додатне або найменше від‘ємне число можуть бути або .

Нехай - найбільше додатне число, тоді крок виберемо настільки малим щоб:

, тоді з (5):

, що неможливо з Лемою.

Щоб отримати таблицю розв‘язків диференціального рівняння при всіх розглянутих крайових умовах слід розв’язати одним з відомих способів систему (3) і (4) або (4‘) або (4‘‘).

Якщо позначити наближений а - точний розв‘язки, тоді величина: вказує похибку в вузлі . Можна довести, що похибка оцінюється таким співвідношенням:

(7)

Заваження1: Розглянутий метод можна застосувати також для рівнянь виду:

(8)

з умовами:

(9)

Якщо другу похідну замінити так само як в попередніх випадках, а першу похідну так:

, отримаємо:

(10)

Зауваження2: Похідні замінюють таким співвідношенням:

(11)

В усіх розглянутих випадках похибка оцінюється нерівністю (7).

§6. Метод прогонки розв’язування крайових задач для звичайних диференціальних рівнянь другого порядку.

Нехай задане диференціальне рівняння:

(1)

(2)

де:

Де функції неперервна на

і розіб‘ємо точками:

Домовимось, що , ,

Замінивши похідні в (1) різницевими співвідношеннями, маємо:

(3)

(4)

Тоді: (5)

З крайових умов:

(6)

Розв‘яжемо (5) відносно змінної

(7)

Припустимо, що за допомогою повної системи рівнянь (5) вдалося з останнього рівняння виключити змінну , тоді з (7):

(8)

де: і - невідомі коефіцієнти.

Знайдемо формули для їх обчислення.

При з (7) маємо:

і з першого рівняння системи (6):

Після арифметичних перетворень маємо:

З другого боку з (8):

Таким чином отримаємо:

(9)

При з (8):

Підставивши знайдене значення в (7) маємо:

З (8) маємо:

(10)

З (9) і (10) можна послідовно обчислити і де - на цьому завершується прямий хід.

Підставивши в (8) замість , маємо:

(11)

Використовуючи формули (8) і першу крайову умову з (6) ми можемо знайти

(12)

52