2301

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

Лекция

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

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

Украинкский

2013-01-06

106.36 KB

114 чел.

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

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


 

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

79774. Разработка устройства, позволяющего контролировать и управлять бытовым оборудованием через сеть GSM 388 KB
  При этом помимо необходимости охраны существует ещё необходимость отслеживать состояние бытового оборудования, прежде всего системы отопления, поскольку на значительной части территории России оставлять жилой доме без работающего отопления нельзя. При этом на рынке нет систем, предназначенных только для контроля работы системы отопления. Некоторые охранные сигнализации позволяют это делать, однако они их цена значительно выше, чем у специализированного устройства. Поэтому разработка устройства для удалённого контроля бытового оборудования целесообразна.
79775. Исследование взаимосвязи гендерных ролей и физического пола 214.5 KB
  Если век назад ребенок мог даже по внешнему виду определиться к какому полу принадлежит другой человек взрослый подросток длинные волосы прическа юбка брюки туфли ботинки косметика материал одежды то теперь его окружают мужчины с длинными волосами маникюром женщины коротко стриженные с бритым затылком в брюках рубашке не говоря уже о сигарете это картины нашего мира сейчас. Упомяну и последнии коллекции Готье Ив Сен Лорана Версачче и еще ряда модельеров одевших мужчин в юбки. Гипотезой исследования было: современные...
79776. Моделирование на уроках информатики в начальных классах как способ развития логического мышления 250.5 KB
  Современное образование требует преодоления разрозненности учебных предметов. Каждой научной дисциплине свойственно свое особое сочетание формализованных и неформализованных методов моделирования явлений, процедур доказательства и объяснения....
79777. Исследование принципов и методов рабочих процедур видеопроизводства на основе систем управления контентом 9.25 MB
  В данной дипломной работе проведено исследование принципов и методов рабочих процедур видеопроизводства на основе систем управления контентом. Разработана система процедур рабочего процесса по работе с медиаданными, создана возможность автоматизации рабочих процедур при работе в видеомонтажных программах в видеостудии кафедры ИКТ. Разработан фирменный стиль экспериментального Интернет-телеканал МИЭМ-ТВ.
79779. Анализ формирования портфеля ценных бумаг в коммерческом банке 876 KB
  Теоретические основы формирования портфеля ценных бумаг. Анализ политики формирования портфеля ценных бумаг в ЗАО АБ Газпромбанк. Состав портфеля ценных бумаг банка и принципы его формирования. Оценка эффективности портфеля ценных бумаг банка.
79781. Разработка проекта по совершенствованию управления финансовой устойчивостью на предприятии «Вятка» 1.68 MB
  Теоретические аспекты управления финансовой устойчивостью предприятия. Краткая экономическая характеристика предприятия. Мероприятия по укреплению финансовой устойчивости предприятия. Разработка модели оптимизации финансовой устойчивости предприятия.
79782. Особенности исчисления налога на добавленную стоимость на предприятии сферы услуг (на примере ООО «НетСервис.ру») 567 KB
  Несмотря на то что НДС является относительно новым налогом он используется в качестве основной формы косвенного налогообложения в большинстве государств в разных частях света и находящихся на различных этапах экономического развития. Из стран с наиболее развитыми...