Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Моделювання економічної динаміки (методичні вка...doc
Скачиваний:
0
Добавлен:
01.04.2025
Размер:
502.78 Кб
Скачать

2.2. Лінійні диференціальні рівняння вищих порядків

Розглянемо узагальнення класу рівнянь першого порядку на випадок рівнянь вищих порядків.

Визначення. Диференціальне рівняння n-го порядку називається лінійним, якщо воно має вигляд:

у(n) + a1(t)у(n-1) + a2(t)у(n-2) +.+ an(t) у = f(t) (5)

де а1(t), а2(t),..., аn(t), f(t) - неперервні функції.

Позначимо через L(y):

L(y)= у(n) + a1(t)у(n-1) + a2(t)у(n-2) +.+ an(t) у (6)

Відмітимо, що вираз такого вигляду називається лінійним диференціальним оператором n-го порядку.

З урахуванням наших позначень рівняння (5) може бути записано у вигляді

L(y)= f(t) (7)

Рівняння

L(y)= 0 (8)

називається лінійним однорідним рівнянням, відповідним рівнянню (7). У протилежність цьому рівнянню (8) (при f(x)0) називається неоднорідним. Наступне твердження зв'язує розв’язки рівнянь (7) і (8).

Теорема. Загальний розв’язок неоднорідного рівняння (6) є сума часткового розв’язку (х) цього рівняння і загального розв’язку відповідного йому однорідного рівняння (7).

Загальний розв’язок однорідного диференціального рівняння:

f(t C1, C2,..., Cn)= C1y1(t) +C2y2(t) +... + Cnyn(t).

Загальний розв’язок неоднорідного диференціального рівняння

у(t)=y(х) + f(t C1, C2,..., Cn).

Для пошуку загального розв’язку однорідного рівняння розв’язують характеристичне рівняння:

a0n + а1 n-1 +... + а n-1. +аn = 0

Тоді у(t)= С1 exp(1t)+ С2 exp(2 t) +... + Сn exp(n t).

Якщо * є кратним коренем кратності т n, то і

texp(*t), t2 exp(*t), ..., tmexp(*t)

також є корені однорідного рівняння.

При отриманні комплексно-зв'язаного кореня 1,2 = .

у(t)= 1 cost +С2 sint ).

Частковий розв’язок неоднорідного рівняння отримують підстановкою загального розв’язку однорідного рівняння в початкове (метод невизначених коефіцієнтів).

2.3. Поняття про різницеві рівняння

Різниця першого порядку представлена:

yt=f(t+1)-f(t)=yt+1 - yt.

Різниця другого порядку представлена:

2yt=yt+1 - yt= (уt+2 – yt+1) - (yt+1 – уt) = уt+2 – 2yt+1t.

Різниця n – порядку представлена:

nyt=n-1yt+1 + n-1yt

Визначення. Рівняння вигляду:

F(t, yt yt+1,.yt+k)=0 (9)

де k – фіксоване число, а t - довільний момент часу y1, y2, ., yn+1, yn+k члени деякої числової послідовності, називається різницевим рівнянням k-го порядку.

Розв’язати різницеве рівняння означає знайти всі послідовності { хе}, що задовольняють рівнянню (9). Різницеві рівняння часто використовуються в моделях економічної динаміки з дискретним часом, а також для наближеного розв’язку диференціальних рівнянь.

Визначення. Різницеве рівняння вигляду

В0yn + В1yt+1 + Вkyt+k=K (10)

де В0, В1, ...,Вk - деякі функції від t, називається лінійним різницевим рівнянням k-го порядку.

При K = 0 рівняння називають однорідним.

Частковий розв’язок неоднорідного рівняння отримують підстановкою загального розв’язку однорідного рівняння в початкове. У разі, коли коефіцієнти В0, В1 ...,Вk є константами, методи розв’язку даного класу рівнянь багато в чому аналогічні розв’язку лінійних диференціальних рівнянь з постійними коефіцієнтами. Покажемо це для різницевих рівнянь другого порядку:

yt +2 + pyt+1 + qyt=K (11)

Так само як і для лінійних диференціальних рівнянь, загальний розв’язок рівняння (11) визначається по формулі

yt = у(t )+ Y(t) (12)

де у(t) - деякий частковий розв’язок рівняння (12), Y(t) - загальний розв’язок відповідного однорідного рівняння (випадок K = 0). Для знаходження загального розв’язку однорідного рівняння необхідно спочатку вирішити характеристичне рівняння

2 + p + q=0 (13)

Після цього можуть виникнути три варіанти.

1. Обидва корені 1 і 2 - дійсні і різні. Тоді загальний розв’язок знаходиться по формулі:

Y(n)=C1 1t + C22t (14)

де С1, і С2 - довільні константи.

2. Обидва корені дійсні і рівні (1 =2=), тоді

Y(n)=(C1 + nC2)t (15)

3. У разі комплексно зв'язаного кореня 12 = r (cos ± i sin)

Y(n)= rt (C1 cos t+ C2 sin t) (16)