Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
высшая математика теоретический материал часть 3.pdf
Скачиваний:
102
Добавлен:
26.07.2016
Размер:
1.43 Mб
Скачать

c2

(x) = −

 

e x

 

dx = − ln(1 + e x ).

1

+ e

x

 

 

 

 

70. Запишем частное решение у*:

y* = x ln(1 + ex )ex ln(1 + ex ), y* = x (1 + ex )ln(1 + ex ).

80. Запишем ответ – общее решение уравнения: y = y + y*,

y = c1 + c2ex + x (1 + ex )ln(1 + ex ).

Пример 5.2.2. Найти общее решение дифференциального уравнения y''+ y = cos1 x .

 

x2

ln x

Ответ:

y =

 

 

 

 

 

 

2

 

 

 

3

 

 

 

x2 + c + c

2

x e2 x .

 

 

4

1

 

 

 

 

 

6. Линейные неоднородные уравнения ІІ порядка с постоянными коэффициентами. Метод Лагранжа неопределенных коэффициентов

Продолжаем рассматривать методы решения уравнения. y''+ py + gy = f (x)

Как выяснилось на предыдущем занятии, метод Эйлера вариации произвольных постоянных, с помощью которого отыскиваются частные решения у*, связан с интегрированием функций c1 '(x) и c2 '(x), что представляет определенные практические

трудности. Имеются случаи, когда частное решение у* можно найти проще, не прибегая к интегрированию. Речь пойдет о широко применяемых в науке дифференциальных уравнениях, у которых правая часть имеет вид: f (x) = edx [Px (x)cos βx + Qm (x)sin β (x)] , где Px (x), Qm (x) - заданные многочлены одной или разных степеней.

Поставим в соответствие уравнению (6.1) с правой частью (6.2) число α ± β i и назовем его основным параметром уравнения.

Сконструируем функцию вида y* = edx [M n (x)cos βx + Nn (x)sin βx]xr ,

гдеMn (x), Nn (x) - многочлены степени n = max{k, m}, записанные пока с неопределенными коэффициентами (отсюда название метода);

r - кратность корня характеристического уравнения, равного параметру α ± βi .

31

Таблица 4

f (x)) и соответствующие решения уравнения у*

 

 

 

 

 

 

Формы правой части

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Правая часть

 

Основной

Сравнение параметра

Конструкция

 

 

 

с корнями

 

 

N

 

 

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

уравнения f (x)

 

параметр α ± βi

характеристического

 

 

 

 

 

уравнения

у*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α = β = 0

 

0 не является корнем

B

 

 

 

 

 

1

А

 

 

0 однократный корень

Bx

 

 

 

 

 

α ± β i = 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 двукратный корень

Bx2

 

 

 

 

 

 

 

 

 

0 не является корнем

M n (x)

 

 

 

 

 

α = β = 0

 

 

 

 

 

 

 

 

 

2

Pn (x)

 

0 однократный корень

M n (x) x

 

 

α ± β i = 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 двукратный корень

M n (x) x 2

 

 

 

 

 

 

α не является корнем

Beαx

 

 

 

 

 

 

β = 0

 

 

 

 

 

 

 

 

 

 

3

Aeαx

 

 

α однократный

Beαx x

 

 

 

 

α ± β i

= α

 

 

 

 

 

 

корень

 

 

 

 

 

 

 

 

 

 

 

 

 

α двукратный корень

Beαx x2

 

 

 

 

 

 

 

α не является корнем

M n (x)eαx

 

 

P (x)eαx

 

β = 0

 

 

 

 

 

 

4

 

 

α однократный

M

n

(x)eαx x

 

 

n

 

α ± β i

= α

корень

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α двукратный корень

M n (x)eαx x 2

 

 

A cos βx +

 

α = 0;

 

± β i

не являются

C cos βx + D sin βx

5

 

 

корнями

 

 

+ B sin βx

 

β 0

 

 

 

 

 

 

 

 

 

 

 

 

± β i

- корни

C cos βx

+

D sin βx

x

 

 

 

 

 

 

 

(

 

 

)

 

 

Pk (x)cos βx +

 

 

 

± β i

не являются

M n (x)cos βx + Nn (x)sin βx

6

 

α = 0;

 

корнями

 

 

n = max{k;m}

 

+ Qm (x)sin βx

 

β 0

 

 

 

 

 

 

 

 

 

 

 

 

± β i

- корни

( Mn

(x)cos βx +

 

 

 

 

 

 

+ Nn

(x)sin βx )x

 

 

 

 

 

 

 

 

 

7

