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

Практикум по спецглавам высшей математики (ТФКП, ОИ, ТП). Учебное пособие

.pdf
Скачиваний:
0
Добавлен:
12.08.2026
Размер:
1 Мб
Скачать

 

 

 

 

 

 

 

i 32t

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2i

 

2

2

i 274

 

 

2

 

 

 

i6

 

2

 

 

 

 

 

27

i sin

27

 

3i e

 

 

 

 

 

 

 

 

 

 

3 e

 

 

 

 

 

3 e

 

 

 

 

3 cos

 

 

 

4

 

 

 

 

4

 

2

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

2

 

1

 

 

i 2

 

1

 

 

2

 

 

 

 

 

 

2

 

 

 

2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cos6 isin 6

 

 

 

 

 

 

 

 

 

 

 

i

.

 

 

 

 

 

 

2

 

 

 

 

3

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

2

 

 

 

 

3

 

3

 

 

 

 

 

 

3

 

 

 

3

 

 

 

 

 

 

Решение 2. Т.к. подынтегральная функция аналитиче-

ская, то

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

z2 dz 2 4

z

6

 

 

2

 

4 i

6

 

4 1

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

1

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

cos6

2

 

 

 

 

 

 

 

isin 6

2

 

 

 

 

 

 

 

1 6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

27

 

 

 

 

 

27

 

 

 

 

 

 

 

 

2

 

2

 

 

 

 

1

 

 

 

 

 

2

 

 

 

1

 

 

 

 

 

 

 

 

 

3

cos

 

 

isin

 

 

 

 

 

 

1 3

3

 

 

 

 

 

i

3

 

 

 

 

 

 

 

 

 

4

 

 

 

4

 

 

 

 

 

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

i

 

2

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Задача 14. Вычислить интеграл с помощью вычетов

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

a b 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a bcos x

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение. Применяя подстановку eix

 

z , имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

ix

 

 

 

 

 

ix

 

 

 

 

1

 

 

 

 

 

 

 

1

 

 

 

 

z2 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cos x 2 e

 

 

 

 

e

 

 

 

 

 

 

 

2

z

 

z

 

 

 

 

 

 

 

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dz d eix ieix

izdx dx dz

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

iz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

eix

 

 

 

 

 

 

 

cos2 x sin2 x 1,

x 0,2 ;

 

z

 

 

 

 

 

cos x isin x

 

 

 

 

 

 

 

 

т.е. отрезок 0,2 изменения

x

переходит в окружность

z 1 в комплексной плоскости z . Тогда после преобразования, используя формулу (32), получим

41

 

 

 

 

 

I

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b z

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

2

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 iz a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

zdz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

4i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4i 2 i res F zk

 

,

 

 

z

bz

2

2az b

2

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

k 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

zk

- полюсы подынтегральной функции F z , заклю-

ченные внутри области

 

z

 

1 .

 

 

 

z

 

1 при условии a b 0

 

 

 

 

 

 

 

Внутри единичного круга

 

 

 

 

 

 

 

 

 

находится только один двукратный полюс

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z a a2 b2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вычет функции F z

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

относительно

 

bz2

 

2az b 2

 

этого полюса

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

 

 

 

 

 

 

 

z z z1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

res F z1 zlimz dz

 

 

2

z

z1

 

2

z

z2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

z z2

 

 

 

 

 

 

 

 

 

a a2 b2 23

,

 

 

 

 

2

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

b

 

 

z z

2

 

 

 

z z1

4

 

 

 

 

 

 

 

 

 

где

z

a

 

a2 b2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Окончательно

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I

 

 

 

 

 

 

 

2 a

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b2 2

 

 

 

 

 

 

 

 

42

ЧАСТЬ 2

ОПЕРАЦИОННОЕ ИСЧИСЛЕНИЕ

2.1. КРАТКИЕ СВЕДЕНИЯ ИЗ ТЕОРИИ

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

Определение. Действительная или комплексная кусоч- но-интегрируемая функция f(t) действительной переменной называется оригиналом, если

b

– для любого отрезка [a, b] | f (t) | dt

a

(локальная интегрируемость f(t));

f(t) = 0 при t < 0;

– найдутся такие M > 0 и s R, что | f(t) | M est

