Тема: Розв’язування диференціальних рівнянь при заданих крайових умовах та з частинними похідними


Скачати 143.02 Kb.
Назва Тема: Розв’язування диференціальних рівнянь при заданих крайових умовах та з частинними похідними
Дата 04.04.2013
Розмір 143.02 Kb.
Тип Документи
bibl.com.ua > Астрономія > Документи
ПРАКТИЧНА РОБОТА № 8

Тема: Розв’язування диференціальних рівнянь при заданих крайових умовах та з частинними похідними.

Мета: Придбання навичок роботи з методами вирозв’язку диференціальних рівнянь і систем диференціальних рівнянь в системі MathCad.

Теоретичні відомості

Розв’язку диференціальних рівнянь і систем диференціальних рівнянь. Нелінійні диференціальні рівняння та системи з такими рівняннями, як правило, не мають аналітичних методів розв’язку. Для їх розв’язку застосовуються чисельні методи, реалізація яких в системі MathCad здійснюється низкою функцій. Так для знаходження розв’язків звичайних диференціальних рівнянь, представлених в звичайній формі Коші, в MathCad передбачені наступні функції:

• rkfixed (y, x1, x2, n, F) - повертає матрицю рішень методом Рунге-Кута системи звичайних диференціальних рівнянь з початковими умовами у векторі y, праві частини яких записані у символьному векторі F на інтервалі від x1 до x2 при фіксованому числі кроків n.

• Rkadapt (y, x1, x2, n, F) - повертає матрицю рішень методом Рунге-Кута із змінним кроком для системи звичайних диференціальних рівнянь з початковими умовами у векторі y, праві частини яких записані в символьному векторі F на інтервалі x1 до x2 при фіксованому числі кроків n;

Функція rkfixed повертає матрицю з р +1 стовпцями (р - кількість рівнянь або порядок рівняння) і n рядками: перший стовпець (в MathCad - нульовий стовпець) - це n значень аргументу x від x1 до x2, наступні стовпчики - значення ординат розв’язку.

Можна розв’язати завдання більш точно (швидше), якщо зменшувати крок там, де похідна змінюється швидко, і збільшувати крок там, де вона веде себе більш спокійно. Для цього передбачена функція Rkadapt. Але, незважаючи на те, що вона при вирішенні диференціального рівняння використовує непостійний крок, функція Rkadapt надає відповідь для точок, що знаходяться на однаковій відстані, заданому користувачем. Аргументи і матриця, яка повертається функцією Rkadapt, такі ж, як при rkfixed. Диференціальні рівняння пов'язують аргументи, функції і похідні. Результат - це функція, яка при підстановці перетворює рівняння в тотожність. Дифрівняння має нескінченну безліч рішень. Для одержання частинного розв’язку із загального використовуються початкові умови. Для цього розв’язку в MathCad потрібно задати:

  • початкові умови,

  • список точок, в яких потрібно знайти розв’язку,

  • записати дифрівняння в спеціальній для MathCad формі,

  • застосувати одну з вбудованих функцій

Приклад пошуку розв’язку диференціального рівняння першого порядку:

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

.

Приклад пошуку рішення даного рівняння в системі MathCad має вигляд:

Крайова задача формулюється так: нехай на відрізку [a,b] потрібно знайти рішення диференціального рівняння (для простоти виклад будемо вести на прикладі дифрівняння другого порядку): image864, при граничних умовах у(а) = А, у(b) = В.

У цьому випадку Mathcad пропонує використовувати функцію sbval, щоб знайти відсутні початкові умови в точці а.

Sbval (v, а, b, D, load, score) - повертає вектор, що містить відсутні початкові умови в точці а.

Вектор v задає початкові наближення, а, b - граничні точки інтервалу рішень, D (x, y) - функція-вектор з першими похідними невідомих функцій. load (а, v) - функція-вектор, що повертає значення початкових умов у точці а. score (b, y) - функція-вектор, кожен елемент якого містить різницю між початковою умовою заданою в точці b, і значенням шуканого рішення в цій точці.



Після того, як ці відсутні початкові умови будуть отримані, можна вирішувати звичайну задачу з початковими умовами - завдання Коші, використовуючи будь-яку з функцій, описаних вище. Приклад розв'язання крайової задачі:

Порядок виконання практичної роботи:

Завдання 1

  1. Розв'язати задачу Коші: image868, y (0) = 1 з кроком h = 0.1 на відрізку [0, 1]:

  • методом Ейлера;

  • методом Рунге-Кутта

  • використовуючи функцію rkfixed

  • використовуючи функцію Rkadapt..

  1. Побудувати графіки функцій в одній системі координат, отриманих методом Ейлера та за допомогою функцій rkfixed та Rkadapt..

 №

f(x, y)



f(x, y)



f(x, y)

1

x + y

6

2 y - cos 2 x

11

2 y + 3 e - x

2

2 x 2 + 2 y

7

y - e x / 2 + 2

12

y / 2 - e - x

3

e x - 3 y

8

3 y - 2 sin x

13

y + (cos x) /3

4

y - sin x

9

e 2 x - y

14

y - 4 x + 5

5

y / 3 - x 2

10

2 sin x + y

15

2 x - y / 3 - e x

Завдання 2. Розв'язати задачу Коші для системи ОДУ при заданих початкових умовах на відрізку [0, 2] c кроком h = 0.2 за допомогою функції rkfixed. Побудувати графіки функцій u(t) і v(t). 

 №


Система

Початкові умови




Система

Початкові умови

 

u(0)

