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

№179.11.Ястребов

.pdf
Скачиваний:
96
Добавлен:
15.02.2015
Размер:
335.55 Кб
Скачать

31

Наложим на функции C1, C2 условие:

 

 

 

 

 

 

 

 

 

 

 

 

 

C y + C

y

2

= 0;

 

 

 

(23)

 

 

 

 

1

1

2

 

 

 

 

 

 

 

 

тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

yчн

= C1 y1

+ C2 y2

 

 

 

 

 

 

Дифференцируем повторно:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

′′

′ ′

 

′′

 

 

′ ′

 

′′ .

 

 

 

 

yчн = C1 y1

+ C1 y1 + C2 y2

+ C2 y2

 

 

 

 

Подставим

выражения

для

yчн, yчн, yчн

в исходное

уравнение и

 

 

 

 

′′

сгруппируем по отдельности слагаемые с C1 и с C2 :

 

 

 

 

′′

 

 

′′

 

 

 

 

 

′ ′

f ( x)

.

C1( y1 + py1 + qy1) + C2 ( y2 + py2

+ qy2 ) + C1 y1

+ C2 y2 =

 

Множители при C1 и C2 равны тождественно нулю, поскольку функции

y1 и y2 являются решениями однородного уравнения, так что

 

 

 

 

 

C1 y1

+ C2 y2 = f (x)

.

 

 

(24)

 

 

 

 

′ ′

 

 

 

 

 

 

Соотношения (23) и (24) дают систему двух уравнений относительно неиз-

вестных C, C:

1 2

 

 

 

 

 

+ C2y2

= 0,

C1y1

 

 

 

 

 

 

(25)

C

y

+ C

 

y

= f (x).

 

 

2

 

2

 

1

1

 

 

 

 

Решая эту систему, получим выражения для C, Cчерез уже известные

1 2

функции

y1, y2

, y1

, y2

C1(x), C2

(x)

на-

 

, после чего сами функции

 

 

ходятся интегрированием.

В качестве итога сформулируем алгоритм решения неоднородного уравнения методом вариации произвольных постоянных:

1. Решение соответствующего однородного уравнения, получение функций y1, y2 .

32

2. Запись общего решения соответствующего однородного уравнения в виде yоо = C1 y1(x) + C2 y2 ( x) .

3.Вычисление производных y, y.

12

4.Запись системы (25) для отыскания C, C.

12

5.Решение системы, получение функций C, C.

12

6.Нахождение каких-либо первообразных C1(x), C2 (x) интегри-

рованием функций

C1

, C2

 

.

7. Запись частного

 

решения неоднородного уравнения в виде

yчн = C1(x) y1(x) + C2 (x) y2 (x) .

8. Запись общего решения неоднородного уравнения: yон = yоо + yчн.

Пример. Решим методом вариации произвольных постоянных урав-

нение y′′ + 2 y′ = x. Следуя алгоритму, последовательно получаем:

1. y′′ + 2 y′ = 0; z2 + 2z = 0; z1 = −2, z2 = 0 ;

y1 = e2 x ; y2 = 1.

2.

y = C e2 x

+ C

2

.

 

 

 

 

оо

 

1

 

 

 

 

 

 

 

 

 

 

3.

y1 = −2e

2 x

; y2 = 0

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 x

 

 

 

 

1 = 0

 

 

 

4.

C1e

 

 

+ C2

 

 

 

 

 

 

 

 

 

 

 

 

+ 0 = x.

 

 

 

C(2e2 x )

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

5.

C1

= −

1

 

xe

2 x

;

C2

=

1

x

.

 

2

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6. Интегрируя по частям, имеем:

 

 

 

 

33

 

 

 

 

C (x) = − 1

xe2 xdx = −

1 1 xe2 x

1 e2 x

=

1

2

 

 

2

4

 

 

 

 

2

 

 

= −

1

xe

2 x

+

1

e

2 x

= e

2 x

1

x +

1

 

 

 

 

 

 

 

 

 

4

 

8

 

 

 

 

 

4

8

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C2 (x) =

2

xdx =

4

x .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7. yчн = e2 x

 

1

x +

1

 

 

 

 

1

x2

1

=

1

x2

1

x +

1

 

 

4

8

e2 x +

4

4

4

8

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8. yон = yоо + yчн = C1e2 x + C2 + 14 x2 14 x + 18 .

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

Это метод отыскания частного решения линейного неоднородного уравнения с постоянными коэффициентами

y′′ + py′ + qy = f (x)

в случае, когда правая часть f (x) имеет специальный вид, самая общая запись которого:

f (x) = eα x (P (x) cos β x + Q (x)sin β x),

(26)

 

m

 

n

 

 

 

 

где Pm (x), Qn (x) — многочлены степени m и n соответственно:

 

P (x) = p xm

+ p xm1

+ ... + p

 

 

