книги из ГПНТБ / Видаль П. Нелинейные импульсные системы
.pdfт. е.
z2x (0) — 2zx (0) + zx ( 1) — Tzy (0) z 2— 2z -I- 1
(2-82)
Если начальные условия нулевые, то получаем снова уравне ние (2-75).
2-7. АНАЛИЗ УСТОЙЧИВОСТИ НЕЛИНЕЙНОЙ СИСТЕМЫ
Метод свертки 2-:преобразований позволяет найти ре шение класса нелинейных дискретных систем, описывае мых нелинейными разностными уравнениями вида
YikPx(n |
р) A- F \х{п), х{п-\~ 1),... ,х { п + г — 1)J = 0, |
||||
|
|
|
|
|
(2-83) |
где п — дискретная |
переменная |
времени; |
Е[х(п), |
||
х ( п + 1), |
. . х{п + г— 1)] — нелинейный |
член, |
содержа |
||
щий по |
меньшей мере сомножители |
х(п + р) |
во второй |
||
степени. |
|
|
|
|
|
В этом параграфе, основанном на работах Э. И. Джу- |
|||||
ри и М. |
А. Пая (Л. |
2-6, 2-9, 2-14, 2-15], |
приближенное |
||
решение уравнения |
(2-83) получено на основании неко |
торых допущений с помощью свертки 2-преобразований.
а) ДОПУЩЕНИЯ
1. Состояние равновесия х0 нелинейного разностного уравнения (2-83) определяется из решения уравнения
Г
2 M o + F (-V х 0, . . . , х 0) = 0. |
(2-84) |
2. Нелинейная функция F[x(n), ..., х(п + г— 1)] отно сительно своих переменных х(п), ..., х(п+г— 1) аналитична в окрестностях точки равновесия (например, поли номиальная функция, часто встречаемая при исследова нии импульсных систем, содержащих элемент
снасыщением).
3.Точка равновесия асимптотически устойчива в том
смысле, что х(п) стремится к точке равновесия х0, если
п — »-оо.
6— 352 |
81 |
Главная особенность при решении разностного урав нения (2-83) заключается в предположении, что F(z) = = 'Z{x(n)\ — мероморфная функция вида
00
|
А |
(2-85) |
|
а |
1 |
||
|
|||
/=0 ■ |
|
||
где а; — различны и по |
модулю меньше единицы, |
а й0=1.
Задача тогда состоит в нахождении ^-преобразования функций типа Z[f(n + k)1} или Z[f (п+ z) hf (п+ q)1] и т. д.
Для решения этой задачи воспользуемся теоремой о свертке ^-преобразований.
б) СВЕРТКА В z-ОБЛАСТИ НЕКОТОРЫХ ФУНКЦИЙ
Предположим, что f(n) — функция |
времени, |
извест |
|||
ная для дискретных величин /г— 0, |
1, |
2, |
а ее г-преобра- |
||
зование имеет следующий вид: |
|
|
|
|
|
|
00 |
|
|
|
|
Z |
= Р (г) = |
при |
= 1 • |
(2-86) |
|
|
1=0 |
|
|
|
|
Выражение (2-20), являющееся интегралом свертки, |
|||||
позволяет записать: |
|
|
|
|
Z ff m = ~ \ p ~ lp{-p)F ( f ) dp- |
(2-®7) |
г |
|
где Г — контур, охватывающий все особенности p~iF(p). Согласно (2-20)
00
1=о
Таким же образом
Z\f(ny] =
и т. д.
Гi~ О
|
|
00 |
00 |
|
|
|
|
AiA, |
( 2-88) |
|
|
|
ataiz- |
|
|
|
|
|
|
|
|
i—0 fcO |
|
|
СО |
00 |
00 |
|
|
Ц V |
р |
\— aialamz - ' |
(2-89) |
|
|
1 |
|
||
(= 0 |
/= 0 |
m = О |
|
|
82
Так же легко получить (преобразование других функ ций, встречающихся в уравнении (2-83). Согласно фор
муле (2-6), являющейся теоремой о смещении г-преобра- зований:
1 [ f(nAr k y \ = z k [ Z [ f { n y ] - Y l |
[f(< 7 )j**- 9 ]. (2-90) |
||||
Таким же образом |
|
|
я—о |
|
|
|
|
|
|
|
|
Z[f(n)f(n + k)\ |
2ЩJ |
р ■X |
|
||
|
k-\ |
|
Г |
|
|
|
|
|
|
|
|
X F ( p ) - ^ f ( q ) p~q F |
dp, |
(2-91) |
|||
|
<7-0 |
|
|
|
|
где J — контур, охватывающий осооенности |
|
||||
q-'pb |
F (р) - 2 |
f {q)p~4 |
|
||
огкуда |
я-о |
|
|
|
|
ео |
оо |
|
h , , |
|
|
|
|
|
|||
Z[f(n)f(n + |
k)]: |
|
ai AtAi |
(2-92) |
|
|
1—ata,z~ 1 |
||||
|
|
|
|
i=o г=о
в) РЕШЕНИЕ НЕЛИНЕЙНЫХ РАЗНОСТНЫХ УРАВНЕНИИ
Рассмотрим следующий тип нелинейных разностных уравнений:
х(п-г2) +bx(n+ 1) + x(ti)[c-\-kx2{n)]=$u(n), (2-93)
где и(п) — единичная ступенчатая функция; b, с, k, § — постоянные коэффициенты.
Предположим, что
|
А, |
А й |
Аг |
|
i = 0 |
|
|
— а |
|
|
|
|
|
|
А, |
|
|
|
|
,2__1 |
|
3 , - j |
2, - , |
|
—а |
1—aU |
- a:,z |
|
|
|
|
1— ЛТЛа2Г |
■^12 |
+ . . . (2-94) |
■axa2z- |
|
<7,^2’ |
||
|
|
|||
и, с другой стороны, |
1 > |а, |> |
|а2 |
|
|
6* |
|
|
|
83 |
Члены Oi{i^ 3) |
имеют вид а“, |
а| (а |
и |
[3 — положи |
||||||
тельные целые числа, |
значения |
которых |
больше |
или |
||||||
равны 2). |
|
|
|
|
|
|
|
|
|
|
^-преобразование уравнения (2-93) запишется: |
|
|||||||||
X (г) (z2- f bz -\-c) + |
kZ lx3 (n)\= |
|
+ |
|
||||||
|
-f- (z2-f- bz)x{0) + |
гл; (1). |
|
|
||||||
Если Zi и г2 — корни уравнения гг+Ьг + с = 0, то |
по |
|||||||||
лучаем: |
|
|
|
|
|
|
|
|
|
|
Х(г): |
|
kZ[x3(n)\ |
| |
|
1 |
|
|
|||
(г — г,) (г — г2) |
1 (г — 2,)(г — г2) X |
|
||||||||
X |
(г _ 17 + |
(г3 + |
fe) |
|
(0) + |
^ ( 1 ) |
(2-95) |
|||
или |
|
|
|
|
|
|
|
|
|
|
Z [а-3 ( л )] : |
Л0 |
ЗА^А, |
|
3AqA2 |
|
|||||
|
|
1 — tfjZ' 1 |
1 |
1 — a2z~ |
|
|||||
3 A qA 2 0 - j - 3 A 0A j |
З А ^ А 30 + A ^ -)- 6Д 0A 1A 20 |
|
||||||||
|
-a,z |
1 |
|
|
1— a\z-' |
|
|
|
||
ЗЧдЛ40 -f- 3A20Ao+ 6/1 |
30 -f- 3a \a 23 |
(2-96) |
||||||||
|
|
|
1— a\z-' |
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
||
Заменяем |
X(z) |
его |
выражением, |
взятым априорно |
||||||
[уравнение (2-94)], |
a Z(x3(n)] — уравнением |
(2-96); |
что |
бы исключить неизвестные коэффициенты при z~h, доста точно приравнять члены с равными степенями z~h. Это
осуществляется |
путем |
последовательного умножения |
|
обоих членов |
на |
(z— 1), |
(z—at), (z—a2) ..., если пола |
гать z = l , Ci, |
..., |
fli. Коэффициенты A0, Au ..., A{ опре |
делены. Тогда возможно, учитывая нелинейность и на чальные условия, определить условия сходимости отно сительно х(п).
Отметим, что, вообще говоря, довольствуются первы ми членами разложения, что обусловлено, с одной сто роны, большой сложностью расчетов, а с другой сторо ны, приближенностью решения разностного уравнения. Возможности метода, теоретически пригодного для ре-
84
шения любых разностных уравнений (2-94), в действи тельности ограничены исследованием некоторых нелиней ных разностных уравнений первого или второго порядка.
г) ПРИМЕРЫ ПРИМЕНЕНИЯ
Пример 2-12. [Л. 2-12]. Рассмотрим нелинейное разностное урав нение первого порядка
|
|
|
|
&Уп + 1 А'^УпнУ^ — 8у^п — |
3. |
(2-97) |
|||
|
Предположим, |
что |
|
|
|
|
|
||
|
|
|
/ \ |
А и |
t |
А * |
I |
* |
(2-98) |
|
|
|
(г) |
1 — 2_i 4 "^ — агг ~ |
1 — a2z~ |
||||
|
|
|
|
||||||
|
z -преобразование уравнения (2-97) запишется: |
|
|||||||
|
|
8г \у (г) — г/0] + 8Z [у„+,у,г ] + |
8Z \угп ] = -J Z T |
(2-99) |
|||||
|
С учетом (2-88), (2-89) и (2-98) вышеприведенное уравнение |
||||||||
запишется так: |
|
|
|
|
|
|
|
||
|
|
|
оо |
|
г Е |
AjAiAmcim |
|
||
|
|
|
|
|
8 |
|
|||
|
|
|
/—о |
I—aialamz~x |
|
||||
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|||
|
|
|
|
оо |
|
,l,m —О |
|
|
|
|
|
|
|
А А |
|
3z |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
- |
1— aya^1 |
z — \ + 8zy°- |
( 2- 100) |
|||
|
|
|
Е1,1=0 |
|
|
|
|
|
|
на |
Сравнение |
результатов |
после |
последовательного деления |
на z, |
||||
(z— 1), .. ., |
(z—at) позволяет получить: |
|
|
||||||
|
|
|
|
А + А + А |
-(- ... = у0', |
|
|||
|
|
|
8А |
+ 8Л^— 8^ |
= 3, |
где А |
= 0,5; |
|
|
|
|
8А + 8а- 1 (А* А,а, + 2А* Аг) — 8« f 1 (2 А А ) = 0. |
|
||||||
|
|
А , (д, 4 - Лд а, + |
2Ад — 2А ) |
= 0 -> а, = 0,4; |
|
||||
|
А |
8 (A2Aq<z^ 4" 2АгАд 4" 2Л] |
А0а, |
4 - aiA ) — 8 (-Л|4-2АА) = 0 , |
|||||
т. |
е. А |
|
1 |
2 |
|
откуда |
|
|
|
= — - д - А\ и так далее, |
|
|
|||||||
F (z) = |
|
А _________ 1 |
^ |
(2- 101) |
|||||
|
0 ,4 2 -' |
3 1 — 0, Юг- 1 + •" |
|||||||
|
|
|
|
|
85
и |
|
уп = 0,5 + Л, (0,4)" - 1/3А\ (0,16)" + |
|
Уо — 0.5 + |
( 2- 102) |
Аг — 1/ЗЛ^ + ... |
|
Нужно отметить, что состояние равновесия будет при х=0,5 и |
|
что А, можно определить из |
уравнения (2-102) как функцию у(0). |
В зависимости от величины у(0) у (и) получаем из уравнения (2-102). Однако точное решение получить практически почти невозможно, что
придает методике скорее теоретический, чем |
практический характер. |
||
З а м е ч а н и е . Очевидно, что уравнение |
положения равновесия |
||
позволяет определить Ац. Действительно, принимая |
yn +i= yn, |
полу |
|
чаем уравнение. |
|
|
|
8 - - 8 у 2п Ч - % п = 3 , |
|
|
|
единственное решение которого уп = 0,5. |
у {г) |
уравнения |
(2-98) |
После разложения на простые элементы |
|||
конечное значение y(z) выражается как |
|
|
|
z — 1
Уоо = Iim — ; — У {г) = А = 0 ,5.
г-*1 6
Рассмотрим нелинейную импульсную систему. Метод применим для исследования импульсных систем со сла бонелинейными элементами. Он позволяет найти выход ной сигнал в дискретные моменты времени как функцию
Рис. 2-10. Нелинейная импульсная система.
любого входного сигнала. Недостатком метода является представление нелинейности в виде разложения в конеч ный ряд. Это приводит к линеаризации нелинейностей типа реле. Метод довольно удобен для характеристик типа насыщения. Локальная устойчивость системы сле дует из уравнения линеаризованной системы.
Рассмотрим пример, предложенный Джури [Л. 2-16]. Структурная схема представлена на рис. 2-10. Фиксатор нулевого порядка характеризуется тем, что выходной сигнал в промежутках между моментами пТ и (п +1)Г постоянен и равен г(пТ). Передаточная функция фикса тора нулевого порядка будет: •
1—е~тР
• >
Р
86
откуда
м |
1 |
1 |
1. о -'1'Р ' |
(2-103) |
Е(г) |
р + |
|
||
Разностное уравнение системы будет: |
|
|||
Sn+i + s„[a—ае~т—е^т]+ 6(1 —e-T)s3n= |
|
|||
|
= |
е |
„1 —( е - т ) - |
(2-104) |
Рассмотрим случай, когда входной сигнал является единичной ступенчатой функцией, причем линеаризиро ванная система описывается уравнением
(* + 1 - 2e~r ) S (г) |
|
z (1 — ё |
7 \ |
S{z)-. |
(2-105) |
(z — 1) (г + 1— 2е~г)
при а=1, решение которого получаем при помощи таб лицы обратных преобразований:
s(n7')=0,5[l — (—0,264)"] при Т= 1. (2-106)
Линеаризированная система устойчива, если при п— v
— >-оо s(tit)— >-0,5. Полагаем 6 достаточно малым по сравнению с единицей, таким, чтобы нелинейная система была устойчивой, например 6= 0,1. В этом случае по ана логии с методом свертки имеем:
|
ОО |
S <z) = |
1-Д г -1’ а ° = 1; Й2 —-ai ; ^ = "a‘i • (2-107) |
|
1 = 0 |
Подставляя это уравнение, а также выражение (2-89) в разностное уравнение (2-104), получаем:
|
|
|
|
00 |
|
1 |
|
|
|
|
V |
А( |
|
|
|
|
|
1 |
||
|
|
|
|
7 j |
1— atz |
‘ |
|
00 |
00 |
00 |
:=о |
|
|
х |
У 1 |
|
|
|||
У |
1 |
У ! |
|
|
||
|
' Zj Li Li |
1— axaxamz - ' |
||||
|
i=.Q / " О |
т = ;0 |
|
|
:^ |
1 |
'S |
|
Чi! |
X |
||
z-\-a — ae 1—е ‘ |
z(1 —ё
’_ |) (z + а_ ае- т— е~т) '
(2- 108)
87
предполагаем, что начальные условия пулевые.
Сравним коэффициенты в обеих частях при z(z— I ) - 1,
z(z—а)-1, ..., z(z—а*)-1. ... |
и положим z — 1, |
ait |
..., +. |
|||
В частном случае, когда 7=1, а= 1, получаем: |
|
|||||
|
, |
0,6326 |
лз — 0,5; |
|
|
(2-109) |
|
т- |
1,264 |
Ло |
|
|
|
, , |
|
0-6326 |
п.2 л — q. |
|
(2- 110) |
|
Л « + |
а , + 0,264 й л о Л 1 — |
и’ |
|
|||
|
|
|||||
А, |
|
(З-^о 4 + 3 4 |
Ло) = |
0. |
( 2- 111) |
|
a f + |
0,264 |
|
|
|
|
Уравнения позволяют рассчитать методом последо вательных приближений все неизвестные; например, если Ь = 0,1, получаем:
4 — 0.494; я, = -0 ,3 1 ; |
= |
(2-112) |
Решение s(tiT) устойчиво, так как при a i< l e~ait — исчезающая во времени функция, и аналогично этому то же происходит и для так как <+ = а*<1. Тогда для нулевых начальных условий из уравнения (2-107) получаем:
|
2 |
4 = 0. |
|
|
(2-113) |
|
i= 0 |
|
|
|
|
Подставляя соотношения (2-112) в равенство |
(2-113) |
||||
и ограничиваясь первым |
членом |
разложения, |
можно |
||
определить (в основном численно) |
Аи откуда |
|
|
||
0,494 |
|
0,461 |
0,049 |
.... |
(2-114) |
S(z) = 1—2 -' |
1+ 0,312~1 |
|
|||
0,0962- 1 |
|
||||
т. е. |
|
|
|
|
|
s (п7) =0,494—0,461 (—0,31) " — |
|
|
|||
—0,049 (0,096) ” + 0,007 (—0,03)"+ ... |
(2-115) |
88
Сохраняя шесть первых членов s(tiT), Джури пока зал на этом примере, что полученный результат близок к точному, полученному вычислением рекуррентного со отношения с точностью до третьего знака.
2-8. ПРИЛОЖЕНИЯ МЕТОДА z-ПРЕОБРАЗ ОБАИИЯ
а) РЕШЕНИЕ РАЗНОСТНЫХ УРАВНЕНИЙ, КОЭФФИЦИЕНТЫ КОТОРЫХ — ПЕРИОДИЧЕСКИЕ ФУНКЦИИ
Рассмотрим уравнение
■*(0 = 2 U(0 ■* (« + 0 + |
2 8i (0 У(я + 0. (2-116) |
t=0 |
1=0 |
где я — целое число; fi и gi — периодические функции, характеризующие изменения соответствующих перемен ных с периодом, равным единице:
fi(t)=U(t + 1);
gi(t)=gi(t+l ) .
Джури и Муллин [Л. 2-17] предложили следующий метод.
Принимаем:
t= n+ p + nr,
р |
ч |
|
х (я + р + т) = Yi fi {т) х (я + г) + £ gi (яг) у (я + г), |
||
i= о |
г=о |
|
тогда |
|
|
/Дя + р + яг) =/г(/я) |
и g ;(я + р + яг) = gy(яг). |
|
Но по определению |
|
|
X(z, т) = Z m[x(n— 1+яг)]. |
|
|
Возьмем ^-преобразование от обоих членов |
и по |
|
лучим: |
|
|
Zm[х (я + р -}- яг)] = £ U(яг) Z [я (я + г] - f |
|
|
|
г=о |
|
+ 2 f t ( m ) Z [ y ( / i + 1)]. |
(2-117) |
|
(=0 |
|
|
89
Согласно уравнению (2-6)
Z [ x { n - \ - i)\ = ziX (г) — X I x (/)
/=о
p
Z \ x ( n p +m)] — zP+1X(z, tn) — '£ix(i, m)zr>+l-i.
/=о
Подставим эти два уравнения в уравнение (2-117). Получаем:
р
zp+'X (г, т) = |
h i,n) ziX (z) — |
i= 1 |
i=0 |
|
|
h (m)IjX(i) zi-j |
Ё |
ft |
(« ) zlY (г) |
|
|
|
1 = 0 |
/•=,0 |
;=o |
|
|
|
|
ч |
1=1 |
p |
|
|
|
■ |
s gi (m) Ц у (/) z*~ i |
+ S |
* |
(/> m ) z p +1~ >, |
( 2 - 1 1 8 ) |
|
i= 0 |
/=o |
j=o |
|
|
||
где X (z, m) |
и X (г) — неизвестны. |
не имеет разрывов в моменты |
||||
Предполагая, что функция x(t) |
||||||
квантования, |
возможно принять т= 1 и X(z, |
1) = Х ( г ) : |
|
|||
|
Х { г ) zP+‘ - E |
f t ( |
|
|
|
|
|
|
г=о |
|
|
|
|
|
<7 |
Р |
|
*’= 1 |
|
|
= ^ ( г ) Ц ^ ( ! ) г ‘ — S f t ( l ) z‘ l ] x ( / ) z - J — |
|
|||||
|
i=0 |
1=0 |
|
/=о |
|
|
|
q |
£ = 1 |
|
р |
|
|
- S f t ( l ) |
z‘ S £ /( /) z - i |
+ 2P + ‘ I]x (/)2 - M |
(2-119) |
|||
|
1=0 |
/=о |
|
/=0 |
|
Это уравнение позволяет получить Аг(г), и , исполь зуя выражение (2-118), можно получить, наконец, выра жение для X (z, т). Отыскание оригинала х(п, т) осу ществляется просто с помощью обратного модифициро ванного 2-преобразования.
б) РЕШЕНИЕ ЛИНЕЙНЫХ РАЗНОСТНЫХ УРАВНЕНИИ С ПЕРИОДИЧЕСКИМИ КОЭФФИЦИЕНТАМИ
Рассмотрим уравнение
х ( п + 1) — а(п)х{п) +у{п), |
( 2- 120) |
90