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

Математика_6

.pdf
Скачиваний:
0
Добавлен:
28.12.2025
Размер:
6.14 Mб
Скачать

где SD – площадь круга D:

x2 y2

R2 , равная

R 2 . В итоге:

A 2 R 2 – искомая работа силы.

 

 

Пример 2. Вычислить интеграл J

x 2 y3dx

dy zdz , если

 

 

Г

 

Г есть окружность x2 2

1 в плоскости z=2, обходимая против

часовой стрелки.

Решение. По формуле Стокса (2.23) исходный интеграл сведем к поверхностному интегралу по кругу Т:

 

 

 

 

T :

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Итак, учитывая, что P(x, y, z)

 

 

x2 y3,

Q(x, y, z)

1,

R(x, y, z) 2 ,

имеем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J

 

 

(1)

 

 

(x 2 y

3

dxdy

 

(2)

 

 

 

(1)

dydz

 

T

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(x 2 y3 )

(2)

 

 

 

 

2

2

 

 

 

 

+

 

 

 

 

 

 

 

 

 

dzdx

3

x

 

y

dxdy .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D

 

 

 

 

 

 

Последний интеграл есть двойной интеграл по кругу D Oxy,

на который проектировался круг Т; D:

x2

2

 

1.

Перейдем к

полярным координатам: x=rcos , y=rsin , [0;2 ], r [0;1]. В итоге:

2

1

 

 

 

 

 

2

 

 

 

 

 

1

 

 

 

 

J 3

0 d

0 r2 cos2 sin2 rdr

3

0 cos2

sin2 d

0 r5dr

 

 

.

8

Пример 3. Найти поток П векторного поля

&

(x2; y2;z2)

F

через полную поверхность Т пирамиды

W:

 

x

2y

3z

6,

x

0, y

0, z

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(рис. 2.19) в направлении внешней нормали к поверхности.

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

3 y

 

 

 

 

 

 

x 6

Рис. 2.19

51

Решение. Поток равен П

x 2dydz 2dzdx z 2dxdy .

 

T

Применяя формулу Остроградского-Гаусса (2.24), сводим задачу к вычислению тройного интеграла по фигуре W -пирамиде:

 

 

П

 

 

 

 

(x 2 )

( y2 )

 

 

 

 

(z2 )

 

dxdydz

 

 

 

 

 

 

 

 

W

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

x

2

x

 

2 y

 

 

 

 

2

(x

 

 

 

 

 

z)dxdydz

6

 

 

 

2

3

 

3

 

 

z)dz

 

 

 

 

 

 

 

 

 

 

2 dx

dy

 

(x

 

 

 

 

W

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

0

 

 

 

 

0

 

 

 

 

 

 

 

 

 

3

 

x

 

 

 

 

8

 

 

 

 

4

 

 

 

5

 

 

 

 

 

14

 

 

 

 

 

8

 

 

 

 

 

 

6

2

 

 

 

 

 

 

 

 

 

 

 

 

x 2

 

 

 

 

 

 

 

y2

 

 

 

 

 

dx

 

 

4

 

x

 

y

 

 

 

 

xy

 

 

 

dy

 

 

 

 

 

 

3

3

9

 

 

 

 

9

 

 

 

0

 

0

 

 

 

 

 

 

 

 

 

 

 

 

9

 

 

 

 

 

 

 

 

 

 

 

6

10

 

 

 

 

 

7

x 2

 

13

 

x 3 dx

 

20.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

6

108

 

 

 

 

&

&

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 4. Найти поток П векторного поля

че-

 

 

F x z)k

рез

 

полную

 

 

 

 

 

поверхность

 

T

 

 

 

 

пирамиды

W:

x

2 y 3z 6; x

 

 

0; y 0 ;

 

z

 

 

 

 

(рис.

 

 

2.20), в

направлении

внешней нормали к поверхности.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

O

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T1

 

6

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 2.20

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение. Применим формулу Остроградского-Гаусса (2.24)

П

 

(x

3 z)

dxdydz

 

3 dxdydz

3V 18 , где V – объем пи-

 

z

 

 

 

W

 

 

 

 

 

 

 

 

 

 

 

W

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

