Допомога у написанні освітніх робіт...
Допоможемо швидко та з гарантією якості!

Лінійні однорідні рівняння другого порядку із сталими коефіцієнтами (реферат)

РефератДопомога в написанніДізнатися вартістьмоєї роботи

Зауваження 3. Використаний метод підбору окремого частинного розв’язку рівняння (91) можна застосовувати лише для певних диференціальних рівнянь, а саме для лінійних рівнянь із сталими коефіцієнтами і з спеціальною правою частиною виду (92) або (98). В інших випадках частинний розв’язок треба шукати методом варіа­ції довільних сталих. Насамперед слід зазначити, що частинний розв’язок… Читати ще >

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

Лінійні однорідні рівняння другого порядку із сталими коефіцієнтами Розглянемо лінійне однорідне рівняння другого порядку.

у ' ' + р у ' + gy = 0 (85).

де р, qдійсні числа.

Ейлер запропонував шукати частинні розв’язки цього рівняння у вигляді.

у = e kx (86).

де k — стала (дійсна чи комплексна), яку треба знайти. Підставивши функцію (86) в рівняння (85), дістанемо.

e kx (k2 + pkq) = 0.

Оскільки e kx /= 0, то.

k2 + pk + q = 0 (87).

Отже, якщо k буде коренем рівняння 87), то функція (86) буде роз­в'язком рівняння (85). Квадратне рівняння (87) називається характеристичним рівнянням диференціального рівняння (85).

Позначимо корені характеристичного рівняння через k1 і k2. Можливі три випадки:

І. k 1 і k 2  — дійсні і різні числа ( k 1 /= k 2 );

II. k 1 і k 2  — комплексні числа ( k 1 . 2 = ± );

III. k 1 і k 2  — дійсні і рівні числа ( k 1 = k 2 ).

Розглянемо кожен випадок окремо.

І. Корені характеристичного рівняння дійсні і різні: k 1 /= k 2 . У цьому випадку частинними розв’язками рівняння (85) є функції.

у 1 = е k 1 x , y 2 = e k 2 x .

Ці розв’язки лінійно незалежні, тому що при k 1 /= k 2 .

y 1 y 2 = e ( k 1 - k 2 ) x /= const .

Згідно з теоремою 4 (п. 3.2) загальний розв’язок рівняння (85) знаходять за формулою.

y = C 1 e k 1 x + C 2 e k 2 x (88).

II. Корені характеристичного рівняння комплексно-спряжені:

k 1 = ± , k 2 = - .

Підставивши значення k 1 та k 2 у формулу (86), знайдемо розв’язки.

у 1 = е ( + ) x , y 2 = e ( - ) x .

За формулою Ейлера.

e i = cos + i sin .

маємо.

e ( + ) x = e e i = e ( cos + i sin ) ;

e ( - ) x = e e - i = e ( cos - i sin ) .

Зауважимо, що коли функція z (х) = u (х) + iє розв’язком рівняння (85), то розв’язками будуть також функції u (х) та Дійсно, підставивши функцію z (x) в рівняння (85), дістанемо:

u ' ' + ' ' i + p ( u ' + ' i ) + g ( u + ) 0 ,.

або.

( u ' ' + p u ' + gu ) + i ( ' ' + p ' + g ) = 0 .

Остання тотожність можлива, коли вирази в дужках дорівнюють нулю (гл. 7, п. 1.4). Це означає, що функції u та розв’язки рівняння (85). Згідно з цим зауваженням частинними розв’язками рівняння (85) є функції.

y 1 = e sin , y 2 = e cos .

Ці розв’язки лінійно незалежні, оскільки.

y 1 y 2 = e sin e cos = tg /= const ,.

тому загальний розв’язок рівняння (85) запишеться у вигляді.

y = e ( C 1 sin + C 2 cos ) (89).

ІІІ. Корені характеристичного рівняння дійсні і рівні: k 1 = k 2 = k . За формулою (86) дістанемо один з розв’язків:

у = e kx .

Другий розв’язок шукатимемо у вигляді y 2 = ue kx , де u — невідома функція від х. Знайшовши y 2 ' і y 2 ' ' та підставивши їх у рівняння (85), дістанемо.

( u ' ' + 2 u ' k + uk 2 ) e kx + p ( u ' + uk ) e kx + gue kx = 0 .

або.

u ' ' + ( 2 k + p ) u ' + ( k 2 + pk + g ) ) u = 0 .

