книги из ГПНТБ / Видаль П. Нелинейные импульсные системы
.pdfпозволяет вычислить в моменты квантования величины, необходимые для решения задачи.
В этой главе мы вспомним методы графов и переменных состоя
ния и рассмотрим их обобщение на случай нелинейных импульсных систем.
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) -Е* tр); |
(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 являются |
преобразованиями Лапласа |
переменных |
||
в узлах 1—3. |
|
получаем: |
|
|
|
|
Исходя из уравнений (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) |
|
|
X» |
|
|
||
|
р + т - 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