x + p

 

m

0

1

m1

m

 

Q (x) = q xn + q xn1

+ ... + q

 

x + q .

 

n

0

1

n1

 

n

 

Укажем характерные частные случаи функции (26).

 

 

 

 

1. α = 0; β = 0 . В этом случае правая часть

f (x) = Pm (x) является

многочленом (поскольку e0 = 1; cos0 = 1, sin 0 = 0). Например:

34

1) y′′ + py′ + qy = 3x2 4

y′′ + py′ + qy = e0 x [(3x2 4) cos(0x) + 0sin(0x)]; 2) y′′ + py′ + qy = −3

y′′ + py′ + qy = e0 x [3cos(0x) + 0sin(0x)].

2. α = 0. В этом случае f (x) = Pm (x) cos β x + Qn ( x)sin β x .

Например:

y′′ + py′ + qy = 4 cos x;

y′′ + py′ + qy = −3cos x + 7x3 sin x.

3. β = 0; Pm (x) = 1. В этом случае f (x) = eα x . Например:

y′′ + py′ + qy = ex ; y′′ + py′ + qy = eπ x .

4. Qn (x) = 0. В этом случае f (x) = eα x Pm (x) cos β x . В част-

ности, если степень многочлена m = 0, то многочлен является постоянным числом: P0 (x) = p0 0

Например:

y′′ + py′ + qy = −8e9 x cos x;

y′′ + py′ + qy = e2 x (5x + 1) cos12x.

5. Pm (x) = 0 . В этом случае f (x) = eα xQn (x)sin β x . Напри-

мер:

y′′ + py′ + qy = sin 2x;

y′′ + py′ + qy = e2 x (5x + 1)sin12x.

35

Рассмотрим комплексное число r = α + β i , где α и β берутся из записи (26) правой части f (x) неоднородного уравнения. Корни характе-

ристического уравнения z2 + pz + q = 0 записываются в виде

 

z

 

= −

p

±

 

 

D

 

 

i ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1,2

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

где D — дискриминант уравнения.

 

 

 

 

 

Кратность k числа r

 

как корня

характеристического уравнения

может иметь одно из трех значений:

 

 

 

 

 

1) k = 0; это означает,

что

 

r не является корнем характеристического

уравнения;

 

 

 

 

 

 

 

 

 

 

 

2) k = 1; это означает,

что

r является простым корнем характеристиче-

ского уравнения; при этом дискриминант уравнения D отличен от нуля;

3) k = 2 , то есть r — кратный корень; последнее возможно в том и толь-

ко том случае, когда D = 0; корень r при этом вещественный, так как его

мнимая часть β =

 

 

D

 

 

= 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

На практике для определения кратности k нужно решить (квадрат-

ное) характеристическое уравнение и сравнить его корни с числом r .

 

Пусть l

— максимальная из степеней m и n многочленов Pm и Qn

в записи правой части (26).

 

 

 

 

 

 

Можно убедиться непосредственной подстановкой в дифференци-

альное уравнение, что функция

 

 

 

 

 

 

y

чн

= xk eα x ( A (x) cos β x + B (x)sin β x) ,

(27)

 

 

 

 

 

 

 

l

 

 

 

l

 

 

где Al и Bl — многочлены степени l :

 

 

 

A = a xl + a xl 1

+ a

2

xl2

+ ... + a

x + a ,

 

 

l

0

1

 

 

l1

l

 

36

Bl = b0 xl + b1xl1 + b2 xl2 + ... + bl 1x + bl ,

при некоторых значениях коэффициентов ai , bi этих многочленов явля-

ется частным решением неоднородного уравнения. Для отыскания этих (неопределенных вначале) коэффициентов у функции (27) вычисляются

первая и вторая производные; затем

yчн, yчн,

yчн

′′

подставляются в ис-

ходное уравнение, после чего приравниваются коэффициенты при одинаковых функциях-слагаемых обеих частей равенства. Это дает систему уравнений для отыскания коэффициентов ai , bi .

В итоге отыскание частного решения проводится по следующему ал-

горитму:

 

 

 

1.

Нахождение корней характеристического уравнения.

2.

Определение величин α , β ,

m, n, l .

3.

Выявление кратности k числа α + β i как корня характеристиче-

ского уравнения.

 

 

 

4.

Запись частного решения

yчн в виде (27) с неопределенными ко-

эффициентами многочленов Al , Bl .

 

5.

Вычисление производных

yчн, yчн

′′ .

6.

Подстановка

yчн, yчн, yчн

в дифференциальное уравнение.

′′

7.Приравнивание коэффициентов при одинаковых функциях в обеих частях равенства.

8.Решение получившейся системы уравнений, получение коэффици-

ентов многочленов Al , Bl .

9. Запись частного решения с найденными коэффициентами.

После этого можно записать общее решение в виде