( Acos βx +

 

α ± β i

α ± β i

не являются

(C cos βx + D sin βx)eαx

+ B sin βx )eαx

 

корнями

 

 

 

 

 

 

 

 

 

 

 

 

α ± β i

- корни

C cos βx

+

D sin βx eαx x

 

 

 

 

 

 

(

 

 

)

 

 

( Pk (x)cos βx +

 

 

 

α ± β i

не являются

( M n cos βx +

 

8

 

α ± β i

корнями

 

 

+ N n sin βx )eαx

 

+ Qm (x)sin βx )eαx

 

 

 

 

 

 

 

 

 

 

 

α ± β i

- корни

(M n cos βx + N n sin βx)eαx x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Как видим, конструкция функции у* определяется как формой правой части уравнения – функцией f (x), так и видом левой его части – корнями характеристического уравнения. Доказано, что при соответствующем выборе значений коэффициентов для многочленов

Mn (x), Nn (x) функция у* является частным решением уравнения (6.1).

Втаблице 4 приведены различные формы правой части f (x) (частные случаи

α= 0, β = 0, α = β = 0 ) и соответствующие решения уравнения у*.

Пример 6.1. Для данных неоднородных уравнений подобрать частное решение у* в форме, соответствующей форме из таблицы 4.

32

Решения неоднородных уравнений примера 6.1

 

Правая часть

Основной

Сравнение параметра с

корнями

уравнения

параметр

характеристического

 

f(x)

α ± βi

 

уравнения

 

 

 

 

 

 

 

1

2

3

4

 

 

 

 

 

 

α = βi=0

0 не является корнем

1

А

0 однократный корень

α ± βi=0

 

 

0 двукратный корень

 

 

 

 

 

 

 

 

 

α = βi=0

0 не является корнем

2

Pn(x)

0 однократный

α ± βi=0

корень

 

 

 

 

 

0 двукратный корень

 

 

 

 

 

 

β=0

α не является корнем

3

Aeαx

α – однократный

α ± βi

 

 

корень

 

 

 

α – двукратный корень

 

 

 

 

 

 

β=0

α не является корнем

4

Pn(x) eαx

α – однократный

α ± βi

 

 

корень

 

 

 

α – двукратный корень

 

 

 

 

 

 

α=0

± βi не являются

5

Acosβx + Bsinβx

корнями

α ± βi

 

 

± βi – корни

 

 

 

 

Px(x)cosβx +

α=0

± βi не являются

6

корнями

Qm(x)sinβx

α ± βi

 

± βi – корни

 

 

 

 

 

 

 

 

(Acosβx + B

 

α ± βi не являются

7

α ± βi

корнями

sinβx)eαx

α ± βi – корни

 

 

 

 

(Px(x)cosβx +

 

α ± βi не являются

8

α ± βi

корнями

Qm(x)sinβx) eαx

α ± βi – корни

 

 

 

 

 

 

 

Таблица 5

Конструкция частного решения

у*

5

B

Bx Bx2

Mn(x)

Mn(x) · x

Mn(x) · x2

Beαx

Beαx · x

Beαx · x2

Mn(x)eαx

Mn(x)eαx · x

Mn(x)eαx · x2

Ccosβx + Dsinβx

(Ccosβx + Dsinβx) · x

Mn(x)cosβx + Nn(x)sinβx)

(Mn(x)cosβx + Nn(x)sinβx)

· x

n = max{k, m}

(Ccosβx + Dsinβx)eαx

(Ccosβx + Dsinβx)eαx· x

(Mn(x)cosβx + Nn(x)sinβx)eαx

(Mn(x)cosβx + Nn(x)sinβx)eαx· x

n = max{k, m}

Пример 6.2. Подобрать для данных неоднородных уравнений частное решение у* в форме, соответствующей форме из таблицы 4.

а) y''+2y = 2x ;

г)

y''+4y = 3sin 2x ;

б ) y + 2y'= x3 + x ;

 

д) y''8y'+16y = xe4 x ;

в) y''+2y'= 3e2 x ;

е)

y''6y'+34y = e3x cos5x ..

Ответ: а) y* = Ax + B ;

 

 

33

б) y* = Ax4 + Bx3 + Cx2 + Dx ; в) y* = Axe2 x ;

г) y* = (Acos 2x + B sin 2x)x ;

д) y* = (Ax3 + Bx2 )e4 x ;

е) y* = (Ax cos5x + Bx sin 5x)e3x .

Пример 6.3. Найти общее решение дифференциального уравнения y′′ + 4 y' = x3 + 1.

