Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Басов4

.pdf
Скачиваний:
9
Добавлен:
16.04.2015
Размер:
330.45 Кб
Скачать

В. В. Басов Курс лекций по ОДУ

Разделив теперь это уравнение на W (x); получим ЛОУ (4.3) порядка n; в котором, в частности, p1(x) = r1(x)=W (x):

Остается проверить, что '1(x); : : : ; 'n(x) – это его решения. Это очевидно, так как при подстановке любого 'j(x) в левую

часть уравнения (4.7) она будет обращаться в нуль, поскольку в определителе появятся два одинаковых столбца.

60: Формула Лиувилля.

Оказывается, что вычислить определитель Вронского удается, даже не зная ФСР, по которой он строится. Достаточно только задать начальные данные для решений, входящих в ФСР, или W (x0):

Для вывода требуемой формулы напомним правило дифферен-

цирования определителя. Если матрица (x) = f'ij(x)gi;jn

=1; то

 

 

 

 

 

'110

 

: : : '10 n

+

'11

 

: : : '1n

+: : :+

'11 : : : '1n

:

 

(x) 0 = '21 : : : '2n

'210

 

: : : '20 n

'21 : : : '2n

j

 

j

 

..

 

 

 

..

 

 

 

 

..

 

 

 

 

 

 

 

..

 

 

 

 

 

 

 

 

 

..

 

 

 

 

 

 

..

 

 

 

.

 

 

: : : .

 

 

 

 

.

 

 

: : : .

 

 

 

 

 

 

 

 

 

.

 

 

 

 

: : : .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

'n1

 

: : : 'nn

 

'n1

 

: : : 'nn

 

 

 

 

 

 

 

 

'n0

1

 

: : : 'nn0

 

 

Поэтому

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

'1

 

: : :

 

 

'n

 

0

 

 

 

 

 

'10

 

: : :

 

 

 

'n0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

'1

 

: : :

 

 

'n

 

 

 

 

 

 

 

 

'10

 

: : :

 

 

 

'n0

 

 

 

 

 

 

 

 

W 0(x) = ...

 

: : :

 

 

 

...

 

 

 

 

=

 

 

 

 

 

..

 

: : :

 

 

 

 

..

 

 

 

+ : : :

 

 

 

 

 

 

 

'

(n

 

2)

: : : '

(n

 

2)

 

 

 

 

 

 

'

(n. 2)

: : : '

(n.

2)

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

n

 

 

 

 

 

 

 

 

 

1

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(n

 

1)

 

 

 

(n

 

1)

 

 

 

 

 

 

 

 

(n 1)

 

 

 

 

 

(n

 

1)

 

 

 

 

 

 

 

 

 

 

 

 

 

'1

 

: : : 'n

 

 

 

 

 

 

 

 

 

'1

 

 

: : : 'n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

'1

 

 

 

: : :

'n

 

 

 

 

 

 

 

 

 

'1

 

 

 

 

: : :

 

'n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

'10

 

 

 

: : :

'n0

 

 

 

 

 

 

 

 

'10

 

 

 

 

: : :

 

'n0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

: : : +

..

 

 

 

 

: : :

..

 

 

+

..

 

 

 

 

: : :

 

..

 

 

 

 

= r

(x);

 

 

 

 

 

'

(n.

 

1)

 

: : : '

(n.

1)

 

 

 

 

 

'

(n.

 

2)

 

: : : '

(n.

2)

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(n

 

1)

 

 

 

(n

1)

 

 

 

 

 

 

(n)

 

 

: : :

 

(n)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

'1

 

 

 

: : : 'n

 

 

 

 

 

 

 

'1

 

 

 

 

'n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

так как

в сумме все определители,

кроме последнего,

 

имеют по две

одинаковые

строки.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В итоге получили ЛОУ первого порядка W 0(x) = p1(x)W (x):

Разделяя

переменные и интегрируя по s от x

 

 

до x

(x

; x

2

(a; b));

 

 

x

 

 

 

 

 

 

 

x dW s

 

 

 

 

 

W x

0

 

 

 

 

 

 

 

 

 

0

 

 

 

имеем: Zx0

p1(s) ds = Zx0

 

( )

 

 

= ln

( )

;

 

откуда

 

 

 

 

 

 

W (s)

 

 

W (x0)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W (x) = W (x0) exp Zx0

p1(s) ds

формула Лиувилля.

 

11

§ 3. ЛИНЕЙНОЕ НЕОДНОРОДНОЕ УРАВНЕНИЕ

