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

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

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

(рис. 1.16,а). Здесь уместно сделать следующее замечание. Пе­ редаточной функции вида (1.41) соответствует также оператор, зада ваемый соотношением

аяУ{п1(*) +------1- «іУ 0 it) + а0у (t) = bQx (t) -1-----+ bmxW(t).

(1.4 4)

Этот оператор отличается от оператора (1.42) —(1.43) тем, что сначала осуществляется m-кратное дифференцирование входного сигнала x(t), а далее определяется решение дифференциального уравнения п-го -порядка с постоянными коэффициентами а0, а\,.., ап. Этот оператор определен лишь для входных сигналов, имею­ щих т производных, которые удовлетворяют в нулевой точке ус­ ловию х(0) =0, ... , л("‘-Ч(0)=0. Переходной функции этот опера­ тор не имеет, так как входной сигнал вида x(t) — l(t) не являет­ ся дифференцируемым. Для входных сигналов, дифференцируе­

мых пг раз, оба оператора совпадают, т. е.

 

 

 

A t X (t ) = A2x(t),

(1.45)

где А,

— оператор,

задаваемый соотношениями

(1.42) —

 

(1-43);

.

 

Ач — оператор, задаваемый соотношениями (1.44);

 

x(t)

— функция, имеющая пг производных.

 

а)

Р и с . 1.16. Переходные функции

Докажем (1.45). Применяя преобразование Лапласа к левой и правой частям соотношений (1.42) — (1.43), с учетом того, что z (0) =0, . .. , *1 (0) = 0, получим

2 (р) = ---- —---- ------

; X (р),

апРп + - ' ■Jr aiPJr ao

Уі [Р) = ФтРт1-------

'rbxp + bü)z{p).

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

 

 

w -

 

х\ р ) .

апр п + - - - - г а хр

+ а0

47

Далее, применяя преобразование Лапласа к левой и правой частям соотношения (1-44), с учетом того, что *(0) =0 ..., _\;(m-l)(0) =0, получаем

У2 (Р) = КгРт + - - - +

ъіР + bp д р)

ап Р" + • • • + Oj р + Q-0

В таком случае У\(р) = у*(р)

и, следовательно, Уі(Ь)=Уі (і).

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

W (р)= ^ ^ в основном соответствует оператор вида (1.43) —

А (р)

(1.44), Примером этого является оператор (1.27), описывающий

передающие свойства 7?С-цепочки

с передаточной

функцией

W(p) —— ——

И, наконец, отметим

следующее. Если

ввести

Тр + 1

 

 

anD + ■• • +

D +

а0 и

дифференциальные операторы A (D) =

В (D) = bmDm+

