Материал: Михайлов ТФКП Практикум 2013

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

 

 

 

 

Y ( p) =

1

 

; y (x) = e−x sinx.

 

 

 

 

p2 + 2 p + 2

 

 

 

 

1

1

Решение задачи Коши имеет вид:

 

y(x) = y

0

y

0

(x) + y y (x) = A e−x (cosx −sinx) + B e−x sinx.

 

 

1

1

 

 

Перейдем к нахождению решения задачи Коши с ненулевой правой частью f (x) и нулевыми начальными условиями, т.е.

a

0

y(n)

(x) + a y(n−1)

(x) +.... + a

n

y(x) = f (x) ;

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

y(0) = y′(0) =... = y(n−1) (0) = 0 .

Умножив обе части уравнения на e− px

и проинтегрировав по х

от 0 до +∞, перейдем к изображениям:

 

 

 

 

 

 

(a

0

pn + a pn−1 +... + a

n

) ·Y ( p)

= F ( p),

 

 

 

 

1

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

F( p)

 

 

 

 

 

 

 

 

Y ( p) =

 

 

 

 

 

,

 

 

 

 

 

a pn +a pn−1

+... + a

 

 

 

 

 

0

 

 

1

 

 

n

из этого выражения для Y ( p)

обратным преобразованием находим

искомую функцию

y( х) (обычно это делается с помощью выче-

тов).

Пример 15.2 ([4], № 50 (1)). Решить задачу Коши y′′+ y = x3 +6x, y(0) = y′(0) = 0.

Решение. Сразу переходя к изображениям, получим

 

 

 

p2Y ( p)

+ Y ( p) =

6

+

 

6

 

;

 

 

 

p4

p2

 

 

 

 

 

 

 

 

 

 

Y ( p) = 6

 

1

 

 

1

 

 

 

 

6

 

; y ( x) = x3.

 

 

+

 

 

 

=

 

 

 

( p2 +1)

p2

 

 

 

p4

 

p4

 

