2301

Чисельне вирішення задачі Коші для звичайних диференціальних рівнянь І-го порядку

Лекция

Математика и математический анализ

Основні типи рівнянь інженерної практики. Методи вирішення диференціальних рівнянь. Постановка задач для звичайних диференціальних рівнянь (ЗДР). Метод Ейлера. Модифіковані методи Ейлера та Ейлера-Коші. Метод Рунге-Кутта. Приклад вирішення задачі Коші для ЗДР І-го порядку в середовищі системи Mathcad.

Украинкский

2013-01-06

106.36 KB

110 чел.

Лекція 14. Чисельне розв’язування задачі Коші для звичайних диференціальних рівнянь І-го порядку

План

1. Основні типи рівнянь інженерної практики.

2. Методи розв’язування диференціальних рівнянь.

3. Постановка задач для звичайних диференціальних рівнянь (ЗДР)

4. Методи розв’язування задачі Коші для ЗДР.

4.1. Метод Ейлера.

4.2. Модифіковані методи Ейлера та Ейлера-Коші.

4.3. Метод Рунге-Кутта.

5. Приклад розв’язування задачі Коші для ЗДР І-го порядку в середовищі системи Mathcad.

1. Основні типи рівнянь інженерної практики

Всі рівняння інженерної практики можна розділити на два класи:

  1.  числові рівняння – описують стаціонарні статичні процеси і їх розв’язком є числа;
  2.  функціональні рівняння – описують нестаціонарні динамічні процеси і їх розв’язком є функції.

Функціональні рівняння діляться на два класи:

  1.  диференціальні рівняння (ДР) – функція під знаком диференціала;
  2.  інтегральні рівняння – функція під знаком інтеграла.

Диференціальні рівняння можна розділити на два класи:

  1.  звичайні диференціальні рівняння (ЗДР) – функція залежить від однієї змінної;
  2.  диференціальні рівняння в частинних похідних – функція залежить від двох або більше змінних.

Порядком ДР називається найвищий порядок похідної або диференціала, який входить в це рівняння. Розв’язком ДР називається функція, підстановка якої в рівняння перетворює його в тотожність.

Існує три типи розв’язків ЗДР:

  1.  загальний розв’язок – сімейство розв’язків ЗДР, яке лежить від довільних постійних, кількість яких співпадає з порядком цього рівняння;
  2.  частковий розв’язок – може бути отриманий із загального розв’язку при певних числових значеннях довільних постійних, які в нього входять;
  3.  особливий розв’язок – розв’язок, який у всіх своїх точках не задовільняє умову єдиності, тобто в будь-якому околі особливого розв’язку існує хоча б дві інтегральні криві, які проходять через цю точку.

2. Методи розв’язування диференціальних рівнянь

Серед методів розв’язування ДР можна виділити такі:

  1.  класичні (точні) методи – дають можливість отримати розв’язки у вигляді формул шляхом аналітичних перетворень і табличного інтегрування елементарних функцій, але ці методи не завжди можна застосувати, а у випадку застосування часто дають складні та громіздкі розв’язки;
  2.  наближені методи:

а) графічні методи – наближений розв’язок у вигляді графіків;

б) аналітичні методи – базуються на спрощенні заданих рівнянь так, щоб більш просте рівняння можна було б розв’язати класичними методами, але при цьому виникає питання про достовірність отриманого розв’язку;

  1.  числові методи – наближені розв’язки дають у вигляді таблиць.

Числові методи, на даний час, отримали дуже широке застосування у зв’язку з інтенсивним розвитком обчислювальної техніки. Ідея цих методів базується на заміні диференціальних рівнянь і додаткових умов системою алгебраїчних рівнянь, розв’язок якої отримуємо у табличному вигляді.

3. Постановка задач для ЗДР 

Для всіх фізичних задач характерним є наявність межі Γ області G, в якій вивчається той чи інший процес. Ці межі можуть бути скінченними або нескінченними. Оскільки математична модель повинна адекватно описувати певне фізичне явище в даному виділеному середовищі, то вона включає в себе не тільки диференціальне рівняння або систему диференціальних рівнянь, але й додаткові або крайові умови, які задаються у вигляді значень шуканої функції або її похідних для деяких значень незалежних змінних, тобто в окремих точках, або у вигляді залежностей шуканої функції та її похідних на деякій області.

Додаткові (крайові) умови – це сукупність граничних і початкових умов. Граничні умови задають режим фізичного процесу на межі Γ області G, а початкові умови накладають обмеження на функцію u та її похідні по часу до (n-1)-го порядку в деякий початковий момент t=t0.

Для знаходження часткового розв’язку диференціального рівняння необхідно задати додаткові (крайові) умови, кількість яких повинна бути не меншою, ніж порядок рівняння.