yон = yоо + yчн

37

Примеры. 1. Рассмотрим уравнение y′′ + 2 y′ = 20x + 2 .

Выполняем последовательно инструкции алгоритма:

1. z2 + 2z = 0; z1 = −2; z2 = 0 . Общее решение соответствующего однородного уравнения: yоо = C1e2 x + C2e0 x = C1e2 x + C2 .

2. Находим параметры правой части:

20x + 2 = e0 x ((20x + 2) cos(0x) + 0sin(0x)); α = 0; β = 0; m = 1; n = 0; l = 1.

3.α + β i = 0 + 0i = 0; кратность k = 1.

4.yчн = xA1(x) = x(a0 x + a1) = a0 x2 + a1x .

5.

′′

.

 

yчн = 2a0 x + a1; yчн = 2a0

 

6.Подставляем в уравнение: 2a0 + 2(2a0 x + a1) = 20x + 2 .

7.4a0 x + 2(a0 + a1) = 20x + 2;

4a0

= 20

; a0

= 5; a1 = −4 .

8.

 

2(a0 + a1) = 2

 

 

9. yчн = x(5x 4) = 5x2 4x.

Общее решение неоднородного уравнения:

yоо = C1e2 x + C2 + 5x2 4x .

2. Рассмотрим уравнение

y′′ − 10 y′ + 21y = e6 x (3x2 + 4x + 2).

Выполняем последовательно инструкции алгоритма:

38

1. z2 10z + 21 = 0; z1 = 3; z2 = 7 . Общее решение соответствую-

щего однородного уравнения: yоо = C1e3x + C2e7 x . 2. Находим параметры правой части:

e6x (3x2 + 4x + 2) = e6x ((3x2 + 4x + 2)cos(0x) + 0sin(0x));

α = 6; β = 0; m = 2; n = 0; l = 2.

3.α + β i = 6 + 0i = 6; кратность k = 0.

4.yчн = x0e6 x A2 (x) = e6 x (a0 x2 + a1x + a2 );

5.yчн′ = 6e6 x (a0 x2 + a1x + a2 )+ e6 x (2a0 x + a1 ) =

=e6 x (6a0 x2 + (2a0 + 6a1)x + (a1 + 6a2 ));

′′

 

6 x

2

.

 

 

(36a0 x + (24a0 + 36a1)x + (2a0 + 12a1 + 36a2 ))

yчн = e

 

 

6. Подставляем в уравнение:

 

 

e6 x (

36a0 x2 + (24a0 + 36a1)x + (2a0 + 12a1 + 36a2 ))

10e6 x (6a0 x2 + (2a0 + 6a1)x + (a1 + 6a2 ))+ +21e6 x (a0 x2 + a1x + a2 ) = e6 x (3x2 + 4x + 2).

7. Приводим подобные члены и приравниваем коэффициенты при x2e6 x , при xe6 x и при e6 x :

3a0 = −3

4a0 3a1 = 4

2a0 + 2a1 3a2 = 2.

39

8.

Решаем систему: a0 = 1;

 

a1 = 0;

a2 = 0.

 

9.

yчн = e6 x (x2 + 0x + 0) = x2e6 x .

 

Общее решение неоднородного уравнения:

 

y

оо

= C e3x + C e7 x + x2e6 x .

 

 

 

1

2

 

 

 

 

 

 

 

3. Рассмотрим уравнение y′′ + 4 y = 2cos x 5sin x .

1.

z2 + 4 = 0;

z = ±

4 = ±

4i = 0 ± 2i .

Общее решение соот-

 

 

 

1,2

 

 

 

 

 

ветствующего однородного уравнения:

 

y

оо

= e0 x (C cos 2x + C

2

sin 2x) = C cos 2x + C sin 2x .

 

1

 

 

 

1

2

2. Находим параметры правой части:

2 cos x 5sin x = e0 x (P0 (x) cos x + Q0 (x)sin x); α = 0; β = 1; m = 0; n = 0; l = 0 .

3.α + β i = 0 + 1i ; кратность k = 0.

4.yчн = x0e0 x (a0 cos x + b0 cos x) = a0 cos x + b0 cos x ;

5.yчн′ = −a0 sin x + b0 cos x ;

′′

.

yчн = −a0 cos x b0 sin x

 

6. Подставляем в уравнение:

 

a0 cos x b0 sin x + 4(a0 cos x + b0 sin x) = 2 cos x 5sin x ; 3a0 cos x + 3b0 sin x = 2 cos x 5sin x .

7. Приравниваем коэффициенты при cos x и при sin x :

3a0 = 23b0 = −5.

8. a0 = 23 ; b0 = − 53 .

40

9. yчн = 23 cos x 53 sin x .

Общее решение неоднородного уравнения:

yоо = C1 cos 2x + C2 sin 2x + 23 cos x 53 sin x .