10. Определим тип уравнения:

y′′ + 4 y' = x3 + 1 - линейное, неоднородное, ІІ порядка, с постоянными коэффициентами, со специальной правой частью.

20. Запишем формулу общего решения: y = y + y * .

30. Найдём общее решение однородного уравнения - y : y''+4 y' = 0,

k 2 + 4k = 0, k1 = 0, k2 = −4,

y = c1 + c2e4 x

40. Проведём анализ правой части уравнения: x3 + 1 = eox ((x3 + 1)cos ox + o sin 0x),

α = 0, β = 0,

.

f (x) = P3 (x).

 

50. Вычислим основной параметр уравнения:

α ± βi = 0 .

60. Определим параметр r :

Основной параметр α ± βi = 0 является однократным корнем характеристического уравнения, следовательно r = 1.

70. Сконструируем частное решение – у*: y* = M 3 (x)x = (Ax3 + Bx2 + CX + D)x .

80. Вычислим коэффициенты функции у*: 8.1. Найдём производные от функции у*: y* = Ax4 + Bx3 + Cx2 + Dx,

(y *)'= 4Ax3 + 3Bx2 + 2Cx + D, (y *)''= 12Ax2 + 6Bx + 2C.

8.2. Поставим функцию у* и её производные в данное уравнение:

16 Ax3 + (12A + 12B)x 2 + (6B + 8C)x(2C + 4D) = x3 + 1 .

8.3. Приравняем коэффициенты при подобных членах левой и правой части равенства:

34

16 A = 1,

12A + 12B = 0,

6B + 8C = 0,2C + 4D = 1.

8.4. Решим систему:

A = 161 , B = − 161 , C = 143 , D = 12829 . 90. Запишем частное решение у*: y* = 161 x4 161 x3 + 643 x2 + 12829 x .

100. Запишем ответ – общее решение уравнения: y = c1 + c2e4x + 161 x4 161 x3 + 643 x2 + 12829 x .

Пример 6.4. Найти решение задачи Коши для дифференциального уравнения y''+ y = 8sin x , удовлетворяющее условиям y(0) = 1, y' (0) = 0 .

10. y''+ y = 8sin x - линейное, неоднородное, ІІ порядка, с постоянными коэффициентами, со специальной правой частью.

20. Запишем формулу общего решения y = y + y * .

30. Найдём общее решение однородного уравнения - y : y''+ y = 0,

k 2 + 1 = 0, k1,2 = ±i,

y = c1 cos x + c2 sin x .

40. Проведём анализ правой части уравнения: 8sin x = eox (o cos x + 8sin x),

α = 0, β = 1,

f (x) = p0 cos x + Q0 sin x.

50. Вычислим основной параметр уравнения:

α ± βi = ±i .

60. Определим параметр r :

Значения основного параметра ± i являются однократными корнями характеристического уравнениями, следовательно r = 1.

70. Сконструируем частное решение – у*: y* = (Acos x + B sin x)x .

80. Вычислим коэффициенты функции у*: 8.1. Найдём производные от функции у*: y* = Ax cos x + Bxsin x ;

(y *)'= (A + Bx)cos x + (B Ax)sin x;

(y *)''= (2B Ax)cos x + (2A 8x)sin x

35

8.2. Поставим функцию у* и её производные в данное уравнение: 2B cos x 2Asin x = 8sin x

8.3. Приравняем коэффициенты при подобных членах левой и правой части равенства:

2B = 0

2A = 8

8.4. Решим систему:

A = −4, B = 0

90. Запишем частное решение у*: y* = −4x cos x .

100. Запишем общее решение уравнения: y = (c1 4x)cos x + c2 sin x .

110. Найдём значения произвольных постоянных c1 и c2 :

y = (c1 4x)cos x + c2 sin x,

y'= (c2 4)cos x + (4x c1 )sin x

При х=0, у=1, у’=0 имеем c1 = 1,c2 = 4 .

120. Запишем ответ – частное решение уравнения: y = (1 4x)cos x + 4sin x .

Пример 6.5. Найти общее решение дифференциального уравнения y''2y'+2y = 2x .

Ответ: y = (c1 cos x + c2 sin x)ex + x + 1.

Пример 6.6. Найти решение задачи Коши для дифференциального уравнения y''+4y'+4y = xe2 x ,

удовлетворяющее условиям y(0) = 0, y'(0) = 1.

 

 

 

 

 

 

 

1

 

3

 

2 x

 

Ответ: y = x +

 

x

 

e

 

.

6

 

 

 

 

 

 

 

 

36