10: Структура общего решения ЛНУ.

Рассмотрим линейное неоднородное уравнение (4.1) порядка n : y(n) + p1(x)y(n 1) + : : : + pn 1(x)y0 + pn(x)y = q(x) или Ly = q:

Пусть y = (x) решение ЛНУ (4.1) на (a; b); т. е. L q: Используя его, сведем (4.1) к линейному однородному уравнению (4.3).

Для этого сделаем "сдвигающую" замену y = z + (x):

Имеем: L(z + ) = q , Lz + L = q , Lz = 0:

Пусть z = '(x; C1; : : : ; Cn) – общее решение ЛОУ Lz = 0: Тогда, подставляя его в замену, устанавливаем, что функция

y = '(x; C1; : : : ; Cn) + (x)

(4:8)

является общим решением ЛНУ (4.1).

Иными словами, общее решение ЛНУ есть сумма общего решения ЛОУ и частного решения ЛНУ.

20: Метод вариации произвольной постоянной.

Теорема (о нахождении частного решения ЛНУ). Пусть набор '1(x); : : : ; 'n(x) является фундаментальной системой решений ЛОУ (4.3) на (a; b): Тогда частное решение y = (x) ЛНУ (4.1) может быть найдено в виде квадратур от '1(x); : : : ; 'n(x); коэффициентов p1(x); : : : ; pn(x) и неоднородности q(x):

Д о к а з а т е л ь с т в о . Частное решение будем искать в виде

(x) = C1(x)'1(x) + : : : + Cn(x)'n(x);

(4:9)

где функции C1(x); : : : ; Cn(x) 2 C1((a; b)); т. е. варьировать произвольные константы из формулы общего решения '(x; C1; : : : ; Cn) =

C1'1(x) + : : : + Cn'n(x) ЛОУ (4.3).

Подставляя решение y = (x) в виде (4.9) в ЛНУ (4.1), получаем только одну связь на n неизвестных функций C1(x); : : : ; Cn(x): Остается возможность разумным образом наложить еще (n 1) -у связь на эти функции, рассчитывая затем их однозначно определить из появившейся системы из n уравнений.

Идея выбора связей заключается в недопущении роста числа слагаемых в ходе нахождения производных функции (x); необходимых для подстановки их в ЛНУ.

12

для 8 x0; x 2 (a; b):

В. В. Басов

 

 

 

 

 

 

 

 

 

 

 

Курс лекций по ОДУ

Итак,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(x) = C1'1 + : : : + Cn'n;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

}|

 

 

 

 

{

 

 

 

 

 

0(x) = C1'10 + : : : + Cn'n0 + C10 '1 + : : : + Cn0 'n;

 

 

 

 

 

 

 

 

 

 

 

 

=0

 

 

 

 

 

 

 

 

 

 

 

...

 

z

 

 

 

 

}|

 

 

{

 

 

 

 

 

 

1

1

 

 

n

 

 

 

 

 

00(x) = C1'100 + : : : + Cn'n00 + C0

'0 + : : : + C0

 

'n0 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=0

 

 

 

 

 

( ) =

1 1

+ : : : + Cn'n

 

 

+ z

 

 

 

 

 

 

}|

 

 