(|f(t)| может расти не быстрее некоторой показательной функции). Нижняя грань s0 чисел s, удовлетворяющих указанному неравенству, называется пока-

зателем роста оригинала f(t)

Замечание. Для локальной интегрируемости f(t) достаточно, чтобы f(t) была кусочно-непрерывной, то есть на каждом отрезке имела не более чем конечное число разрывов I рода.

Примером оригинала является функция Хевисайда

1,

t 0;

 

(t)

 

Рис. 1

0,

t 0.

Ее показателем роста является s0 = 0. Произведения

eat (t), sin at (t),

tn (t) ( n N) также являются ори-

гиналами. Их показатели роста равны соответственно a, 0, 0 (напомним, что степенные функции при достаточно больших значениях переменной растут медленнее любой показатель-

ной функции). Функции et2 (t),

tg t (t),

1

(t)

ориги-

налами не являются.

 

t

 

 

 

 

g(t)

 

Замечание. Оригиналы вида g(t) (t), где

элемен-

43

тарная функция, принято записывать как g(t) , подразумевая, что множитель (t) всегда присутствует. Тогда, в частности, функцию Хевисайда можно записать как 1, 1= (t) . Аналогично sin t sin t (t) и т. д.

Определение. Если f(t) – оригинал, то комплексная

функция F(p) комплексной переменной p = s + i , определяемая равенством

F( p) f (t)e pt dt ,

называется преобразованием0 Лапласа для f(t) или ее изображением. Соответствие оригинала и изображения принято

обозначать так: f(t) F(p);

f(t), t ( , ) F(p), p C.

2. Свойства преобразования Лапласа

Изображения функций удобно обозначать той же буквой, что и функцию-оригинал, только заглавной. Например, f(t)

F(p); g(t) G(p), h(t) H(p) и т.д.

Теорема. Пусть f(t) – оригинал и s0 – его показатель роста. Тогда

a)его изображение F(p) является аналитической функцией в полуплоскости Re p > s0, то есть в точках p = s + i , для которых s > s0;

б) lim F( p) 0.

Re p

Рис.2

Свойства изображений.

 

 

 

1.

Линейность ( , - числа):

 

 

 

 

f (t) g(t) F( p) G( p) .

2.

Теорема подобия (a > 0):

f

 

aF(ap).

 

f at a1 F

p

,

t

3.

a

a

Дифференцирование оригинала: если f(t) непрерывна

44

F( p) . p

на [0, ) и f (t) – оригинал, то

f (t) pF(p) – f(0).

Аналогично, если f (t) непрерывна на [0, ) и f (t) – оригинал, то

f (t) p(pF(p) – f(0)) – f (0),

т.е.

и т.д.

f (t) p2F(p) – pf(0) – f (0)

4. Дифференцирование изображения:

t f (t) dpd F( p). 5. Интегрирование оригинала:

t

f (t)dt

0