Оскільки k  — корінь рівняння (87), то k 2 + pk + g ) = 0 і за теоремою Вієта 2 k = - р, тому 2k + і u" =0, звідки u = С1х + С2, де С1, С2 — довільні сталі. Поклавши С1 = 1, С2 = 0 (нас цікавить який-небудь розв’язок u (х) /= 0), знайдемо другий частинний розв’я­зок рівняння (85):

у2 = х хe kx .

Розв’язки y1 та у2 — лінійно незалежні, тому загальний розв’язок рівняння (85) має вигляд.

у = хe kx (С1 +С2х) (90).

Приклади.

1. Знайти загальний розв’язок рівняння.

y ' ' - 5 y ' + 6 y = 0 .

Складемо характеристичне рівняння k2 — 5 k + 6 = 0 і знайдемо його корені k 1 =2, k 2 =2. За формулою (88) шуканий розв’язок має вигляд:

y = C 1 e 2 x + C 2 e 3 x .

2. Розв’язати рівняння.

y ' ' + 4 y ' + 13 y = 0 .

Характеристичне рівняння k 2 + 4 k + 13 = 0 має комплексні корені k 1.2 = - 2 ± 3і. Загальний розв’язок дістанемо за формулою (89):

y = e - 2 [ ( C 1 sin 3 x + C 2 cos 3 x ) .

Неоднорідні диференціальні рівняння другого порядку із сталими коефіцієнтами. Рівняння із спеціальною правою частиною Розглянемо неоднорідне диференціальне рівняння.

y ' ' + p y ' + gy = f ( x ) (91).

де р, q — задані дійсні числа, f ( x ) 0 — задана функція, неперерв­на на деякому проміжку (а, b).

Згідно з теоремою п. 3.3, загальний розв’язок такого рівняння являє собою суму частинного розв’язку рівняння (91) і загального розв’язку відповідного однорідного рівняння. Загальний розв’язок однорідного рівняння ми вже знаходити вміємо, тому розглянемо детальніше питання про знаходження частинного розв’язку неоднорідного рівняння.

Насамперед слід зазначити, що частинний розв’язок неоднорідного диференціального рівняння (91) можна знайти в квадратурах методом варіації довільних сталих (п. 3.4). Проте для рівнянь із спеціальною правою частиною частинний розв’язок можна знайти значно простіше, не вдаючись до операції інтегрування.

Розглянемо деякі з таких рівнянь.

І. Нехай права частина в рівнянні (91) має вигляд.

f ( x ) = e P n ( x ) (92).

де дійсне число, Рn (х) — многочлен степеня n.

Можливі такі випадки:

а) число е є коренем характеристичного рівняння.

k2 + pk + q = 0 (93).

Тоді диференціальне рівняння (91) має частинний розв’язок виду.

y Q n ( x ) e = ( A 0 + A 1 x + . . . + A n x n ) e , (94).

де А0, А1,…, Аn — невизначені коефіцієнти.

Справді, підставляючи функцію (94) в рівняння (91), після скорочення на е дістанемо.

Q n ' ' ( x ) + ( 2 + p ) Q n ' ( x ) + ( 2 + p + g ) Q n ( x ) P n ( x ) , (95).

де Q n ' ' ( x )  — многочлен степеня n — 2t Q n ' ( x )  — многочлен степеня n — 1, a Q n ( x ) і P n ( x )  — многочлени степеня n. Таким чином, зліва і справа в тотожності (95) стоять многочлени степеня n. Прирівнюючи коефіцієнти при однакових степенях n, дістанемо систему n + 1 лі­нійних алгебраїчних рівнянь, з якої визначимо n + 1 невідомих коефіцієнтів аі многочлена Q n ( x ) .

Не зупиняючись далі на доведеннях (див. [24]), вкажемо форму, в якій потрібно шукати частинний розв’язок рівняння (91), залежно від виду правої частини f ( x ) цього рівняння;

б) якщо число бігається з одним коренем характеристичного рівняння (93), тобто є простим коренем цього рівняння, то частинний розв’язок рівняння (91) треба шукати у вигляді.

y xQ n ( x ) e , (96).

в) якщо число двократним коренем рівняння (93), то частинний розв’язок рівняння (91) шукають у вигляді.

y x 2 Q n ( x ) e (97).

Об'єднаємо випадки а) — в): якщо права частина рівняння (91) має вигляд (92), то частинний розв’язок цього рівняння треба шукати у вигляді.

