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

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

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

позволяет вычислить в моменты квантования величины, необходимые для решения задачи.

В этой главе мы вспомним методы графов и переменных состоя­

ния и рассмотрим их обобщение на случай нелинейных импульсных систем.

5-1. ГРАФ ЛИНЕЙНОЙ ИМПУЛЬСНОЙ СИСТЕМЫ

а) КРАТКОЕ ОПРЕДЕЛЕНИЕ ПОНЯТИЯ ГРАФА

Граф передачи сигнала изображается в классической форме гра­

фа цепи в виде ряда ветвей, соединяющих узлы. Следует при этом

различать два основных составных элемента: узлы (точки ввода или

суммирования сигналов), ветви (направленные в сторону передачи

данного сигнала и характеризующие его передаточную функцию).

Более

подробное и строгое описание метода графов можно найти

в [Л.

5-1 и 5-2].

Однако

напомним

важное прави­

ло, которое

иллюстрируется рис. 5-1:

 

п

 

 

= X

(5-1)

1

Граф передачи исследуемой системы получается из уравнений, характеризую­ щих линейную зависимость между раз­ ными переменными. После нескольких элементарных преобразований, основан­

ных на вычетах, можно определить передаточную функцию. Метод Мезона (Л. 5-3] позволяет определить передаточную функцию, при­ менив общее правило, которое гласит: передаточная функция между двумя узлами графа определяется соотношением

ZThAh

 

==_^

 

(5-2)

где

(— l^SG i-f- (— l ) 2SG iG j+ (— l)32GiG)Gk+

. . ( 5 - 3 )

Д= 1+

SGj — сумма

передаточных функций всех контуров;

S G iG j— сумма

двойных произведений передаточных

функций

всех

контуров;

SGjGj-Gft — сумма тройных произведений передаточных функций всех

контуров;

Гк — передаточная функция k-ro каскада, соединяющего

входной

узел

с рассматриваемым; Дк — величина Д,

взятая

для

части графа,

отделенного k-м каскадом. «Отделенный»

каскад

сле-

% h i

Х8

о -

х 1

л

Рис. 5-2.

141

дует понимать как не имеющий ни одной общей точки с рассматри­ ваемым каскадом.

Пример 5-1. Ищем проводимость Т= хе1хх графа рис. 5-2. Мы показали на рис. 5-3 различные контуры

 

2Gi = bg+ch+di+fghi

(5-4)

 

ZGiGj-bgdi.

(5-5)

 

a

h i

Z6,

 

 

 

Я

f

 

 

U i

d

 

 

 

 

Рис. 5-3.

 

Каскады приведены на рис. 5-4. Первая передаточная функция

Tt = abcde; Ai = l

(5-6)

не содержит разомкнутых контуров. Вторая передаточная функция

откуда получаем выражение для общей передаточной функции

х 6 _

 

abcde + afe (1 ch)

(5-8)

х j

1

(bg -f- ch -f- di + fghi)

-f- bgdi

 

б) ГРАФ ИМПУЛЬСНОЙ СИСТЕМЫ

Мы видели, что передаточная функция непрерывной линейной системы может быть получена с помощью метода графов. Следую­ щий этап состоит в распространении этого метода на анализ им­ пульсных систем. Однако очевидно, что применение метода Мезона не всегда так просто, как в рассмотренном выше случае.

142

Действительно, нельзя заменить импульсный элемент с помощью эквивалентной передаточной функции, и метод Мезона не может применяться ни к каким другим графам, как только к таким, кото­ рые содержат узлы, характеризующие либо непрерывные сигналы, либо дискретные сигналы. Однако большинство импульсных систем в зависимости от того, какая из их координат нас интересует, будет нести информацию как в непрерывной, так и в дискретной формах. Поэтому необходимо модифицировать классический метод графов. Это было проделано в различных публикациях Р. Ашем и др. [Л. 5-4], М. Зальцером [Л. 5-5] и В. С. Ку [Л. 5-6]. Рассмотрим одно­ контурную импульсную систему, представленную на рис. 5-5, ошибка которой е — квантованная величина. Тогда имеем:

5* (p) = G*(p)

G*(P) -Е*);

(5-9)

 

+ 773* (р)

 

5(2)

G{z)

(5-10)

■Е(г).

 

1+773 (z)

 

Построение графа производится последовательно.

1. Построение предварительного графа

Исходя из структурной схемы системы, строим первый граф, эквивалентный системе. Так, структурная схема рис. 5-5 представ­ ляется в виде графа рис. 5-6.