(Напомним, что смысл интеграла в левой части таков:

t

f ( )d . Это интеграл с переменным верхним преде-

0

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

6.Интегрирование изображения: если f t(t) - оригинал, то

f (t)

 

F( p)dp.

t

 

p

7.Теорема запаздывания (запаздывающий оригинал равен нулю при t < a; a > 0; его график получается смещением графика оригинала f(t) на a единиц вправо):

f(t a) e ap F( p).

8.Теорема смещения (изображения) (a С):

eat f (t) F( p a) ,

e at f (t) F( p a).

9.Изображение свертки оригиналов f и g, то есть оригинала со значениями

t t

( f g)(t) f ( ) g(t )d g( ) f (t )d

0 0

45

(напомним, что f g = g f):

t

f ( ) g(t )d F( p) G( p).

0

10. Изображение произведения:

 

1

s i

 

f (t) g(t)

 

 

F(z)G( p z)dz,

2 i

s i

Здесь s0 – показатель роста произведения; рется по прямой s = const.

3. Примеры изображений

s s0.

интеграл бе-

Таблица 1

 

 

 

f(t)

 

 

 

F(p)

 

 

 

 

f(t)

 

 

 

F(p)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1.

 

 

1

1

 

 

 

 

 

 

2.

 

 

t

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p2

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3.

tn , n N

 

 

 

 

n!

 

 

 

 

 

4.

 

 

eat

 

1

 

 

 

 

 

 

 

 

 

pn 1

 

 

 

 

 

 

p a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5.

 

sin t

 

 

 

 

 

 

 

 

 

 

6.

 

cos t

 

 

 

 

p

 

 

 

 

 

 

 

p2 2

 

 

 

 

p2 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7.

 

sh t

 

 

 

 

 

 

 

 

 

 

8.

 

ch t

 

 

 

 

p

 

 

 

 

 

 

 

p2 2

 

 

 

 

p2 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

9.

e

at

sin t

 

 

 

 

 

 

 

 

 

 

10.

e

at

cos t

 

 

 

p a

 

 

p a 2

2

 

 

p a 2 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 2

 

 

 

 

 

 

 

 

 

 

 

f(t)

 

 

 

 

 

 

F(p)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

По свойству 7

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(запаздывание)

 

 

11.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

из формулы 1

 

 

 

 

 

 

f(t) =

(t a)

 

 

 

 

таблицы

 

 

 

 

 

 

 

 

 

 

 

 

F(p) = e ap

1

.

 

 

 

 

 

 

 

 

 

 

Рис.3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

46

Таблица 3

Заметим, что умножением функции g(t) на (t a) мы превращаем ее в нуль при t < a, оставляя ее неизменной в других точках (не путать с запаздыванием, то есть с параллельным переносом графика).

f(t)

12.

f(t) = t t (t a) f(t) = t – (t – a) (t a) –

a (t a).

Рис.4

F(p)

По свойствам 1 (линейность) и 7 (запаздывание) из формулы 2 и 1 таб- лицы F(p) =

=p12 e ap p12

a e ap 1p .

Таблица 4

 

f(t)

 

 

 

F(p)

 

 

13.

 

Применим свойство 7 (за-

 

 

паздывание) к оригиналу из

 

 

примера 10 при a = 2:

 

 

 

f(t) рассмотрим как запаз-

F( p) e p

1

e 2 p

1

 

 

 

p2

 

 

 

p2

 

 

дывающий на 1 оригинал

 

1

 

.

 

 

 

f1(t) из примера 12 (a = 2),

a e 2 p

 

 

 

 

f1(t) = t t (t – 2).

 

p

 

 

 

 

Рис.5

 

 

 

 

 

 

 

4.Нахождение оригинала по изображению

4.1.Теорема обращения и ее следствие. Нахождение оригинала с помощью теории вычетов.

Теорема обращения преобразования Лапласа. Если f(t) – оригинал, являющийся дополнительно на каждом отрезке кусочно-непрерывной и кусочно-монотонной функцией, и s0 – его показатель роста, то для всех точек t, в которых f(t) непрерывна,

47

 

 

f(t) =

1

 

s i F( p)e pt dp, s s ;

 

 

 

 

 

 

2 i

 

0

 

 

 

 

 

s i

 

интеграл берется по прямой s = const.

 

С помощью теоремы обращения можно доказать следу-

ющую теорему.

 

 

 

Теорема. Любая правильная рациональная дробь F(p) =

 

Mm ( p)

, где Mm ( p) и Mn ( p)

– многочлены степеней m и n

 

Mn ( p)

 

 

 

 

 

 

соответственно (m < n), является изображением некоторого

оригинала, который можно найти по формуле

 

 

f (t)

n

 

 

 

Res [F( p)e pt ] ;

 

 

 

 

 

p p

 

 

 

 

 

k 1

k

вычеты берутся по всем особым точкам функции F( p)e pt .

4.2. Нахождение оригинала с помощью таблицы и свойств преобразования Лапласа.

В основе этого метода нахождения оригинала лежит применение свойств изображений к известному изображению более простого оригинала в процессе последовательного его преобразования в заданный оригинал. Особо отметим, что в случае, когда изображение является правильной раци-

ональной дробью, т.е. F(p) = Mm ( p) , где Mm ( p) и Mn ( p) –

Mn ( p)

многочлены степеней m и n соответственно (m < n), можно использовать возможность разложения F(p) в сумму более простых слагаемых. Приведем соответствующую теорему.

Теорема. Каждая правильная рациональная дробь может быть единственным образом представлена в виде суммы элементарных дробей четырех видов:

I.

 

A

.

II.

 

 

A

,

k N, k 2.

p a

 

p a k

 

 

 

 

 

 

 

III.

Bp C

, b2 4c 0.

 

 

p2 bp c

 

 

 

 

 

 

 

 

 

 

 

IV.

 

Bp C

 

,

b2 4c 0,

k N, k 2.

p2 bp c k

 

 

 

 

 

 

 

 

Найдем оригиналы для дробей первых трех видов, используя свойства линейности и смещения.

48

 

 

A

Aeat ;

 

 

 

 

 

A

 

 

 

Aeat

 

 

tk 1

;

 

 

 

 

 

 

 

 

 

p a k

 

 

 

 

 

 

p a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(k 1)!

в дроби III сначала выполним преобразования

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Bp C

 

 

 

 

 

Bp C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p2 bp c

p b2 2 c b42

 

 

 

 

 

 

 

 

и введем обозначения

 

 

 

 

 

b

 

 

b2

 

 

2

 

 

 

 

 

 

 

2 a,

 

c

4 .

Получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Bp C

 

 

 

Bp C

 

 

 

 

B( p a) C Ba

 

p2 bp c

 

 

 

 

p a 2 2

 

 

 

p a 2 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

p a

 

 

 

(C Ba)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

p a 2 2

 

 

 

p a 2 2

Отсюда

 

 

 

 

Bp C

Be at cos t

C Ba

e at sin t.

 

 

p2 bp c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5. Решение линейных дифференциальных уравнений с постоянными коэффициентами с помощью преобразо-

вания Лапласа

 

 

 

 

Уравнение

 

 

 

 

a0 y(n) a1 y(n 1)

... an 1 y an y f (t),

y y(t)

запишем кратко

L[ y] f (t) ,

 

 

 

 

 

где L[ y] a0 y(n) a1 y(n 1)

... an 1 y an y

– линейный

дифференциальный оператор с постоянными коэффициен-

тами, а функция f(t) – оригинал. Обозначим

 

y(t) Y ( p),

f (t) F( p).

 

Тогда по свойству дифференцирования оригинала

 

 

pY ( p) y(0),

 

 

 

y (t)

 

 

 

2

 

 

 

y (t)

p Y ( p) py(0) y (0),

 

---------------------------------

 

y(n) (t) pnY ( p) pn 1 y(0) ...

y(n 1) (0).

49

Применяя преобразование Лапласа к обеим частям уравнения, получим

a0[ pnY ( p) pn 1 y(0) ... y(n 1) (0)] ... an 1[ pY ( p) y(0)]

илиanY ( p) F( p)

a0 pn a1 pn 1 ... an 1 p an Y ( p) G( p) F ( p).

Здесь G(p) – многочлен, определяемый начальными

условиями y(0),

y (0), …, y(n-1)(0). Если все эти условия нуле-

вые, G(p) 0. Обозначим

 

 

 

 

 

a pn a pn 1 ... a

n 1

p a

L( p) ,

0

1

 

n

тогда

 

 

 

 

F( p) G( p)

 

L( p)Y ( p) G( p) F( p);

 

Y ( p)

.

 

 

 

 

 

 

 

L( p)

Уравнение решено в изображениях. Остается перейти к оригиналу y(t).

Замечание. В случае нулевых начальных условий можно избежать нахождения изображения для правой части f(t). В

этом случае G(p) = 0 и Y ( p) FL((pp)) . Если найти оригинал

h(t) для L(1p) , то оригинал дляY(p) можно найти как свертку оригиналов h(t) и f(t) (используем свойство 9 изображений):

t

t

y(t) h( ) f (t )d

или y(t) f ( )h(t )d .

0

0

2.2. ПРИМЕРЫ РЕШЕНИЯ ЗАДАЧ ПО ОПЕРАЦИОННОМУ ИСЧИСЛЕНИЮ

Задача 1. Найти изображение оригиналов

а) f1(t) cos 2t sin 5t, б) f2 (t) sh 3t sh t .

Решение.

а) Оригинал f1(t) преобразуем по формуле

sin cos 12 sin( ) sin( )

и применим свойство линейности изображений, взяв изображение синуса из таблицы.

50

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]