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

КР№3 НАША ПО МАТЕМАТИКЕ

.pdf
Скачиваний:
8
Добавлен:
31.05.2015
Размер:
2.37 Mб
Скачать

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

x 2 y2

R 2 , равная

R 2 . В итоге:

A

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

 

 

 

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

x 2 y3dx

dy zdz , если

 

 

 

Г

 

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

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

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

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

T :

x 2

y2 1

z

2

 

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

J

 

(1)

 

 

(x 2 y3

 

x

 

y

T

 

 

+

(x 2 y3 )

 

 

(2)

 

 

z

 

x

 

 

 

x 2 y3, Q(x, y, z) 1, R(x, y, z) 2 ,

dx dy

(2)

(1)

dydz

 

y

 

z

 

 

 

 

dzdx

3 x 2 y2dx dy .

 

 

D

 

 

 

Последний интеграл есть двойной интеграл по кругу D Oxy, на который проектировался круг Т; D: x 2 y2 1. Перейдем к

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

 

2

1

2

1

 

 

J

3 d

r2 cos2 r2 sin2 rdr

3 cos2 sin2 d

r5dr

 

.

8

 

0

0

0

0

 

 

 

 

 

 

 

 

 

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

F (x 2; y2;z2 )

через полную поверхность Т пирамиды 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 y2dzdx z2dx dy .

T

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

П

 

 

 

(x 2 )

 

( y2 )

 

 

 

 

(z2 )

 

dx dydz

W

 

 

x

 

y

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

x

 

2

x

 

 

 

2 y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

2

 

 

3

 

3

 

 

2

(x y

 

 

z)dx dydz

 

 

2 dx

 

 

dy

 

 

 

(x

 

W

 

 

 

 

 

 

 

 

 

 

 

 

0

 

0

 

 

 

0

 

 

 

 

 

 

3

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

2

 

 

 

8

 

 

4

 

 

 

5

x 2

 

14

 

 

 

 

 

 

8

y2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

4

 

x

 

y

 

 

x y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0

 

 

 

3

 

3

 

 

9

 

 

9

 

 

 

 

 

9

 

 

6

 

 

 

 

7

 

 

 

 

13

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

10

 

x

x 2

 

x 3 dx

 

20.

 

 

 

 

 

 

 

 

0

 

6

108

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

y z)dz

dy

(x 3z)k че-

рез

полную

 

 

поверхность

 

T

 

пирамиды

W:

x 2y

3z 6; x

0; y 0 ;

z 0

(рис.

2.20), в направлении

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

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

T3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T4

 

 

 

 

 

 

T2

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

O

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T1

 

6

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 2.20

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

П

(x 3z)

dxdydz

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

 

W

z

W

 

 

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

52

П

(x 3z)dxdydz

 

 

 

;

T

 

T1

T2

T3

T4

,

 

 

 

 

 

 

(x

3z)dxdy

(x 3z)dxdy

0;

 

 

T3

 

T4

 

 

 

 

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

(x

3z)dxdy

 

xdxdy;

 

 

 

 

 

T1

 

 

G

 

 

 

 

 

 

(x

3z)dxdy

 

(x 6

x

2 y)dxdy;

T2

 

 

G

 

 

 

 

 

 

 

 

 

 

3

6 2 y

 

 

П

(6 x

2 y)dxdy

dy

 

(6 x

2 y)dx

G

 

 

 

0

0

 

 

 

3

2 y)2

 

(6 2 y)2

 

1

 

3

2 y)2 dy 18.

(6

 

 

dy

 

 

 

(6

2

2

 

0

 

 

 

0

 

y

3 x+2y=6

G

O

6

x

 

Рис. 2.21

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

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

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

1)

 

f (t)

 

0 при t 0 ;

2)

существуют

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

 

 

f (t)

 

MeS0 t

для всех t;

 

 

 

 

 

 

 

 

53

3) при t 0 функция 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)

f (t) .

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 3.1

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f (t) при t

0

 

 

 

 

