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

книги из ГПНТБ / Цейтлин Я.М. Проектирование оптимальных линейных систем

.pdf
Скачиваний:
2
Добавлен:
24.10.2023
Размер:
9.35 Mб
Скачать

Это отношение не зависит от

как следует из (1.32), тогда и

только тогда, если

 

 

. д$

дф , .

 

Эф

дх

дх

ду

ду

.

-

-.

,

или

т. е. если соблюдаются условия Коши—Римана (1.31). Когда эти условия выполнены, функция W = W (г) называется аналитиче­ ской функцией или просто функцией комплексного переменного Z,

а отношение

=

W

(z) — производной этой функции

ПО Z.

Эта производная,

которая, как следует из (1.32) и (1.33),

равна

dW _

w ,

, .

_ _Эф_ ,

• _£Ф =

JN>_ _

_5ф_

 

dz

 

^ '

дх '

дх

ду

ду '

 

находится путем обычного дифференцирования по аргументу г аналитической функции комплексной переменной W (z). Функция W (z) — аналитическая в данной области R; в каждой точке этой области она однозначна, непрерывна и имеет конечную произ­ водную.

Точки, в которых функция W (z) не аналитическая, называются особыми. Особенности функции комплексного переменного очень важны, так как они определяют ее поведение в комплексной пло­ скости. В особой точке производная функции W (г) или не суще­ ствует или зависит от способа (пути) приближения к этой точке.

Различают три типа особых точек: полюсы, существенно осо­ бые точки и точки ветвления. Нас в дальнейшем будут интересо­ вать полюсы.

5. ИНТЕГРАЛ от Ф У Н К Ц И И КОМПЛЕКСНОГО

ПЕРЕМЕННОГО

 

 

Пусть функция / (г) = ср (х,

у) + ity (х, у)

и ее производная

/' (г) однозначны и непрерывны

во всех точках

некоторой части

плоскости. Положим, что Г есть дуга некоторой кривой, проведен­ ной в этой части плоскости.

Интегралом функции /

(z), взятым по дуге Г, называется

кри­

волинейный интеграл

 

 

 

 

J f (z) dz =

J (ф +

Щ (dx + і dy) =

 

 

г

г

 

 

-

} (ф dx -

tydy) +

/ j (i|> dx + ф dy),

(1.34)

 

г

 

г

 

взятый по этой

кривой.

 

 

 

Заметим, что оба подынтегральные выражения в правой части (1.34) представляют собой полные дифференциалы. Действи­ тельно, если первое подынтегральное выражение — полный диф­ ференциал некоторой функции Р (х, у), а второе — полный диф­ ференциал некоторой функции Q (х, у), то

Ф dx 1|5 dy

tydx •

дР . , дР ,

дх

(1.35)

dQ

 

Из первого выражения (1.35) следует, что

9 =

дР

 

дР

 

 

дх

 

ду

 

 

откуда

 

 

 

 

 

д*Р

_

ду

dty

(1.36)

дхду

 

ду

дх '

Из второго выражения

(1.35) вытекает, что

 

 

,

dQ

 

dQ

 

 

откуда

 

 

 

 

 

 

=

а ф =

_ ^ _

 

(1.37)

дхду

дх

ду

v

'

Выражения (1.36) и (1.37) представляют собой условия

Коши—

Римана (1.31).

 

 

 

 

 

Таким образом, подынтегральные выражения (1.34) действи­ тельно являются полными дифференциалами, в результате чего интеграл (1.34) не зависит от пути интегрирования, т. е. от кри­ вой Г, а зависит только от положения начальной и конечной то­ чек z0 и гг этой дуги.

Таким образом, условия независимости контурного интеграла (1.34) от пути интегрирования совпадают с условиями Коши— Римана (1.31), т. е. с условием аналитичности (регулярности) функции. Отсюда следует основная теорема Коши, состоящая в том, что интеграл, взятый по произвольному простому замкнутому контуру от любой аналитической функции комплексного пере­ менного, как для всех точек плоскости внутри контура, так и на самом контуре равен нулю.

Рассмотрим случай, когда функция F (z) при некотором г = а обращается в бесконечность, причем точка z = а расположена внутри замкнутого-контура Г.

В точке z = а функция F (г) имеет полюс. В этой точке она не аналитическая 1 . Во всех же других точках внутри контура Г и на контуре функция F (г) является аналитической.

Замкнутый контур Г изображен на рис. 1, а. Точке М соот­ ветствует г = а. Опишем вокруг точки М круг радиусом р, вы­ берем на этом круге две точки G и F и соединим их произволь­ ными линиями с произвольными точками А и В контура. В ре­ зультате получатся два простых замкнутых контура I (ABFCGA)