• • ■+ b jD + é (/ где D —

 

то соотношения

(1.42) к (1.43) и соотношение (1.44)

 

могут

быть

формально

записаны

в виде:

 

 

 

 

 

 

A(D)z{t) = x(t)-

 

 

 

 

 

 

 

A (D)y (t) = В (D)x(i).

 

 

 

y(t) = B(D)z(t);

 

 

 

 

 

Применяя формальное деление на A(D), для обоих случаев по­

лучим

 

 

 

 

В(Р)

 

 

y { t ) =

^ )x { t ) = W ( D ) x { t y ,

W (D) =

 

(1.46)

 

 

 

 

 

A (D)

 

 

Отчетливого математического смысла выражения (1.46) не имеют и являются по существу символической записью опера­ торов вида (1.42) —(1.43) и (1.44). Однако такая запись приме­ няется в литературе, причем символ W(D) принято называть пе­ редаточной функцией системы в операторной форме.

Весовая функция g(t) линейной стационарной системы

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

В силу данного определения изображение g(p) весовой функ­ ции g(t) линейной стационарной системы с передаточной функ­ цией W(p) будет равно:

g ( p ) - W (p ) L [ b { t) \ = W(p).

(1.47)

48

В таком случае весовая функция является обратным преобра зованием Лапласа от передаточной функции системы

g{t) = L-'[W{p)).

(1-48)

Если левую и правую части соотношения (1.48) разделить на р и взять от них обратные преобразования Лапласа, то получим

Ь~1 g ( p ) — =

L

W(p) —

[

р

 

L

р J

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

L- 1

 

 

g (■=) <*•

 

 

 

О

 

Учитывая ■-1 Щ р ) —

= h(t)t получаем

 

р

 

/

 

 

 

 

 

 

А (0 =

J g[*) di,

(1.49)

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

Дифференцируя (1.49) по верхнему пределу, получаем, что весовая функция является производной от переходной функ­ ции

g ( t ) - d- ^ P .

(1.50)

at

 

Весовая функция g(t), как уже было указано, может быть найдена как оригинал от передаточной функции g(t) =L~l[W(p)].

П р и м е р . Определить весовые функции g\(t) и g2(t) ей-

Trt

стем с передаточными функциями

 

W x{p) =

----- — и W 2(p) =

__

к

 

 

 

Тр + 1

 

 

 

 

~

тр + Г

 

 

 

t__

 

Г тр 1

 

 

 

 

I fl

 

1

т .

 

Тр +1

 

 

Тр + 1.

 

 

к

k

г

 

 

Тр + 1

— е

 

 

Т

 

 

 

Графики этих функций изображены на рис. 1.17,а и б.

4. Изд. № 5312

49

Кроме того, весовая функция может быть найдена как про­ изводная от переходной функции h(t), определенной методом ре­ шения дифференциального уравнения. Найдем этим методом

весовые функции gi(t) и gz(t) данных в последнем примере си­ стем. В соответствии с (1.40) имеем

 

 

 

A

t

 

 

gi Ѵ) =

~ і г Аі [t)

I (Oe“

= 8(^)e

1(0 =

 

at

 

dt

 

t

 

 

 

 

 

 

 

 

 

 

 

 

T ;

 

g2( * ) - - £ - M O

 

 

k

г

 

 

 

При вычислении g\(t)

путем дифференцирования следует учиты­

вать разрыв, который терпит функция h{(t)

в точке 0=0, так как

,

t <

0;

 

 

 

А і(0—

1 е

 

(рис. 1.16,а).

 

 

 

т £>0

 

 

 

Р и с . 1.17. Весовые функции

Отметим, что реакция g(t) всех реальных систем на 8-функ­ цию 8(£), подаваемую на вход в момент времени t=0, отлична от нуля лишь при значениях времени t, больших нуля,

О

t < 0;

Ф 0

(1.51)

t ]> 0.

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

В силу стационарности реакция линейной стационарной сис­ темы на входную функцию, равную 8-функции, поступающей на вход в момент времени t, будет равна весовой функции, смещенной на время -с, g (t —т).

50

Ясно, что для физически осуществимых систем эта функция должна удовлетворять соотношению g{t — т) = 0 при t «£т.

 

 

 

 

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

 

 

1. Определите

 

весовые

функции систем

с передаточными

функциями

(2р + 1Р

2

 

 

5р + I

 

 

 

------- --------------- — ,

— ——— и нарисуйте их графики.

 

(Зр+1)2Д (/>2+4)/;

10/7+1

 

 

 

Ответы:

 

 

 

 

 

 

 

г(0 -

1

1

te

-сг

5

...

1

cos2^J;

 

 

 

g(t) = — [! -

27

t

io

20

2. Определите переходные функции этих же систем и нарисуйте их графики

§ 1.6. РЕАКЦИЯ ЛИНЕЙНОЙ СТАЦИОНАРНОЙ СИСТЕМЫ НА ВХОДНОЙ СИГНАЛ ПРОИЗВОЛЬНОГО ВИДА.

ИНТЕГРАЛ ДЮАМЕЛЯ

Пусть на вход линейной стационарной системы с передаточ­ ной функцией W(p') с момента времени t=0 поступает входной

сигнал x(t)

произвольного

вида, удовлетворяющий условиям

|х(()І < С

и x (t) = 0

при t<^0. В таком случае изображе­

ние у(р) выходного сигнала y(t) существует и связано с изобра­

жением х(р)

входного сигнала соотношением

 

 

y(p) = W(p) х(р).

(j.52)

Учитывая,

что L~'[W(p)]=g(i), по теореме о свертке из (1.52)

получаем

s

 

 

 

 

У(0 = $ g (f - Х) * ( Х)Ж- '

(1-53)

 

о

 

Соотношение (1.53), называемое интегралом Дюамеля, задает связь между входным и выходным сигналами с помощью весовой функции g(t), которая выступает здесь в роли весовой функции интегрального оператора. Таким образом, оператор любой ли­ нейной стационарной системы может быть задан в виде интег­ рального оператора.

П р и м е р . Оператор линейной стационарной системы с ле-

к

редаточной функцией W (р)=——— - может быть задан как диф­

ференциальным уравнением

Ту(() + у ( 0 = hx{t) у (0) = 0,

4*

51

гак и интегральным оператором

t

/ — X

е т x(t)dx.

о

Последнее соотношение вытекает из того, что

k - V

g { t ) = L -

Тр + 1 J

Следует заметить, что интеграл Дюамеля иногда представля­ ют в несколько иной форме, чем (1.53). Соотношение (1.52) мо­ жет быть представлено в виде:

y(P) = —

W ІР) \рх (Р) — *о] + — w (Р)х (0).

Р

Р

Применяя к последнему выражению теорему о свертке и учи­

тывая, что L-1 і Щ /?)— = h [ i ) и L~l \ рх ( р) х 0] = X (t)i

L

Р

получаем

 

 

(1.54)

Последняя форма записи интеграла Дюамеля удобна в тех случаях, когда известна переходная функция h(t) системы. Функ­ ция h(t) может быть, например, непосредственно получена из эксперимента путем измерения реакции системы на единичный входной сигнал. Весовая функция g(t) непосредственно из экс­ перимента получена быть не может, так как о-функцию нельзя представить в виде реального сигнала любой физической приро­ ды.

В дальнейшем будем пользоваться лишь представлением интеграла Дюамеля в форме (1.53).

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

X (t) = 0 при t<C_ t0.

Если to>0, то условие х(і) =0 при /< 0 выполняется, и в силу формулы (1.53) имеем

t

і

У ( 0 = j g V — т)-*(х) dz =

j g(t - т) * (x) dx-

0

(0

52

Изменение нижнего предела интегрирования здесь следует из того, что на отрезке интегрирования [0, /0] функция х(х) =0. Ес­ ли ^о<[0, то функция х(1) не является оригиналом и не имеет изо­ бражения по Лапласу. Для того, чтобы обойти эту трудность, рассмотрим вспомогательную функцию x(t + t0), задержа;нную по отношению к входному сигналу на время (10) (рис. 1.18). Функ­ ция x(t+t0), рассматриваемая как функция переменного t, удовлетворяет условию

x(t + tü) =0 при /!<0.

Рис. 1.18. К выводу формулы интегра­ ла Дюамеля

Выходной сигнал y\(t), являющийся реакцией на входной сигнал x(t+t0), в соответствии с (1.53). будет равен:

t

Уі(*) = J g [t x)x (t + t0) dx.

0

В силу стационарности системы выходной сигнал y(t), соответ­ ствующий сигналу x(t), будет равен сигналу y t (t), взятому в момент времени .т = t — t0 (см. рис. 1.18).

і- Іо

y ( t ) = y i ( t — to) = f g (t - t0 - *)x (T + tQ)di-

6

Путем замены переменной интегрирования х -)-£0= т ' в интеграле

получаем y(t) = J

dx'.

Опуская штрихи в по-

^0

 

 

следнем равенстве, окончательно имеемУ

 

У (*) =

J g [t — *)x(x)dx.

(1.55)

Итак, входной x(t) и выходной y(t) сигналы связаны соотноше­ нием (1.55), справедливым для произвольного момента времени to начала поступления входного сигнала. Это соотношение, обоб­ щающее формулу (1.53), также будем называть интегралом Дюамеля.

53

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

1. Весовая функция системы g ( t ) = 2 е~ 2^І (/). Найти значение выходного сигнала этой системы в момент времени t=2, если на вход системы подается

сигнал x(t) = 2е4'. При решении задачи передаточной функцией не пользо­ ваться.

Ответ:

[es ~ e - 4l.

 

 

3

1

 

 

 

3р +

1

2. Передаточная функция системы равна W(p) =

о. Найти опера­

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

§ 1.7. ВЫХОДНОЙ СИГНАЛ ОДНОМЕРНОЙ ЛИНЕЙНОЙ СТАЦИОНАРНОЙ СИСТЕМЫ С НЕНУЛЕВЫМИ НАЧАЛЬНЫМИ

УСЛОВИЯМИ

Пусть одномерная линейная стационарная система имеет пе­

редаточную функцию W \р)= ^

и задается оператором вида:

 

А [р\

Г а)

a„zW{t)\- an_xz^-V{t) +

■■■± a lzM{t)+a0z(i)-=x(t)m

I б)

y{t) = bmz ('mHt) +

-------1-bxz^(t)-\-bQz{t)

 

 

(1.56)

Пусть далее начальные условия дифференциального уравнения не обязательно являются нулевыми.

2(ч-і) (0) = *«-»; . . . . 2W(0) = г>0, z (0) = z0°-

Начальное условие по выходному сигналу у ( 0) в таком случае имеет вид:

y(°) = bmz™+ £m_ i2 '" - 4 -------

г M o 1 + bnz 0°.

Преобразовав по Лапласу с учетом начальных условий диффе­ ренциальное уравнение (1.56,а), получим

zip)

, A~l {p)

о J _

[ а - ^ - ' Ч р ) т П—\

>

zip)--

+

 

А (р)

0

А(р)

А (р)

 

 

 

где

 

 

 

 

 

А ~ 1(р) = апрп-'

Н-------і-агр + щ;

 

 

А ' Чр ) ~ апРп~'

Н-------

\-*і+іР + а»

 

 

А -(П-Х) {р) = дя.

 

 

 

 

54

Из дифференциального соотношения (1.56,6) получаем

 

У (Р ) = В ( р ) г (р) - 5 - 1(р) z0° ---------

(^) zm-it

где

В-'[р) =

Ътря-1-\-------

\-Ь2р + Ь »

 

 

ß-(m-l) _ _

£ - m =

ß-{m+ ] ) =

^ - ( n - l ) = Q .

В таком случае

У(Р) = ~ ^ ^ ( Р } +

В(Р)А-Ңр) - В - Ч р )

А(р)

 

А(р)

+

 

В(р)Ап~1(р) 3 ,я_

гох + ■■■+

. >1 (/»

 

А{р)

*о° +

4P) гП—1

(1.57)

На основании теоремы об изображении свертки имеем

t

П— 1

 

У (0 = J s ( t -

dl + £ф,(*)го0’

(1.58)

О

/='о■

 

где

В(р)А-ѵ+» - В - ^ Ц р )

А(р)

Структурная схема системы с учетом ненулевых начальных ус­ ловий имеет вид (рис. 1.19).

Р и с. 1.19. Структурная схема системы с учетом начальных условий

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

55

П р и м е р . Найти выходной сигнал ^С-цепочки (см. рис. 1.6) при условии, что на вход системы поступает сигнал х(і) = U\(t), а конденсатор С в начальный момент времени заряжен до на­ пряжения ди­

электрические процессы в этой цепочке описываются соотно­ шениями

u2{t) =

Ri{t);

iil [t) = ivc{t) + u2{t);

 

 

t

 

t

uc(t)= uc0 +

 

г(т) d^ = ucQ+ £ ^ - |и2(т)</т =

 

и

 

 

О

 

t

 

 

= а

+

u2(x)dX ,

где а = 1/CR.

а

о

 

 

 

 

 

 

 

It

t

 

 

 

р

u„J)d^=z{t)

то эти соотношения могут

Если обозначить— _j_ J

а

о

 

 

 

быть записаны в виде:

z( l)+a z(t)=ul (t); z (0) = z0= иJ a ;

u2(t)=z(t).

 

 

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

 

pz(p)—z( 0) + az(p) = и+р);

 

u2(p)=pz(p)—z(0).

 

Далее имеем

 

 

и2(р) = р '«1 (P) 1 zo

-Z o = —^— ihiP) + [

P

p + a p + a \

P + a

P +0-

ll «о- J

(1.59)

Мы видим, что выражение (1.59) полностью соответствует об­ щей формуле (1.57), так как в данном случае

В(р)=р, В~'(р) = \, А( р)=р+а, А~Ңр) = 1.

Окончательно получаем

 

 

- z.-1

а

Ue~at azQe~at—

ІР+ a\

P +a + (p) J

 

= ІУ -

ил ) е~а1.

 

56

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