
- •Тема 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).
Оцінка похибки
Для даної задачі:
.
Отже,
,