2301

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

Лекция

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

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

Украинкский

2013-01-06

106.36 KB

109 чел.

Лекція 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. Розвязати задачу Коші для ЗДР першого порядку   на проміжку .

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


 

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

39861. Модернизация конструкции передвижного стола пресса ковочного гидравлического пресса усилием 150 МН 3.45 MB
  В контексте всего вышесказанного в данном дипломном проекте рассмотрены вопросы производительности гидравлического пресса и на основании исследований были произведены изменения в цилиндре стола с целью увеличения производительности также проведены мероприятия по охране труда и гражданской обороне. В качестве аккумуляторов в кривошипных и винтовых прессах применяют маховики в гидропрессахгидроаккумуляторы и в гидровинтовых прессах маховики и гидроаккумуляторы. Однако особый интерес представляют пресса большой мощности к ним относятся...
39862. Экономическое обоснование создания нового предприятия (на примере ООО «Blue bird») 847.5 KB
  В дипломном проекте дана характеристика разрабатываемой станции технического обслуживания, проанализирован рынок сбыта и конкуренция. Также было подсчитано количество необходимого рабочего персонала для оптимальной работы предприятия, разработан план маркетинга и план производства.
39863. УСТРОЙСТВА ОПТОЭЛЕКТРОНИКИ 469.5 KB
  Созданы волоконные световоды с малыми потерями: затухание сигнала = 1 дБ км в ближней ИК области спектра. Наиболее широкополосны одномодовые световоды в области длин волн 126 132 мкм где материальная дисперсия кварцевых стёкол ближе к 0; полоса пропускания составляет 1011 Гцкм. Важными свойствами такого перехода является наличие обедненной носителями области перехода концентрирующей относительно сильное поле и области поглощения где поглощается падающий свет захватываются фотоны. Структура рn перехода: 1 обедненная область; 2 ...
39864. Обществознание. Учебник для 11 класа 2.83 MB
  азделение общества на группы называют социальной дифференциацией. Слово «дифференциация» происходит от латинского корня, означающего «различие». Многие исследователи считают, что дифференциация свойственна любому обществу. Даже в первобытных племенах выделялись группы в соответствии с полом и возрастом, с присущими им привилегиями и обязанностями.
39865. АВТОМОБИЛЬНЫЙ РЕЧЕВОЙ ИНФОРМАТОР 675 KB
  Курсовой проект выполнен на листах содержит чертежей.5 В каждый тарный ящик должен быть вложен упаковочный лист составленный по форме принятой предприятием изготовителем и согласованной с представителем заказчика. Упаковочный лист и сопроводительная документация должны находиться в пакете из полиэтиленовой плёнки ГОСТ 1035482. Сам корпус изготовляется из листового алюминиевого сплава толщиной 1 мм.
39866. Автомобильный речевой информатор 320.5 KB
  Были произведены необходимые инженерные расчеты. РАСЧЕТНОТЕОРЕТИЧЕСКАЯ ЧАСТЬ Анализ технического задания описания работы информатора речевого автомобильного 1. Проверочный схемотехнический расчет 1. Компоновочный расчет 1.
39867. Разработка стенда для снятия фазо-токовых характеристик ферритовых фазовращателей 1.11 MB
  Нормы времени на разработку определяем из Нормативы времени на разработку КД в НИИ и КБ по всем направлениям техники на ОАО ММЗ; п – количество листов документации ед; К комплексный коэффициент учитывающий различные условия организации труда и вид разрабатываемой аппаратуры. К1 = 10 – тип производства – единичное; К2 = 10 – условия применения РЭС – стационарная; К3 = 14 – количество чертежей; К4 = 06 – одновариантная разработка чертежа; К5 = 12 – перевод на иностранный язык; К6 = 085 – метод размножения документов; К7 = 09 –...
39868. АЦП с буферной памятью 1.45 MB
  Схема АЦП с буферной памятью состоит из следующих блоков: генератор тактовых импульсов, счётчик формирователь адресов, буферную память составляет динамическое ОЗУ, мультиплексор, регистр последовательного приближения, буферный регистр, компаратор, ЦАП и три логических элемента.
39869. Разработка электронного блока для ИИС мониторинга тепловых полей растений 1.38 MB
  В данном дипломном проекте представлен электронный блок для мониторинга температуры деревьев предназначенный для мониторинга окружающей среды. Блок обеспечивает длительные измерения и регистрацию температуру стволов деревьев в атомном режиме. Между состоянием деревьев их водным режимом и температурой стволов существует определенная связь. Температуре деревьев свойственен свой суточный ход который коррелирует с суточным ходом солнечной радиации и температуры воздуха и оказывается смещен по сравнению с ними во времени в зависимости от...