{

 

 

 

 

10

 

 

1

+

+ n0

n

(n

 

1) x

C '(n 1)

(n 1)

C

'(n 2)

 

 

 

: : : C

'(n 2);

 

(n)(x) = C1'(n)

+ : : : + Cn'n(n)

 

+ C0

'(n 1)

+ : : : + C0

'n(n 1);

 

 

 

1

 

 

 

 

 

 

1

 

 

1

 

 

 

 

 

n

 

 

т. е. в ходе дифференцирования

на

C10 ; : : : ; Cn0

была наложена

ровно (n 1) -а связь.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теперь для подстановки решения в уравнение домножим левую

иправую части первого равенства на pn(x); второго – на pn 1(x)

итак далее, предпоследнего – на p1(x) и последнего – на единицу,

после чего все равенства сложим и приравняем к q(x): В результате получим: L = C1L'1 + : : : + CnL'n + C10 '(1n 1) + : : : + Cn0 '(nn 1) = q:

Но L'j 0 (j = 1; n); поскольку 'j(x) – решения ЛОУ (4.3). Собрав вместе все связи на C10 ; : : : ; Cn0 ; получим линейную неоднородную алгебраическую систему n уравнений с n неизвестными

C

0

'1(x) + : : : + C0 'n(x) = 0

 

 

8 ...

1

(n

2)

n

(n

2)

(x) = 0

:

(4:10)

>C10

'1

 

(x) + : : : + Cn0

'n

 

 

 

>

 

 

 

 

 

 

 

 

 

<

 

(n

1)

(x) + : : : + Cn0

(n

1)

(x) = q(x)

 

 

>C10

'1

 

'n

 

 

 

>

 

 

 

 

 

 

 

 

 

:

 

 

 

 

 

 

 

 

 

Определитель ее матрицы (x); выписанный в (4.5), – это ОВ W (x); который отличен от нуля при 8 x 2 (a; b); так как по условию '1(x); : : : ; 'n(x) является фундаментальной системой решений.

Следовательно система (4.10) по формуле Крамера имеет единственное решение Cj0 = Anj(x) q(x)=W (x) (j = 1; n); где Anj – алгебраическое дополнение n -й строки и j -о столбца матрицы :

 

n

x Anj(s)

 

X

Согласно (4.9) частное решение (x) =

j=1 'j(x) Zx0

 

q(s) ds

W (s)

13

§ 3. ЛИНЕЙНОЕ УРАВНЕНИЕ С ПОСТОЯННЫМИ КОЭФФИЦИЕНТАМИ

10: ЛОУ с постоянными коэффициентами.

Пусть в ЛОУ (4.3) коэффициенты p1(x); : : : ; pn(x) постоянны, тогда оно имеет вид

Lcy = y(n)+a1y(n 1)+: : :+an 1y0+any = 0 (a1; : : : ; an 2 R1): (4:3c)

При этом, очевидно, что любое решение уравнения (4:3c) определено на всей вещественной оси.

Общее решение ЛОУ (4:3c) с постоянными коэффициентами можно всегда найти в явном виде, в отличие ЛОУ (4.3) с произвольными непрерывными коэффициентами, структура общего решения которых исследована в § 2, ФСР существует, но отсутствуют способы гарантированного нахождения в явном виде решений из ФСР.

Для того чтобы решить ЛОУ (4:3c); предварительно потребуются

два технических результата.

 

Утверждение 1. Функции xk1 e 1x; : : : ; xkn e nx; где kj 2 Z+;

(kj; j) 6= (kl; l) (j; l =

 

 

j 6= l); линейно независимы на R1:

1; n;

Д о к а з а т е л ь с т в о .

Следуя определению, покажем, что

xk1 e 1x + : : : + xkn e nx R1 0 верно только при ; : : : ; = 0:

1 n 1 n

Приводя в этом тождестве подобные члены при одинаковых показателях экспоненты, получаем

 

m

 

 

где 1; : : : ; m

Xs=1 P (s)(x)e sx 0

(1 m n);

(4(:s11))

различны и любой коэффициент многочленов P –

какое-либо j

из набора 1; : : : ; n;

так как пары (kj; j) различны.

Докажем методом математической индукции по m; что в тождестве (4.11) многочлены P (1)(x); : : : ; P (m)(x) 0; а значит, равны нулю все их коэффициенты, и набор 1; : : : ; n тривиален.

При m = 1 в (4.11) имеется одно слагаемое: P

(1)

 

x

0;

 

(x)e 1

поэтому P (1) 0:

 

 

 

 

 

 

 

 

 

 

 

 

R

1

 

Предположим, что из

(4.11) следует, что P (1)(x); : : : ; P (m)(x)

 

0:

 

 

 

m+1

 

(s)

x

0:

 

 

 

 

Рассмотрим тождество

 

s=1

P

 

(x)e s

 

 

 

 

 

 

После деления на

e

 

 

оно равносильно тождеству

 

 

 

 

 

 

m+1

P

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

( s m

+1 6= 0):

 

 

Xs=1 P (s)(x)e( s m+1)x + P (m+1)(x) 0

 

 

14

В. В. Басов Курс лекций по ОДУ

Пусть многочлен P (m+1)(x) имеет степень m: Продифференцировав последнее тождество (m + 2) раза, имеем:

sm=1 Q(s)(x)e( s m

+1)x 0; где для 8 s =

 

многочлен Q(s)(x)

1; m

P

 

 

 

имеют ту же степень, что и P (s)(x); поскольку в результате дифференцирований коэффициент при старшей степени Q(s) равен коэф-

фициенту при старшей степени P (s) (какому-то js ); умноженному

на ( s m+1)m+2:

По индукционному предположению все Q(s) 0; следовательно все P (s) 0; а тогда и все P (m+1) 0:

При n = 1 ЛОУ (4:3c) имеет вид y0 + a1y = 0 и общее решение y = e a1x: Поэтому, чтобы выяснить какие решения может иметь уравнение (4:3c) при n 2; подставим в него функцию y = e x с пока что произвольным показателем :

Чтобы e x оказалось решением, должно выполняться тождество

Lce x = ( n + a1 n 1 + : : : + an 1 + an)e x 0

(4:12):

Df. Функция g( ) = n +a1 n 1 +: : :+an 1 +an называется характеристическим многочленом ЛОУ (4:3c); а его нули или корни характеристического уравнения g( ) = 0 называются характеристическими числами.

По основной теореме высшей алгебры характеристическое уравнение имеет n корней 1; : : : ; n; причем комплексные корни могут появляться только комплексно сопряженными парами из-за предположения о вещественности коэффициентов уравнения.

Благодаря тождеству (4.12) ясно, что функции e 1x; : : : ; e nx являются решениями ЛОУ (4:3c) тогда и только тогда, когда показатели 1; : : : ; n – характеристические числа.

Кроме того, если характеристические показатели различны (у них всех кратность – единица), то набор e 1x; : : : ; e nx является ФСР, правда, не обязательно вещественной ЛОУ (4:3c):

А что делать, если кратность k0 какого-либо характеристического числа, обозначим его 0; окажется больше единицы? Как тогда набрать n линейно независимых решений?

15

Утверждение 2. Пусть 0 – характеристическое число кратности k0; тогда ЛОУ (4:3c) имеет следующие линейно независимые решения: e 0x; xe 0x; : : : ; xk0 1e 0x:

Д о к а з а т е л ь с т в о . Согласно (4.12) Lce x = g( )e x: Продифференцируем это равенство k раз по ; тогда в левой

части с учетом линейности оператора Lc получим: dk(Lce x)=d k = Lc(dke x)=d k) = Lc(xke x); а в правой части k -я производная про-