y x r e Q n ( x ) .

де Q n ( x )  — многочлен з невизначеними коефіцієнтами того самого степеня, що й многочлен P n ( x ) , а r — число коренів характеристич­ного рівняння, які дорівнюють . Якщо не є коренем характерис­тичного рівняння, то приймаємо г = 0.

II. Нехай права частина в рівнянні (91) має вигляд.

f ( x ) = e ( P n ( x ) cos + R m ( x ) sin ) , (98).

де P n ( x )  — многочлен степеня n, Rm (x) — многочлен степеня mа дійсні числа. (Функція (92) є окремим випадком функції (98) і утворюється з неї при).

Частинний розв’язок рівняння (91) треба шукати у вигляді.

y x r e ( Q s ( x ) cos + L s ( x ) sin ) (99).

де Q s ( x ) та L s ( x )  — многочлени степеня s з невизначеними коефі­цієнтамиs- найвищий степінь многочленів Rm (х) та P n ( x ) , тобто s = max (n, m) — r — число коренів характеристичного рівняння, які дорівнюють.

Зокрема, якщо права частина рівняння (91) має вигляд.

f ( x ) = A cos Вsin (100).

де А, В — відомі дійсні числа, то частинний розв’язок цього рівняння треба шукати у вигляді.

y x r ( a cos + b sin ) , (101).

де a, b — невідомі коефіцієнтиг — число коренів характеристичного рівняння (93), які дорівнюють /p>

Зауваження 1. Шукані многочлени Q n ( x ) у формулах (94), (96) і (97) мають бути повними, тобто містити всі степені х від 0 до n, незалежно від того, чи повним є заданий многочлен P n ( x ) . Те саме стосується многочленів Q s ( x ) та L s ( x ) у формулі (99), причому невизначені коефіцієнти при одних і тих же степенях х у цих многочленах повинні бути, взагалі кажучи, різними.

Зауваження 2. Якщо права частина рівняння (91) є сумою декількох різних за структурою функцій виду (92) або (98), то для відшукання частинного розв’язку потрібно використати теорему про накладання розв’язків (п. 3.4).

Зауваження 3. Використаний метод підбору окремого частинного розв’язку рівняння (91) можна застосовувати лише для певних диференціальних рівнянь, а саме для лінійних рівнянь із сталими коефіцієнтами і з спеціальною правою частиною виду (92) або (98). В інших випадках частинний розв’язок треба шукати методом варіа­ції довільних сталих.

Приклади.

1. Розв’язати рівняння у" - 2у' + у = 2х + 3.

Характеристичне рівняння k2 — 2k + 1 = 0 має корені k1= k2 = 1, тому загальний розв’язок однорідного рівняння має вигляд у ( ч ) = у х ( С 1 + С 2 х ) . Оскільки правою частиною даного рівняння є функція виду Р 1 ( х ) е 0 х , причому, k 1 , /= k 2 , то за формулою (94) частинний розв’язок шукаємо у вигляді y Q 1 ( x ) e 0 x , тобто у* = А + Вх, де, А і В — невідомі коефіцієнти. Знайшовши похідні y ' = В, у*" = 0 і підставивши їх у рівняння, дістанемо.

  • -.2 В + А + Вх = 2х + 3.

Прирівнюючи коефіцієнти при однакових степенях, дістаємо систему рівнянь.

{ B = 2 - - 2 B + A = 3, .

звідки В = 2, А = 7. Отже, частинний розв’язок даного рівняння має вигляд у* = 7 + 2х, тому.

y = y ( x ) + y ( x ) = e x ( C 1 + C 2 x ) + 2 x + 7 .

шуканий загальний розв’язок.

2. Розв’язати рівняння у" - 3у' + 2у = 8е3х.

Характеристичне рівняння k2 — 3k + 2 = 0 має корені k 1 = 1 і k 2 = 2, тому загальний розв’язок однорідного рівняння має вигляд у = С 1 е х + С 2 е 2 х ) . Оскільки правою частиною даного рівняння є функція виду Р0(х)е3х, причому, k 1 , /= k 2 , то частинний розв’язок шукаємо у вигляді.

y Q 0 ( x ) e 3 x , тобто y Ae 3 x .

де, А — невідомий коефіцієнт.

Знайшовши похідні (у*)' = 3Ае3х, (у*)" = 9Ае3х і підставивши їх у рівняння, дістанемо.

9 Ae 3 x - 9 Ae 3 s + 2 Ae 3 x = 8 e 3 x ,.