В залежності від способу задання додаткових (крайових) умов для відшукання часткового розв’язку диференціального рівняння розглядають задачу Коші та крайову задачу.

Якщо ці умови відносяться до однієї точки, то задача називається задачею Коші, додаткові умови – початковими умовами, а сама точка – початковою. Якщо ж ці умови відносяться до більш, ніж однієї точки, то така задача називається крайовою задачею, а відповідні умови – граничними умовами.

  1.  Задача Коші. Знайти функцію y=y(x), яка задовольняє диференціальне рівняння   (1)

і початкові умови:  (2)

При постановці задачі Коші для ЗДР n-го порядку задається значення функції y(x) в деякій точці, а також значення похідних цієї функції в цій точці до (n-1)-го порядку включно.

  1.  Крайова задача. Знайти на деякому проміжку [a, b] функцію y=y(x), яка задовольняє всередині відрізка диференціальному рівнянню    (3)

а на кінцях відрізка – граничним умовам:

  (4)

4. Методи розв’язування задачі Коші для ЗДР

Нехай на проміжку [x0, b]задана задача Коші для ЗДР І-го порядку

Необхідно знайти функцію y=y(x), яка є розв’язком рівняння (5) і задовольняє початкову умову (6).

На проміжку [x0, b] з допомогою точок побудуємо рівномірну різницеву сітку з (k+1) вузлом, відстань між якими дорівнює h:    (7)

Задача полягає у знаходженні значень у вузлах сітки (7) шуканої функції y(x), яка є розв’язком задачі Коші (5)-(6)

y(x0)= y0, y(x1)= y1, y(x2)= y2, …, y(xk-1)= yk-1, y(xk)= yk. (8)

Одним з найбільш простих та універсальних числових методів розв’язування диференціальних рівнянь є метод скінченних різниць, який базується на розкладі шуканої функції y(x) в ряд Тейлора в h-околі точки

(9)

4.1. Метод Ейлера

З розкладу (9) беремо перші два члени

(10)

що, з використанням (8), можна записати у вигляді

(11)

У відповідності з (5) це рівносильно

(12)

або (13)

де – залишковий член.

З (13) отримаємо різницеву схему, яка апроксимує в точках (7) розв’язок задачі Коші та має вигляд рекурентних формул

 (14)

де (15)

Коли розв’яжемо систему числових рівнянь, то отримаємо значення

(16)

що відповідають шуканим значенням (8).

З (12) випливає, що метод Ейлера має перший порядок точності і з віддаленням від початкової точки похибка накопичується. На практиці частіше використовують різні модифікації методу Ейлера, які дають більш високу точність.

4.2. Модифіковані методи Ейлера та Ейлера-Коші

Якщо з розкладу (9) взяти перші три члени, то отримаємо розрахункові формули для модифікованих методів Ейлера та Ейлера-Коші, які мають другий порядок точності , що при малих h забезпечує більш високу точність, ніж метод Ейлера. Ці методи можна представити такими рекурентними формулами:

модифікований метод Ейлера

 (17)

модифікований метод Ейлера -Коші

 (18)

4.3. Метод Рунге-Кутта

Якщо з розкладу (9) взяти перші п’ять членів, то отримаємо розрахункові формули для методу Рунге-Кутта, який має четвертий порядок точності , що забезпечує більш високу точність, ніж розглянуті раніше, які фактично є методами Рунге-Кутта відповідно І-го і ІІ-го порядків.

Цей метод можна представити такими рекурентними формулами:

(19)

5. Приклад розв’язування задачі Коші для ЗДР І-го порядку в середовищі системи Mathcad

Диференціальні рівняння першого порядку можуть, за означенням, містити , крім шуканої функції , тільки її першу похідну . В більшості випадків диференціальне рівняння можна записати в стандартній формі

і тільки з такою формою може працювати обчислювальний процесор MathCad.

Крім диференціального рівняння , потрібно задати початкову умову – значення функції в деякій точці . Необхідно визначити функцію на проміжку від до .

Для чисельного інтегрування звичайного диференціального рівняння (ЗДР) у користувача системи MathCad є вибір – або використати обчислювальний блок Given/Odesolve, або одну із вбудованих функцій Rkfixed, Rkadapt, Bulstoer. Розглянемо послідовно обидва варіанти розв’язування.

Обчислювальний блок Given/Odesolve

Обчислювальний блок для розв’язування одного диференціального рівняння, який реалізує чисельний метод Рунге-Кутта, складається із трьох частин:

  1.  Given ключове слово;
  2.  звичайне диференціальне рівняння та початкова умова, записані за допомогою логічних операторів, причому початкова умова повинна бути записана у формі ;
  3.  Odesolve() – вбудована функція для розв’язання звичайного диференціального рівняння відносно змінної на проміжку .

