2301

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

Лекция

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

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

Украинкский

2013-01-06

106.36 KB

88 чел.

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

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


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

 

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

39031. Базовые технологии доступа к БД в Borland C++ Builder. Лабораторная работа 125.5 KB
  Указания к выполнению лабораторной работы Ранее уже говорилось что наборы данных представляют собой группы записей переданных из базы данных в приложения для просмотра и редактирования. Каждый набор данных инкапсулирован в специальном компоненте доступа к данным. Основные свойства и методы базового класса наборов данных TDtSet уже были рассмотрены нами ранее.
39032. Современные технологии доступа к БД 117 KB
  Технология ODBC все еще признается в качестве отраслевого стандарта доступа к базам данных однако также не развивается поэтому разработчики программного обеспечения все чаще обращают свое внимание на современные и возможно более эффективные технологии. С конца 90ых годов и до настоящего времени фирма Microsoft пытается создать универсальную платформу доступа к разнородным хранилищам данных для семейства своих операционных систем Windows. Источниками данных в данном случае могут являться не только SQLсервера Баз Данных но и иные...
39033. Технологии разработки распределенных информационных систем 121.5 KB
  Получить представление о компонентной объектной модели COM. Развитием трехуровневой архитектуры является так называемая многоуровневая nуровневая организация вычислений Multitier computing когда информационная система состоит из большого количества удаленных друг от друга объектов серверов каждый из которых может предоставлять другим объектам клиентам разнообразные информационные услуги. В настоящее время наиболее популярны следующие технологии организации распределенных информационных систем: Технология COM и ее развитие...
39034. Разработка технического задания на разработку ИС 77.5 KB
  Техническое задание ТЗ – это завершающий предпроектную стадию документ который содержит цели и обоснование проектирования а также определяет основные требования к ИС и исходные данные необходимые при разработке. В настоящее время при составлении технического задания обычно руководствуются требованиями следующих ГОСТов: 34.60289 Техническое задание на создание автоматизированной системы – описывает состав и содержание ТЗ которые распространяются на автоматизированную информационную систему в целом в том числе: общесистемные...
39035. Базовые технологии доступа к БД в Borland C++ Builder и их принципы 156 KB
  Указания к выполнению лабораторной работы Современные информационные системы не могут существовать без Баз Данных. По этой причине современные средства разработки приложений должны обеспечивать программиста средствами которые бы: Обеспечивали универсальный механизм доступа к базам данных построенных с использованием различных СУБД; Обеспечивали приемлемый уровень эффективности; Позволяли быстро разрабатывать полнофункциональные приложения для работы с БД любого размера. В Borlnd C Builder предлагается большое количество компонентов...
39036. Использование Borland C++ Builder в качестве средства быстрой разработки приложений 46 KB
  Компоненты для изучения: BitBtn StringGrid вкладка dditionl DteTimePicker вкладка Win32. Компоненты для изучения: StringGrid вкладка dditionl PgeControl вкладка Win32 Timer вкладка System. Компоненты для изучения: ColorBox вкладка dditionl RichEdit вкладка Win32 Timer вкладка System. Компоненты для изучения: MskEdit вкладка dditionl PgeControl вкладка Win32 Timer вкладка System.
39037. Разработка многооконных приложений с использованием Borland C++ Builder 88.5 KB
  Для добавления новой формы в проект приложения необходимо вызвать команду File New Form главного меню или нажать соответствующую кнопку на панели инструментов. Для каждой формы приложения вызывается метод CreteForm создать форму объекта приложения ppliction. Главная форма автоматически отображается на экране при запуске приложения в то время как остальные формы которые иногда называют вторичными будут созданы но останутся невидимыми для пользователя до тех пор пока не будут явно...
39038. Количество информации. Мера Хартли и мера Шеннона 80.5 KB
  Рассмотрение предложенных способов измерения количества информации удобно начать с примера. Тем не менее только на основе априорной информации мы не можем точно сказать какое именно число очков выпало в результате конкретного подбрасывания. С поступлением новой информации о результате подбрасывания эта неопределенность может уменьшаться.
39039. Понятие информационной системы 98.5 KB
  Сейчас пришло время дать этому понятию более точное определение: Информационная система ИС это взаимосвязанная совокупность средств методов и персонала используемых для хранения обработки и выдачи информации в интересах достижения поставленной цели. Как видно из определения информационные системы обеспечивают сбор хранение обработку поиск выдачу информации необходимой в процессе принятия решений задач из любой области. Основными элементами работы информационных систем являются: ввод новой информации и выдача текущей информации по...