изведения дает: dk(g( )e x)=d k =

k

C g( )( )xk e x:

 

В итоге

 

k

P =0

k

 

 

(k0)

 

 

 

 

 

 

X(k0 1)

 

 

 

 

 

 

 

 

 

 

Lc(xke x) =

Ck g( )( )xk e x:

 

(4:13)

 

 

 

=0

 

 

 

 

 

 

 

 

 

Но

g(

) = g0(

) = : : : = g

(

) = 0;

а

g

 

(

) = 0;

так как

0

0

 

0

 

 

 

 

0

6

корень 0 имеет кратность k0:

Поэтому при = 0 и k = 0; k0 1 в (4.13) Lc(xke 0x) = 0: Теорема (о ФСР ЛОУ с постоянными коэффициентами). Пусть

характеристическое уравнение ЛОУ (4:3c) имеет корни j

крат-

ности kj (j =

 

): Тогда ФСР уравнения (4:3c) имеет вид:

1; m

e 1x; xe 1x; : : : ; xk1 1e 1x; : : : ; e mx; xe mx; : : : ; xkm 1e mx:

(4:14)

Д о к а з а т е л ь с т в о . Предложенный набор содержит n функций (k1 + : : : + kn = n); они линейно независимы на R1 по утверждению 1 и являются решениями (4:3c) по утверждению 2.

Замечание 2. Если ЛОУ (4:3c) имеет вещественные коэффициенты, то комплексные решения из ФСР (4.14) можно и нужно овеществить, как это сделано в лемме об овеществлении. В частности, если = + i ; то Re (xke x) = xk cos x; Im (xke x) = xk sin x:

Пример 1. Рассмотрим ЛОУ второго порядка

y00 + 2y = 0 ( > 0);

описывающее свободные колебания математического маятника без трения, так как отсутствует слагаемое с y0:

Характеристический многочлен этого уравнения g( ) = 2 + 2; характеристические числа 1;2 = i ; ФСР образуют ei x; e i x; а вещественную ФСР – cos x; sin x:

В результате общее вещественное решение имеет вид: y = C1 cos x + C2 sin x (C1; C2 2 R1):

16

В. В. Басов

Курс лекций по ОДУ

20: Метод неопределенных коэффициентов

 

для ЛНУ с постоянными коэффициентами.

 

Рассмотрим ЛНУ с постоянными коэффициентами

 

Lcy = y(n) + a1y(n 1) + : : : + an 1y0 + any = q(x);