Приклад 1. Розвязати задачу Коші для ЗДР першого порядку   на проміжку .

Розв’язання:

Вбудована функція Rkfixed

Так як розвязування за допомогою вбудованої функції Rkfixed мало чим відрізняється від попередньго способу ( за допомогою обчислювального блоку), то приведемо приклад його використання. Звернемо лише Вашу увагу на необхідність явного задання кількості точок інтегрування ЗДР M=100, а також на отримання результату, на відміну від обчислювального блоку, не у вигляді функції, а у вигляді матриці розмірності . Ця матриця складається із двох стовпців: в одному знаходяться значення аргументу , а в другому відповідні значення шуканої функції .

Приклад 2. Розвязати задачу Коші для ЗДР першого порядку   на проміжку .

Розв'язання:


 

А также другие работы, которые могут Вас заинтересовать

49318. Расчет АЧХ, ФЧХ и переходной характеристики полосового фильтра 129.03 KB
  Техническое задание 3 АЧХ и ФЧХ Переходная характеристика Заключение Техническое задание Построить АЧХ ФЧХ и переходную характеристику цепи. Схема расчета АЧХ и ФЧХ Операторное передаточное сопротивление: где G =1 R Заменим p на jw АЧХ: ФЧХ: Подставим значения Найдем собственную частоту контура Определим характеристическое сопротивление Рассчитаем добротность Полоса пропускания Таблица 1. ФЧХ 0.
49319. Проект очистных работ для лавы 111.6 KB
  Поэтому проверим возможность применения механизированного комбайнового комплекса в условиях заданной лавы: При выборе механизированного комплекса следует учесть что вынимаемая мощность в данном случае будет складываться из мощности угля и мощности ложной кровли: где мощность пласта угля м; мощность ложной кровли м. В условиях данной лавы возможно применение следующих комбайновых комплексов: 1КМ88; 2КМ87УМА; 1КМ87УМВ; 2КМ87УМВ. при работе по простиранию падению восстанию 15 10 10 20 10 10 20 10 10 20 10 10 Устойчивость...
49320. Разработка программы на алгоритмическом языке программирования Си 195.01 KB
  Программа на языке СИ. Необходимо ответить на вопросы: “Что заданКакой должен быть получен результат†“Как получить результат†Задача моего варианта курсовой работы заключается в проверке истинности высказывания: Все цифры данного числа различныЭто значит что мне нужно используя знания полученные на курсах информатики а так же при необходимости используя дополнительную литературу составить программу на языке си которая могла бы определить – все ли цифры различны в заданном трехзначном...
49321. Пароходик догоняет большой пароход 49.92 KB
  Целью работы также является проведение сравнительного анализа языков программирования. Сравнить языки потребуется как в общем, так и применительно поставленной задаче. Одной из подзадач является создание описание алгоритма программы и составление математической модели. Пароходы будут состоять из линий и окружностей. Для написания выбран язык программирования, среда Delphi 7.
49322. Моделирование логических игровых программ средствами Delphi 747.5 KB
  Объект исследования – применение среды программирования Borland Delphi с целью изучения возможности отображения графической информации, построения фракталов.
49323. СИНТЕЗ СХЕМЫ ГЕНЕРАТОРА ЧИСЕЛ СО СТРУКТУРОЙ АВТОМАТА МУРА 6.4 MB
  Синтезировать схему генератора чисел 0-15-2-1-5-6-10-9 0-13-1-7-5-2-11-6-12 со структурой автомата Мура и Мили на RS и D триггерах в базисе ИЛИ-НЕ, определить схему с минимальным количеством входов, проверить правильность синтеза в MicroCap.
49324. ПРОЕКТИРОВАНИЕ И РЕАЛИЗАЦИЯ МОДУЛЯ ПРОХОЖДЕНИЯ ТЕСТИРОВАНИЯ СИСТЕМЫ ПРОВЕРКИ УЧАЩИХСЯ НА ЗНАНИЕ АЛГОРИТМОВ 1.58 MB
  Обзор систем тестирования Приложение Визуальная студия тестирования Система тестирования INDIGO
49325. Методы локализации неисправностей на аппаратуре СВ и РМ 1.63 MB
  После записи числа Х в ячейку памяти У при наличии свободных оперативных регистров контролируем содержимое ячейки ЗУ: на информационном поле оперативного пульта управления набираем адрес У; нажимаем клавиши НУ ЗАП ССП ПУСК; на поле индикации при переключателе режимов установленном на значении ОР число Х не отображается. Вычислительное устройство ВчУ является основным операционным устройством СВ предназначенным для обработки цифровой и логической информации реагирования на сигналы прерывания внешних устройств и управления...