загрузка...

2301

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

Лекция

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

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

Украинкский

2013-01-06

106.36 KB

85 чел.

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

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


Данной работой Вы можете всегда поделиться с другими людьми, они вам буду только благодарны!!!
Кнопки "поделиться работой":

 

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

42408. Работа с простыми объектами в Corel Draw 66 KB
  Порядок выполнения работы Для того чтобы нарисовать линию необходимо воспользоваться горячей клавишей [F5] или иконкой на панели инструментов. Для того чтобы вставить прямоугольник можно воспользоваться иконкой на панели инструментов или горячей клавишей [F6]. Для того чтобы вставить окружность можно воспользоваться горячей клавишей [F7] или воспользоваться иконкой на панели задач. Как и с прямоугольником можно воспользоваться дополнительными клавишами [Ctrl и или Shift].
42409. Работа с текстом и заливкой в Corel Draw 383.5 KB
  Для того чтобы залить объект заливкой можно воспользоваться одним из инструментов: Fill Color Dilog [ShiftF11] Fountin Fill Dilog [F11] Pttern Fill Dilog Texture Fill Dilog PostScript Fill Dilog Not Fill. Fill Color Dilog Fountin Fill Dilog Pttern Fill Dilog Texture Fill Dilog PostScript Fill Dilog Not Fill – убирает заливку Содержание отчета Тема и цель лабораторной работы Отчет о проделанной работе Вывод Контрольные вопросы Как с помощью горячих клавиш вставить текст Какие параметры текста можно выставить на...
42410. Создаем иллюзию стекла в Corel Draw 356.5 KB
  Создайте продолговатый треугольник. Выделите созданный вами треугольник откройте панель Effects Extrude выдавливание и нажмите кнопку Edit редактирование. Под вашим треугольником пунктирной линией будет отображаться копия. Если вы точно следовали инструкциям ваш осколок должен быть похож на тот что изображен на рисунке: Теперь выделите исходный треугольник на созданном вами осколке.
42411. КОМПОЗИЦИЯ. Симметричные и ассиметричные композиции 81 KB
  Характеристики композиции Вы можете грамотно построить тени и перспективу подобрать цветовые сочетания шрифтовое оформление но работа не станет своеобразной и выразительной. Это и есть характеристики композиции. Формы элементов композиции и форма композиции в целом могут оказывать огромное эмоциональное воздействие.
42412. Композиция в технике 149 KB
  Например литая несущая конструкция должна быть такой формы чтобы не возникало сомнений это именно литье а не сварная или какаялибо иная конструкция. Поэтому можно говорить о тектонике литой формы тектонике легких штампованных несущих элементов и тектонике пластмассовых конструкций. Образно говоря тектоника это искренность формы в отношении конструкции и материала. Объемнопространственную структуру можно определить как эстетически осмысленную взаимосвязь формы предмета с его внутренним строением и внешним пространством рассмотрите...
42413. Построение изображения на плоскости 183.5 KB
  Точка съемки определятся следующими параметрами координатами: а удаленностью от объекта т. расстоянием с которого ведется съемка; б высотой установки фото или видеокамеры; в смещением фото или видеокамеры в сторону от ее центрального положения относительно снимаемого объекта определяющем направление съемки. Удаленность от объекта определяет масштаб изображения который увеличивается с приближением точки съемки к объекту и уменьшается с увеличением расстояния между точкой установки камеры и снимаемым объектом.
42414. Компьютерная дискретная математика 180.5 KB
  Высказывание  повествовательное утверждение которое имеет значение истинности т. Простое высказывание называется атомом сложное молекулой. Например: не Р это высказывание земля не плоская; Р или Q земля плоская или Маша доктор; Р и Q земля плоская и Маша доктор. Обозначим через Р высказывание логика забава а через Q сегодня пятница.
42415. Логика и доказательство. Доказательство: прямое, обратное, от противного. Метод математической индукции 73 KB
  Метод математической индукции. Рассмотреть метод математической индукции. Метод математической индукции можно сравнить с прогрессом. Принцип математической индукции  это следующая теорема: Пусть мы имеем бесконечную последовательность утверждений P1 P2 .
42416. Теория множеств. Операции над множествами. Диаграммы Венна 758 KB
  Тип данных представляет собой множество объектов со списком стандартных операций над ними. Множество  это совокупность объектов называемых элементами множества. Объекты которые образуют множество называются элементами этого множества. Пример: Множество S = {3 2 11 5 7}  элементы множества записывают в фигурных скобках.