книги из ГПНТБ / Видаль П. Нелинейные импульсные системы
.pdfРиг. 2-4. Оригинал мо дифицированного г-пре- образования.
|
|
|
___I |
I I L— |
||
|
|
|
о |
т гт зт 4т |
||
откуда |
|
|
1_ р-пт |
|||
f (пТ, т) |
= е ~ |
|||||
тТ |
|
(2-44) |
||||
На этом примере |
видно, |
что |
f(nT, |
tri) |
характеризует разрывы |
|
в моменты квантования |
(рис. |
2-4). |
|
|
|
2-4. ПРИЛОЖЕНИЕ К ИМПУЛЬСНЫМ СИСТЕМАМ
Широкое использование метод г-преобразования на ходит при исследовании линейных импульсных систем.
а) ИМПУЛЬСНАЯ ПЕРЕДАТОЧНАЯ ФУНКЦИЯ
Рассмотрим систему с передаточной функцией G(p) и импульсными элементами, включенными так, как это показано на рис. 2-5. Определяем импульсную переда точную функцию как отношение G(z)=S(z)/E(z), при чем предполагается, что выходной сигнал также кван-
Рис. 2-5. Импульсная си- |
Рис. 2-6. Замкнутая импульсная си |
стема. |
стема. |
туется с той же частотой некоторым фиктивным импульсным элементом, работающим синхронно с дей ствительным. Найдем импульсную передаточную функ цию замкнутой системы, показанной на рис. 2-6.
71
Получаем:
G,e (z) = GtE (2) — G, s (2) GZHG1(2);
S (2) = Gl e (2) G2 (2),
или
G2 (z)G,£(z)
1+ G2HG1 (г)
G i = 1 и G2= G, то
S(z) _ |
G,(z) |
E G) |
1 4 . GH (z) - |
(2-45)
(2-46)
Символ над GH(z) означает, что ^преобразование берется для произведения функций GH. Именно так и следует поступать, а не использовать произведение z-преобразований G и Я.
б) МОДИФИЦИРОВАННАЯ ИМПУЛЬСНАЯ ПЕРЕДАТОЧНАЯ ФУНКЦИЯ
Модифицированная импульсная передаточная функ ция, соответствующая системе, изображенной на рис. 2-5, определяется как отношение G(z, m )=S(z, m)/E(z). Предполагается, что выходной сигнал как бы пропу скается через элемент запаздывания (1—т)Т и затем уже квантуется фиктивным импульсным элементом.
Таким же образом находится импульсная передаточ ная функция замкнутой системы, изображенной на рис. 2-6:
5 (г> т) = |
М ?. |
m)G,F(z) ^ |
(2-47) |
|
|
1 + |
G2//G , (г) |
|
|
В частном случае, когда G i= l, G2= G , |
это соотноше |
|||
ние запишется так: |
|
|
|
|
S (г, т) __ G (г, т) |
(2-48) |
|||
Е (г) |
1+ |
GTI (г) |
||
|
2-5. АНАЛИЗ УСТОЙЧИВОСТИ ЛИНЕЙНОЙ ИМПУЛЬСНОЙ СИСТЕМЫ
Главной задачей при исследовании линейных систем является отыскание условий устойчивости. Для линей ных непрерывных систем задача сводится к тому, чтобы корни характеристического уравнения F (р) были распо-
72
ложены в левой полуплоскости р. Эти необходимые и до статочные условия были развиты для импульсных систем, описываемых с помощью метода ^-преобразования. Шур [Л. 2-10], Кон {Л. 2-11] и Марден [Л. 2-12] доказали, что корни устойчивого характеристического дискретного уравнения F(z) на плоскости z должны располагаться внутри окружности единичного радиуса, Э. И. Джури [Л. 2-13] предложил также упрощенный метод. Кроме того, были развиты также и геометрические критерии, являющиеся обобщением критерия Найквиста.
Мы укажем только лишь на критерии, нашедшие наи большее практическое применение.
а) АЛГЕБРАИЧЕСКИЕ КРИТЕРИИ УСТОЙЧИВОСТИ
1. Критерий Шура — Кона
Дискретная линейная система устойчива, если лю бому ограниченному входному сигналу соответствует ограниченный выходной. Импульсная система, изобра женная на рис. 2-6, устойчива, если корни ее характе ристического уравнения
F{z) = \+ G ^ G l{z) |
(2-49) |
по модулю меньше единицы (за исключением случая скрытых колебаний ‘ ) .
Критерий Шура — Кона позволяет узнать, будут ли
у полинома |
(2-50) |
F (z) = йо4"UiZ-j- ... + anzn |
корни по модулю меньшими единицы.
Для этого образуют определитель порядка 2/г
а0 |
0 |
. |
. |
a t |
а0 . |
. |
. |
0 я» ап - 1 •
0 0 ап . * * ап - ft “ 2
я*- 1 |
йЛ-2 |
• |
• а0 0 Я г |
• ■ « « |
(2-51) |
“ п |
0 |
. |
. 0 «0 « 1 |
• - “ f c - 1 |
|
в „- h+1 f i n - f c + Z • |
• |
• « л 0 0 |
• ■ в» |
|
k= 1, 2, ..., п; ап сопряжено с ап.
1 Под устойчивостью системы понимается устойчивость дискрет ной последовательности sn.
73
Критерий. Полином F(z) имеет все корни, меньшие единицы, если Д*<0 при к нечетном и А/г> 0 при к чет ном.
Пример 2-8. Рассмотрим импульсную систему, представленную на рис. 1-9. Она может быть представлена в виде классической структурной схемы рис. 2-7. Импульсная часть системы характери зуется тем, что выходной сигнал в промежутках между моментами п'Г и (п+ 1)Т постоянен и равен nkren (рис. 1-10). Система состоит из интегратора с последовательно включенным фиксатором нулевого порядка, а ее передаточная функция будет:
|
|
|
|
|
1 |
<?-? _ |
z — 1 |
1 |
|
|
|
|
|
|
__ |
_ _ |
_ |
_ , |
|
|
|
|
|
е. |
|
|
|
|
|
Рис. 2-7. Структурная схема |
|
Х(Р) |
kr |
Z — |
I |
1 |
|||
|
к (Р) |
г |
р2' |
||||||
линейной импульсной системы. |
|
|
|
|
|
|
|||
По табл. 2-1 определяем импульсную передаточную функцию |
|||||||||
IX (2) _ |
2 — |
2 |
|
k. |
|
|
(2-52) |
||
К |
(2 ) |
|
2 |
( 2 — |
1)2 |
Z — |
|
|
|
При импульсной передаточной функции объекта, равной к0, по |
|||||||||
лучаем: |
|
|
|
|
|
|
|
|
|
|
|
Y (г) |
k0kr |
|
|
|
|
(2-53) |
|
|
|
^ W = |
(7=TT |
|
|
|
|||
Характеристическое уравнение рассмотренной системы запи |
|||||||||
шется: |
|
|
Mr |
|
|
|
|
|
|
|
F(2) = 1+ |
+ 2 — |
|
|
|
(2-54) |
|||
Применение критерия Шура — Кона к полиному |
|
|
|
||||||
дает: |
|
|
z+ k — 1= 0 |
|
|
|
|
|
|
|
ао=й—1; а( = 1; |
|
|
|
(2-55) |
||||
|
|
|
|
|
|||||
|
Д = |
а0 |
ах |
= а02 — af < 0, |
|
|
|
||
|
о , |
б 0 |
|
|
|
||||
|
|
|
|
|
|
|
|
||
|
или |
(к— 1) 2< |
1, т. е. k<2. |
|
|
|
(2-56) |
||
Это соотношение было найдено при исследовании системы с по |
|||||||||
мощью разностных уравнений. |
|
|
|
|
|
порядка |
|||
З а м е ч а н и е . |
Если |
бы |
объект был звеном первого |
||||||
с постоянной времени т, то уравнение (2-53) |
было |
бы |
2-преобразо |
||||||
ванием от ( l + t p ) -1. |
|
|
|
|
|
|
|
|
2. Модифицированный критерий Рауса
Критерий Рауса позволяет узнать, являются ли все действительные части корней полинома отрицательными.
74
Эта теорема применяется к /преобразованному с помо щью билинейного преобразования уравнению (2-50)
1+ w 1—w
Уравнение (2-50) запишется тогда так:
F ( w ) = w n + b, w n~ ... + b n- iw + bn.
Образуем определитель
ь3 i |
о |
о |
о |
. . . |
о |
ь3 ь2 ь, 1 о |
. . . о |
||||
Ъь bt |
b3 |
Ь2 |
Ьг |
1 0 |
0 |
о ................... |
|
|
|
о |
ьп |
(2-57)
(2-58)
(2-59)
Все действительные части корней F(w) отрицатель ны, если все миноры А; положительны:
|
ь, |
1 |
Д |
Ь, |
1 |
о |
|
&i = bx; А3 = |
ь3 ьг |
ъх ... |
|||||
|
|
||||||
|
|
, |
64 |
= |
|||
|
|
|
|
Ьь |
Ъъ |
Пример 2-7. Вернемся к примеру, приведенному в предыдущем
параграфе: |
и (2 —k) -\-k |
|
1 -f- w |
(2-60) |
|
F (да) = 1—ш + k— |
1 — w |
Следовательно, нужно, чтобы й<2.
3. Упрощенный критерий Э. Джури
Э. И. Джури [Л. 2-13] предложил в 1962 г. упрощен ный критерий, позволяющий получить условия устойчи вости для систем второго, третьего, четвертого и пятого порядков.
1. Пусть
F(z) = a 0 + aiz + a2z2-, а2>0. |
(2-61) |
Условия устойчивости записываются в виде нера венств
\а01< а2 |
или а0 — а2 < 0; |
'j |
|
|
«о + я. + |
аа> 0 ; |
/=■(!)>0; |
|
(2‘62) |
ай— ал+ « 2> 0 ; |
F (— 1) > 0. |
J |
|
|
2. Пусть |
|
|
|
(2-63) |
F(z) = а 0 + aiZ-\-azz2 + a3z3-, а3> 0, |
75
тогда |
|
|
|
|
|
F( 1) > |
0; |
|
|
|
|
F |
( — 1) < |
0 ; |
|
(2-64) |
|
I |
I <~C ^s> |
|
|
||
|
|
|
|||
a\ — a\<C. a0a2~ |
aias- |
, |
|
||
3. Пусть |
|
|
|
|
|
|
|
|
|
(2-65) |
|
F(z) = a 0 + aiz + a22:2 + a32:3-rfl424; a4>0, |
|||||
тогда |
F (1 )> 0 ; |
|
|
|
|
|
|
|
|
||
|
F (— 1) > |
0; |
|
|
|
al — a4 — аоаз“b aia4 |
0; |
|
|||
<1,2 |
“f" <70йз |
|
Tliai |
0) |
( 2-66) |
a\ + 2a0a2a4 + a,asa4 — a0a4 — a2a4 — a0a| — |
|||||
— a?0a4 — Лд a2 — |
+ |
a4 + |
a0a,as > |
0. |
|
|
|
|
|
|
) |
б) ГЕОМЕТРИЧЕСКИЕ КРИТЕРИИ УСТОЙЧИВОСТИ
1. Обобщение критерия Найквиста
Критерий Найквиста, предложенный как метод иссле дования устойчивости линейных непрерывных систем, применим и для дискретных систем. Задача состоит в построении годографа функции
|
G(z) = GiHG2(z) |
(2-67) |
в плоскости z, если полагать, что |
|
|
z = e ia; 0 < ш < ^ . |
(2-68) |
|
Для того чтобы система была устойчивой, нужно,что |
||
бы число оборотов |
1+ G(z) вокруг начала |
координат |
[или число охватов |
годографом G(z) точки — 1] было |
равно количеству неустойчивых полюсов функции G(z).
Пример 2-8. Рассмотрим еще раз систему, изображенную на рис. 1-12, для которой
k |
_______k_______ _____ k_____ kj sin <f |
||
^ = z — 1 — |
cos <p-(- / sin <p— 1 |
z |
2(1 — cos <p) ’ ' |
76
Достаточно нескольких точек, чтобы |
|
|
|||
определить кривую |
Найквиста |
(рис. 2-8): |
|
||
г = 1; |
G (г) == оо |
|
|
|
|
2 = — /; G (г) = ~- / -1 |
|
|
|
||
|
|
k |
|
|
|
■ Г , |
G (г) = 7/ - 1 |
2 ( ! + / ) ; |
|
|
|
|
G (г) = |
"o'- |
|
Рис. 2-8. |
Кривая |
|
|
|
|
||
Кривая Найквиста является |
прямой и |
Найквиста |
для им |
||
система |
устойчива, |
если/5:<2. |
|
пульсной |
системы |
|
|
|
|
рис. 2-7. |
|
|
|
2. |
Метод |
корневых годографов |
Этот важный метод, используемый при синтезе непре рывных систем, находит применение и для дискретных линейных систем. Для того чтобы узнать, устойчива ли система, строим геометрическое место корней, т. е. гео метрическое место корней 1+ G(г) = F {z) =0, при k, изменяющемся от 0 до оо, и определяем значения к, для которых это геометрическое место расположено внутри окружности единичного радиуса. Правила те же, что и
для непрерывных систем, а именно: |
уравнения |
||||||
1. Каждому |
корню характеристического |
||||||
F(z) = 1 + G (z) |
соответствует ветвь годографа. |
|
|||||
2. Годографы корней для d—п асимптотических на |
|||||||
правлений {исходя |
из степени знаменателя |
G(z), п — |
|||||
числитель G(z)] задаются следующим образом: |
|
||||||
argA = |
^ |
( l + |
2 Я), Я = 0, 1, 2... |
(2-70) |
|||
Асимптоты |
пересекают |
действительную |
ось в |
точке |
|||
с абсциссой |
|
|
|
|
|
|
|
ър - |
ъг |
|
I Р — полюсы G(г); |
\ |
(2-71) |
||
d — n |
’ |
Где \ |
г — нули G(г). |
| |
|||
|
3. d ветвей годографа исходят из d полюсов G(z), и
вэтой точке имеем k— 0.
4.d ветвей годографа стремятся к п нулям G(z),
a d—п — к бесконечности, и для этих точек £ = оо.
77
5. Годографы, состоящие из участков действительной оси, соединяющей полюсы или нули G(z), начинаются
вполюсе или нуле, лежащем правее.
6.Исходные точки на оси — это кратные корни k, получающиеся при решении
d |
' N (а) |
—[О, где G(z) |
N (г) |
(2-72) |
da |
D (а) |
D(z) |
7.Пересечение с мнимой осью соответствует чисто мнимым корням характеристического уравнения.
8.Годограф определяется и вычисляется с помощью двух соотношений:
|
(Z--г,) (г —za) ... |
1 |
||
|
(г —Р\) |
(z — Pi) ••• |
k |
|
|
|
|
|
(2-73) |
и |
|
|
|
|
П |
|
|
d |
|
£ |
afg (г ~~ гз) |
£ arg ( z - |
Pi) = |
|
/=1 |
|
|
/=1 |
|
Рис. 2-9. Корневые годо- |
= |
7t —f—2 Ян. |
|
|
графы G(z)=kl(z—1). |
k |
|
|
|
Пример 2-9. Пусть G (г) = |
|
|
|
|
-р. |
|
|
|
|
Корневой годограф состоит |
из участка |
действительной оси |
||
(рис. 2-9). Зависимость от k определяется условием |
|
|||
(г — г ,)(.г — г2) ... |
_1_ |
' |
|
|
(г — р,)(г — Рг) ... |
k |
|
Система устойчива, если k<2.
2-6. СООТНОШЕНИЯ МЕЖДУ г-ПРЕОБРАЗОВАНИЯМИ И ЛИНЕЙНЫМИ РАЗНОСТНЫМИ УРАВНЕНИЯМИ
Этот параграф содержит краткое исследование раз ностных уравнений с помощью z-преобразования. Мы ограничимся несколько поверхностным изложением зада чи, так как целый ряд публикаций [Л. 2-3, 2-4] был по священ этой проблеме.
78
а) МЕТОД РАЗНОСТНЫХ УРАВНЕНИИ
Этот метод очень полезен для практического исполь зования, если располагают вычислительной машиной. Он позволяет выявить связь, существующую между 2-преоб разованием и разностными уравнениями. Эта связь ана логична связи, существующей между преобразованием Лапласа и дифференциальными уравнениями.
Найдем оригинал
(2-74)
Задача заключается в разложении N (z) и D(z) в ряд по степеням z~l.
Определение коэффициентов при соответствующих членах с определенной степенью z~l позволяет получить разностное уравнение, соответствующее F(z). Это урав
нение решается так же, как и |
разностное |
уравнение. |
Оно позволяет последовательно определить f(nT). |
||
Пример 2-10. Найдем оригинал |
|
|
Tz |
А/ (г) |
(2-75) |
' ( г> = 1 1 = 1 )Г = т ф г |
Это уравнение можно записать следующим образом:
г2Л/ (г) — 2zN (г) + |
N (г) = |
TzD (г), |
|
|
т. е. |
|
|
|
|
г2 [N (0) + N (Т) z - ' + N (27) z ~ 2 + |
...] — 2z[iV (0) +iV (7 ) z - >+ |
...] + |
||
+ W(0) + A A ( 7 ) z ~ ‘ - f ... = |
Tz [D(0) + D ( 7 ) z - 4 - ...], |
|
||
или |
|
0; |
|
|
H(0) = |
|
|
||
JV(7) — 21V (0) = TD (0); |
|
|||
N (27) — 2N (7) + |
2N (0) = |
TD (7); |
|
|
N(nT) — 2N[(n — 1) 7] + |
N(nT) = |
TD[(n + 1 ) 7 ] . |
|
|
Выражению (2-75) соответствует разностное уравнение |
|
|||
N[ (п+2) 7]—2Л7[(я +1) 7] ■+ N (пТ) = TD[ (я + 1) Г]. |
(2-76) |
В зависимости от величины D{t) последовательное решение этого уравнения дает решение задачи. Отметим, однако, что соотношение (2-76) годится для любого значения я, если только 7<0, D(t) и N(t) равны нулю.
Интересно отметить, что разностное уравнение (2-76) очень про сто получить, исходя из уравнения (2-75), если рассматривать z как оператор смещения (т. е. z{E(z)— >-£[(«+()7]).
79
б) МЕТОД z-ПРЕОБРАЗОВАНИЯ
Этот метод состоит в решении линейного разностного уравнения с помощью z-преобразования. Рассмотрим разностное уравнение
У! Щ,Хп+1 ~~ У] ЬгУп+и |
(2-77) |
|
1=0 |
1=0 |
|
где а{ и bi — константы; |
уп — возмущающая |
функция, |
известная для всех значений я^гО.
Задача состоит в отыскании хп как функции параме тра пТ (я-положительное целое).
Применим z-преобразование к уравнению (2-77); с помощью соотношения (2-6) получим:
p |
t = l |
|
я |
|
|
i=l |
2 CLiZ1 |
X(z)-- 2 ^ 2_j |
— 2 |
h zi |
Y(z) -T > yiz ~j |
||
/=0 |
1=0 |
|
0 |
|
/=0 |
|
|
|
|
|
|
|
(2-78) |
где X(z) |
и Y(z)— z-преобразования |
х, у\ |
Xi и г/,- — на |
|||
чальные условия. |
|
и получаем: |
|
|||
Решаем уравнение (2-78) |
|
|||||
|
Я |
|
р |
/ 1 = 1 |
|
\ |
|
Y (г) У] btz! + |
У] atzi ( У] |
Xjz~1 |
— |
||
X(z) = |
|
|
/—о |
|
||
|
|
|
|
|
||
|
я |
|
i=0 |
|
|
|
|
|
г=1 |
\ |
|
||
|
■У] Ьгхг |
2 |
Уз*-’ |
|
|
|
|
t = 0 |
|
/=о |
|
|
(2-79) |
|
|
|
|
|
|
Тогда обратное z-преобразование позволяет опреде лить Хп- Этот метод обладает наряду с многими другими тем преимуществом, что начальные условия вводятся автоматически и нет необходимости искать постоянные интегрирования.
Пример 2-11. Возьмем приведенное выше разностное уравнение
х{п + 2) —2х ( я + 1) +х(п) =Ту(п+1). |
(2-80) |
Применив метод z-преобразования, получим: |
|
ггХ (z) —z2x (0) —zx (1) —2zX (z) + |
|
+ 4zx(Q) +X(z) ^TzY(z)— Tzy(0), |
(2-81) |
80