- •Тема 1 рівняння математичної фізики 6
- •Тема 2 зведення рівнянь другого порядку до канонічного вигляду за допомогою заміни змінних 28
- •Тема 3 метод фур'є 55
- •Тема 4 метод сіток для рівняння параболічного типу 111
- •Тема 5 спеціальні функції математичної фізики 121
- •Передмова
- •Тема 1 рівняння математичної фізики
- •1.1 Рівняння малих поперечних коливань струни
- •1.2 Рівняння малих поздовжніх коливань стержня
- •1.3 Рівняння малих поперечних коливань мембрани
- •1.4 Телеграфне рівняння
- •1.5 Рівняння теплопровідності
- •1.6 Рівняння поширення тепла в стержні
- •1.7 Основні рівняння математичної фізики
- •Питання для самоперевірки
- •Завдання для самостійної роботи
- •Тема 2 зведення рівнянь другого порядку до канонічного вигляду за допомогою заміни змінних
- •Розв’язування
- •Розв’язування
- •2.1 Рівняння гіперболічного типу
- •Розв’язування
- •2.2 Рівняння еліптичного типу
- •Розв’язування
- •2.3 Рівняння параболічного типу
- •Розв’язування
- •Розв’язування
- •Розв’язування
- •Розв’язування
- •Розв’язування
- •Розв’язування
- •2.4 Розв’язування задачі Коші для рівняння коливання струни методом характеристик (формула д’Аламбера)
- •Розв’язування
- •Питання для самоперевірки
- •Завдання для самостійної роботи
- •Тема 3 метод фур'є
- •3.1 Розв’язання методом Фур’є першої крайової задачі для рівняння малих поперечних коливань струни
- •3.2 Розв’язання методом Фур’є першої крайової задачі для рівняння теплопровідності
- •3.3 Розв’язання методом Фур’є першої крайової задачі для рівняння поширення тепла у нескінченному стержні
- •3.4 Приклади розв’язання задачі Коші для рівняння теплопровідності
- •Розв’язування
- •Розв’язування
- •Розв’язування
- •Розв’язування
- •Розв’язування
- •Розв’язування
- •Питання для самоперевірки
- •Завдання для самостійної роботи
- •Тема 4 метод сіток для рівняння параболічного типу
- •Розв’язування
- •Питання для самоперевірки
- •Завдання для самостійної роботи
- •Тема 5 спеціальні функції математичної фізики
- •5.1 Інтеграл Ейлера першого роду
- •5.2 Інтеграл Ейлера другого роду
- •5.3 Функція Бесселя
- •5.4 Рекурентні формули для функції Бесселя
- •5.5 Інтегральне представлення Пуассона функції Бесселя та його використання
- •5.6 Сферичні функції. Поліноми Лежандра
- •5.7 Виробнича функція для поліномів Лежандра
- •5.8 Рекурентні формули для поліномів Лежандра
- •Питання для самоперевірки
- •Завдання для самостійної роботи
- •Відповіді
- •Література
- •Предметний покажчик
Розв’язування
Виберемо для аргументу х крок h = 0,1. Оскільки , отримаємо для аргументу t такий крок :
.
Помічаємо, що початкова і крайові умови задачі симетричні відносно прямої . Тому і розв’язок буде також симетричним відносно цієї прямої.
Спочатку обчислимо значення функції на нульовому шарі , тобто для та ; 0,1; 0,2; 0,3; 0,4; 0,5.
Результати обчислень зручно подати у вигляді таблиці.
Початковий рядок цієї таблиці заповнюється на основі заданої початкової умови
, .
Матимемо такі значення функції :
…
У перший та у останній стовпці записуємо дані крайових умов
,
тобто,
,
.
Решта рядків таблиці послідовно заповнюється на основі розрахункової формули (5.9):
.
Таблиця 4.1 – Розв’язання рівняння теплопровідності методом сіток
j |
і |
0 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
|
0 |
0,1 |
0,2 |
0,3 |
0,4 |
0,5 |
0,6 |
0,7 |
0,8 |
0,9 |
1,0 |
|
0
|
0 |
0,0000 |
0,3090 |
0,5878 |
0,8090 |
0,9511 |
1,0000 |
0,9511 |
0,8090 |
0,5878 |
0,3090 |
0 |
1
|
0,005 |
0 |
0,2939 |
0,5590 |
0,7695 |
0,9045 |
0,9511 |
0,9045 |
0,7695 |
0,5590 |
0,2939 |
0 |
2
|
0,010 |
0 |
0,2795 |
0,5317 |
0,7318 |
0,8603 |
0,9045 |
0,8603 |
0,7318 |
0,5317 |
0,2795 |
0 |
3
|
0,015 |
0 |
0,2659 |
0,5057 |
0,6960 |
0,8182 |
0,8603 |
0,8182 |
0,6960 |
0,5057 |
0,2659 |
0 |
4
|
0,020 |
0 |
0,2529 |
0,4610 |
0,6620 |
0,7782 |
0,8182 |
0,7782 |
0,6620 |
0,4610 |
0,2529 |
0 |
5
|
0,025 |
0 |
0,2405 |
0,4575 |
0,6296 |
0,7401 |
0,7782 |
0,7401 |
0,6296 |
0,4575 |
0,2405 |
0 |
При цьому, звичайно, доцільно враховувати симетрію шуканої функції.
Заповнимо перший рядок таблиці, обчислюючи значення шуканої функції за формулою:
,
використовуючи отримані значення функції з початкового рядка і крайові умови для .
Таким чином, матимемо:
;
;
;
;
.
Враховуючи те, що шукана функція симетрична відносно прямої , то наступні значення можна записати, не обчислюючи їх:
;
;
;
;
.
Записуємо отримані значення у перший рядок таблиці. Після цього переходимо до обчислення значень функції на другому шарі за формулою:
,
використовуючи отримані значення функції з першого рядка (з першого шару) і відповідні крайові умови.
Аналогічно заповнюємо решту таблиці, обчислюючи значення шуканої функції за формулою (4.9).
Оцінка похибки
Для даної задачі:
.
Отже, ,