L( f (t))

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t n

 

 

 

 

 

 

 

 

n!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

pn

1

 

 

 

 

 

 

 

e t

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p 2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cos t

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p 2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sh

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p 2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ch

t

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p 2

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(t n f (t))

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

 

1,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 0 .

 

 

9.

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

L( f1 * f2 )

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

f1 * f2

 

f1(s) f2 (t

 

s)ds

 

– свертка

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

функций 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

b2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

55

 

 

 

 

 

 

 

 

 

shat sin bt

1

e

at

sin bt

 

 

1

e

at

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

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

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

1

L(eat sin bt

1

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

2

2

 

 

 

 

 

 

 

 

 

 

 

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

1

F( p

a)

1

 

F( p a)

1

 

 

 

 

 

b

 

1

 

 

 

 

b

 

 

 

.Это и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

2

 

 

2 ( p

a)2

b2

 

 

2 ( p

a)2

 

b2

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

t

 

 

 

 

t

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

Решение.

t * et

set

s ds

 

set e

 

s ds

 

 

et

se

s ds . Вычис-

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

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

 

et (1

te 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)

F ( p)

 

1

 

 

 

 

 

 

 

2 p

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p2

1

 

 

( p2

1)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

G( p

2)

 

 

 

 

2( p

 

2)

 

 

.

 

 

 

 

 

 

 

 

 

 

 

(( p2

1)2

 

1)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

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

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

функцией p :

F( p)

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

– многочлены

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

ние

v( p)

на

простейшие

множители

имеет

вид

v( p) ( p p )k1

( p p )k2

( p

p )kr , k k

2

k

r

n ,

то,

1

2

 

r

1

 

 

 

56

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

 

 

 

 

 

 

Ajs

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

;

j 1, r; s

1, k j . Итак,

 

( p

p j )k j

s 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

k j

Ajs

 

 

 

 

 

 

 

 

 

 

 

 

F ( p)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(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

 

 

 

 

 

 

 

 

 

 

 

 

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

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

ся: v( p) ( p

p1)(p

p2 ) ( p

 

 

pn );

 

 

( p j

 

 

 

pk

при j

k) ;

 

 

 

n

 

 

Aj

 

 

 

 

 

 

 

 

 

u( p j )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F ( p)

 

 

 

 

 

 

, где

 

Aj

 

 

 

 

 

 

.

(3.3)

 

j 1 p

p j

 

 

 

v ( p j )

 

 

 

 

 

 

 

 

 

 

 

 

 

После отыскания тем или иным способом разложения F( p)

на простейшие дроби оригинал

 

 

f (t) находится так:

 

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

 

 

 

 

 

 

 

r

 

k j

 

 

t k j

 

s

 

 

p

j

t

 

 

 

 

 

 

f (t)

 

 

 

Ajs

 

 

 

 

 

 

e

 

 

;

 

 

(3.4)

 

 

 

 

 

(k j

 

 

 

 

 

 

 

 

 

 

 

 

j

1 s

1

 

 

 

s) !

 

 

 

 

 

 

 

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

 

 

 

 

f (t)

r

 

 

u( p

j

)

e

p j t

.

 

 

 

(3.5)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

1 v ( p j )

 

 

 

 

 

 

 

 

Пример

1. Найти

 

оригинал

 

 

 

f (t) ,

 

 

если

известно, что

F ( p) L( f (t))

 

 

 

p

1

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p( p

1)( p

2)( p 3)

 

 

 

 

 

 

 

 

Решение.

У изображения

 

F( p)

 

 

в данном случае все корни

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

u( p)

p 1;

v( p)

p( p 1)( p 2)( p

3)

p 4

6 p3

11 p 2

6 p; v ( p) 4 p3

18 p 2 22 p 6.

57

Корни

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v( p) : p1

0;

 

 

p2

 

 

 

 

1;

 

 

 

 

p3

2;

 

 

 

 

p4

3;

 

u( p1 )

 

 

1

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v ( p1 )

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u( p2 )

 

1;

 

 

 

 

u( p3 )

 

 

 

 

 

 

3

 

;

 

 

u( p4 )

 

 

2

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v ( p2 )

 

 

 

 

v ( p3 )

 

 

 

 

2

 

 

 

v ( p4 )

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отсюда

 

 

 

 

по

 

 

 

 

 

 

 

 

 

формуле

 

 

 

 

 

(3.5)

 

 

 

 

 

находим

f (t) :

f (t)

 

1

 

 

 

e

t

 

 

3

e

2t

 

2

e

3t

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

 

 

 

 

 

2

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример

 

2.

 

 

Найти

 

оригинал

 

 

f (t)

 

по

 

его

изображению

F( p)

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( p 1)3 ( p 2)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

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

F( p)

 

 

 

 

 

 

A11

 

 

 

 

 

 

 

 

 

 

 

 

A12

 

 

 

 

 

A13

 

 

 

 

 

 

 

 

A21

 

 

 

 

 

 

A22

 

 

. (3.6)

 

 

 

 

( p 1)3

 

 

 

 

 

( p 1)2

 

p 1 ( p 2)2

 

 

 

 

 

p 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Находим коэффициенты Aij

 

 

по формуле (3.2)

 

 

 

 

 

 

A

 

 

 

 

 

 

 

1

 

 

lim (( p

 

 

1)3 F ( p))

 

lim

 

 

 

 

 

 

 

 

p

 

 

 

 

1

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

11

 

 

 

 

0! p

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p 1

( p

2)2 9

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

1

 

lim

 

d

(( p 1)3 F ( p))

 

 

lim

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

12

 

 

 

 

1!

 

 

 

 

 

dp

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

1

 

 

( p

2)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

lim

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 p

 

 

 

 

 

 

1

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( p 2)3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p 1 ( p 2)2

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 p 1

 

 

 

 

 

 

 

 

 

 

 

 

 

27

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Аналогично получим 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

2

e

t

te

t

e

t

2te

2t

2

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) a y(n 1)

(t) a y(t) f (t) ,

1

 

n

правая часть которого

f (t) является оригиналом. Тогда и реше-

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

y(0) y

0

, y (0)

y

0

,

, y(n

1) (0) y (n

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

 

 

 

 

 

0

 

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

Обозначим

изображение

искомого

решения 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)

задачи Коши.

Аналогично

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

 

 

 

 

 

 

 

 

Пример

1.

Решить

 

ЛДУ

y'' (t)

2y (t)

3y(t)

e3t , если

y(0)

y (0) 0 .

 

 

 

 

 

 

 

 

 

 

Решение. Обозначим L( y(t))

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

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

 

оригинала

имеем

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

L( y (t)) p2S( p)

py

y

0

p2S( p) .

Тогда

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

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

59

 

 

 

 

уравнение таково:

p2 S ( p) 2 pS( p) 3S ( p)

 

1

 

. Отсюда

 

 

 

p

 

3

 

 

 

 

S( p)

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( p

3)( p2

2 p

3)

 

 

( p

 

1)( p

3)2

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

.

 

 

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

S ( p)

на про-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

A

 

 

B

 

C

 

 

. Ищем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( p 1)( p

3)2

 

 

 

( p

 

 

3)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 0 ,

получа-

ем 1

A

3B

 

9C ,

откуда

B

 

1

( A

9C

1)

 

 

1

 

,

 

A

 

1

.

3

 

16

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S ( p)

 

1

 

 

 

1

 

 

 

1

 

 

1

 

 

1

1

 

 

 

 

 

y(t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4 ( p 3)2

 

 

 

16 ( p

3)

 

16 ( p

1) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

te3t

 

1

e3t

 

1

 

e t .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

16

16

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

x

2 y;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 2.

 

 

Решить

систему

ЛДУ

 

 

dt

 

 

 

если

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2x

 

y

1,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x(0)

1,

 

y(0)

5 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

L(x(0))

 

T ( p),

 

 

L( y(t))

 

S( 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