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

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

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

т. е.

z2x (0) — 2zx (0) + zx ( 1) — Tzy (0) z 22z -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 ^ -)-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), (za2) ..., если пола­

гать 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^

= 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

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