книги из ГПНТБ / Цейтлин Я.М. Проектирование оптимальных линейных систем
.pdfВторое слагаемое этой функции •
|
|
|
Г_00 = r ( s ) - r + ( s ) |
||
имеет все полюсы в правой полуплоскости. |
|||||
Теперь |
определение |
функции Г+ (s) очевидно. |
|||
Если |
Г (s) = |
|
|
— рациональная функция, |
|
|
Фхг |
|
XX |
|
|
то Г+ (s) = |
(S) |
в |
соответствии |
с (II 1.68) представляет со- |
|
®7х |
|
||||
|
(S) |
+ |
|
полюсы которых находятся |
|
бой сумму тех элементарных дробей, |
в левой полуплоскости. Функция Г_ (s) представляет собой сумму остальных элементарных дробей, полюсы которых находятся в пра вой полуплоскости.
Если же Г (s) = Ф * г ^ не является рациональной функцией,
Фхх (s)
то прежде всего следует найти функцию у (t) как обратное пре образование Фурье функции Г (s), т. е.
Теперь функцию Г + |
(&) можно найти |
как |
прямое |
преобразова |
|||||||
ние Лапласа функции |
у |
(t), |
т. е. |
|
|
|
|
|
|
||
|
|
|
со |
|
|
|
|
|
со |
|
|
Г+ (s) = L |
[у (01 = |
J Y (0 е - 5 ' |
^ |
= |
j |
V* it) e-s ' |
Л , |
||||
|
|
|
0 |
|
|
|
|
—со |
|
|
|
где |
|
|
|
|
|
|
|
|
|
|
|
|
Y* (0 |
= |
о, t |
<0 |
|
|
|
||||
Точнее, |
у |
(0 |
t |
> о. |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
Y* (0 |
= |
0, |
t |
< а |
|
|
|
|||
|
|
|
|
|
V (0, |
* > а , |
|
|
|||
где а — абсцисса |
абсолютной |
сходимости. |
|
|
|||||||
Функция Г+ (s), являющаяся прямым преобразованием Фурье |
|||||||||||
функции времени v* (t), |
равной нулю при t |
<30, имеет все полюсы |
|||||||||
в левой полуплоскости. |
|
|
|
|
|
|
|
|
|
||
Таким образом, искомая передаточная функция оптимальной |
|||||||||||
системы F (s) всегда может быть определена |
по формуле (II 1.65). |
Функция веса оптимальной системы найдется как обратное пре образование Лапласа передаточной функции (III.65), т. е.
2m'J r j F(s)estds l(0 = L-> [F(s)\.
—їсо
23.СИНТЕЗ ЛИНЕЙНЫХ ОПТИМАЛЬНЫХ СИСТЕМ
СОБРАТНОЙ СВЯЗЬЮ
Вподавляющем большинстве случаев системы автоматического управления представляют собой замкнутые динамические системы (системы с обратной связью).
Такие системы |
могут |
быть |
представлены |
в виде |
структурной |
||
схемы, изображенной на |
рис. |
17. Здесь т (t) |
— полезный |
сигнал |
|||
на входе; п |
(t) — аддитивная |
помеха на входе; у (t) |
— фактиче |
||||
ский сигнал |
на |
выходе; |
z (t) |
— желаемый |
Сигнал |
на |
выходе; |
d (t) — помеха на выходе |
объекта |
(«возмущение нагрузки»); |
||
т |
W |
ъ |
, |
Q |
|
|
|
I
I
|
Рис. |
17 |
a (t), |
b (t) — сигналы на входе |
регулятора (корректирующего |
звена) |
и объекта (фиксированной |
части системы); є (t) — ошибка |
системы; W (р), G (р), Н (р) — операторы регулятора (коррек тирующего звена), объекта (фиксированной части системы) и желае мого (идеального) преобразования.
Из рис. 17 непосредственно следуют уравнения: |
|
|||||||||
|
WG (т + |
|
п — у) |
+ d = |
у, |
|
||||
откуда |
|
|
|
|
|
|
|
|
|
|
|
|
|
/ |
I |
\ |
I |
1 |
|
А |
|
|
|
|
|
|
|
|
|
|
|
|
|
y = TTwG^m |
|
+ n^ + |
|
TTwGd' |
|
||||
и |
|
|
|
|
|
|
|
|
|
|
|
а = т-\-п |
— у = Y^fWG (m |
+ n |
~" d)> |
|
|||||
|
b=Wa= |
|
W |
(т + п — d). |
|
|||||
Ошибка |
l |
+ W G |
|
|||||||
системы |
|
|
|
|
|
|
|
|
|
|
_ |
—(я |
W |
G |
\ |
|
|
W G |
n |
1 |
Иa |
e — z — y—yn |
j + |
WGjm |
|
l + |
W |
G |
i + |
wG - |
9* |
131 |
В этих уравнениях т, п, у, z, а, Ь, є представляют собой пре образования Лапласа функций т (t), п (t), у (t), г (t), a (t), Ь (t),
є (0 и
W = W(s) = W (p) U ;
|
G |
= G(s) |
= |
G (p) |
| p = s ; |
|
|
|
Н |
= H(s) |
= |
H (p) | p = s . |
|
||
Обозначим |
передаточные |
функции |
|
|
|||
F = |
1 |
F |
- |
W |
• |
WG |
(111.70) |
|
|
|
|
|
|
||
a |
l + W G ' - b |
|
l + |
WG' |
WG ' |
|
Между функциями (III.70) существуют простые соотношения, позволяющие выразить каждую из них через любую другую:
Fa= 1 - FbG = 1 |
|
F„ |
||
Fh = I - |
F |
a |
G |
(III.71) |
|
G |
~ |
|
|
- |
Fa |
= |
FbG. |
|
Эти уравнения можно переписать в виде
у = Fy |
(т + п) + Fad; |
||
а = |
Fa(m + |
п — |
d)\ |
b = |
Fb (т + |
п — |
d); |
є = (Я — Fy) т — Fyn — Fad.
Очевидно, что рассматриваемая система автоматического управ
ления |
будет |
устойчива |
только |
в том случае, |
если |
все три переда |
|||
точные |
функции { Fa, |
Fb, |
Fy \ |
не имеют полюсов |
в правой полу |
||||
плоскости. |
|
|
|
|
|
|
у (t), |
a (t) |
|
Действительно, только |
при |
этом условии |
функции |
||||||
и b (t) |
по |
истечении |
времени |
переходного |
процесса |
будут |
на |
ходиться в установившемся режиме, поскольку при отсутствии полюсов передаточных функций в правой полуплоскости пере ходные процессы, обусловленные любыми ненулевыми началь ными условиями, «затухают», т. е. соответствующие возмущения с течением времени t становятся меньше (по абсолютному значе нию) любого сколь угодно малого наперед заданного числа В > 0.
Задача синтеза — найти такую передаточную функцию регу лятора (корректирующего звена) W = W (s), при которой система удовлетворяла бы требуемому критерию качества, например ми нимуму дисперсии ошибки є (і) в установившемся режиме. Так как
Fb |
_ |
F y |
(III.72) |
I-Ftf |
~ |
G(l-Fy) |
' |
то сформулированная задача сводится к определению передаточ ных функций { Fa, Fb, Fy \, удовлетворяющих требуемому крите рию качества, после чего передаточная функция регулятора (кор ректирующего звена) находится по одной из формул (III.72). Поскольку передаточные функции связаны между собой соотно шениями (III.71), может показаться, что достаточно определить одну (любую) из них, а затем найти остальные две на основании этих соотношений.
Однако поступать таким образом можно только в одном част ном случае, когда объект (фиксированная, заданная часть си стемы) является устойчивым и минимально-фазовым, т. е. когда передаточная функция G (s) не имеет ни полюсов, ни нулей в пра вой полуплоскости.
Пусть
G(s) ~ B«(s)
где Ag (s), Bg (s) — полиномы. Рассмотрим все возможные случаи.
1. В результате решения задачи синтеза определена переда точная функция Fa, соответствующая устойчивой системе. Не обходимость устойчивости очевидна. Тогда
F - |
1 ~ F a |
- |
V~F')B8. |
ь |
G |
~ |
ЧAs |
|
|
|
|
|
F y ^ \ - F a . |
||
Передаточная функция Fy не |
будет иметь полюсов в правой |
||
полуплоскости. Что же |
касается |
передаточной функции Fb, то |
она не будет иметь полюсов в правой полуплоскости только в том
случае, если полином Ag не имеет нулей |
в правой |
полуплоскости, |
т. е. если объект является минимально |
-фазовым. |
Поэтому пере |
даточная функция Fa может быть принята за основу только в том случае, если объект является минимально-фазовым.
2. Определена передаточная функция Fb, соответствующая устойчивой системе. Тогда
Fa=*\-FbG=\—^\
Fy = F b G = - ^ .
Передаточные функции Fa и Fy не будут иметь полюсов в пра вой полуплоскости только в том случае, если полином Bg не имеет нулей в правой полуплоскости, т. е. если объект устойчив. По этому передаточная функция Fb может быть принята за основу только'.в том случае, если объект устойчив.
3. Определена передаточная функция Fy, соответствующая устойчивой системе. Тогда
Fa=\ |
— Fy, |
||
F |
— -fjL |
— |
Р У В Е |
Передаточная функция |
Fa |
не |
будет иметь полюсов в правой |
полуплоскости. Что же касается передаточной функции Fb, то она не будет иметь полюсов в правой полуплоскости только в том
случае, |
если полином |
Ag не |
имеет полюсов в правой полупло |
|||
скости, |
т. е. если объект является минимально-фазовым. |
Только |
||||
в |
этом |
случае передаточная |
функция |
Fy может быть |
принята |
|
за |
основу. |
|
|
|
|
|
|
Таким образом, если объект является устойчивым и мини |
|||||
мально-фазовым (оба |
полинома Ag и Bg |
не имеют нулей в |
правой |
полуплоскости), за основу может быть принята любая из пере
даточных функций | Fa, Fb, Fy\. |
|
Если объект является устойчивым, но не минимально-фазовым |
|
(полином Ag |
имеет нули в правой полуплоскости), за основу может |
быть принята только передаточная функция Fb. |
|
Если объект является неустойчивым, но минимально-фазовым |
|
(полином Ag |
не имеет нулей в правой полуплоскости, а полином Bg |
может их иметь), за основу могут быть приняты только передаточ
ные |
функции Fa |
или |
Fy |
[6, 9 ] . |
Если объект |
является неустойчивым и не минимально-фазовым |
|||
(оба |
полинома |
Ag и |
Bg |
имеют к<уши в правой полуплоскости), |
за основу не может быть принята ни одна из передаточных функ
ций { Fa, Fb, |
Fy\. |
|
оптимальных си |
|
Получим конкретные алгоритмы определения |
||||
стем. Будем считать, что сигналы т (t), |
п (t) и d |
(t) являются |
ста |
|
ционарными |
случайными функциями, |
причем положим для |
про |
стоты, что эти функции взаимно независимы, т. е.
Ф«п |
= ®ппг = |
®md |
= ®dm = ®nd = |
®dn = |
0; |
|
|
Фтт = Фт, |
Ф„п = Ф«, |
Фdd |
= Фd. |
|
|
Рассмотрим все возможные случаи синтеза системы. |
||||||
1. Объект устойчивый и минимально-фазовый. |
|
|||||
Определим |
передаточную |
функцию |
Fy: |
|
|
|
|
є = (Я — Fy)m |
— Fyn — (1 — Fy)d; |
|
|||
Фее = Фт(Н- |
Fy)(Я- |
Fy) |
+ 0„FyFy |
+ <J)d(l-Fy) |
(1 - Fy). |
|
Уравнение |
Винера—Хопфа |
|
|
|
^ і - = ~ Ф т ( Я - ^ ) + Ф ^ - Ф л і - ^ ) = і,
at у
где I — неизвестная функция, аналитическая в левой полупло скости, т. е. все полюсы функции I расположены в правой полу плоскости, или
откуда
|
|
|
|
|
F, |
|
1 |
г в |
|
|
|
|
(111.73) |
|
|
|
|
|
У |
АЛ |
A- |
J + |
|
|
|
||
|
|
|
|
|
|
|
|
|
|
||||
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
А = Фт |
+ Ф„ + Ф/, |
В = ФтН |
+ |
Фй. |
|
||||||
По |
формулам*(Ш.71) |
|
|
|
|
|
|
|
|
||||
|
|
|
Fa=l-Fy=l |
_ |
_ |
1 |
|
|
|
|
|
||
|
|
|
У ~ |
1 |
А + |
|
|
+ |
|
(111.74) |
|||
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
В |
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
F B |
~ |
G - |
A G A + |
A- |
|
|
|
|
|
Определим |
передаточную функцию Fb. |
|
|
|
|||||||||
|
|
|
є = |
( Я — FbG)m |
— FbGn |
— {\ |
—FbG)d- |
|
|||||
|
Фее |
— Фт(Н — FbG) (Я - а д + <UnFbGFbG |
|
||||||||||
|
|
|
|
|
+ Ф Л і - а д ( і - З Д ; |
|
|
||||||
дФее |
|
— |
|
|
|
|
~ |
|
- |
|
|||
dFb |
|
-Фто ( я - ад - f ®nFbGG - Фао (і - ад - % |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||
или |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
FbAGG |
— B G = |, |
|
|
|
|
|||
откуда |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
BG |
|
|
1 |
|
В 1 |
(IH.75) |
|
По |
формулам |
(III.71) |
Л" (GG)- |
_ + |
|
A+G |
. Л" J+ |
|
|||||
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
Г |
В |
' |
|
|
|
|
|
|
|
|
|
^ G = |
4 L |
А- |
+; |
|
|
|
|
|
|
|
|
Fa |
= |
|
|
|
1 |
- в |
- |
|
|
|
|
|
|
|
|
|
|
А' |
- |
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
||
Так |
как |
для |
устойчивого и минимального фазового объекта |
||||||||||
|
|
|
|
|
(GG)+ = G и (GG)~ = G, |
|
|
|
|||||
выражение |
для |
|
принимает |
вид |
|
|
|
|
|
||||
|
|
|
|
|
1 |
Г В |
1 |
_ |
fig |
Г В |
] |
|
|
|
|
|
F„ = A+G |
L А- \ + |
|
|
|
|
|||||
|
|
|
|
AGA- L А + |
J |
|
Определим передаточную |
функцию |
Fa: |
|
e = |
(H—l+Fa)m-(l-Fa)n |
— Fad; |
|
Фее = Ф т ( Я - |
1 + Fa) ( Я - |
1 + Fa) + |
Ф « (1 ~ Fa) ( ~ Fa) + |
+Ф * № ,
|
д |
Ф ™ |
- |
Ф,„ ( Я — 1 4 - Fa) — |
Ф„ (1 — Fa) + |
<bdFa = g, |
||||||
или |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
F a 4 |
- |
[ Ф т |
(1 - |
Я ) + |
Ф„1 |
= g, |
||
откуда |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
F |
= |
— |
Ф |
т ( 1 - Я ) |
+ Ф„ |
|
||
|
|
|
|
|
|
|
|
J + |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
Числитель выражения в квадратных скобках можно пред |
||||||||||||
ставить |
в |
виде |
|
|
|
|
|
|
|
|
|
|
|
|
|
Ф т ( 1 - Я ) + Ф„ = Ф т + Ф „ - Ф т Я = |
|||||||||
|
|
|
|
= A — Od |
— B + Od |
= A — B. |
||||||
Тогда |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
Г А |
— в |
1 |
|
і |
/1+ |
В |
1 |
|
|
|
|
|
L л- J + - л+ /1 |
л- J + |
||||||
|
|
- |
1 |
л+ |
|
|
|
|
= 1 |
|
1 |
|
|
|
— ^ г л |
|
|
А- |
+ |
А+ |
А- |
||||
По |
формулам (II 1.71) |
|
|
|
|
|
|
|||||
|
|
|
|
F„=\ |
|
-F |
|
—— |
В |
|
|
|
|
|
|
|
|
|
А- |
+ |
|
||||
|
|
|
|
|
|
|
г а |
— |
А + |
|
||
|
|
|
|
Fh |
= |
1 |
|
|
|
В |
} + |
|
|
|
|
|
|
|
|
|
|
|
. А- |
|
Таким образом, .в случае устойчивого и минимально-фазового объекта можно синтезировать оптимальную систему, беря за ос нову любую из трех передаточных функций. Иными словами, устой чивый и минимально-фазовый объект не накладывает никаких ограничений на синтез системы.
2. Объект устойчивый, но не минимально-фазовый.
В этом случае можно определить только передаточную функ
цию F b |
. Общее выражение для передаточной |
функции Fb опти |
|||
мальной |
системы дано в |
формуле (II 1.75). |
|
||
Так как объект |
устойчивый и не |
минимально-фазовый, |
|||
|
(BgBgy |
^Bg,^ |
(BeBt)j= |
Ве, (AgAgf |
= |
= Ag, (AgAg) = Ag. |
(111.76) |
Поэтому (III.75) с учетом (III.76) принимает вид
Fh |
= Js- |
АаВ |
|
AgA- J+ • |
|||
|
AgA+ |
||
|
|
По формулам (II 1.71)
••FhG=-=± |
АеВ |
|
АаА- |
||
АеА+ |
3. Объект минимально-фазовый, но неустойчивый
с,
I
т + |
I |
w ъ.
Рис. 18
(III.77)
(III.78)
В |
этом случае следует определять передаточную функцию Fy |
|
или |
Fа — 1 — Fy. Выражения для этих функций |
уже получены |
в формулах (III.73), (III.74). |
|
|
Таким образом, в зависимости от динамических |
характеристик |
объекта (нахождение нулей полиномов Ag и Bg в левой или правой полуплоскости) приходится применять различные методы синтеза, изложенные выше.
В случае же неустойчивого объекта рекомендуется сначала стабилизировать его путем введения дополнительной обратной связи, после чего решать задачу синтеза системы с «новым» (уже устойчивым) объектом1 [13].
Рассмотрим этот |
алгоритм. |
|
||
Структурная |
схема |
синтезируемой |
системы приведена на |
|
рис. 18, где К (s) |
= |
a (s |
передаточная |
функция дополнитель |
ной обратной связи, осуществляющей стабилизацию объекта.
1 В работе [9] предлагается весьма громоздкий метод синтеза, не требующий предварительной стабилизации и сводящийся к введению сложной системы изо метрических ограничений.
Передаточную функцию К (т. е. полиномы а и J3) следует вы брать таким образом, чтобы «новый» объект, обведенный на рис. 18 штриховой линией, оказался устойчивым.
Из схемы непосредственно следует
|
|
|
G [Ьх — |
К (у — |
d)] |
= |
у —• d, |
откуда |
|
|
(1 + |
KG) (y—d) |
= |
Gbu |
|
|
|
|
|||||
или |
|
|
|
|
|
|
|
|
|
|
У— d— l + K G - Ь г |
= G A , |
|||
где Gx |
= |
д _^ |
— передаточная |
функция «нового» объекта, |
|||
а также |
|
|
|
|
|
|
|
|
|
|
|
і = ~Q~ ІУ ~ |
d), |
||
откуда |
Fb |
= ~ |
Fy. |
|
|
|
|
Подставляя в выражение для Gx G = - ф - и К = — , получим tig а
Для устойчивости такого объекта требуется, чтобы знамена тель (II 1.79) удовлетворял условию
|
Q = p\4g + |
аВе = QH, |
|
|
(III.80) |
|
где |
<3H — полином, не имеющий |
корней в правой |
полуплоскости. |
|||
Задавшись видом полинома Q, не имеющего корней в правой |
||||||
полуплоскости (например, Q = |
1, Q = а + |
bs, |
где а > 0, 6 > 0 |
|||
и т. |
д.), всегда можно найти |
полиномы |
а |
и |
(3, |
удовлетворяю |
щие |
(111.80). |
|
|
|
|
|
Разумеется, эта задача не имеет однозначного |
решения, так |
|||||
как |
условию (III.80) при заданном Q удовлетворяет неограничен |
|||||
ное количество полиномов а и р . Следует |
иметь в виду, что и сам |
|||||
полином Q (не имеющий корней в правой |
полуплоскости) может |
|||||
быть |
задан бесчисленным количеством способов. |
|
Важно то, что поставленная задача всегда имеет решение, которое не является однозначным.
На практике задание полинома Q и вид полиномов а и Р сле дует выбирать, исходя из простоты реализации полученных за висимостей. В частности, полином а следует искать в классе по линомов, не имеющих корней в правой полуплоскости. При этом дополнительная обратная связь, осуществляющая стабилизацию неустойчивого объекта, сама будет устойчива.
После стабилизации объекта рассмотренным выше способом решается задача оптимизации системы с устойчивым объектом.
Так как «новый» объект может оказаться не минимально-фазо вым, поскольку полином Ag может иметь корни в правой полу плоскости, целесообразно определять передаточную функцию Fbv играющую роль функции Fb в предыдущей задаче.
В соответствии с формулой (II 1.75), заменив Fb на Fbl и G на G J ( получим
1 |
|
BGt |
(Ш.81) |
А+ (G&Y |
A- |
|
|
|
|
Теперь
ааАцАо
так как полином Q не имеет корней в правой полуплоскости, так же как и полином а.
Поэтому (III.81) примет вид
|
|
|
|
|
|
Q |
|
|
BaAgQ |
|
|
|
|
|
|
А+а (AgAg)+ |
QA-a |
(AgAg) |
|
||||
|
|
|
|
|
|
Q |
|
|
B~Aa |
(Ш.82) |
|
|
|
|
|
|
|
|
A- (AgAgY |
||||
|
|
|
|
А+а (AgAg)+ |
|
||||||
По формулам |
(III.71) |
|
|
|
|
|
|
|
|||
|
|
Fy |
= |
Fbfii |
|
= |
|
|
|
AaB |
|
|
|
|
|
|
|
(AgAgYA- |
(Ш.83) |
||||
|
|
|
|
|
|
(AgAg)+A+ |
|
||||
|
|
|
|
|
|
Fn |
1 - |
Fy, |
|
||
Кроме |
того, |
|
|
|
|
|
|
|
|
|
|
|
|
|
Fb^-Q |
= (AgAgYA+ |
|
AaB |
(ІП.84) |
||||
|
|
|
|
(AgAgYA- |
|||||||
В частности, если неустойчивый объект минимально-фазовый, |
|||||||||||
(AgAg)+ = |
Ag |
и |
(AgAg)' |
=Ag. |
Тогда |
(111.82) — (111.84) |
прини |
||||
мают вид: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ГF |
|
= |
—+ |
|
|
|
|
|
|
|
|
У |
|
А |
|
|
|
|
|
|
|
|
|
F |
|
—-ЁМ- |
|
А- |
(ІП.85) |
||
|
|
|
|
b |
|
~ |
AgA+ |
|
Fn=\ А+