u’(0)

v(0)

v’(0)

 

u(0)

u’(0)

v(0)

v’(0)

1

image869

1.5

1.5

1

1

9

image870

2

0

-1

1

2

image871

-1

1

-1.5

3

10

image872

-1

2

-1.5

0

3

image873

1.5

1.5

1

1

11

image874

1.5

1.5

-1

-1

4

image875

1

1.5

0

2

12

image876

-1

1.5

0

-2

5

image877

0.5

1.5

-1

2

13

image878

-0.5

1

-1

2

6

image879

0.5

2

1

2

14

image880

0

-2

0

2

7

image881

5

5

-1

1

15

image882

3

3

-1

1

8

image883

1.5

1

3

1

 

 

 

 

 

 

Завдання 3. На відрізку [a, b] с використанням функцій load, score и sbval перетворити краєву задачу: image884= f(x, y, y) при граничних умовах y(a) = А, y(b) = В до задачі Коші і знайти розв’язок заданного дифрівняння в 10 проміжних точках за допомогою функції rkfixed.



варіанта

f(x, y, y’)

Граничні

умови

 

a

b

y(a)

y(b)

1

ex y + cos x

1

2

0

0

2

y sin x + e -x

2

3

1

0

3

y cos x + tg x

0

1

0

0.45

4

x3 y + cos x

0

1

1

0

5

x + ex y/(1 - x)

2

4

1

0.14

6

x2 y + 1/(1 + x)

1

3

0

0.17

7

y cos x + cos 2x

1

2

0

0

8

(2 + x) y + arctg x

0

3

0

0.22

9

(5 - x) y + x

2

4

0

-1.2

10

e -x y + 2 e -x

0

1.5

2.4

0

11

e -x y/x + x

-3

-2

3

0

12

(x2 + 1/x) y + 1/x 2

2

3

0

0

13

(10 - x) y + x

-1

0

2

0

14

y/x2 + x

1

3

1.5

0

15

y ln x + 1 + x

7

8

0

0

Завдання 4

Розв'язати диференціальне рівняння вищого порядку, використовуючи функцію rkfixed і функцію Rkadapt. Побудувати графіки знайденої функції та її похідних в одних осях.



Вар.


Завдання

Початкові умови


X1


X2


Y(x1)

Y'(x1)

Y''(x1)





2

5

1

-2








0

2.0

4.0

-3.0








1.0

3.0

-1.0

0








0

1.0

0

0








0

2.0

3.0

-1.0

1.0





0

2.0

0

-3.0

3.0





1.0

2.0

1.0

-2.0








2.0

3.0

1.0

-2.0








2.0

3.0

1.0

-2.0








2.0

10.0

1.0

-2.0








8.0

10.0

1.0

-2.0








2.0

3.0

1.0

-2.0








1.0

3.0

1.0

-2.0








1.0

3.0

1.0

-2.0








1.0

3.0

1.0

-2.0




Контрольні питання:

  1. Які чисельні методи розв’язання дифрівнянь вам відомі?

  2. Які існують типи диференціальних рівнянь другого порядку з частинними похідними?

  3. Як класифікують крайові задачі?

Зразок виконання роботи:





Графіки?


Схожі:

Тема. РОЗВ'ЯЗУВАННЯ РІВНЯНЬ ВИЩИХ СТЕПЕНІВ Заняття 1
Розв'язування рівнянь виду (х+а)(х+b)(х+с)(х+d) = А за умови, що а + b = с + d, або а + с = b + d, або а + d = b + с, де А Розв'язування...
Тема уроку: Розв’язування тригонометричних рівнянь
Навчальна: ознайомити учнів з іншими способами розв'язування тригонометричних рівнянь; навчити раціонально вибирати метод їх розв'язування;...
Урок №105 Тема. Розв'язування задач за допомогою рівнянь
Раціональні числа і дії над ними Тема Рівняння. Розв’язування рівнянь з однією змінною
Урок №63 Тема
Тема. Підсумковий урок з теми «Квадратний тричлен. Розв'язування рівнянь, що зводяться до квадратних рівнянь та їх використання для...
Розв’язування найпростіших тригонометричних рівнянь
Мета: Повторити властивості тригонометричних функцій та загальні розв’язки найпростіших тригонометричних рівнянь, розвивати вміння...
Урок №73 Тема. Системи двох лінійних рівнянь із двома змінними та...
Ня щодо залежності кількості розв'язків системи лінійних рівнянь від співвідношення коефіцієнтів a, b, c цих рівнянь; ви­роблення...
Урок №80 Тема. Розв'язування задач за допомогою систем лінійних рівнянь
Мета: відпрацювати навички застосування схеми розв'язання текстових задач на складання системи лінійних рівнянь із двома змінними...
Тема. Розв’язування тригонометричних рівнянь
Мета: узагальнити і систематизувати матеріал за темою “Розв'язування тригонометричних рівнянь ”, розвивати логічне мислення, уяву,...
Тема. Рівняння. Мета: освітня
Вироблення навичок розв’язування рівнянь за правилами знаходження невідомих компонентів додавання, віднімання і множення; формувати...
УРОК №71 Тема уроку. Системи рівнянь
Мета уроку: формування понять: «система рівнянь з двома змінними»; «розв'язки системи лінійних рівнянь з двома змінними»; «ознайомлення...
Додайте кнопку на своєму сайті:
Портал навчання


При копіюванні матеріалу обов'язкове зазначення активного посилання © 2013
звернутися до адміністрації
bibl.com.ua
Головна сторінка