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

книги из ГПНТБ / Видаль П. Нелинейные импульсные системы

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

Риг. 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 <p1

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(z1).

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

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