звідки, А = 4, тому у* = 4e3х — частинний розв’язок даного рівняння, а у = С 1 е х + С 2 е 2 х + 4 е 3 х  — його загальний розв’язок.

3. Розв’язати рівняння у" + у = tg x.

Характеристичне рівняння k2 + 1 = 0 має корені k 1 . 2 = ±і, тому y = C 1 sin x + C 2 cos x  — загальний розв’язок відповідного однорідного рівняння. Права частина рівняння f ( x ) = tg x не є функцією спеціального виду (92) або (98), тому частинний розв’язок даного рівняння методом підбору шукати не можна.

Знайдемо цей розв’язок методом Лагранжа. Складемо систему виду (84) і розв’я­жемо її:

{ C 1 sin x + C 2 cos x = 0 - C 1 cos x - C 2 sin x = tgx , .

С 1 = sin x , C 2 = - sin 2 x cos x .

Інтегруючи, дістанемо.

C 1 = sin xdx = - cos x + C 1 .

C 2 = - sin 2 x cos x dx = - 1 - cos 2 x cos x dx = - dx cos x + cos xdx = - ln | tg ( x 2 + 4 ) | + sin x + C 2 .

де C 1 , C 2  — довільні сталі. При C 1 = C 2 = 0 дістанемо частинний розв’язок:

y - sin x cos x + ( sin x - ln | tg ( x 2 + 4 ) | ) cos x ,.

тоді.

y = C 1 sin x + C 2 cos x - cos x ln | tg ( x 2 + 4 ) | .

— загальний розв’язок даного рівняння.

Лінійні диференціальні рівняння n-гo порядку Застосуємо методи знаходження розв’язків диференціальних рівнянь другого порядку до рівнянь вищих порядків. Не зупиняючись детально на теорії (див. [26]), сформулюємо необхідні твердження і розглянемо приклади.

Нехай маємо лінійне диференціальне рівняння nго порядку.

y ( n ) + a 1 y ( n - 1 ) + . . . + a n y = 0 (102).

де а1, а2, … аn — сталі дійсні числа.

Характеристичним для рівняння (102) називається алгебраїчне рівняння nго степеня виду.

k n + a 1 k n - 1 + . . . + a n - 1 k + a n = 0 (103).

де k — невідоме дійсне чи комплексне число.

Як відомо (гл. 7, п. 1.5), рівняння (103) має n коренів. Позначимо ці корені через k 1 , k 2 , …, k n .

Теорема. Кожному простому кореню k рівняння (103) відповідає частинний розв’язок, e kx рівняння (102), а кожному кореню k кратнос­ті m > 1 відповідає m частинних розв’язків вuду ekx, xekx, …, xm-1 ekx.

Кожній парі простих комплексно-спряжених коренів рівняння (103) відповідає два частинних розв’язки e sin та e cos рівняння (102), ожній парі а ± комплексно-спряжених коренів кратності - 1 відповідає 2р частинних розв’язків виду.

e sin , xe sin , . . . x p - 1 e sin x .

e cos , xe cos , . . . x p - 1 e cos x .

Загальна сума кратностей всіх коренів рівняння (103) дорівнює п, тому кількість всіх частинних розв’язків рівняння (102), складених згідно з цією теоремою, дорівнює п, тобто збігається з порядком рівняння (102). Позначимо ці частинні розв’язки через y1, y2, …, уn. Мож­на показати, що знайдені частинні розв’язки є лінійно незалежними, і загальний розв’язок рівняння (102) знаходиться за формулою.

y = C 1 y 1 + C 2 y 2 + . . . + C n y n (104).

Нехай тепер задано неоднорідне рівняння n-го порядку.

y ( n ) + a 1 y ( n - 1 ) + . . . + a n - 1 y ' + a n y = f ( x ) .

де а1, а2, …, аn — сталі дійсні числа, f ( x ) /= 0 — неперервна на деяко­му проміжку функція.

Як і для рівнянь другого порядку, загальним розв’язком рівняння (105) є функція.

у = у ( х ) + у ( х ) .

де у (х) — загальний розв’язок відповідного однорідного рівняння (102), а у*(х) — частинний розв’язок рівняння (105).

Побудову загального розв’язку у (х) рівняння (102) з’ясовано. Проаналізуємо знаходження частинного розв’язку у*(х). Якщо права частина f ( x ) рівняння (105) є функцією спеціального виду (98), то час­тинний розв’язок цього рівняння треба шукати за формулою (99). Як­що права частина f ( x ) не є функцією виду (98), то для знаходження у*(х) застосовують метод варіації довільних сталих. Стосовно рів­няння (105) суть цього методу така.

Нехай функція (104) є загальним розв’язком відповідного однорідного рівняння (102). Знаходимо частинний розв’язок рівняння (105) за тією ж формулою (104), вважаючи, що величини C1, С2, … Сn — функції від х, тобто покладемо.

y ( x ) = C 1 ( x ) y 1 + C 2 ( x ) y 2 + . . . + C n ( x ) y n (106).

де C 1 ( x ) , C 2 ( x ) , . . . , C n ( x ) невідомі функції.

Складемо систему рівнянь.

{ C 1 y 1 + C 2 + . . . + C n y n = 0, C 1 y 1 + C 2 y 2 + . . . + C n y n = 0 . . . . . . . . . . C 1 y 1 ( n - 2 ) + C 2 y 2 ( n - 2 ) + . . . + C n y n ( n - 2 ) = 0 C 1 y 1 ( n - 1 ) + C 2 y 2 ( n - 1 ) + . . . + C n y n ( n - 1 ) = f ( x ) .

Розв’язуючи цю систему, знаходимо похідні С і ( х ) , і = 1, 2, … n, а потім інтегруванням і самі функції Сі(х). Якщо взяти всі сталі інтегрування рівними нулю і підставити функції Сі(х) в рівність (106), то матимемо частинний розв’язок рівняння (105) — якщо у рівність (106) підставити функції Сі(х) + С і де С і  — довільні сталі, то зразу ді­станемо загальний розв’язок.

Приклади.

1. Розв’язати рівняння y ( 5 ) + 4 y ( 3 ) = 0 .

Характеристичне рівняння k 5 + 4 k 3 = 0 має корені k 1 = k 2 = k3 = 0, k 4 = 2 = - 2і. Згідно з теоремою маємо частинні розв’язки: y 1 = 1, y 2 = x , y 3 = x 2 , y 4 = cos 2 x , y 5 = sin 2 x . Загальний розв’язок даного рівняння знаходимо за форму­лою (104):

y = C 1 + C 2 x + C 3 x 2 + C 4 cos 2 x + C 5 sin 2 x .

2. Знайти розв’язок рівняння.

y ' ' ' = C 1 + C 2 x + C 3 x 2 + C 4 cos 2 x + C 5 sin 2 x ), який задоволь­няє початкові умови у (0) = 2, у'(0) = 2, у" (0) = - 1.

Маємо неоднорідне лінійне рівняння третього порядку. Характеристичне рівняння k 3 - 2 k 2 + k = 0 має корені k1 = 0, k 2 = 1, k 3 = 1. Загальний розв’язок відповідного однорідного рівняння має вигляд y = C 1 + ( C 2 + C 3 x ) e x . Правою частиною даного рівняння є функція виду (100), де, А = В = 4,. Оскільки число 1 · і = і не є коренем характеристичного рівняння (г = 0), то окремий роз­в'язок шукаємо у вигляді (101):

y cos x + b sin x .

де а, b — невідомі коефіцієнти. Знайшовши похідні y ' , y n і підставивши їх у дане рівняння, після спрощень дістанемо.

2 х + 2 b sin х = 4 cos х + 4 sin x,.

звідки, а = b = 2, тому у* = 2 cos х + 2 sin х — частинний розв’язок неоднорідного рівняння, а.

y = y + y C 1 + ( C 2 + C 3 x ) e x + 2 ( sin x + cos x ) .

— загальний розв’язок. Продиференціювавши його двічі, знайдемо.

y ' = ( C 2 + C 3 ( 1 + x ) ) e x + 2 ( - sin x + cos x ) .

y n = ( C 2 + C 3 ( 2 + x ) ) e x - 2 ( cos x + sin x ) .

Скориставшись початковими умовами у (0) = 2, y ' (0) = 2, у" (0) = -1, дістанемо систему рівнянь:

{ C 1 + C 2 + 2 = 2 C 2 + C 3 + 2 = 2 C 2 + 2 C 3 - 2 = - 1 .

звідки С1 = 1, С2 = -1, С3 = 1. Отже, шуканий розв’язок має вигляд.

y = 1 + ( x - 1 ) e x + 2 ( cos x + sin x ) .

Показати весь текст
Заповнити форму поточною роботою