Рис. 5-5.

Обозначим все источники сигналов Е, а зависимые переменные и переменные на выходе импульсного элемента как х2—е*. Ветви, исходящие от узлов, обозначающих зависимые переменные, описы­

ваются передаточными функциями, вычисляемыми по правилу Мезо­ на (G, Н).

2. Построение импульсного графа

Этот этап требует более внимательного отношения. Рассмотрим рис. 5-6. Так как граф имеет импульсный элемент, то нельзя непо­ средственно применить правило Мезона. Однако можно записать:

X=E = E—HXS;

(5-11)

Х2= Е *= А *,;

(5-12)

Xs= S= GX2,

(5-13)

143

или для большей простоты обозначим,

например, Е*(р)

через Е*;

Xi,

Х2,Х3 являются

преобразованиями Лапласа

переменных

в узлах 13.

 

получаем:

 

 

 

Исходя из уравнений (5-11) — (5-13),

 

 

 

 

Х, = — HGX^ + E;

 

(5-14)

 

 

Х2= Х \ ;

 

 

(5-15)

 

 

X3=GX*, = GX2.

 

(5-16)

 

Рассмотрение этих

уравнений

показывает,

что

можно

исключить импульсный элемент из рис. 5-6 при условии, что к узлу

Х2 прикладывается

сигнал X2= X*i= E*. Если применить

дискретное

преобразование вышеприведенных уравнений, то получим:

 

X*, =

(HG)* А'*, + Е*, (HG)* = HG*;

(5-17)

 

Х*2= Х*и

(5-18)

 

X*3=G*X*i.

(5-19)

Эти три уравнения содержат только дискретные переменные и позволяют построить дискретный граф рис. 5-7. Использование пра­ вила Мезона приводит тогда к выражению

G* (в)

5* {р) = ----- 1 = ^ ~ - Е* (р), (5-20) 1 + HG*(p)

1

■£* (р).

Е* (р) :

\+ Ш ( р )

Рис. 5-7.

(5-21)

Если система содержит несколько контуров, метод может быть

применен без существенных изменений.

 

3.

Полный граф

Показано, что правило Мезона позволяет выразить квантован­ ные величины графа как функции входных сигналов.

Но граф содержит не только квантованные сигналы. В ряде слу­

чаев желательно получить выражение

для непрерывной переменной

 

как

 

функцию

входных

сигналов.

нй*

Нужно найти

граф,

отличающий­

ся

от предыдущего

тем,

что он

 

является

комбинацией

предвари­

 

тельного

и импульсного

графов.

 

Поэтому мы называем его полным

 

графом

(composite

signal flow

 

graph).

 

 

 

 

 

 

 

Импульсный элемент предвари­

 

тельного

графа

заменяется

искус­

 

ственным источником (входом).

 

Например, если импульсный эле­

 

мент располагается между

узлами

 

1 и 2

(рис. 5-6),

то операций кван­

 

тования

представляется с помощью

 

ветви

с

единичной

передаточной

144

функцией между узлом x*i импульсного графа и узлом х2 предвари­ тельного графа. Совокупность двух графов, соединенных таким обра­ зом, составляет полный граф и позволяет рассчитать неквантованпые переменные. Такой граф изображен на рис. 5-8 и имеет X2=X*t.

Величина непрерывной выходной величины графа получена с по­ мощью правила Мезона

G(p)

S(P)

i + т * (р) Е* (Р),

(5-22)

а также

 

 

 

Е (р) = £(/>)■

С (р) Н (р) Е* (р)

 

Е(Р)

1 + HG* (р)

1 +

(5-23)

 

HG* (р)

5-2. ИСПОЛЬЗОВАНИЕ

ПЕРЕМЕННЫХ

 

 

 

состояния

Импульсная система благодаря своему дискретному характеру

хорошо поддается исследованию с помощью

матричных уравнений.

В ряде исследований показано [Л. 5-7, 5-8], что использование ма­ триц весьма удобно при исследовании и расчете импульсных систем с переменными во времени элементами, а также с нелинейными эле­ ментами.

Метод переменных состояния является мощным средством иссле­ дования импульсных линейных систем. С одной стороны, матричное описание приводит систему к системе дифференциальных или разно­ стных уравнений первого порядка, решение которых легко найти. С другой стороны, начальные условия могут быть легко учтены, и исследование разных типов линейных импульсных систем можно рассматривать в общем стандартном виде.