(4:1c)

вкотором a1; : : : ; an 2 R1; а неоднородность q(x) 2 C((a; b)): Общее решение соответствующего (4:1c) ЛОУ (4:3c) известно –

это линейная комбинация решений из ФСР (4.14). Поэтому частное решение ЛНУ (4:1c) всегда может быть найдено в квадратурах методом вариации произвольной постоянной и явный вид общего решения ЛНУ (4:1c) даст формула (4.8).

Однако метод вариации весьма громоздок. Согласно (4.9) частное решение приходится раскладывать в сумму из n произведений и в каждое произведение, помимо решений ЛОУ, входят функции Cj(x); для вычисления которых требуется решать систему (4.10) и затем n раз интегрировать. При этом на практике после приведения подобных членов в формуле (4.9) формула частного решения существенно упрощается.

Существует еще один метод нахождения частного решения ЛНУ (4:1c); не связанный со знанием ФСР ЛОУ и интегрированием, что является его несомненным достоинством. Это – метод неопределенных коэффициентов, при котором, зная только корни характеристического уравнения, удается указать структуру частного решения с точностью до тех самых неопределенных коэффициентов, задающих многочлены, входящие в искомое решение.

"Расплатой" за сравнительную простоту метода неопределенных коэффициентов является его не универсальность.

Метод можно применять только в случае, когда неоднородность q имеет специальный вид: q(x) = p(x)e x; где 2 C; а p(x) – многочлен, возможно, с комплексными коэффициентами.

Итак, рассмотрим ЛНУ с постоянными коэффициентами

Lcy = p(x)e x;

(4:15)

в котором j – нуль характеристического многочлена кратности kj

(j = 1; m); многочлен p(x) = Pl pjxj:

j=0

17

Df: Говорят, что в уравнении (4.15) имеет место резонанс порядка k; если = k кратности k (k = 1; m); и резонанс отсутствует (его порядок равен нулю), если 6= j для 8 j = 1; m:

Теорема (о построении решения ЛНУ методом неопределенных коэффициентов). Пусть из правой части (4.15) совпадает с корнем характеристического уравнения кратности k: Тогда существует и единственно частное решение ЛНУ (4.15) вида

(x) = xkr(x)e x;

где многочлен r(x) = Pl rsxs:

s=0

Д о к а з а т е л ь с т в о . Будем искать пока еще неопределенные коэффициенты rs многочлена r(x): Для этого подставим решение y = (x) в уравнение (4.15), получая тождество по x на R1

 

 

 

 

 

 

 

l

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

 

 

 

Lc

rsxk+se x

pjxje x:

(4:16)

 

 

 

 

 

 

 

s=0

 

 

j=0

 

 

 

Используя линейность Lc

и равенство (4.13), преобразуем левую

часть тождества (4.16): Lc(

 

l

 

P

l

 

 

s=0 rsxk+se x) =

s=0 rsLc(xk+se x) =

 

l

r

s

k+s C

 

g( )( )xk+s e x:

 

 

 

P

s=0

 

=0 k+s

 

 

P

 

 

 

 

 

 

 

P

 

 

 

 

 

 

 

 

 

Сокращая теперь в (4.16) обе части на e x и учитывая, что g( ) =

g0( ) = : : : = g(k 1)( ) = 0;

а

g(k)( ) = 0;

получаем тождество

 

 

 

 

 

 

6

l

 

k+s

 

 

 

 

 

 

l

 

 

X X

 

 

g( )

 

xk+s

 

X

xj;

 

r

s

C

 

 

p

:

s=0

 

=k

 

 

 

 

 

 

j=0

 

 

в котором, как легко заметить, степень переменной x в обеих частях не превосходит l:

Приравняем в (4.17) коэффициенты при различных степенях x; начиная со старшей степени l:

Пусть в левой части степень икса k + s = l или k = l s: Но l s; поэтому k ; а значит, k = ; откуда l = s:

В результате при xl в (4.17) имеем равенство rlCkk+l g(k)( ) = pl; из которого однозначно находим коэффициент rl:

Предположим теперь, что однозначно найдены коэффициенты rl 1; : : : ; ri+1 (0 i l 1); и будем искать коэффициент ri; стоящий при xi:

18

В. В. Басов Курс лекций по ОДУ

Итак, пусть в левой части (4.17) k+s = i или k = i s 0; а значит, индекс s может принимать значения i; i + 1; : : : ; l:

Если s = i; то = k и в левой части (4.17) при xl имеется только одно слагаемое с таким s – это riCkk+i g(k)( ): А если s > i; то в любое слагаемое множителем входит коэффициент rs; который по индукционному предположению уже однозначно найден.

Следовательно при xi из тождества (4.17) получаем равенство riCkk+i g(k)( )+h(ri+1; : : : ; rl) = pi с известной функцией h и из него однозначно находим коэффициент ri:

Приведем два соображения, которые полезны при практическом нахождении частного решения ЛНУ (4:1c):

Первое соображение обычно называют "принцип суперпозиции".

 

 

 

 

 

 

 

 

 

 

 

 

 

Он применяется к уравнениям (4:1c) вида

 

 

 

 

Lcy = q0(x) + q1(x) + : : : + qm(x)

 

(m 0)

(4:18)

и заключается в следующем: если

j

– частное решение уравнения

Lcy = qj(x) (j =

 

); то функция

 

 

0; m

(x) =

0(x) + : : : + m(x)

будет частным решением уравнения (4.18).

 

 

 

 

Действительно, Lc( ) = Lc(

 

m

j) =

m

 

m

 

 

j=0

j=0 Lc( j) =

j=0 qj:

 

Чтобы использовать этот

принцип, произвольную неоднородность

 

 

P

 

 

P

 

 

P

q(x) разбивают в сумму неоднородностей q0(x); q1(x); : : : ; qm(x); в которой для 8 i = 1; m функции qi(x) имею вид неоднородности из правой части уравнения (4.15), т. е. qi(x) = p(i)(x)e (i)x; и частное решение y = i(x) каждого из уравнений (4.15) ищется методом неопределенных коэффициентов, а частное решение уравнения Lcy = q0(x) ищется методом вариации произвольной постоянной. После этого все найденные решения складываются.

Разумеется, при m 1 возможно, что q0(x) 0:

Второе соображение связано с решением вещественных линейных

неоднородных уравнений (4:1c): Метод неопределенных коэффициентов в них применим для q(x) = p1(x)e x cos x + p2(x)e x sin x:

~

 

~

 

Введем = + i ;

q~(x) = (p1(x) ip2(x))e

 

:

Тогда q = Re q~ и, если y = ~(x) – частное (комплексное) решение уравнения Lcy = q;~ то функция y = Re ~ будет частным решением вещественного уравнения Lcy = p1(x)e x cos x + p2(x)e x sin x:

19

Пример 2. Рассмотрим ЛНУ второго порядка

 

y00 + 2y = b cos !x ( ; ! > 0; b 2 R1);

(4:19)

описывающее вынужденные колебания математического маятника под воздействием 2 =! -периодической возмущающей силы.

В примере 1 было найдено

2 = -периодическое общее решение

y = A sin( x + ) с A = (C2

+ C2)1=2; = arctg (C1=C2) ЛОУ

1

2

y00 + 2y = 0; имеющего характеристические числа 1;2 = i : Следуя изложенному выше соображениям, рассмотрим ЛНУ

y00 + 2y = bei!x;

(4:20)

вещественная часть неоднородности которого совпадает с правой частью уравнения (4.19).

Вид частного решения уравнения (4.20) зависит от величины !:

1) ! 6= – частота вынужденных колебаний не совпадает с частотой свободных колебаний, поэтому резонанс отсутствует.

Тогда частное решение уравнения (4.20) согласно теореме о построении решения ЛНУ методом неопределенных коэффициентов ищем в виде ~ = aei!x:

Подставляя ~ и ~00 = a!2ei!x в (4.20), имеем a!2 + a 2 = b; откуда ~(x) = b( 2 !2) 1ei!x:

Выделяя вещественную часть ~; находим общее решение (4.19)

y= A sin( x + ) + b( 2 !2) 1 cos !x – ограниченная функция x:

2)! = – частота вынужденных колебаний совпадает с частотой свободных колебаний, поэтому имеет место резонанс.

Тогда частное решение ЛНУ (4.20) надо искать в виде ~ = axei!x: Подставляя ~ и ~00 = (2ia! a!2)ei!x в (4.20), имеем 2ia! = b;

откуда ~(x) = b(2i!) 1xei!x:

Выделяя вещественную часть ~; находим общее решение (4.19) y = A sin( x + ) + b(2!2) 1x sin !x – не ограничено с ростом x:

20

Тут вы можете оставить комментарий к выбранному абзацу или сообщить об ошибке.

Оставленные комментарии видны всем.