и I I (AGDFBA), обходимых

в положительном направлении, т. е.

так, чтобы ограничиваемая

площадь оставалась слева.

Рис. 1

Для каждого из этих двух контуров

\F(z)dz = 0;

\F{z)dz = 0,

і

и

так как точка М лежит вне этих контуров. Поэтому

\F(z)dz+

\F(z)dz = Q.

іи

В состав этих двух интегралов вошли

J F (z) dz = — f F (z) dz

AG

GA

И

 

\F(z)dz =

JF(z)dz,

BF -

FB

 

которые можно не учитывать,

так

как в сумме они исчезают.

Тогда из остальных частей образуются два контура — первона­ чальный контур Г, обходимый в положительном направлении, и

1

Если в некоторой точке z —. b функция F (г) обращается в нуль, то функ­

ция

F (г) имеет нуль в этой точке.

круг

радиусом

р,

обходимый

в

отрицательном

направлении

(рис.

1, а). Поэтому

 

 

 

 

 

\F(z)dz

+

\ F ( Z ) dz =

j F(z)d(z) \F{z)dz

= 0,

 

і

 

а

г

p

 

откуда следует, что

 

 

 

 

 

 

 

\F(z)dz

= \F

(2) dz.

 

 

 

 

Г

p

 

 

Таким' образом, интересующий нас интеграл по замкнутому кон­ туру равен интегралу по окружности радиусом р, описанной

вокруг

полюса, при неограниченном убывании радиуса

р, так

как величина этого интеграла не зависит

от него.

 

 

Аналогичное положение

имеет место

и

в том случае,

когда

внутри

замкнутого контура

Г функция F (г), имеет более

одного

полюса,

т. е. обращается в

бесконечность

в точках

z = at (і =

~- I , 2,

. . ., п) (рис. 1, б).

Описав вокруг

каждого

из полюсов

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

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

полюса,

получим

 

 

 

 

 

 

\F{z)dz=

£

\F{z)dz.

 

 

 

 

i = i

р( -

 

 

Таким образом,

J F (z) dz

равен

сумме интегралов

j

F (г) dz

по окружностям

радиусами

р,-, описанным, вокруг полюсов, при

неограниченном

убывании этих радиусов.

 

 

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

контуру,

который описан, как из центра, из той точки, где подынтегральная функция обращается в бесконечность (т. е. из полюса), и разде­ ленный на 2яі (і = У1), называется вычетом и обозначается Resa [F (z)]. [Символ Res— сокращение термина «residu» (вычет),

предложенного Коши]. Таким

образом,

 

 

п

j

F (г) dz = 2ш 2 Resa . [F (г)].

Предположим,

что функция

F (z) может быть представлена

в виде

 

 

F ( г )

где / (z) — функция, аналитическая во всех точках внутри зам­ кнутого контура Г (в том числе и в точке z = а) и на контуре. Точка z — а является простым полюсом функции F (г). Этот полюс имеет кратность, равную единице.

Рассмотрим

функцию

 

 

 

 

 

 

 

 

 

 

 

 

 

1 w

2 а

 

 

 

 

 

 

 

Эта функция [в отличие от функции F (г) ] является

аналитической

во всех точках внутри контура

Г

(в том

числе

и в точке г

= а)

и на контуре. Действительно, в точке z =

а функция

Fx

(z)

обра­

щается не в бесконечность, а в производную аналитической

функ­

 

 

 

 

 

ции,

вычисленную

в

точке

 

 

 

 

 

z

= а.

как функция F\

(z)—

 

 

 

 

 

 

Так

 

 

 

 

 

аналитическая,

 

 

 

 

 

 

 

 

Г

 

 

 

Г

 

 

 

 

 

 

 

 

 

Отсюда

следует,

что

 

 

 

 

 

 

 

 

 

 

 

fl')

 

dz-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

г

 

 

г

 

 

 

Рис. 2

 

 

или, в

соответствии

с

пока­

 

 

 

занным ранее,

 

 

 

 

\F{z)dz

=

\F

(г) dz = J

-±&L

dz =

/ (a)

J

dz

 

 

 

 

 

 

где I интеграл

по

окружности

бесконечно

малого

радиуса,

р

 

 

 

 

 

 

 

 

 

 

 

 

 

описанной вокруг точки г а.

 

 

 

 

 

 

 

 

 

 

Для всех точек

окружности

радиусом

р

(рис.

2)

 

 

 

 

 

 

z — а

 

ре'и .

 

 

 

 

 

 

 

Отсюда

 

 

 

 

 

 

 

 

 

 

 

 

 

dz — г'ре1'9 dQ.

Поэтому

JF(z)dZ=:j.

f{z)

dz

=

 

 

 

toe*9 dQ =

/(a)t j * d0

2nif (a),

pe'e

 

 

 

так как обходу точки г = а по окружности соответствует измене­ ние угла 0 от 0 до 2л, или

 

Resa F(z) =

f (a).

 

Отсюда

становится понятным смысл

множителя

, вводимого

Коши

в определение вычета.

 

 

Если функция F (z), имеющая простой полюс в точке z = а, не может быть представлена в виде F (г) = / І ! д» вычет в полюсе

г = а определяется по формуле

 

 

 

 

Resa [F(z)] = lim

 

 

l(z-a)F(z)].

