- •В. Н. Бєловодський Вісім лекцій по чисельних методах
- •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. Завдання
Завдання
1. Знайти методом Пикара три перших наближення рішення заданої (відповідно до варіанта завдання) задачі Коші для звичайного диференціального рівняння першого порядку. Привести всі необхідні розрахунки. Побудувати графіки відповідних функцій.
2. Для тої ж задачі Коші за допомогою ручного прорахунку знайти інтегральну криву методом Рунге-Кутта 4-го порядку на відрізку одиничної довжини (починаючи від X0) із кроком 0,2. Привести всі проміжні розрахунки. Остаточні результати розрахунку звести в таблицю. За результатами побудувати графік.
3. Побудувати Simulink-Модель для заданої задачі Коші за умовами попереднього пункту. Зрівняти результат (представлений у вигляді графіка) з результатом п.2.
Варіанти завдань
Лекція 8. Основи спектрального аналізу
У якій нагадуються основні положення теорії рядів Фур'є, виробляється висновок комплексної форми й приводиться їхній дискретний аналог для рівновіддалених вузлів.
Елементи загальної теорії
Нагадаємо основні положення теорії рядів Фур'є.
Нехай є деяка функція , задана на проміжкуй розглядається нескінченна система функцій
.
Ставиться задача про подання даної функції у вигляді тригонометричного ряду
(1)
Таке подання за певних умов можливо і його коефіцієнти обчислюються по наступних формулах
(2)
Тригонометричний ряд (1) з коефіцієнтами (2) називається поруч Фур'є функції .
Справедлива наступна теорема про розкладність (т. Дирихле):
Якщо кусочно-монотонна функція й має не більш, ніж кінцеве число крапок розриву першого роду, то її ряд Фур'є (1), (2) сходиться до значенняв крапках її безперервності й до середнього арифметичного її однобічних меж у крапках розриву.
Так, наприклад, якщо крапка розриву, то сума ряду Фур'є в цій крапці дорівнює
- (Малюнок 1),
де
, .
Малюнок 1. Пояснення до теореми Дирихле
Т.о. за винятком, бути може, кінцевого числа крапок сума ряду Фур'є (1), (2) дорівнює .
Зауваження 1. У крапках ряд Фур'є сходиться до середньої арифметичної правої й лівої меж функціїв крапкахвідповідно.
Оборотний увага на особливості розкладання в ряд Фур'є парних і непарних функцій.
Якщо парна на відрізку, те, також парна, а, - непарна. Тому коефіцієнти
і розкладання (1) приймає вид
Якщо ж функція непарна, тотакож непарна, а, - парна.
Тому
,
і розкладання (1) приймає вид
Зауваження 2. Вираження в розкладанні (1), уводячи допоміжний кут, можна представити у вигляді, дей ряд у цілому
Тоді доданки називаютьсягармонійними складовими або гармоніками, коефіцієнти -амплітудами гармонік, частотами, -початковими фазами. Іноді гармоніка називаєтьсяосновний, гармоніки , -побічними.
Зауваження 3. Іноді тригонометричні розкладання заданої функції будуються на проміжку . У цьому випадку розглядається система тригонометричних функцій
1, ,
і розкладання має вигляд
(3)
де
(4)
.
Теорема Дирихле для проміжку формулюється відповідно аналогічним образом.
Комплексна форма рядів Фур'є
У додатках досить популярної є комплексна форма рядів Фур'є. Їхнє одержання базується на формулах Эйлера
звідки
, . (5)
Приведемо комплексну форму для розкладання (3). З (4) з використанням (5) маємо
або
,
де
.
Аналогічним образом,
,
або
Тоді гармонійна складова ряду (3)
і ряд, у цілому,
(6)
де
, ,,. (7)
Розкладання (6), (7) і є шуканим.
Оборотний увага, що амплітуда n-ой гармоніки в цьому випадку дорівнює
.