- •В. Н. Бєловодський Вісім лекцій по чисельних методах
- •7.080407 «Комп'ютерний еколого-економічний моніторинг»)
- •Лекція 1.Елементи теорії погрішностей
- •.Типи й джерела погрішностей
- •1.2.Абсолютні й відносні погрішності наближених чисел
- •1.3. Погрішності виконання арифметичних операцій
- •.Погрішність обчислення функції
- •1.5.Запис наближених чисел
- •1.6.Правила дій над наближеними числами
- •1.7.Погрішності при машинному поданні чисел
- •1.8. Варіанти індивідуальних завдань
- •Лекція 2. Інтерполяція функцій
- •2.1. Постановка задачі
- •2.2. Алгебраїчна інтерполяція, існування й одиничність інтерполяційного багаточлена
- •2.3. Інтерполяційний багаточлен Лагранжа
- •2.4. Кінцеві й розділені різниці
- •2.5. Інтерполяційний багаточлен Ньютона
- •2.6. Порівняльний аналіз інтерполяційних багаточленів
- •2.7. Погрішності інтерполяційних формул
- •2.8. Інтерполяційні формули для рівновіддалених вузлів
- •2.9. Сплайн - інтерполяція
- •2.10. Варіанти індивідуальних завдань
- •Лекція 3. Методи рішення систем лінійних алгебраїчних рівнянь
- •3.1. Попередні зауваження
- •3.2. Точні методи рішення
- •3.3. Наближені методи рішення
- •3.4. Збіжність і погрішність наближених методів
- •3.6. Варіанти індивідуальних завдань
- •Лекція 4. Вирішення нелінійних рівнянь
- •4.1. Попередні зауваження
- •4.2. Методи, засновані на алгебраїчній інтерполяції
- •4.3. Метод послідовних наближень
- •4.4. Завдання
- •Лекція 5. Рішення систем нелінійних рівнянь
- •Метод ітерацій
- •Метод Ньютона
- •5.3. Порівняльний аналіз методів
- •5.4. Завдання
- •Вимоги до програми
- •Вимоги до звіту
- •6.1. Вступні зауваження
- •6.2. Формули Ньютона-Котеса
- •6.3. Найпростіші квадратурні правила
- •6.4. Погрішності квадратурних формул
- •6.5. Поняття про методи Монте-Карло
- •6.6. Завдання
- •Варіанти завдань
- •Лекція 7. Методи рішення звичайних диференціальних рівнянь
- •Вступні зауваження
- •Аналітичні методи
- •Чисельні методи рішення. Правило Рунге
- •Завдання
- •Лекція 8. Основи спектрального аналізу
- •Елементи загальної теорії
- •Дискретна форма рядів Фур'є
- •8.4. Завдання
Варіанти завдань
Лекція 7. Методи рішення звичайних диференціальних рівнянь
У якій розглядаються найпростіші аналітичні послідовних наближень і рядів Тейлора) і чисельні (Эйлера, трапецій і Рунге-Кутта) методи рішення задачі Коші звичайних диференціальних рівнянь.
Вступні зауваження
Нижче розглядаються деякі, що одержали поширення аналітичні й чисельні методи рішення задачі Коші для звичайного диференціального рівняння першого порядку
(1)
. (2)
Умова (2) називається початковим.
Достатні умови можливості розв'язання таких задач формулюються теоремою існування й одиничності, один з найбільш простих варіантів якої, полягає в наступному.
Нехай у крапці (x0, y0) і деякої її околиці функція f(x, y) і її похідна безперервні. Тоді існує єдина функціяy=y(x) безупинно дифференцируемая в крапці x0 і деякої її околиці задовольняючому диференціальному рівнянню (1) і початковій умові (2).
Надалі ці умови передбачаються виконаними.
Аналітичними називають методи, що дозволяють одержати вираження, що описують наближені значення шуканої функції в будь-якій крапці заданого відрізка [a, b]. Чисельними називаються методи, що дозволяють одержати наближені значення лише в окремих крапках заданого відрізка.
Аналітичні методи
Розглянемо два з них.
Метод послідовних наближень. Представимо в рівнянні (1) у вигляді відносини диференціалів, тоді
Далі, інтегруючи обидві частини отриманого співвідношення на проміжку [x0, x] одержуємо
,
або
(3)
Співвідношення (3) являє собою інтегральне рівняння, еквівалентне задачі (1), (2). На його основі будується наступний обчислювальний процес
(4)
де , що і називаєтьсяметодом послідовних наближень. За певних умов послідовність функції сходиться до точного рішення задачі (1), (2). А саме,
Нехай у прямокутнику виконані умови теореми існування й одиничності й. Тоді в проміжку, депослідовність (4) сходиться до точного рішення. Причому справедливо наступну оцінку
де .
Для ілюстрації цього методу розглянемо наступний
Приклад. Знайти рішення задачі Коші
, .
Знайти два перших наближення до рішення, оцінити погрішність.
Рішення. Візьмемо як область квадрат [-1, 1; -1, 1]. Тут,,. Тоді
т.е. М=2,
т.е. N=2,
.
Погрішність n-го наближення
,
звідси погрішність другого,- .
Знайдемо наближення Думаючи, що, з (4) маємо
,
т.е. .
Тоді
т.е.
Зауваження. У тому випадку, коли дана погрішність наближеного рішення число ітерацій, тобто послідовних наближень, можна знайти, зажадавши
Так, наприклад, у розглянутому прикладі
,
звідки й значеннявстановлюється послідовним перебором.
Метод рядів Тейлора. У цьому випадку рішення задачі (1), (2) шукається у вигляді ряду
.
Значення похідних, необхідних для побудови рішення перебувають шляхом послідовного диференціювання рівняння (1). Так, безпосередньо з нього треба
Далі,
,
тоді
і т.д.
Очевидно, що дана процедура дозволяє одержати рішення з як завгодно високою точністю лише в тому випадку, коли функція є нескінченно дифференцируемой у крапціВ противному випадку даний метод може бути взагалі не застосуємо. Див., наприклад, задачу
де вже не існує.