( p2 +1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Возвращаясь к решению задачи Коши для однородного уравнения, выделим из фундаментальной системы решений ψk (x) реше-

ние с номером k = n −1 , т.е. у которого все начальные данные рав-

ны 0, кроме ψ(nn−−11) (0) :

 

 

 

 

ψ

n−1

(0) = ψ'

(0) = ... = ψ(n−2)

(0) = 0;

ψ( n−1)

(0) =1 .

 

n−1

n−1

 

n−1

 

В соответствии с полученным выражением для yk (x) получаем

231

Yn−1 ( p) =

 

a0

 

 

,

a pn + a pn−1

+... + a

 

0

1

 

n

что позволяет решить задачу Коши для неоднородного уравнения, но с нулевыми начальными условиями с помощью этой функции Yn−1 ( p). Действительно, так как

Y ( p)

=

 

F ( p)

 

=

1

·F ( p)·Y

( p),

a pn +a pn−1 +... +a

 

a

 

 

n

n−1

 

 

 

0

1

 

0

 

 

то по теореме о свертке

 

 

 

 

 

 

 

 

 

y(x) =

1

∫x f (τ)·yn−1 ( x − τ)dτ,

 

 

 

a

 

 

 

 

0

0

 

 

 

 

 

Обычно, функцию ψn−1(x) обозначают как g(x) и называют

функцией источника для уравнения

a0 y(n) (x) + a1 y(n−1) (x) + ... + an y(x) = 0.

Таким образом, решение задачи Коши для неоднородного уравнения с нулевыми начальными данными можно получить с помощью функции источника g(x) = ψn−1 (x) в виде интеграла

y(x) = 1 ∫x g(x − τ)· f (τ)dτ ,

a0 0

называемого интегралом Дюамеля. Решим эту же задачу:

y′′+ y = x3 +6x, y(0) = y′(0) = 0

с помощью функции источника. Найдем функцию источника g(x), которая является решением задачи: g ''+ g = 0; g (0) = 0; g '(0) =1.

Её изображение

G( p) =

 

1

 

 

g(x) = sinx .

p2 +

1

 

 

Искомое решение

 

 

 

 

 

y(x) = ∫x (t3 +6t)·sin ( x −t )dt = x3.

0

 

 

 

 

 

Пример 15.3 ([4], №49 (8)).

Решить задачу Коши

y ''+4y = 2cos 2x ;

 

y(0) = 0; y '(0) = 4 .

232

Решение. Переходим к уравнению для изображений:

 

 

 

 

 

 

p2 [Y ( p)

−

4

 

] + 4Y ( p)

= 2

 

 

 

p

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

p2

p2

+ 4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Y ( p) =

 

 

 

 

2 p

 

 

 

+

 

 

 

 

4

 

=

−

 

1

 

 

 

'

+

 

 

4

;

 

 

 

 

( p

2

 

2

 

 

 

p

2

+ 4

 

 

2

 

 

 

 

p

2

 

 

 

 

 

 

 

 

 

+ 4)

 

 

 

 

 

 

 

 

p

 

+ 4

 

 

 

 

 

+ 4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y(x)

=

 

1 x sin 2x + 2sin 2x .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 15.4 ([4], №50 (7)). Решить задачу Коши

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y ''+9 y = f (x); y(0)= y '(0)= 0.

 

 

 

 

 

 

 

 

Решение.

 

p2 Y ( p) +9Y ( p) = F( p);

Y ( p) =

 

F( p)

. По теореме о

 

 

p2 +9

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

свертке получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y(x) = 1

∫x sin 3τ

f (x − τ)dτ .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 15.5 ([4], №50 (9)). Решить задачу Коши

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y ''+ 4 y '+ 4 y = f (x) .

 

 

 

 

 

 

 

 

 

 

 

Решение.

p2 Y ( p) + 4 p Y ( p) + 4Y ( p) = F( p);

Y ( p) =

 

F( p)

 

;

 

( p + 2)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

так как

1

 

 

 

 

 

 

 

1

 

 

 

'

 

 

то

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= −

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( p + 2)

2

p +2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y(x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xe

−2 x

= ∫0τe

−2τ

·f (x −τ)dτ .

 

 

 

 

 

 

( p +2)2

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 15.6 ([4], № 51 (2)). Найти общее решение уравнения

 

 

 

 

 

 

 

y '''−6 y ''+11y '−6 y

= 12x2e3x − e2 x .

 

 

 

 

 

 

Решение. Найдем общее решение однородного уравнения: y '''−6 y ''+11y '−6 y = 0 .

Перейдем к изображениям:

p3[Y ( p) − Cp1 − Cp22 − Cp33 ] −6 p2[Y ( p) − Cp1 − Cp22 ] +

233

 

 

 

 

 

 

+11p[Y ( p) − C1

] −6Y ( p) = 0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C p2

+(C

2

−6C ) p +C −6C

2

+11C

 

 

 

 

Y ( p) =

1

 

 

 

 

 

 

1

 

 

 

3

 

 

 

 

1

=

 

 

 

 

 

 

 

 

p3 −6 p2 +11p −6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

C p2

+(C

2

−6C ) p +C −6C

2

+11C

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

1

 

3

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

( p −3)( p −2)( p −1)

 

 

 

 

 

 

 

 

 

 

y

 

 

(x)

= ex

(C +C −6C +C −6C +

11C ) −

 

 

 

однор

 

 

 

2

 

 

1

 

2

 

 

1

 

3

 

 

2

 

 

1

 

 

 

 

−e2 x (4C +2C

2

−12C +C −6C +11C ) +

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

1

 

3

 

 

2

 

1

 

 

 

 

 

+e3x (9C +3C −

18C +C −

6C +11C ) =

 

 

 

 

 

2

 

 

1

 

 

2

 

 

1

3

 

 

2

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=ex

(6C −5C

2

+C ) −e2 x (3C −4C

2

+C )+e3x

(2C −3C

2

+C )

2

1

 

 

 

3

 

 

 

 

 

 

1

 

 

 

3

 

2

 

 

1

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

однор

(x) =Сех

+С

е2 х +С е3х .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2

 

 

 

 

3

 

 

 

 

 

 

 

Частное решение неоднородного уравнения найдем, например, с помощью функции источника:

G( p) =

1

 

 

 

 

 

 

 

 

 

 

 

1

x

 

 

2 x

 

1

3x

 

 

;

 

g(x) =

2 e

 

−e

 

 

+

2 e

 

.

( p −3)( p −2)( p −1)

 

 

 

 

 

Откуда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

3t

 

2t

)

1

 

x−t

 

2( x−t)

 

1

 

 

3( x−t)

 

 

 

2

e

−e

e

−e

+

e

 

 

yчаст. (x) = ∫(12t

 

 

 

2

 

 

 

2

 

 

dt.

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В результате интегрирования получим:

yчаст (x) = xe2 x + (2x3 −9x2 + 21x)e3 x .

Пример 15.7 ([4], № 52 (2)). Найти частное решение уравнения, удовлетворяющее заданным начальным условиям:

 

0

(0 < x < a),

y ''+ y ' =

 

(a < x < b), y (0) = 0; y '(0) = 2.

f ( x), где f ( x) = 1

 

 

0 (x > b),

 

 

234

Решение. Запишем

f (x) через σ0 (х):

 

 

 

 

 

 

 

 

 

 

 

f (x)=σ0 (х−a) −σ0 (х−b).

 

 

 

 

Далее (по теореме запаздывания)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F( p) =

eap

−

ebp

.

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Переходим к операторному уравнению:

 

 

 

 

 

 

 

p2[Y ( p) −

 

2

 

] + p Y ( p) = eap −ebp ,

 

p2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

Y ( p)

=

 

 

eap −ebp

 

+

 

 

 

2

 

=

 

 

 

 

 

 

 

 

p2 ( p +1)

 

p( p +1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= (e

ap

−e

bp

 

1

 

 

 

1

 

1

 

 

 

 

 

1

 

 

 

1

 

 

 

 

)

 

 

 

 

 

−

 

+

 

 

 

 

 

 

+ 2

 

−

 

 

 

 

;

 

 

p

2

 

 

p

 

p +1

p

 

p

+1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y(x) = (x −a −1 +e−( x−a) ) σ0 (x −a) −

 

−(x −b −1 +e−( x−b) ) σ0 (x −b) + 2 −2e−x .

 

Пример 15.8 ([4], № 53 (2)).

Решить систему уравнений:

y′−2 y −4z = cos x,

y (0) = 0; z (0) = 0.

 

 

 

+ 2z = sin x ,

z′+ y

 

 

 

 

 

 

 

 

 

 

 

 

Решение. Переходим к операторным уравнениям:

 

 

 

 

p

 

 

 

 

 

 

 

p·Y ( p) −2·Y ( p) −4Z ( p) =

 

 

 

 

,

 

 

 

 

 

 

p

2

+1

 

 

2

 

 

 

 

 

 

Z ( p) = −

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

p

2

( p

2

+1)

 

p·Z ( p) +Y ( p) + 2Z ( p) =

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

p

 

+1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 2 ( p21+1 − p12 ) z(x) = 2sinx −2x.

Из второго уравнения исходной системы

y = sin x − z '−2z = sin x −2cos x +2 −4sin x +4x = = 2 + 4x −3sin x −2cos x.

Пример 15.9 ([4], № 54 (2)). Найти общее решение уравнения с переменными коэффициентами: xy ''−(1+ x) y '+ y = 0 .

235

Источник: https://studfile.net/preview/16708793/