рамиды. Сравним с решением непосредственного вычисления потока ( T1,T2 ,T3,T4 – грани пирамиды).

52

П

x

z)dxdydz

 

 

 

 

;

T

 

 

 

T1

T2

T3

T4

,

x

z)dxdy

( x

z)dxdy

;

 

 

 

T3

 

 

T4

 

 

 

 

 

так как проекция граней T3, T4 на плоскость Oxy имеет нулевую площадь (рис. 2.21),

 

x

z)dxdy

 

xdxdy;

 

 

 

T1

 

 

 

G

 

 

 

 

 

x

z)dxdy

 

( x

x

y)dxdy;

T2

 

 

G

 

 

 

 

 

 

 

 

 

 

3

6

2 y

 

 

П

 

6 x

y)dxdy

dy

6

x

y)dx

 

G

 

 

0

 

0

 

 

3

(6

y)2

(6

y)2

dy

1 3

6

2 y)2 dy 18.

0

 

2

 

2 0

 

 

 

 

 

 

 

y

3 x+2y=6

G

O

 

x

6

Рис. 2.21

3.ЭЛЕМЕНТЫ ОПЕРАЦИОННОГО ИСЧИСЛЕНИЯ

3.1.Оригинал и его изображения

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

1)

 

f (t)

 

0 при t

0 ;

2)

существуют

такие постоянные M 0 и S0 0 , что

 

 

f (t)

 

MeS0t для всех t;

 

 

 

53

3) при t функция f (t) непрерывна или имеет конечное чис-

ло точек разрыва первого рода на каждом конечном интервале оси Ot.

Изображением функции f (t) по Лапласу называется функ-

ция F( p) комплексного переменного p i, определяемая

равенством F( p)

e pt f (t)dt . Если f (t) – оригинал, интеграл

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

в правой части последнего равенства сходится при Re p

S0 .

Тот факт, что F( p)

является изображением оригинала

f (t) , бу-

дем обозначать так:

.

 

 

 

 

 

 

 

 

 

.

F( p) L( f (t)) или F( p) f (t), F( p)

(t) .

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 3.1

Изображение основных элементарных функций

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f (t) при t

 

 

 

 

 

L( f (t))

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t n

 

 

 

 

 

 

 

 

n!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

pn

1

 

 

 

 

 

 

 

e t

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cos

t

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

p2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sh

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p 2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ch

t

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

p2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t sin

t

 

 

 

 

 

 

 

2 p

 

 

 

 

 

 

 

 

 

 

 

( p2

2 )2

 

 

 

 

 

 

 

 

 

 

 

t cos

t

 

 

 

 

 

p2

2

 

 

 

 

 

 

 

( p2

2 )2

 

 

 

 

 

 

 

 

 

 

54

3.2. Основные теоремы операционного исчисления

1. Теорема линейного изображения.

Для любых оригиналов f (t) и g(t) и любых чисел a, b L(af (t) bg(t)) aL( f (t)) bL(g(t)) .

 

Пусть всюду в дальнейшем L( f (t))

F( p) .

 

 

 

 

2.

Теорема подобия (изменения масштаба). Для любого

постоянного C

0

 

L( f (Ct))

1

 

F

p

.

 

 

 

 

 

 

 

C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C

 

 

 

 

 

 

 

3.

Теорема смещения. Для любого числа

 

 

 

 

:

L(e t f (t)) F( p

) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4.

Теорема о дифференцировании оригинала. Если функции

f (t), f

(t),

, f (n) (t)

являются оригиналами, то

 

 

 

 

 

L( f (t))

pF( p)

f (0);

 

 

 

 

 

 

 

 

 

 

 

 

 

L( f (t))

p2 F( p)

pf (0)

f (0);

 

 

 

 

 

 

 

L( f (n) (t))

pn F( p) pn

1 f (0)

pn 2 f

(0)

 

 

f (n 1) (0).

 

5.

Теорема о дифференцировании изображения.

 

L(tn f (t))

