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

книги из ГПНТБ / Основы авиационной автоматики учеб. пособие

.pdf
Скачиваний:
19
Добавлен:
22.10.2023
Размер:
13.63 Mб
Скачать

Такие часто встречающиеся в автоматике функции времени, как постоянная функция, линейно возрастающая функция, экс­ понента, степенные функции и др., не имеют преобразования Фурье. По этой причине постараемся расширить класс функций, к которому можно применить преобразование, подобное преоб­

разованию Фурье. Рассмотрим

функции ф(^),

удовлетворяю­

щие условиям:

 

 

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

значений аргу­

1. ф(£)

равняется нулю для

мента: ф {t) = 0 при £<Т).

 

 

 

 

 

| ф(/?)| <се^т

2. ф(0

имеет экспоненциальную мажоранту

т. е. существуют такие числа

с и

ß, при которых указанное не­

равенство выполняется для всех значений ^ >

О-

 

т

,

,

,,,

( 1 при t

>

О

, которую при-

Так, например, функция

1

(t) =

I 0 при t

<

О

 

 

 

 

 

пято называть «единичным скачком» или «ступенчатой функци­

ей», в этот класс входит. Функции

в

 

е'3 ’ W ' T z r ' W *

этот класс не входят, так как не имеют экспоненциальных мажо­ рант. Функция sin t в этот класс не входит, так как не выполня­

ется первое условие,

а функции

е5' 1 (t),

eSi ts 1(t),

t 31 (t)

входят, так как выполняются оба условия.

 

 

 

 

Минимальное значение чисел ß (точнее говоря, нижняя грань

inf ß),

при

котором

выполняется второе

 

условие, называется

степенью роста функции

ф(£).

Степень роста

функции

будем

обозначать буквой В. Так, для функции

5e3/sin 21 1 (t)

степень

роста В = 3, а для функции sin 2t \ (t) степень роста 5 = 0.

 

Далее введем в рассмотрение функцию

 

f(t) = ф(^)е—т/,

где

т[ — произвольное число, превосходящее

степень

роста

функции т > 5 .

 

 

 

является абсолютно интегриру­

В таком случае функция f(t)

емой

 

 

 

 

 

 

 

 

 

 

 

 

J |/ ( 0

\ d t <

fce-O-^v dt =

 

<С оо.

 

 

---- СО

 

ö

 

 

 

 

 

К функции f(t)

применимо преобразование

Фурье /

(/<о) =

со

 

 

оо

 

 

 

 

 

 

 

= J

f{t) e~imt dt — j" ф (^) е _(т+-'шН dt.

 

Обозначим комплекс-

— со

 

7 + уш

0

 

 

 

В таком случае инте-

ное число

буквой/?: 7 + у'ш = р.

грал

Jсо ф(Д e~pl dt,

зависящий от параметра р, задает некото-

о

рую функцию от этого параметра, которую мы будем обозначать

ф(/?)= [ 'b(t)e-Pt clt.

(1.3)

О

 

27

Применяя к функции/(/ш) =

ф (/7) обратное преобразование

Фурье, получим

1

со

^

 

 

 

 

 

f (t) = è (t) е - і ‘ = ----

Г Ф(р) eiast dm,

р = у-|-/ш;

Ф(^)= ет<—

2тс

<>

 

f ф(/7) &р‘ d<o.

Гф(р)еіи>‘ da> =

1

00

 

I

00

2тг

J

 

J

 

 

— со

 

— со

 

Проведя в последнем интеграле замену переменной интегрирова­

ния р =

Ч + у'ш, получим

 

 

 

 

Ф(0 =

1

Tt ;co -

(1.4)

 

-----г

[ Ф{p)&P‘ dp.

 

 

 

7 - у 'с о

 

Итак,

функция ф (0

представляется в виде интеграла

(1.4)

но прямой линии, проходящей параллельно оси ординат справа от точки В.

Формулы (1.3) и (1.4) носят название прямого и обратного преобразования Лапласа. В дальнейшем мы будем пользоваться обозначениями

 

 

со

 

 

 

L [ф ^)) =

Ф(р) =

I" ф'(0 £~pl d t

 

 

 

6

 

 

 

^-_МФ(/7)1 =

Ф(0 = ^

1

7V “

-

(1.5)

Т

J

Ф(p)eP‘ dp

 

 

 

’v

t - j oo

 

 

и называть функцию ф(£), принадлежащую указанному

выше

классу функций, оригиналом, а функцию ф (р) — ее изображе­ нием.

Приведем здесь примеры преобразований Лапласа от часто встречающихся функций, которые мы в дальнейшем, если нет специальных оговорок, будем считать равными нулю при отри­ цательных значениях аргумента t:

м

=

\ eal e~pt dt =

 

 

( 1.6)

= )

 

 

 

 

Преобразование Лапласа от функции

1 (0 = |

^

полУ"

чим, положив в предыдущей формуле а = 0

L\ 1(г?)] = — . Диффе-

ренцируя равенство

(1.6)

по а оправа и слева п

Р

получим

раз,

Гtnе*e~pt dt

= ------- —-----

 

 

 

J

 

(p — fl)"+I

 

 

 

28

Последнее же означает, что L \іпея/| = -------------

.

Учитывая, что

— а)л+1

L fe(e+^s)/j = L [е"1' cos

9J + /е 0/ sinß/] =

IР CL “I“ yö

=----------- —= -------- -—— и приравнивая в этом выражении p - a - j Q a)2+Q2 * v F

действительные и мнимые части, получим/. [e“* cosö/] =

Р ----- ^

 

г Г

t

Г\Л

 

Й

 

.

 

= ----------------- и

L I

sin Уг

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

(здесь мы ,вос-

( р - а ) 2+ Й 2

1

 

J

( р - а )> + 92

у

Лапласа,

пользовались свойством линейности

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

которое доказывается ниже).

 

 

 

 

 

Сведем полученные результаты в таблицу.

 

 

О ригинал x(t)

^at

1 1)

tneal

іп

sin

cos 0 /

e^ sin Q ^

e 0,,cos Ш

И зобр аж ен и е

1

1

 

п\

п\

Q

P

Q

p—a

х(Р)

р - а

р

(р-а)п+ 1

рп+ 1

р2+

S 2 pi+Ü2 ( p - a ) 2+Q3 ( p - ß ) 2+Q*

§

1.2. СВОЙСТВА ПРЕОБРАЗОВАНИЯ ЛАПЛАСА

 

 

 

1.

Свойство линейности

 

 

Убедимся в том, что принцип суперпозиции

 

 

Цсі х і (0 + сгх ч(0] =

сх L [х (0] + c2L [x2(t)\;

 

L - 1 [с1х 1(р) + c2x2{p)\=c1L~:l [х(р)} + CzL-'X'ip)}

выполняется и для прямого и для обратного преобразования Лап-,

ласа.

 

 

 

 

 

 

 

 

 

Действительно,

из (1.5)

имеем

 

 

 

 

 

 

 

 

со

 

 

 

 

 

L[cxx }(t)-\-c2Xi(t)=

 

c2x2(t)]e~pfdt=

 

 

со

 

 

о

со

 

 

 

 

 

 

 

 

 

 

 

 

=

cl jxrftfrr&dt

+ с2 I" x2(t)e~pt dt =

 

 

 

о

 

 

 

о

 

 

 

 

 

 

= СіХ1(р) +С2х 2(р).

 

 

 

~

 

 

 

1

т+J“

f

[ C i 4

> 1 ( p

H - c 2 4 ' 2 ( Jo

L~l[C\X,(p) + c2x2{p)\ = - —-

 

 

 

 

 

£"-J

1 —j c o

 

 

 

 

=-\-c2x 2[t).

Свойство линейности часто используется для вычисления ориги­

налов от изображений х(р), представляющих собой рациональ­ ную функцию аргумента р-

29

П р и м е р .

Требуется определить

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

по его изоб-

'■/

1

Ьр + 1

 

 

раженшо х(р) =

----------------- .

 

 

 

 

(Р+2)2(Р + 1)

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

можно полу­

Используя свойство линейности,

чить, не обращаясь к общей формуле (1.4), путем сведения х(р) к линейной комбинации табличных выражений. Разложим функ­

цию х(р) на элементарные дроби

5р + 1

 

(Р+2)2(Р + 1) (р + 2)2 +

р + 2 р + 1

(р+2)2 + р + 2

р + 1

Числовые значения с\, сг и Сз вычисляются методом неопределен­ ных коэффициентов.

Далее имеем в соответствии с таблицей изображений и ори­ гиналов

x(t) = L~l

+

1

■9L-1

 

+ 4L-1

1

1

1(р + 2)2J

7 + 2 ]

 

(р+2)Цр + \)

 

— 4L~1

1

=

+ 4 е -2' -

4е~'.

 

 

Р+ 1

2.Изображение производной

Изображение производной L[x'(t)]

может

быть выражено

через изображение функции L[x(t)]=x(p):

 

 

со

с о

 

L[x'(t)\ = j x'(t)e~pt di = x(t)e~pt |

— ^ x(t)(—p)e~pt dt =

О

0

0

 

= pL[x{t)} - х { 0 )= р х { р )

-

x(0).

(1.8)

Применяя эту формулу n раз, получим

 

L[ jc<»> (t)\ = pL [xl*~» (*)] - x

^

(0) =

• • • =

= p nx(p) p n~x X (0) — p n~2x (-l'i(0) — • • ■— рДп-21(0)— ^cC^—2) (0)•

( 1-9)

30

t

3. Изображение интеграла у (t) = j х (т) di

о

Учитывая, что y'(t)=x(t) и у(0)=0, из формулы изображе­ ния производной получаем

L [у '(01 = L\x{t)) = pL \y{t)\.

В таком

случае

изображение интеграла

от функции

имеет

вид:

 

 

 

 

 

 

L[y{t)\ = — х(р).

 

( U 0)

 

 

Р

 

 

В дальнейшем мы будем пользоваться соотношениями

 

L [e~a/x[t)]

= х ( р +

а) и L~l [х (р + а)] =

e~at х (t),

(1.11)

справедливость которых непосредственно следует из (1.5).

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

Сверткой двух функций х\(t) и x2(t) называется функция

y ( i ) = j x x[t i ) x 2(i)di = j x x{i)x2{t — i)di.

(1.12)

Будем считать, что функции Xi(t) и х2(t) являются оригина­ лами, т. е. Xi(t) =x2(t) =0 при t<C_0. В таком случае пределы ин­ тегрирования в интеграле свертки могут быть сужены

у [р} = Г у (^) Q-pl dt =

Р и с . 1.1. Область инте-

J

грирования

СО t

 

x x(t — т)х2(т) die~ptdt.

 

о о

 

Справа здесь стоит двукратный интеграл, распространенный на

сектор 5 в плоскости і, t (рис.

1.1). Изменяя в нем порядок ин­

тегрирования, получим

 

ОО

DO

СО

ОО

 

у(р) = j1х 2 (*) d i j

e~ptx x (t—i) d t —

j" x 2{i)orp''di^xx{tx)&-p^ d t x =

 

=

X\(P)x2{p),

где

tx = t — г.

 

31

Итак, изображение свертки двух функций равно произведению их изображений

 

t

 

 

L\

- t)a'2(t)^ t] = Xx(p)x2(p)

и L~l[xx{p) x^{p)\ =

 

0

 

 

 

t

 

 

 

j1A'j (t

dl.

(1 14)

 

6

 

 

§1.3. СВЯЗЬ МЕЖДУ ИЗОБРАЖЕНИЯМИ ВХОДНОГО

ИВЫХОДНОГО СИГНАЛОВ ЛИНЕЙНЫХ СТАЦИОНАРНЫХ СИСТЕМ.

ПЕРЕДАТОЧНАЯ ФУНКЦИЯ

Из введенных нами ранее операторов линейными и стацио­ нарными, как мы показали, являются операторы дифференциро­ вания, оператор умножения на постоянную величину, оператор, задаваемый дифференциальным уравнением с постоянными ко­ эффициентами при нулевых начальных условиях, интегральный оператор с весовой функцией g (t— т), оператор сдвига на время а, многомерный оператор, задаваемый системой дифференциаль­ ных уравнений при нулевых начальных условиях.

Рассмотрим связь между изображениями входного и выход­ ного сигналов, задаваемую этими операторами.

Для дифференциального оператора из связи между оригина-

.

dx{t)

лами ѵ(г) =

------ вытекает в силу теоремы об изображении про.

 

dt

изводной зависимость между изображениями

у (р)=рх( р)—х(0).

Для оператора, задаваемого дифференциальным уравнением я-го порядка с постоянными коэффициентами

а-п У{п) (t) Н-----+ axyW{t) + а0 у [t) = b0x{t)

(1.15)

при нулевых начальных условиях, связь между

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

х(р) и у(р)

входного и выходного сигналов получается, если пре­

образование Лапласа применить к левой и правой частям

(1.15)

L\any [n){t) -I-------\-axyV'{t)-\-a0y{t\\ = L [bx{t)\.

Отсюда в аилу

свойства

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

получаем

 

 

 

У(Р) =

_______ Ь-х(р)

(1.16)

апРп -і-----+ ахр + а0

 

Из теоремы о свертке получаем зависимость между изобра­

32

жениями х(р)

и у(р)

входного и выходного сигналов интегралы-

 

 

/

 

іного стационарного

оператора у(і) = J g(t —

di, кото-

рая, согласно

 

u

 

(1.14), имеет вид:

 

 

 

y(p)=g(p)x(p),

(1-17)

где g( p ) —L[g(t)] — (преобразование Лапласа функции g(t). Для оператора сдвига y(t)=x(td), применяя к левой и пра­

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

DO

у (/?) = L (і — а)] == jjc(t a)e~ptdl*= j X (t — а) e~pt dt =

со

 

 

 

= j“ x ( t ' ) e - pl'

e~apdt'

= e~apx[p}..

 

о

 

 

 

Итак, изображения х(р)

и у(р)

входного и выходного

сигна­

лов оператора сдвига связаны соотношением

 

У(р) = е - арх ( Р).

(1.18)

Из приведенных выше зависимостей между изображениями

входных и выходных сигналов линейных стационарных

систем

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

Передаточной функцией W(p) линейной стационарной систе­

мы называется

отношение

изображения

выходного сигнала к

изображению входного сигнала

 

 

W(p) = УІР)

(1.19)

 

 

х(р)

 

Передаточные функции

для указанных выше линейных ста­

ционарных систем имеют вид:

 

W (р ) = --------------

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

для систем, задаваемых диффе-

апРпЛ-----

+ ахр + а0

 

ренциальным уравнением с постоянными коэффициентами;

W(p)—g(p) — для систем, задаваемых интегральным уравнени­ ем;

W(p) = e~ap — для систем, задаваемых оператором сдвига во времени.

3 . И зд. № 5312

Дифференциальный оператор в случае произвольных входных сигналов не имеет передаточной функции, так как изображение

выходного сигнала у(р') не представимо в виде y(p) — W(p)x(p) в силу наличия члена х(0) в формуле (1.8). Однако, если рас­ сматривать лишь непрерывные дифференцируемые функции (а только на таких функциях определен этот оператор), то для этих функций X(0) = 0; и в таком случае передаточная функция диф­ ференциального оператора W(p)—p.

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

Покажем на примерах, что у нелинейных и у нестационарных операторов передаточной функции в определенном выше смысле

не существует.

 

 

 

 

Рассмотрим

нелинейный оператор

возведения в квадрат

y(i) =x2(t). Если подать на его вход единичный скачок

1(0» то

отношение изображений выходного сигнала

к входному равно:

L{\-(t)}: L[\(t)]=

•гI п

 

вход x(t)—t,

 

= l , a если подать на

то это

 

1

1

 

 

 

]ІП2

 

отно-

отношение равно: L[t2]: L[t\ — —■— = — .Таким образом,

 

Чр

Р

 

 

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

Рассмотрим теперь нестационарный оператор y(l) = tx(t).

Для входного сигнала x(t) 1(f) имеем отношение L[t]: L[l]=

=■—= — , а для входного сигнала x(t)=t — отношение

Чр

р

2/р3

2

L[t2\ : L[t]= -У— =

— . Таким образом, нестационарные опера-

 

 

1ІР2

Р

торы передаточной функции не имеют.

Обратимся теперь к многомерному линейному стационарному оператору, задаваемому системой линейных дифференциальных уравнений с постоянными коэффициентами

Уі(£)= апУ і (0 Н------

^ аіпУп(^) +

(ОН--------

\rblmx m(t);

_y,(0) = 0;

i = 1 , ___ n

 

(1.20)

при нулевых начальных условиях.

Применив преобразование Лапласа к левым и правым ча-

34

стям системы уравнений (1.20), получим систему линейных алге­ браических уравнений:

І Р - а \ \ ШР ) - а і і У і { р ) --------

аь,Уп{Р)-ЬпМР)-іг- ■+Ьы хт(р)\

- а22Уі(Р)+(Р-а2М Р ) ------

а2пУп(Р) b2lx{(p)+.. +Ь2тх т(р);

-^піУі(Р)-а„аУі[р)-- ’ - M P - а„пЪ’п(Р)=ЬПІх1(р)+ •+Ьптхт (р),

 

 

(1.21)

в-которой У\(р),... ,yk{p),

■. . >Уп(Р) выступают в качестве не-

известных.

 

 

Решая эту систему методом Крамера, получим

УкіР)= - г г т (V (Р)*і (Р) +

V (Р) а'2 ІР) + -----h Ѵ г ІР)хт (р)) =

= W k l ( р) X , (р) + U ^ 2

( / 7 ) х2(р) + -------- Ь W km(p)xm( / ? ) ,

где Д (/Р) — определитель алгебраической системы (1.21)

Р Оц

П]2 —Яіз— • • • — 0-1п

---&21 Р--- ^22--- ^23---- • • • — 2п

Д (Р) =

ап2 апЬ • • ■~\~(р апп)

ап\

к 1к(р) — определитель, получающийся из определителя заменой /г-го столбца на столбец bn, Ь12, . . . , Ъ(п\

Др )

Д (Р)

( 1.22)

Д(р)

Функцию (р) принято называть передаточной функ­ цией многомерной системы от і-го входа к k- щ выходу. Сово­ купность передаточных функций многомерной системы, записан­ ных в виде таблицы

w n ( p ) W M . . . w lm{p)

W M W 22( p ) . . . w am(p)

W(p) =

Wnl{p)Wn2{ p ) . . . W nm(p)

принято называть матрицей передаточных функций многомерной линейной стационарной системы.

3*

35

Из (1.22) следует, что

. Ж к1{р) = Щ - ,

(1.23)

х,(р)

если все входные сигналы, кроме г'-го, равны нулю. В таком случае передаточная функция W ы(р) многомерной стационар­ ной системы от г'-го входа к А-му выходу может быть определена

как отношение изображения у к {р) сигнала на А-ом выходе к

изображению x t (р ) на г-ом входе, при условии, что на осталь­ ные входы поступают сигналы, равные нулю.

ВОПРОСЫ И ЗАДАЧИ

1. Движение самолета вокруг продольной оси х приближенно описывает­ ся дифференциальным уравнением

/т ( 0 = Л Г ^ ( * ) + Ж*«»„

где Т (і) — угол крена; 63(1) — угол отклонения элеронов;

Mwx> М*э — частные производные от продольного момента по скорости вра-

X

X

8Э, которые

в рамках данной задачи

щения ш.ѵ и углу отклонения элеронов

следует считать известными числами.

 

 

 

 

Найти передаточную функцию самолета по каналу крена, считая угол

крена выходным сигналом, а угол поворота элеронов — входным.

 

Ответ:

Ж8

 

 

Щ / > ) = ! ^ - =

т =

J

А =

 

р(Тр + 1) Ж" Ж > '

І( Р)

2.Движение системы описывается дифференциальными уравнениями

Уі (0 =

-

2Уі (t) + Зу2 (0 +

 

2хх [ і ) - х 2(0;

Уо(*) =

У і ( 0 + .УаЮ — *і (0 +

2*2(0-

 

 

— 5

- Р

+ 7

О т

 

?2 + Р + 1 : р 2+ р + 1

 

 

 

 

 

 

 

— P — 4

2/? + 5

 

 

72 + Р + 1 ; р 2+ р + 1

§ 1.4. ПРОСТЕЙШИЕ СОЕДИНЕНИЯ ЛИНЕЙНЫХ СТАЦИОНАРНЫХ СИСТЕМ

Входные сигналы одной системы могут подаваться на вход другой системы. Таким образом, системы могут образовывать до­ вольно сложные соединения. В данном параграфе мы рассмот­ рим лишь самые простые из возможных соединений: .последова­ тельное, параллельное и встречно-параллельное соединения.

36

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