книги из ГПНТБ / Основы авиационной автоматики учеб. пособие
.pdf(рис. 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