Последующее рассуждение распространяется без труда на диф­ ференциальные уравнения, с помощью которых описывается непре­ рывная часть системы регулирования.

Рассмотрим переменные состояния линейной системы. Всякая дискретная система может быть описана системой п разностных уравнений первого порядка

пп

x t \(k + 1) Т] = S atiXj (kT) +

2 Ьи ег (kT),

(5-24)

/=1

е~\

 

 

 

 

где Х{ — переменные состояния, рассмотренные

в

моменты

kT

и

( й - И ) Г ; ei — различные величины входного сигнала

системы.

 

че­

Вообще выходной сигнал sm системы

может

быть выражен

рез значения входного сигнала е и переменных состояния х, т. е. получен из уравнения

[kT] = 2

cmix {kT) 4" £ d-meel {kT).

(5-25)

/ =

i

 

10— 352

145

Два предыдущих уравнения запишутся в матричной векторной форме:

Х [(£+1)7’]= АХ(й7’) + ВЕ(й7’);

(5-26)

S(kT)=CX(kT)+DE(kT). (5-27)

Эти матричные уравнения являются дискретными уравнениями относительно переменных состояния системы (в ряде случаев их на­ зывают уравнениями относительно передаточных функций состоя-

Рис. 5-9.

ния). Они описывают систему только в дискретные моменты времени. Решение этих уравнений с помощью прямого и обратного 2-пре­ образования и теоремы свертки не вызывает затруднений. Примене­ ние этого метода на практике облегчается тем, что за импульсным элементом следует фиксатор. В этом случае выходной сигнал эле­ мента, состоящего из импульсного элемента и фиксатора, является

ступенчатой функцией

yi(kT+)=tn(kT),

(5-28)

где

kT+ = kT+ x; 0 < т «£7\

З а м е ч а н и е . Этот метод весьма полезен при численных рас­ четах. Именно поэтому этот метод до его применения для нелиней­ ных систем был главным образом использован для решения различ­ ных задач, связанных с импульсными системами с дискретным ре­ гулятором (рядом дискретных корректоров).

Пример 5-2. Рассмотрим импульсную систему (рис. 5-9). Здесь Вс является фиксатором нулевого порядка [передаточная функция

Матричное уравнение, аналогичное (5-26) и соответствующее системе дифференциальных уравнений (передаточная функция G(p) связывает непрерывные переменные), запишется с помощью преоб­ разования Лапласа каждого из членов:

рХ(р) =АХ(р) +Ъ\ (р) +x(t+0),

(5-29)

откуда

 

 

 

 

X (р) = (pi -

А) - > х (/ + ) +

(pi - А) ■-* BY (р).

(5-30)

Но

х,

(t) = s ( 0 ;

 

 

 

 

 

 

dx i

 

 

 

dt

* г (0 ;

 

(5-31)

dx 2

=

xt (t)

у (0-

 

 

146

В обозначениях уравнения (5-26)

 

О

 

А

°

1

В:

(5-32)

 

0

—1

 

 

 

x ((q ), У ( t o ) — начальные условия х (<)

и j ( / )

в момент t — t^

Теперь можно записать:

1

 

 

 

 

1

 

Xi.Vo)

(р)

Р

Р { Р +

1)

 

+

 

 

 

 

(р)

О

 

 

* 2 (<0+ )

 

Р +

 

1

 

1

и

 

 

 

 

 

 

 

рр ( р + 1)

 

+

1

Y (р)-

(5-33)

 

 

1

 

 

0

р + 1

 

 

 

 

Учитывая

фиксатор

нуле­

 

хга 0)

вого порядка,

можно записать:

 

y(to) = e(to) — * i (*„)■

(5-34)

Граф передачи сигналов импульсной системы с учетом соотношений (5-32) и (5-33) приведен на рис. 5-10. Уравне­ ния для переменных состояния в этом случае могут быть по­ лучены с помощью правила Мезона для случая рис. 5-10.

Например,

Xi (р) = S(p)

~ l P

P2( P + 1)

Xi (to) +

 

P(P~h 1)

*2 (<c) +

1

e ( t 0).

(5-35)

 

 

 

З а м е ч а н и е . Легко записать уравнения для переменных со­ стояния непрерывной части G(p) в виде

 

Хг __

1 .

 

 

 

м2

 

 

(5-36)

 

 

 

 

р + т - 1

 

 

 

Y

 

 

 

Эти выражения можно представить как

 

 

* 2 (t) =

ж, (/,)

<*-*•> + 1 1 - в -

<*—'»>] у (t„);

 

t

 

 

 

 

1(5-37)

Хг (г) = J

х, (/) dt +

О

= х 2(г'о) [1 -

J +

^0 + [< - <0-1 + *-»-*•>] у (*„) + С»*

10*

147

и

X, (()

1

1

xi (to")

+

хг (t)

0

e - <*-'•>

*гР+)

 

 

i

 

\у(Ф |-

(5-38)

 

 

 

 

5-5. МЕТОД ГРАФОВ

а) ОБЩИН МЕТОД

Метод переменных состояния, применяемый совместно с мето­ дом графов, является мощным средством анализа и синтеза линей­ ных импульсных систем. Ту и Вадханафути [Л. 5-10] использовали их для исследования нелинейных импульсных систем, состоящих из

Рис. 5-11.

фиксатора нулевого порядка и линейной импульсной системы, с включенным в нее дискретным корректирующим устройством—зве­

ном с

переменным коэффициентом усиления. Предполагается, что дис­

 

 

r(tg) р

"1

3(„)

кретное

регулирующее

устрой­

 

 

ство является нелинейным уси­

e(t0)

Д о - *-—о 4

G(р)

1

о

лителем, у которого коэффици­

\ y(tg) Y(p)\_______ |J

 

ент усиления

K(kT) постоянен

 

 

 

 

 

 

в интервале

времени

kT<t<,

 

S ( t g )

S ( t 0 )

 

 

 

<(k + l)T.

 

нелинейную

 

 

 

 

Рассмотрим

 

 

Рис. 5-12.

 

 

систему,

изображенную

на

 

 

 

 

рис. 5-11. Здесь Во— фиксатор

e~Tp)p~l;

G(p) — линейная

часть

нулевого

порядка

Во (р) = (1—

системы регулирования. Нели­

нейный элемент НЭ имеет вход s(t)

и выход

y(t).

Он

может

быть

описан как усилитель с коэффициентом К(кТ)=Кк, постоянным

между двумя моментами квантования

k T < t^ (k + l )T ,

но изменяю­

щимся от одного периода квантования

к другому:

 

 

y(kT +)

 

=

е (kT+)

(5-39)

где К* зависит от нелинейной

характеристики элемента, а также и

от периода квантования.

 

 

 

С учетом замечаний, содержащихся в § 5-1 и 5-2, изобразить

граф нелинейной импульсной системы не представляет

труда. Такой

148

граф изображен на рис. 5-12, Уравнения для переменных состояния запишутся в матричной форме

X{(k+l)T)=A(T, Кк)Х(кТ) + й{Т,

Кк) E(kT).

(5-40)

Величины Кк

в этом уравнении заданы

соотношением

(5-39),

в котором г(кТ+)

определяется с помощью уравнения

 

 

e(kT+)—e(kT)s(kT).

 

(5-41)

Если z{kT) найдено, то по нелинейной характеристике элемента находим y(kT+). Теперь легко могут быть найдены значения Кк для каждого периода квантования, а уравнение (5-40) совместно с гра­ фом рис. 5-12 позволяет определить переменные состояния в момен­ ты квантования.

Пример 5-3. Рассмотрим импульсную систему рис. 5-11, причем передаточная функция непрерывной части будет:

G( P) :

Р( Р +

Нелинейность, рассматриваемая в примере, — это реле с зоной нечувствительности D

 

 

 

(

1,

если е (/) ^ D',

 

'1

 

 

 

«/(<) =

<

0,

если

D < е (0 <

D; I

(5-42)

 

 

 

1,

если s (<) < — D.

J

 

Зная период регулирования (например, Т=

I

сек)

и начальные

условия

 

 

 

х ,( 0) = « ( 0);

 

%

 

 

 

 

 

 

 

 

 

 

 

Г

dXi

Г ds 1

 

\

(5-43)

 

 

х*(°)=

1 Г

dt

 

 

 

определим

 

сигнал

выхода

s(t)

 

 

 

 

в моменты квантования при вход­

 

 

 

 

ном сигнале е (< )= 0.

 

 

на

 

 

 

 

Система.

изображенная

 

 

 

 

рис. 5-11,

идентична

системе рис.

 

 

 

 

5-9, за исключением нелинейного

 

 

 

 

звена, а ее граф строится с уче­

 

 

 

том графа рис. 5-10 и переменной

 

 

 

 

величины

коэффициента

усиления

 

 

 

K(to)=K(kT)=KK. Граф пе;редачи

 

 

 

сигналов

системы

изображен на

 

 

 

 

рис. 5-13. Уравнения для расчета

п 1лучаются

из

вышеприведенного

переходных

процессов в

системе

графа и уравнений

(5-34) и (5-37)

 

 

 

 

 

xtl(k+l)T] = l\ - K K(T - \+e - T)]Xl(kT) +

 

 

 

+ (1- е - т ) х 2(кТ) +Кк(Т-\ + е-т)е(кТ)-,

(5-44)

 

 

xA(k+l)T]=—Кк(1— e~T)xi(kT) -f

(5-45)

 

 

+ e - Tx2(kT)+KK(l—e~T)e(kT).

 

 

 

 

 

 

 

 

 

 

149

Получение решения алгебраических уравнений ,в общем виде затруднено, самый простой путь — применение численных методов [Л. 5-6]. Пусть задано, что

х, (0) = 0,138;

х2 (0) = 0,762;

 

 

 

7

=

1;

 

 

 

 

 

 

 

 

(5 - 4 6 )

 

 

 

0 = 0,1;

 

 

 

 

 

 

 

 

 

 

 

 

 

е (/)

=

0.

 

 

 

 

 

 

 

 

 

Для

6= 0

&(0+) = —х Д 0 )= —0,138,

и,

если

учитывать еоот-ноше-

ния (5-42), tf= (0 + )------

1;

к

 

П0+) =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7,25.

 

 

 

 

 

 

 

 

 

 

 

 

 

е(0+)

 

 

 

 

 

 

 

 

 

Уравнения

(5-44)

и (5-45)

принимают вид:

 

 

 

 

 

 

 

* 1 ( 7 ) = (1 — 0 ,3 6 8 A o) * i ( 0 ) + 0 , 6 3 2 * 2 (0 ) = 0 , 2 5 2 ;

 

(5 - 4 7 )

 

 

( 7 ) = — 0 , 6 3 2 K o * i ( 0 ) + 0 , 3 6 8 * 2 ( 0 ) = — 0 ,3 5 2

 

(5 - 4 8)

 

 

в ( 7 + ) = — x t ( 7 ) = — 0 ,2 5 2 ; у ( Т + ) = — 1,

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

Я ,=3,97.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение

задачи

осущест­

 

 

 

 

 

 

 

вляется шаг за шагам: если

 

 

 

 

 

 

 

*i(7 ), х2(7), Л4

 

известны,

то

 

 

 

 

 

 

 

это

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

значения

 

 

 

 

 

 

 

х\(67) =s(kT) , x2(kT)

в момен­

 

 

 

 

 

 

 

ты 67.

 

 

 

 

 

 

по­

 

 

 

 

 

 

 

 

 

Рассмотренный пример

 

 

 

 

 

 

 

казывает,

что ®

системе

при

 

 

 

 

 

 

 

заданных

начальных

условиях

 

 

 

 

 

 

 

устанавливаются

несимметрич­

 

 

 

 

 

 

 

ные автоколебания с 'периодом

 

 

 

 

 

 

 

7 '= 47

(рис.

5-14).

Естествен­

 

 

 

 

 

 

 

но,

что

в зависимости от

на­

 

 

 

 

 

 

 

чальных

 

условий

и

величины

 

 

 

 

 

 

 

входного

сигнала

автоколеба­

ния в системе могут устанавливаться и более медленно

(переходный

процесс

может

отсутствовать)

 

и период этих колебаний

может быть

разным. Так, например:

 

 

 

 

 

 

 

 

 

 

 

 

 

для * i(0) =

0,1,

*г(0)=0,4

получаем

7 '= 4 7 ;

 

 

 

 

 

»

* i (0 ) =

— 0 ,4,

х 2 (0 )

= 1,2

 

»

 

7 ' = 5

7 ;

 

 

 

 

 

»

jt i (0 ) =

— 0,1

х 2 ( 0 ) = 0 , 6

 

»

 

7 ' = 6

7 .

 

 

 

 

 

Величина D зоны нечувствительности является весьма сущест­

венным параметром. Установлено,

что если 7>= 0,6,

то

автоколебания

полностью исчезают.

 

метод

приводит к

 

решению

разностного

З а м е ч а н и е . Этот

 

уравнения, описывающего

поведение системы.

Периоду

автоколеба­

150

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