Термінова допомога студентам
Дипломи, курсові, реферати, контрольні...

Лінійні неоднорідні системи (реферат)

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

X 1 неод. .. ... x n неод righ e pt (C 0 1 t s + C 1 1 t s — 1 +. .. + C s — 1 1 t + C s 1) cos qt e pt (C 0 2 t s + C 1 2 t s — 1 +. .. + C s — 1 2 t + C s 2) cos qt. .. .. .. .. .. .. .. .. .. .. .. .. .. .. .. .. .. .. .. .. .. .. .. e pt (C 0 n t s + C 1 n t s — 1 +. .. + C s — 1 n t + C s n) cos qt righ () + e pt (D 0 1 t s + D 1 1 t s — 1… Читати ще >

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

Реферат на тему:

Лінійні неоднорідні системи

Система диференціальних рівнянь, що записана у вигляді.

x 1 = a 11 ( t ) x 1 + a 12 ( t ) x 2 + . . . + a 1 n ( t ) x n + f 1 ( t ) x 2 = a 21 ( t ) x 1 + a 22 ( t ) x 2 + . . . + a 2 n ( t ) x n + f 2 ( t ) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . x n = a n 1 ( t ) x 1 + a n 2 ( t ) x 2 + . . . + a nn ( t ) x n + f n ( t ) , { { { .

чи у векторно-матричному вигляді.

x = A ( t ) x + f ( t ) .

називається системою лінійних неоднорідних диференціальних рівнянь.

1. Властивості розв’язків лінійних неоднорідних систем

Властивість 1. Якщо вектор x 1 ( t ) x 2 ( t ) . . . x n ( t ) righ ( ) ( ) ( ) x ( t ) = є.

розв’язком лінійної неоднорідної системи, a x 1 0 ( t ) . . . x n 0 ( t ) righ ( ) ( ) ( ) x 0 ( t ) = розв’язком відповідної лінійної однорідної системи, то сума x ( t ) + x 0 ( t )  — є розв’язком лінійної неоднорідної системи.

Дійсно, за умовою.

x ( t ) - A ( t ) x ( t ) f ( t ) і x 0 ( t ) - A ( t ) x 0 ( t ) 0 .

Але тоді і.

d dt [ x ( t ) + x 0 ( t ) ] - A ( t ) [ x ( t ) + x 0 ( t ) ] = [ d dt x ( t ) - A ( t ) x ( t ) ] + [ d dt x 0 ( t ) - A ( t ) x 0 ( t ) ] f ( t ) , .

тобто x ( t ) + x 0 ( t ) є розв’язком неоднорідної системи.

Властивість 2 (Принцип суперпозиції). Якщо вектори x 1 i ( t ) x 2 i ( t ) . . . x ni ( t ) righ ( ) ( ) ( ) x i = , i = 1, n є розв’язками лінійних неоднорідних систем.

x = A ( t ) x + f i ( t ) , i = 1, n ,.

де f 1 i ( t ) . . . . . . . . f ni ( t ) righ ( ) ( ) f i ( t ) = , то вектор x ( t ) = i = 1 n C i x i ( t ) , де C i  — довільні сталі буде розв’язком лінійної неоднорідної системи.

x = A ( t ) x + i = 1 n C i f i ( t ) .

Дійсно, за умовою виконуються n  — тотожностей.

x i ( t ) - A ( t ) x i ( t ) f i ( t ) .

Склавши лінійну комбінацію з лівих і правих частин, одержимо.

d dt [ i = 1 n C i x i ( t ) ] - A ( t ) [ i = 1 n C i x i ( t ) ] = i = 1 n C i [ x i ( t ) - A ( t ) x i ( t ) ] i = 1 n C i f i ( t ) ,.

тобто лінійна комбінація x ( t ) = i = 1 n C i x i ( t ) буде розв’язком системи.

x = A ( t ) x + i = 1 n C i f i ( t ) .

Властивість 3. Якщо комплексний вектор з дійсними елементами u 1 ( t ) . . . u n ( t ) righ v 1 ( t ) . . . v n ( t ) righ ( ) ( ) x ( t ) = u ( t ) + iv ( t ) = є розв’язком неоднорідної системи x = A ( t ) x + f ( t ) , де f ( t ) = p ( t ) + iq ( t ) , p 1 ( t ) . . . . . . . p n ( t ) righ ( ) ( ) p ( t ) = , q 1 ( t ) . . . . . . . q n ( t ) righ ( ) ( ) q ( t ) = , то окремо дійсна і уявна частини є розв’язками системи.

Дійсно, за умовою.

d dt [ u ( t ) + iv ( t ) ] - A ( t ) [ u ( t ) + iv ( t ) ] p ( t ) + iq ( t ) .

Розкривши дужки і перетворивши, одержимо.

[ u ( t ) - A ( t ) u ( t ) ] + i [ v ( t ) - A ( t ) v ( t ) ] p ( t ) + iq ( t ) .

Але комплексні вирази рівні між собою тоді і тільки тоді, коли рівні дійсні та уявні частини, що і було потрібно довести.

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

Доведення. Нехай x 0 ( t ) = i = 1 n C i x i ( t )  — загальний розв’язок однорідної системи і x ( t )  — частинний розв’язок неоднорідної. Тоді, як випливає з властивості 1, їхня сума x ( t ) + x ( t ) буде розв’язком неоднорідної системи.

Покажемо, що цей розв’язок загальний, тобто підбором сталих C i , i = 1, n можна розв’язати довільну задачу Коші.

x 1 ( t 0 ) = x 1 0 , x 2 ( t 0 ) = x 2 0 , . . . , x n ( t 0 ) = x n 0 .

Оскільки x 0 ( t ) = i = 1 n C i x i ( t )  — загальний розв’язок однорідного рівняння, то вектори x 1 ( t ) , . . . , x n ( t ) лінійно незалежні W [ x 1 ( t ) , . . . , x n ( t ) ] /= 0 і система алгебраїчних рівнянь.

C 1 x 11 ( t 0 ) + C 2 x 12 ( t 0 ) + . . . + C n x 1 n ( t 0 ) = x 1 0 - x 1 ( t 0 ) C 1 x 21 ( t 0 ) + C 2 x 22 ( t 0 ) + . . . + C n x 2 n ( t 0 ) = x 2 0 - x 2 ( t 0 ) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . C 1 x n 1 ( t 0 ) + C 2 x n 2 ( t 0 ) + . . . + C n x nn ( t 0 ) = x n 0 - x n ( t 0 ) { { { .

має єдине розв’язок C i 0 , i = 1, n . І лінійна комбінація x ( t ) = i = 1 n C i 0 x i ( t ) + x ( t ) с отриманими сталими C i 0 , i = 1, n є розв’язком поставленої задачі Коші.

2. Побудова частинного розв’язку неоднорідної системи методом варіації довільних сталих

Як випливає з останньої теореми, для побудови загального розв’язку неоднорідної системи потрібно розв’язати однорідну і яким-небудь засобом знайти частинний розв’язок неоднорідної системи. Розглянемо метод, який називається методом варіації довільної сталої.

Нехай маємо систему.

x = A ( t ) x + f ( t ) .

і x ( t ) = i = 1 n C i x i ( t )  — загальний розв’язок однорідної системи. Розв’язок неоднорідної будемо шукати в такому ж вигляді, але вважати C i не сталими, а невідомими функціями, тобто C i = C i ( t ) і x неодн ( t ) = i = 1 n C i ( t ) x i ( t ) , чи в матричній формі.

x неодн ( t ) = X ( t ) C ( t ) ,.

де X ( t ) -фундаментальна матриця розв’язків, C ( t )  — вектор з невідомих функцій. Підставивши в систему, одержимо.

d dt X ( t ) C ( t ) + X ( t ) dC ( t ) dt = A ( t ) X ( t ) C ( t ) + f ( t ) ,.

чи.

[ d dt X ( t ) - A ( t ) X ( t ) ] C ( t ) + X ( t ) dC ( t ) dt = f ( t ) .

Оскільки X ( t )  — фундаментальна матриця, тобто матриця складена з розв’язків, то.

d dt X ( t ) - A ( t ) X ( t ) 0 .

і залишається система рівнянь X ( t ) C ' ( t ) = f ( t ) .

Розписавши покоординатно, одержимо.

C ' 1 x 11 ( t ) + C ' 2 x 12 ( t ) + . . . + C ' n x 1 n ( t ) = f 1 ( t ) C ' 1 x 21 ( t ) + C ' 2 x 22 ( t ) + . . . + C ' n x 2 n ( t ) = f 2 ( t ) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . C ' 1 x n 1 ( t ) + C ' 2 x n 2 ( t ) + . . . + C ' n x nn ( t ) = f n ( t ) . { { { .

Оскільки визначником системи є визначник Вронського і він не дорівнює нулю, то система має єдиний розв’язок і функції C i ( t ) визначаються в такий спосіб.

f 1 ( t ) x 12 ( t ) x 1 n ( t ) f 2 ( t ) x 22 ( t ) x 2 n ( t ) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . f n ( t ) x n 2 ( t ) x nn ( t ) rli | | | | | | C 1 ( t ) = .

x 11 ( t ) f 1 ( t ) x 1 n ( t ) x 21 ( t ) f 2 ( t ) x 2 n ( t ) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . x n 1 ( t ) f n ( t ) x nn ( t ) rli | | | | | | C 2 ( t ) = .

x 11 ( t ) x 12 ( t ) f 1 ( t ) x 21 ( t ) x 22 ( t ) f 2 ( t ) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . x n 1 ( t ) x n 2 ( t ) f n ( t ) rli W [ x 1 ( t ) , . . . , x n ( t ) ] | | | | | | C n ( t ) = .

Звідси частинний розв’язок неоднорідної системи має вигляд.

x неодн = i = 1 n C i ( t ) x i ( t ) .

Для лінійної неоднорідної системи на площині.

x 1 = a 11 ( t ) x 1 + a 12 ( t ) x 2 + f 1 ( t ) x 2 = a 21 ( t ) x 1 + a 22 ( t ) x 2 + f 2 ( t ) { .

метод варіації довільної сталої реалізується таким чином.

Нехай.

X ( t ) = [ x 11 ( t ) x 12 ( t ) x 21 ( t ) x 22 ( t ) ] .

Фундаментальна матриця розв’язків однорідної системи. Тоді частинний розв’язок неоднорідної шукається у вигляді.

C ' 1 x 11 ( t ) + C ' 2 x 12 ( t ) = f 1 ( t ) C ' 1 x 21 ( t ) + C ' 2 x 22 ( t ) = f 2 ( t ) . { .

Звідси.

f 1 ( t ) x 12 ( t ) f 2 ( t ) x 22 ( t ) rli | | C 1 ( t ) = x 11 ( t ) f 1 ( t ) x 21 ( t ) f 2 ( t ) rli | | C 2 ( t ) = .

І загальний розв’язок має вигляд.

( x 1 ( t ) x 2 ( t ) ) = [ x 11 ( t ) x 12 ( t ) x 21 ( t ) x 22 ( t ) ] ( C 1 C 2 ) + ( x 1 ( t ) x 2 ( t ) ) , ( x 1 ( t ) x 2 ( t ) ) = [ x 11 ( t ) x 12 ( t ) x 21 ( t ) x 22 ( t ) ] ( C 1 ( t ) C 2 ( t ) ) ,.

де C 1 , C 2  — довільні сталі.

4. Метод невизначених коефіцієнтів

Якщо система лінійних диференціальних рівнянь з сталими коефіцієнтами, а векторна функція f ( t ) спеціального виду, то частинний розв’язок можна знайти методом невизначених коефіцієнтів. Доведення існування частинного розв’язку зазначеного виду аналогічно доведенню для лінійних рівнянь вищих порядків.

1) Нехай кожна з компонент вектора f ( x ) є многочленом степеня не більш ніж s , тобто.

f 1 ( t ) . . . . . . f n ( t ) righ A 0 1 t s + A 1 1 t s - 1 + . . . + A s - 1 1 t + A s 1 A 0 2 t s + A 1 2 t s - 1 + . . . + A s - 1 2 t + A s 2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A 0 n t s + A 1 n t s - 1 + . . . + A s - 1 n t + A s n righ ( ) ( ) .

а) Якщо характеристичне рівняння не має нульового кореня, тобто i /= 0 , i = 1, n , то частинний розв’язок шукається в такому ж вигляді, тобто.

x 1 неод . . . . . . x n неод righ B 0 1 t s + B 1 1 t s - 1 + . . . + B s - 1 1 t + B s 1 B 0 2 t s + B 1 2 t s - 1 + . . . + B s - 1 2 t + B s 2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . B 0 n t s + B 1 n t s - 1 + . . . + B s - 1 n t + B s n righ ( ) ( ) .

б) Якщо характеристичне рівняння має нульовий корінь кратності r , тобто 1 = 2 = . . . = r , те частинний розв’язок шукається у вигляді многочлена степеня s + r , тобто.

x 1 неод . . . . . . x n неод righ B 0 1 t s + r + B 1 1 t s + r - 1 + . . . + B s + r 1 B 0 2 t s + r + B 1 2 t s + r - 1 + . . . + B s + r 2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . B 0 n t s + r + B 1 n t s + r - 1 + . . . + B s + r n righ ( ) ( ) .

Причому перші ( s + 1 ) n  — коефіцієнти B i j , i = 0, s , j = 1, n знаходяться точно, а інші з точністю до сталих інтегрування C 1 , . . . , C n , що входять у загальний розв’язок однорідних систем.

2) Нехай f ( t ) має вид.

f 1 ( t ) . . . . . . f n ( t ) righ e pt ( A 0 1 t s + A 1 1 t s - 1 + . . . + A s - 1 1 t + A s 1 ) e pt ( A 0 2 t s + A 1 2 t s - 1 + . . . + A s - 1 2 t + A s 2 ) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . e pt ( A 0 n t s + A 1 n t s - 1 + . . . + A s - 1 n t + A s n ) righ ( ) ( ) .

а) Якщо характеристичне рівняння не має коренем значення p , тобто i /= p , i = 1, n , то частинний розв’язок шукається в такому ж вигляді, тобто.

x 1 неод . . . . . . x n неод righ e pt ( B 0 1 t s + B 1 1 t s - 1 + . . . + B s - 1 1 t + B s 1 ) e pt ( B 0 2 t s + B 1 2 t s - 1 + . . . + B s - 1 2 t + B s 2 ) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . e pt ( B 0 n t s + B 1 n t s - 1 + . . . + B s - 1 n t + B s n ) righ ( ) ( ) .

б) Якщо p є коренем характеристичного рівняння кратності r , тобто 1 = 2 = . . . = r = p , то частинний розв’язок шукається у вигляді.

x 1 неод . . . . . . x n неод righ e pt ( B 0 1 t s + r + B 1 1 t s + r - 1 + . . . + B s + r 1 ) e pt ( B 0 2 t s + r + B 1 2 t s + r - 1 + . . . + B s + r 2 ) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . e pt ( B 0 n t s + r + B 1 n t s + r - 1 + . . . + B s + r n ) righ ( ) ( ) .

І, як у попередньому пункті, перші ( s + 1 ) n  — коефіцієнти B i j , i = 0, s , j = 1, n знаходяться точно, а інші з точністю до сталої інтегрування C 1 , . . . , C n .

3) Нехай f ( t ) має вигляд:

f 1 ( t ) . . . . . . f n ( t ) righ e pt ( A 0 1 t s + A 1 1 t s - 1 + . . . + A s - 1 1 t + A s 1 ) cos qt e pt ( A 0 2 t s + A 1 2 t s - 1 + . . . + A s - 1 2 t + A s 2 ) cos qt . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . e pt ( A 0 n t s + A 1 n t s - 1 + . . . + A s - 1 n t + A s n ) cos qt righ ( ) + e pt ( B 0 1 t s + B 1 1 t s - 1 + . . . + B s - 1 1 t + B s 1 ) sin qt e pt ( B 0 2 t s + B 1 2 t s - 1 + . . . + B s - 1 2 t + B s 2 ) sin qt . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . e pt ( B 0 n t s + B 1 n t s - 1 + . . . + B s - 1 n t + B s n ) sin qt righ ( ) ( ) .

а) Якщо характеристичне рівняння не має коренем значення p ± iq , то частинний розв’язок шукається в такому ж вигляді, тобто.

x 1 неод . . . . . . x n неод righ e pt ( C 0 1 t s + C 1 1 t s - 1 + . . . + C s - 1 1 t + C s 1 ) cos qt e pt ( C 0 2 t s + C 1 2 t s - 1 + . . . + C s - 1 2 t + C s 2 ) cos qt . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . e pt ( C 0 n t s + C 1 n t s - 1 + . . . + C s - 1 n t + C s n ) cos qt righ ( ) + e pt ( D 0 1 t s + D 1 1 t s - 1 + . . . + D s - 1 1 t + D s 1 ) sin qt e pt ( D 0 2 t s + D 1 2 t s - 1 + . . . + D s - 1 2 t + D s 2 ) sin qt . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . e pt ( D 0 n t s + D 1 n t s - 1 + . . . + D s - 1 n t + D s n ) sin qt righ ( ) ( ) .

б) Якщо p ± iq є коренем характеристичного рівняння кратності r , то частинний розв’язок має вигляд.

x 1 неод . . . . . . x n неод righ e pt ( C 0 1 t s + r + C 1 1 t s + r - 1 + . . . + C s + r 1 ) cos qt e pt ( C 0 2 t s + r + C 1 2 t s + r - 1 + . . . + C s + r 2 ) cos qt . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . e pt ( C 0 n t s + r + C 1 n t s + r - 1 + . . . + C s + r n ) cos qt righ ( ) + e pt ( D 0 1 t s + r + D 1 1 t s + r - 1 + . . . + D s + r 1 ) sin qt e pt ( D 0 2 t s + r + D 1 2 t s + r - 1 + . . . + D s + r 2 ) sin qt . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . e pt ( D 0 n t s + r + D 1 n t s + r - 1 + . . . + D s + r n ) sin qt righ ( ) ( ) .

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