( 1 n F(n) ( p), n

 

,2, .

 

 

 

 

 

 

 

 

 

6.

Теорема об интегрировании оригинала.

 

 

 

 

t

 

 

F( p)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L f (s)ds

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

P

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7.

Теорема об интегрировании изображения.

 

 

L

f (t)

 

F( y)dy

(если интеграл сходится).

 

 

 

 

t

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8.

Теорема запаздывания. L( f (t t0 ))

e

pt0 F ( p)

,

t0 .

 

 

 

 

9.

Теорема об изображении свертки двух функций.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

L( f1 * f2 )

F1( p)F2 ( p) , где

f1 * f2

0

f1(s) f2 (t

 

s)ds

 

– свертка

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

функций f1(t) и f2 (t) , F1( p)

 

L( f1(t)),

F2 ( p)

 

L( f2 (t)) .

 

Пример 1. Найти изображения функции shat sin bt .

 

Решение. Известно, что

 

 

 

 

 

 

 

 

 

 

 

 

L(sin bt)

b

 

 

 

F( p);

sh at

 

1

(eat

e at ) . Тогда

p2

2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

55

shat sin bt

1

eat sin bt

 

 

1

e at sin bt . По теореме линейности

2

1

2

 

 

 

1

 

имеем L(sh at sin bt)

L(eat sin bt

L(e at sin bt) . В каждом из

 

 

 

2

 

 

 

2

 

полученных слагаемых применим теорему смещения и получаем

1 F( p

)

1

F( p )

1

 

 

 

 

b

 

 

 

1

 

 

 

 

b

 

.Это и

2

2 ( p )2 2

2 ( p )2 2

2

 

 

 

 

есть искомое изображение.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 2. Найти свертку функций t и et и ее изображение.

 

 

 

t * et

t

s ds

 

t

 

 

s ds

 

 

et

t

s ds .

 

 

Решение.

set

 

set e

 

 

 

 

se

 

Вычис-

 

 

 

 

0

 

 

 

 

0

 

 

 

 

 

 

 

 

 

0

 

 

 

ляя интеграл, имеем tet

et (1

t

 

t e t ) . По теореме об изобра-

жении свертки L(t * et )

L(t)L(et )

 

1

 

 

.

 

 

 

 

p2 ( p

1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 3. Найти L(te 2t sin t) .

 

 

 

 

 

 

 

 

 

 

 

 

Решение. Найдем F( p)

L(sin t)

 

 

 

1

 

 

 

. По теореме о

 

 

p2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

дифференцировании изображения

 

 

 

 

 

 

 

 

 

 

 

 

L(sin t

t)

 

( p)

1

 

 

 

 

 

 

2 p

 

 

 

 

 

G( p) . Наконец, по

 

p2

1

 

 

 

( p2

)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

теореме смещения L(e 2tt sin t)

 

G( p

2 )

 

 

 

 

2( p

2 )

 

.

 

 

(( p2

)2 )2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3.3. Отыскание оригинала по изображению

При отыскании оригинала по изображению в простейших случаях используют таблицу изображений основных элементарных функций и теоремы разложения.

Вторая теорема разложения позволяет найти оригинал по известному изображению, являющемуся дробно-рациональной

функцией p :

F( p)

( p) / v( p) , где u( p) и v( p)

– многочлены

от p соответственно степени m и n, причем m n . Если разложе-

ние

v( p)

на

простейшие

множители

имеет

вид

v( p) ( p p

1

)k1

( p p

2

)k2

( p p

r

)kr , k k

2

k

r

n ,

то,

 

 

 

 

 

1

 

 

 

56

как известно, F( p) может быть разложена на сумму элементар-

 

 

 

 

 

 

Ajs

 

 

 

 

 

 

 

 

 

 

 

 

ных дробей вида

 

 

 

 

 

 

;

j

, r; s

, k j . Итак,

 

( p

p j )k j

s

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F( p)

r

k j

Ajs

 

 

 

 

 

 

 

 

 

 

 

(3.1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j 1s 1 ( p p j )k j

s 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Все коэффициенты могут быть найдены по формуле

 

Ajs

1

 

 

 

lim

 

d s

1

 

[( p

p j )

k j

F( p)]

.

 

(3.2)

(s

 

 

 

 

dps

1

 

 

 

 

1)! p p j

 

 

 

 

 

 

 

 

 

 

 

 

Вместо этой формулы для определения коэффициентов Ajs

можно использовать элементарные приемы, применяемые в математическом анализе при интегрировании рациональных дробей. Если все корни многочлена v( p) простые, разложение упрощает-

ся: v( p) ( p 1)(p p2 )

( p

n );

( pj

pk при j

k ) ;

n

Aj

 

 

u( p j )

 

 

F( p)

 

 

, где

Aj

 

.

(3.3)

 

p j

v ( p j )

j 1 p

 

 

 

 

После отыскания тем или иным способом разложения F( p) на простейшие дроби оригинал f (t) находится так:

а) в случае кратных корней знаменателя

r

k j

tk j

s

 

p j t

 

f (t)

Ajs

 

 

 

e

 

;

(k j

s) !

 

j

1 s 1

 

 

 

б) в случае простых корней знаменателя v( p)

r

 

u( p j )

 

p j t

 

f (t)

 

 

e

 

.

1 v ( p j )

 

j

 

 

 

(3.4)

(3.5)

Пример

1. Найти оригинал

 

f (t) , если известно, что

F( p) L( f (t))

 

p

1

 

.

 

p( p 1)( p

)( p

)

Решение.

У изображения F( p) в данном случае все корни

знаменателя – действительные и простые. Поэтому лучше всего воспользоваться формулой (3.5). Имеем

u( p)

p 1;

v( p)

p( p 1)( p )( p

)

p4

p3

11 p2

p; v ( p) 4 p3

18 p2 22 p .

57

Корни

 

u( p1 )

 

1

 

v( p) : p1 ; p2 1; p3 2; p4

;

;

v ( p1 )

6

 

 

 

u( p2 )

;

u( p3 )

v ( p2 )

v ( p3 )

 

3

;

u( p4 )

2

.

2

v ( p4 )

 

3

Отсюда по формуле (3.5) находим f (t) :

f (t)

1

e t

 

6

 

3 e2t

2 e3t .

2

3

Пример 2. Найти оригинал f (t) по его изображению

p

F( p) ( p 1)3 ( p )2 .

Решение. Разложение F( p)

на простейшие дроби имеет вид

F( p)

A11

A12

A13

A21

A22 .

(3.6)

 

( p 1)3

( p )2

p

( p 2)2

p

 

Находим коэффициенты Aij по формуле (3.2)

A

1

lim (( p 1)3 F( p))

lim

 

p

1 ;

 

 

 

11

0! p

1

 

p

p

2 )2

9

 

 

A

1 lim

 

d

(( p 1)3 F( p))

lim

p

 

 

 

 

 

 

 

12

1!

 

 

dp

p

1

( p

2 )2

 

 

 

 

 

p

lim

 

p

( p

2 )2

p 1

2 p

 

1

;

( p 2 )3

27

 

A

1 lim

d 2

(( p

1)3 F( p))

1 lim

p

13

2!

dp2

 

 

2 p 1 ( p 2 )2

 

 

 

 

p

1

 

 

 

 

 

1 lim

 

4

 

 

 

6 p

 

 

 

 

 

1

.

 

 

 

 

 

 

 

 

( p 2 )3

 

( p 2 )4

 

 

 

 

27

 

 

 

 

 

 

 

2 p 1

 

 

 

 

 

 

 

 

 

 

 

 

 

Аналогично получим A

2

;

 

A

1

. Следовательно,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

21

 

27

 

 

22

27

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F( p)

 

1

 

 

3

 

 

 

1

 

 

1

 

 

2

 

 

 

 

1

. Отсюда

 

27 ( p

 

1)3

 

( p

 

1)2

 

p

1

( p

2)2

 

 

p 2

по таблице изображений и теоремам смещения и линейности изображения имеем

58

f (t)

1

3 t 2et

tet

 

t

te 2t 0 e 2t

 

27

2

 

 

 

 

1

(3t 2

2t 2)et

1

(2t 1)e 2t .

54

 

 

 

27

 

 

Заметим, что коэффициенты разложения (3.6) можно найти и таким способом, который применялся в математическом анализе при интегрировании рациональных дробей.

3.4. Решение дифференциальных уравнений и систем дифференциальных уравнений операционным методом

Пусть дано линейное дифференциальное уравнение (ЛДУ) n- го порядка с постоянными коэффициентами

y(n) (t) 1 y(n 1) (t) n y(t) (t) ,

правая часть которого f (t) является оригиналом. Тогда и решение y(t) этого уравнения, удовлетворяющее начальным условиям

y(0) y0 , y (0)

y0 ,

, y(n 1) (0) y0

(n

1) (то есть решение

задачи Коши для данного ЛДУ), тоже будет оригиналом.

Обозначим

изображение искомого

 

решения y(t) через

S( p) , то есть S( p)

L( y(t)) . Используя теорему о дифференци-

ровании оригинала и свойство линейности, находим изображение левой части исходного ЛДУ и приравниваем его к L( f (t)) . В

итоге вместо ЛДУ с начальными условиями получается так называемое изображающее уравнение, которое является линейным алгебраическим уравнением относительно новой неизвестной

функции S( p)

L( y(t)) . Решая изображающее уравнение, нахо-

дим

S( p) . Определяя затем по

S( p) оригинал

y(t) , мы тем са-

мым найдем искомое решение

y(t)

задачи Коши.

Аналогично

решаются и системы ЛДУ.

 

y'' (t)

 

 

e3t , если

 

Пример

1.

Решить

ЛДУ

y (t)

y(t)

y(0)

y (0)

.

 

 

 

 

 

 

 

 

Решение.

Обозначим L( y(t))

S( p) . По теореме о диффе-

ренцировании

 

оригинала

имеем

L( y (t)) pS( p) y0 ;

L( y (t)) 2S( p)

py y

2S( p) .

Тогда

изображающее

 

 

 

0

0

 

 

 

 

 

59

уравнение

таково:

 

 

p2 S( p)

 

pS( p)

 

3 S( p)

1

 

 

 

.

 

 

Отсюда

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S( p)

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

1

 

 

 

. Восстановим теперь

( p

)( p2

p )

 

 

( p

 

1)( p

)2

оригинал y(t)

 

 

 

( p) . Разложим вначале дробь S( p)

на про-

стейшие дроби:

 

 

 

 

1

 

 

 

 

 

 

 

 

A

 

 

 

 

 

B

 

 

 

 

 

C

 

. Ищем

 

( p

 

1)( p

)2

 

 

 

( p

 

)2

 

 

( p

3 )

 

 

p

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A, B, C: 1 A( p

 

1) B( p 3)(p 1) C( p 3)2 . Полагая

p

1,

получаем 1

16C , то есть C 1 16 ;

полагая p

 

3, p

 

 

 

 

,

получа-

ем 1 A B

9C ,

 

 

откуда

B

1 (A

9 C )

 

1

 

,

A

1 .

16

 

Следовательно,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S( p)

 

1

 

 

 

1

 

1

 

 

1

 

 

1

 

 

1

 

 

 

 

 

y(t)

 

 

 

 

 

 

 

 

4 ( p 3)2

 

16 ( p

3) 16 ( p )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

te3t

 

1

e3t

 

1

e t .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

16

16

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение поставленной задачи Коши найдено.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

x

2 y;

 

 

 

 

 

Пример 2.

 

Решить

систему

ЛДУ

 

 

dt

 

 

 

если

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dy

 

2x

 

y

 

1,

 

x(0)

1

 

 

y(0)

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение. Обозначим

L(x(0))

T ( p),

 

 

L( y(t))

 

 

( p) и най-

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

pT( p) ( 1) T ( p) 2S( p)

 

( p 1)T ( p) 2S( p)

1;

 

 

 

pS( p) 5 2T ( p) S( p)

1

 

 

2T ( p) ( p 1)S( p) S

 

1

.

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

Из последней линейной алгебраической системы уравнений

находим неизвестную T ( p)

(например, по формулам Крамера)

 

T ( p)

 

p2

11 p

2

 

 

 

 

 

p2

11 p 2

.

 

 

 

 

 

 

p(( p 1)2

4 )

 

 

 

p( p

1)( p

3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Разложим T ( p)

на

 

простейшие

 

рациональные

дроби:

T ( p)

 

A

 

 

B

 

C

 

 

 

p2

11 p

2

 

.

 

 

 

 

 

 

p

p 1

p 3

 

 

p( p

1)( p

3 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

60