Предположим, что функция F

(z)

может быть представлена

в виде

 

 

 

 

F(z)

 

/(*)

 

 

— a ) 2

'

 

 

где / (z) — аналитическая функция во всех точках замкнутого контура Г (в том числе и в точке z = а) и на контуре. Точка z = а является кратным (имеющим кратность два) полюсом функции F (z).

Продифференцировав выражение

\ ^ d z = 2mf(a)

г

по-параметру а, получим

 

г

откуда

Resa[F(z)] = f (a) =

 

dFdz(г)

Если функция F (z), имеющая полюс кратности два в точке г — а,

не может быть представлена в виде

F (z) = ^J^ay

, вычет в по­

люсе z = а определяется

по формуле

 

Resa IF (z)] = lim

d

l(z - a? F (z)].

 

Продифференцировав

выражение

 

 

\ 1 i

^

¥

d

z

=

2nir(a)

 

г

 

 

 

 

 

 

 

по параметру а, получим

 

 

 

 

 

 

2 \ 1

^

¥

d

z

=

2nir>(a),

 

г

 

 

 

 

 

 

 

ОТкуДа

 

 

 

 

 

 

 

 

 

 

 

 

 

Следовательно, вычет функции F (z) в точке

z =

а, являющейся

полюсом кратности три,

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Rese [F(2)] =

 

 

 

2

 

[(z-ayF(z)]

 

 

- ^ - l i m -d^ -

 

 

 

 

 

 

 

 

 

dz*

 

 

 

 

 

 

Аналогичным

образом

можно

показать,

что

вычет функции

в точке z = а, являющейся полюсом кратности четыре,

 

 

Rese lF(z)] =

^

r

( a

) =

^

^

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Resa [F (г)] =

± - lim - g f -

