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

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

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

Продолжение табл. 1-2

Структура линейного звена Вид передаточной функции Разностные уравнения

X [я 26 (6 + с) + Ь2с ( а + с) + 6( а - f с) + а с 2 (а+Ь )] 4-

 

Г 3 D2D3

 

3 D, D3 .

 

а с [ а 3 ( а Ь) а )

 

Ь3 ( а Ь) с)

 

 

 

3 - D , - D

2

I .

 

 

 

 

 

с 3 (Ь с) (с а )

J *

 

 

 

Н =

0 ,5Т* (А +

1) + F (D,D2 +

D2D3+

D3D, -

1) +

 

+ (a b c y, \а*Ь + с) + ЬЧ + а ) + а 2с ( а + b ) } +

 

 

, , Г D2D%4- 3D, +

3D3-(- 3

,

 

 

 

+ a b c [

a 3 ( a — b) ( с — a )

+

 

 

D,D3 4- 3D, +

3D3-(- 3

D,D24- 3D, + 3D2 4- 31 .

 

b3 ( a b) (b c)

 

c 3 (b c) ( c a )

J ’

/ =

0,5 T 2 (D,D2 + D2D3 +

D3D, d , — D,

D3) +

 

+ F (D,D2D3 — D, D2— Ds) +

c ( a b c ) ~ 2 \ a 2b (b + c) +

 

 

 

 

Г1 — D2 D3 3D2D3 ,

+ b * c { c + a ) + a c * ( a + b ) l + a b c [

a , { a - b ) ( c - a )

+

.

1 — D , D3 3D,D3

1— D ,— D2 — 3D,D21

 

 

b3 ( a — b) ( b — c )

 

c 3 (b c) ( c — a ) J ’

 

СЛ

Структура линейного звена

 

 

 

 

 

Продолжение табл. 1-2

 

 

 

Разностные уравнения

1 =

0,574 (DXD2 +

D2D3 +

а д + Д . + О . + Д , ) - F (D .D .+

+

D2D3 + D3DX) +

D (аЬс)~г \аЧ (b + с) + Ь2с (с + д) +

 

 

 

 

Г 3D2D3 — D2 — D3 ,

 

+ д С2(д + &)н-д6с [

а,^ -1 -6)(-с _ д; +

 

,

3DlD3 — D2 — D3

,

3DxD2 — Dx— D2

 

^

b3 (a b) (b c)

^

c3 (b — c) (c — a) ]

 

 

. , a*b{b+c) +

ЬЧ (д + с ) + a c 2 (д + 6 )1

= } 0 , 5 T > - F + ----------------------

 

 

-------------------------------- j E +

,, Г _______ D2D3________ .________ D1D3_________1

+

д3(д — b)(c — a)

^

b3 (a — b) (b — c) ^

 

DxD2

 

 

c3 (b — с) (c a)

 

D, = exp ^ — y ~ j ;

D2 = exp ^ —

D3 =

exp ( — yr— ]; д =

7 7 1;

6 = Г - 1; c = 7 7 1

П р и м е

ч а н и е . В таблице приняты

следующие обозначения: Е — идеальный импульсный элемент; £0— фиксатор нулевого порядка;

выходной,

е — входной сигнал; Т — период

регулирования; Dt — exp ( ~ Т / Т г); й г = ехр ( — 77 Т).

ii, и частности,

 

 

 

 

 

 

 

 

 

 

 

 

zn+l = zne

TlJ,,= z nDi

c Di= e r,r' ,

(1-140)

где znT+t

является

входным

сигналом

фильтра

(1 +

+ Г2р)-1 и,

следовательно,

описывается

уравнением

 

 

 

 

 

dy

__

 

 

 

 

 

 

 

 

 

 

 

2

d t ~z.

 

 

 

 

 

 

Учитывая уравнение, связывающее znT+t с znT, реше­

ние вышеприведенного уравнения запишем так:

 

 

УпТ+ t

Ч\

 

t/T,

t/Tа

 

 

-

t/T,

Тг-Т»

 

 

 

) ZnT + У п б

 

 

У п + i

Уп&

“ I-

^

 

 

T^

- e ~

TiT')z n.

 

Полагая D2 = е

_Y/т

получаем:

 

 

 

 

 

 

 

2,

 

 

 

 

 

 

 

Уп+i

Д у п

 

— ^

- ( Д - Д ^ .

(1-141)

Характеристическая матрица фильтра, связывающая

Уп+i с ул и z„, будет:

 

 

 

 

 

 

 

 

 

 

 

 

А =

° 2

T x -

T

t ( Д - Д )

 

 

 

(1-142)

 

 

О

 

D,

 

 

 

 

 

 

 

 

Определитель запишется как

 

 

 

 

 

 

 

Д =

 

 

Я -у- _ ‘ уг- (Д Д )

 

 

 

(1-143)

 

X,

 

Д — Я

 

 

 

 

 

 

 

 

 

 

 

 

 

 

откуда

 

 

 

 

 

 

 

 

 

 

 

 

 

Я" — Я (хх-f- Д ) -{- Д х, — х 2 j,

 

(Д — Д ) — О*

Используя правило,

приведенное выше, запишем:

Уп+2 аиуп+1

 

I

 

 

г, ..

| г, , _ .

|,д

 

 

 

vn+1

Д */п + 1+

Д

(а чУп + Д

)

 

П

■(Dl - D

i)(atlyn +

k2n) =

0.

 

(1-144)

 

Л -

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Следует определить kn и An+i, т. е. влияние входного сигнала а(еп). При наличии постоянного входного сиг-

нала а(е„)

(выходной сигнал

фиксатора

нулевого

по­

 

рядка) уравнение (1-139)

запишется так:

 

 

 

 

 

 

 

 

-t/т,

 

 

 

,-</Л

) а (£„г);

 

 

 

znT+t

z nT е

‘ + (1

n

A

 

((11

 

 

 

 

% n + i

 

Z

 

 

-

f -

-145)

 

Уравнение (1-141)

запишется как

 

 

 

 

 

 

УпТ+t

УпТе

Л

 

 

 

 

е

 

tlTl)zn

 

 

 

•А

 

 

 

 

 

а

( £

п )

;

+ [

 

л - л

 

 

 

 

 

»-Чт* .

Т,

-

f/r,

 

№ )

 

 

 

 

УЛ+1

:D2yr

Т гТг- Т .

 

(

А

(

 

 

A

А

 

 

 

г, — уг.

 

 

-

 

А

 

-

+

 

А

 

Т,

 

 

 

 

а Ы>

 

 

 

 

 

 

 

 

 

 

 

 

 

откуда

 

 

 

 

 

 

 

 

 

 

(1-146)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

*

£

=

(

1

А

fe'n = [

l - D,

тг-Тгт 2 (

 

А

a (s „ ). -

А(1-147)

Разностное

уравнение,

связывающее входные пере­

менные с выходным сигналом фильтра,

имеет вид:

 

 

[а

 

 

Т

 

_7 X(П,

X

 

( s n

7, -1 Г,D2

 

+

i

)

A

+D,[1- D.

Тг

Тг

(

А

А

-

т\

■(г

А

 

-« Ы

= 0 А-

(1-)148)[

Производя преобразования и учитывая, что уп = гп,

можно получить

уравнение

импульсной

системы

рис. 1-30.

 

 

 

 

 

 

 

1-7. МОДИФИЦИРОВАННОЕ ДИСКРЕТНОЕ ПРЕОБРАЗОВАНИЕ КАРСОНА

)

)

1

Метод исследования импульсных систем, основанный на ступенчатых функциях и разностных уравнениях, дает информацию о значениях функции f(i) только в момен-

56

ты замыкания импульсного элемента. Чтобы устранить этот недостаток, возможно использовать модифицирован­ ное дискретное преобразование Карсона [Л. 5], описы­ ваемое уравнением

+00

Cm\f(x)] = F*m(p) = p j J f ( x - A T ) e'Pxdx.

— 00

Это преобразование, предложенное С. Венгжином и А. Буковым, аналогично модифицированному z-преоб- разованию Р. Баркера и Э. И. Джури (см. гл. 2) и при­ водит к сложным расчетам. Принимая во внимание, что использование этого метода ограничивается исключи­ тельно линейными импульсными системами, мы не будем более подробно излагать его.

Глава вторая

МЕТОД z-П РЕОБРАЗОВАНИЯ

Метод z-преобразования был создан в 1730 г. Де Муавром [Л. 2-1]. Позже в 1812 г. Лаплас распространил этот метод на теорию вероятностей. Метод z-преобразо­ вания является одним из способов, позволяющим иссле­ довать линейные системы в частотной и временной об­ ластях. Это преобразование может быть использовано для отыскания решения линейного разностного уравне­ ния. Однако, так же как и дискретное преобразование Лапласа или Карсона, применение z-преобразования по­ чти исключительно ограничено линейными системами. Мы ограничимся общим описанием этого метода, цель которого дать необходимые сведения, которые могут быть применены для исследования дискретных нелиней­ ных систем.

2-1. РАЗЛИЧНЫЕ ОПРЕДЕЛЕНИЯ И СВОЙСТВА

В большинстве дискретных систем сигнал рассматри­ вается в дискретные моменты времени пТ (п— О, 1, ...), где Т— константа, называемая периодом квантования (регулирования). Исследование таких линейных дискрет­ ных систем может быть произведено с помощью z-пре­ образования [Л. 2-2—2-6] *.

1 См. также работу энциклопедического характера [Л. 2-18], в которой в других терминах дается полное изложение методов, со­ ответствующих дискретному преобразованию Лапласа {Прим. ред.).

57

а) ОПРЕДЕЛЕНИЯ

Если / ( / ) — непрерывная функция времени и Т— фиксированное положительное число (которое можно принять за единицу), то считают /(/) определенной для /^гО, причем z-преобразование /(/) есть функция ком­ плексной переменной z

оо

(2-1)

Z\f(n)\ = F*(z) = % f(nT)z

п —О

 

для |z | R= р—*, причем р — радиус сходимости

ряда.

Чтобы упростить обозначение, мы опустим звездочку,

являющуюся символом квантования:

 

2(f) =/(0 Г ) +f(T)z-i + f(2T)z-2+ ...

(2-2)

Если функция f(t) имеет разрыв при t= a, то f(nT) является правым пределом /(/). Часто используют обо­ значения f(tiT), f.(n), fn. z-преобразование совпадает, следовательно, с преобразованием Карсона квантован­ ной функции при замене етР на z, а свойства, на кото­ рые мы укажем, такие же, как и свойства дискретного преобразования Карсона.

Пример 2-1. г-преобразовапие функции f(t)=t запишется при помощи (2-1) как

00

Г ( г ) = 2 (0 - = Х пТг~п;

п—О

г (0 = Г г - 1 +

2Тг~2+ ... +

nTz~n + ... ==

=-- Гг-* (I

+

2г - ' + . . . ) -

(2-3)

-JZZryjr',

для I г I

<

1 z (()

 

Тг

( г - 1 ) 2

 

 

 

J

Если Т= 1, то следует сравнить это соотношение с выражением

ер

С (х) ~~ ( е р - 1 ) 2

б) СВОЙСТВА z-ПРЕОБРАЗОВАНИЯ

1. Линейность

При а и й-константах легко доказать, что

Z [ a f + b g ) = a Z (f ) + b Z { g ) .

(2-4)

58

2. Т е о р е м а с м е щ е н и я

Рассмотрим случай смещения в действительной обла­ сти. Если функция f(t) смещена на величину kT

(рис. 2-1), то

Z[f(t—kT)] = z~kF(z).

(2-5)

Если функция /(/) смещена на kT (рис. 2-2), то учет того, что f(t + kT)=0 при /< 0, приводит к тому, что вы­ ражение (2-5) запишется как

Z [ft * + kT)} = zhF (2) - zkf (0) -

V (Г) -

 

k—\

n

 

- . . . - z f \{k— 1)Г] = 2* F { z ) - ^ F ( i T ) z - i .

(2-6)

/~o

 

 

Рис. 2-1. Запаздывающая Рис.

2-2.

Опережающая

функция.

функция.

При смещении в комплексной области легко дока­

зать, что

 

(2-7)

Z[e-atf(i)] = F(zeaT).

Практически метод состоит в замене членов с г на

zeaT.

 

 

3.

Изменение масштаба

Определение z-преобразования

 

позволяет напи­

сать:

 

 

Z H ( 0 ] ^ ( - i r ) .

(2-8)

4. Теоремы о начальных и конечных значениях

Если ^-преобразование функции /(/) будет F(z) и

если предел F (z) при z— >оо существует, тогда

 

lim /(«7’) =

lim/7(2).

(2-9)

<->0

г-*оо

 

59

а) ОПРЕДЕЛЕНИЯ

Если f ( t ) — непрерывная функция времени и Т— фиксированное положительное число (которое можно принять за единицу), то считают /(/) определенной для /^ 0 , причем ^-преобразование f(t) есть функция ком­ плексной переменной г

00

 

Z\f(n)] = F*(z) = £ \{пТ)г~«

(2-1)

п —О

 

для |2 |>i? = p“ 1, причем р — радиус сходимости

ряда.

Чтобы упростить обозначение, мы опустим звездочку,

являющуюся символом квантования:

 

Z(f) —f(OT) + /(Г )2 _1+ /(27)2”2+ ...

(2-2)

Если функция f(t) имеет разрыв при t= a, то f(tiT) является правым пределом /(/). Часто используют обо­ значения f(nT), f,(n), fn. г-преобразование совпадает, следовательно, с преобразованием Карсона квантован­ ной функции при замене етР на 2, а свойства, на кото­ рые мы укажем, такие же, как и свойства дискретного преобразования Карсона.

Пример 2-1. г-преобразопапие функции / ( / ) = / запишется при помощи (2-1) как

ос

Т (г) =

г (/) — S

пТг~п;

 

 

 

 

/г—О

 

 

 

 

г (t) = Г г - 1+

2Tz~s + ... + пТг~п +

... =

Тг-'(1

+

2г - * +

. . . ) =

 

Гг

(2'3)

( г _

, )2

I

,

 

1 г (t) =

Tz

])2

 

для |г |<

■( г _

 

 

)

Если Г=1, то следует сравнить это соотношение с выражением

<?р

С(ер — I)2 ‘

б) СВОЙСТВА г-ПРЕ0БРА30ВАНИЯ

1. Линейность

При а и 6-константах легко доказать, что

Z(af+bg)=aZ(f)+bZ(g). (2-4)

58

2. Т е о р е м а с м е щ е н и я

Рассмотрим случай смещения в действительной обла­ сти. Если функция f(t) смещена на величину kT

(рис. 2-1), то

Z[f(t—kT)] = z~hF(z).

(2-5)

Если функция /(/) смещена на kT (рис. 2-2), то учет того, что f(t + kT)=0 при К О , приводит к тому, что вы­ ражение (2-5) запишется как

Z[f(t + kT)] = *F(z) - z * f (0) — zk~lf (T)

z/ \(k- 1)7’] = zfc F ( z ) - Z F ( i T ) z - i .

( 2-6)

 

i—o

 

Рис. 2-1. Запаздывающая

Рис. 2-2. Опережающая

функция.

функция.

При смещении в комплексной области легко дока­ зать, что

Z{e~atf(t)] = F(zeaT).

(2-7)

Практически метод состоит в замене членов с z на zeaT.

3. Изменение масштаба

Определение 2-преобразования позволяет напи­ сать:

(2-8)

4. Теоремы о начальных и конечных значениях

Если z-преобразование функции /(/) будет F(z) и если предел F (г) при 2— >-оо существует, тогда

limf («7’) =

lim/;’ (z).

(2-9)

t —>0

oo

 

59

Таким же образом, если функция (1—z~l)F(z) не имеет полюса, лежащего вне круга единичного радиуса на z-плоскости, имеем:

lira / (пТ) = lim i - = l F (г).

(2-10)

л-» оо

г-* 1

2

 

 

5.

Преобразование разности

Если

 

 

 

Д/(/г7’) =

/ [ ( я + 1 ) Л ~ / ( я П

(2-П)

и

 

 

(2-12)

vf{nT) =

f{n + T ) - f { n - \ ) T \ ,

имеем:

 

4-1

 

 

 

 

Z[A*f(n7')] = ( z - l)hF ( z ) - z ' £ ( z - 1)*->~1Д/(07’);

/=о

(2-13) £[V ftf('i7’)] = ( l - z - T ^ z ) . (2-14)

Нетрудно найти при помощи этих формул выраже­ ния:

Z [ ,A fm ] = ( z - l) F ( z ) - z F ( 0 );

(2-15)

Z [V /(«T )]= (1 —z-i)E (z).

(2-16)

Эти соотношения будут применены впоследствии для решения линейных уравнений с помощью z-преобразо- ваний.

6. Преобразование суммы

Z

k = o

= Т Г Г ^ (2 ).

(2-17)

J

В табл. 2-1 даны г-преобразования различных функций.

7. Теорема комплексной свертки

Пусть F*(p) и G*(p) — преобразования Карсона от f(t) и g(t). На основании теоремы о свертке в ком­ плексной р-области, сформулированной в [Л. 2-7, 2-8], получаем:

a+joo

с { [ / ( * ) £ ( * ) ] М О } =

 

 

С+/0О

1

dr. (2-18)

2nj j

W " J F^)G(r-q)dq

 

а—/о си

с—уоо

 

 

60

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