[(г -

af F (г)].

 

 

 

 

 

as г ^ а

иг.

 

 

 

 

 

 

В общем случае, продифференцировав

выражение

 

 

 

 

І 7 = ^ ^

= 2ш7(а)

 

 

 

 

k — 1 раз по параметру а,

получим

 

 

 

 

 

 

( А - 1 ) 1 J

/ ( г )

 

 

 

 

 

 

 

 

. d

 

 

 

 

 

 

 

 

 

 

dzk~l

' v ;

 

 

£ (z-a)k

 

 

 

1

w

 

 

 

 

откуда следует, что, если функция F (г) имеет в точке г = а полюс

кратности

п,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( я - 1 ) !

 

dzn

 

(1.38)

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Resa [F (z)] =

 

 

 

 

лп—1

[(z -

а)" F (г)].

 

 

( w

_ 1

} 1

lun

 

(1.39)

Поскольку

0! =

1! =

1,

формулы (1.38)

и

(1.39)

пригодны

для

всех п.

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, если функция F (г) имеет внутри замкнутого

контура

Г в точках

z —

аи

а 2 ,

. . .,

ат т

полюсов, каждый из

которых

имеет

кратность

 

kt

(j

=

1,2,

. . .,

/я),

то

 

j F (z) dz =

2ni ^

1 g

J -

m

-±Tj-r\{z

-

a / /

F (г)] Ц

=

Resa.[F(z)].

В заключение приведем вычисление интегралов, определяющих обратные преобразования Фурье и Лапласа, т. е.

1 0 0

— ІСО - <

в случае преобразования Фурье и С+ІСО

С— l o o

вслучае преобразования Лапласа. В этих формулах s = a + г'со — комплексная переменная.

6.

ОБРАТНЫЕ ПРЕОБРАЗОВАНИЯ ФУРЬЕ

И

ЛАПЛАСА

Обратное преобразование Фурье определяется формулой

too

ІСО

Интегрирование выполняется по всей бесконечной длине мнимой оси. Вычисление этого интеграла можно упростить, если приме­ нить теорему Коши о вычетах. Для этого интеграл вдоль беско­ нечной мнимой оси представляется в виде разности интеграла по замкнутому контуру, состоящему из мнимой оси и дуги полу­

окружности бесконечного радиуса, и интеграла по этой

дуге

(рис. 3).

 

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

как

в правой, так и в левой полуплоскости. Рассмотрим сначала слу­ чай правой полуплоскости (рис. 3). Очевидно, что

1= J - J .

АВ АВСА ВСА

Найдем интеграл

по полуокружности

|

. Нетрудно

показать

 

 

 

 

 

 

 

 

 

 

 

 

ВСА

 

 

 

 

 

(см.,

например,

[7]), что при / >

0 значение этого интеграла ста­

новится

бесконечным

и решение

не существует. При t

< 0

рас­

сматриваемый

интеграл

равен

 

нулю.

Таким

образом,

интеграл

по

контуру,

охватывающему

 

правую

полуплоскость,

опреде­

ляет

/ (0 при t

< 0 ,

причем

|

=

|

. Поэтому функция f(t) Ua-

 

 

 

 

 

 

ЛВ

 

 

АБСА

 

 

 

 

 

 

 

 

ходится

путем

непосредственного

применения

теоремы

Коши

 

 

 

 

 

 

о

вычетах

к замкнутому

контуру:

 

 

 

 

 

 

 

 

 

 

 

 

т

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f(t)

= ^ R e s a

. [ F ( s ) ]

=

 

 

 

 

 

 

 

 

 

 

 

lim

 

 

 

 

d k

r l

X.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d s k r l

 

 

 

 

 

 

 

 

 

 

 

/

/=1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X [(s-a,)kJF(s)],

 

 

t<0,

(1.40)

 

 

 

 

 

 

где

 

kj — кратность

полюса

c;-.

(1.40)

 

 

 

 

 

 

 

 

Суммирование

в формуле

 

 

 

 

 

 

производится

 

по

 

вычетам

во всех

 

 

Р и с

з

 

 

полюсах правой полуплоскости. Знак

 

 

 

 

 

 

минус обусловлен отрицательным на­

 

 

 

 

 

 

правлением

обхода

контура

(рис. 3).

Аналогичным

образом

можно

показать,

что интеграл

по кон­

туру, охватывающему левую полуплоскость, определяет функцию

времени /

(t) при t > 0. Обратное преобразование

Фурье

опре­

деляется при этом формулой

 

 

 

 

 

f (0 =

/=1

Re s *y [F (s)]

 

 

 

 

 

 

 

 

 

.ft,—1

 

 

=

lim

d 1

— [(s-aj)kiF(s)\,

t>0.

(1.41)

 

s-»a.

ts 1

 

 

 

Здесь суммирование производится по вычетам a-s во всех полюсах левой полуплоскости.

П р и м е р

1. Пусть

1

 

 

1

 

 

 

 

 

 

 

F

( S ) =

~~ (s •

a)(s -

Ь) ~

~(а + s) (b-s) '

Функция f (t)

при t <

0 определится полюсом первого порядка s = b в правой

полуплоскости.

 

 

 

 

 

 

По формуле (1.40) при т =

1 и kx = 1 получим

 

ЇМ-

_

_ _L

 

 

1

t<0.

 

/ 5 = ь J

a.-j- b

 

 

 

oi

s-\-a

 

Функция / (t) при

0 определится полюсом первого порядка

пол'уплоскости.

 

 

 

 

 

 

По формуле (1.41) при п = 1 и fcj =

1 найдем

 

 

.si

 

 

 

 

t>o.

 

О! \ b s J$=—a

 

^ - h r e - < ",

 

 

a +

b

 

Таким образом,

 

 

 

 

 

 

 

1

e",

* < 0 ;

 

 

a +

b

 

 

f(t)

b

t~

,

t>0.

 

 

a +

 

 

1

 

at

 

 

 

График этой функции

изображен на рис. 4.

 

 

 

 

 

 

 

 

С

 

 

 

 

 

 

J)

 

 

 

 

 

V

 

 

 

 

 

 

 

0

0

t

 

 

 

£

 

 

 

 

 

Рис.

4

 

 

 

 

Рис. 5

в левой

J Re

Обратное преобразование Лапласа определяется формулой

С+ІСО

Интегрирование выполняется вдоль линии АВ, параллельной мнимой оси и смещенной так, что все особенности лежат левее этой линии (рис. 5). Интеграл можно представить в виде

j

_

f

- і

-

1 -

J.

АВ

 

ABCDEA

CDE

 

ВС

EA

Если

 

 

 

 

 

 

 

l i m | F ( s ) | < - ^

и

t>0,

 

 

s->m

4

 

 

 

значение интеграла по дуге CDE стремится к нулю при стремлении радиуса R к бесконечности, как и в случае преобразования Фурье. Можно также показать [71, что при R —* оо

ВС ЕА

Соседние файлы в папке книги из ГПНТБ