Vedyakov_articles / Операторный метод анализа и синтеза ЛСУ
.pdfY (s) = L[ y(t)]
к изображению входной переменной
U (s) = L[u(t)]
при нулевых начальных условиях
y(0) = py(0) =K = pn−1 y(0) = u(0) = pu(0) =K= pm−1 y(0) .
Напомним, что преобразование Лапласа
∞
F(s) = L[ f (t)] = ∫ f (t)e−st dt
0
ставит в соответствие каждой функции f (t) (оригиналу), для которой несобствен-
ный интеграл сходится, единственную функцию F(s) (изображение) комплексной переменной s . Использование преобразования Лапласа для изучения дифференциальных уравнений основывается на утверждении:
L[ pk f (t)] = sk F(s) ,
если равны нулю значения f (t) и ее производных вплоть до ( k −1 )–ой при t = 0 .
Применяя преобразование Лапласа к правой и левой частям уравнения (2.1) при
нулевых начальных условиях, получим |
|
||
a(s)Y (s) = b(s)U (s) , |
|
||
и, следовательно, передаточная функция - |
|
||
W (s) = |
b(s) |
. |
(2.3) |
|
|||
|
a(s) |
|
|
Передаточная функция у которой полиномы a(s) и b(s) не являются взаимно простыми называется вырожденной передаточной функцией. При выводе передаточной функции W (s) системы необходимо следить за тем, чтобы не произошло со-
кращения ее числителя и знаменателя. При выполнении последнего условия, знаменатель a(s) передаточной функции называется характеристическим полиномом
системы. Устойчивость объекта управления определяется корнями характеристического уравнения a(s) = 0 . Для асимптотической устойчивости необходимо и доста-
точно отрицательность вещественных частей всех корней. Любой полином, все кор-
11
ни которого имеют строго отрицательную вещественную часть, называется устой-
чивым полиномом.
Корни знаменателя передаточной функции называются полюсами, а корни числителя - нулями. Передаточная функция, у которой вещественные части всех нулей отрицательны называется минимально-фазовой. Для существования решения задачи управления достаточно, чтобы передаточная функция была минимально-фазовой. Если W (s) не является минимально-фазовой передаточной функцией, то решение
задачи управления может существовать только при условии, что наибольшее общее кратное полиномов a(s) и b(s) - устойчивый полином.
Если W (s) является минимально-фазовой передаточной функцией, то решение
задачи будем искать в два этапа. Для этого представим сигнал управления в виде
|
|
|
|
u =W0 ( p)(u1 + u2 ) , |
(2.4) |
|
где |
W0 |
( p) = |
b0 ( p) |
- дробно – рациональная функция. На первом этапе определим |
||
a0 ( p) |
||||||
|
|
|
|
|
||
алгоритм формирования переменной u1 для того, чтобы решить задачу стабилиза-
ции, т.е. обеспечить выполнение предельного соотношения lim y(t) = 0 при любых
t→∞
начальных условиях. На втором этапе, выбором сигнала управления обеспечим ре-
шение задачи слежения, т.е. выполнение неравенства | y∞(t) − y (t) |≤ γ, где y∞(t) -
функция времени такая, что
lim[ y∞(t) − y(t)] = 0 ,
t→∞
y (t) - заранее неизвестная переменная величина (задающее воздействие), а γ > 0 -
некоторое заданное постоянное число. Отметим, что определенная подобным обра-
зом функция y∞(t) называется установившейся реакцией системы.
В заключении первого этапа курсовой работы строятся ЛАХ и ЛФХ передаточной функции W (s) объекта управления. Напомним, что логарифмической амплитудно-
частотной характеристикой, соответствующей передаточной функции W (s) , назы-
вается график функции
L(ω) = 20lg | W ( jω) |
12
от логарифма lgω, где ω - действительная переменная, которая называется часто-
той. Функция W ( jω) , которую получают из передаточной функции W (s) при под-
становке в нее s = jω, называется частотной передаточной функцией. Принято измерять значение ЛАХ в децибелах (дБ), а значение наклона ЛАХ – в децибелах на декаду. Декадой называют интервал на котором частота ω меняется в 10 раз. Логарифмической фазочастотной характеристикой называется график функции
ϕ(ω) = argW ( jω)
от логарифма lgω.
2.2. Решение задачи стабилизации
Для решения задачи стабилизации выберем алгоритм формирования сигнала u1 в
виде
|
|
|
|
|
|
|
u1 = −W1( p) y , |
|
|
(2.5) |
|||
где W |
( p) = |
b1( p) |
. Подставляя (2.4), (2.5) в уравнение (2.2) и разрешая его относи- |
||||||||||
|
|
||||||||||||
1 |
|
a1( p) |
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|||||
тельно выходной переменной найдем уравнение замкнутой системы |
|
|
|||||||||||
|
|
|
|
|
|
|
y =Wy ( p)u2 , |
|
|
(2.6) |
|||
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Wy ( p) = |
|
|
|
W ( p)W0 ( p) |
|
= |
|
b( p)b0 ( p)a1( p) |
|
. |
(2.7) |
|
|
1 |
+W ( p)W0 ( p)W1 |
( p) |
a( p)a0 |
( p)a1( p) +b( p)b0 |
( p)b1( p) |
|||||||
|
|
|
|
|
|
||||||||
Для объектов управления с минимально-фазовыми передаточными функциями, по-
линомы a0 ( p), a1( p), b0 ( p), b1( p) |
выбирают из условия |
|
|
|
|
|||||
|
|
b( p)b0 ( p)a1( p) |
|
= |
|
|
1 |
, |
(2.8) |
|
|
a( p)a ( p)a ( p) +b( p)b ( p)b ( p) |
a |
y |
( p) |
||||||
0 |
1 |
0 |
1 |
|
|
|
|
|
||
где ay ( p) - произвольный устойчивый полином степени n − m . Например, полином ay ( p) можно выбрать следующим образом
ay ( p) = ( p +1)n−m , |
(2.9) |
13
а многочлены a0 ( p), a1( p), b0 ( p), b1( p) -
a0 ( p) = a0,n−1 pn−1 + a0,n−2 pn−2 +K+ a0,0 ,
a1( p) = a1,n−1 pn−1 + a1,n−2 pn−2 +K+ a1,0 ,
b0 ( p) = b0,n−1 pn−1 +b0,n−2 pn−2 +K+b0,0 ,
b1( p) = b1,n−1 pn−1 +b1,n−2 pn−2 +K+b1,0 ,
где 4(n −1) неизвестных коэффициентов этих полиномов находятся из тождеств
|
|
|
|
a ( p) = b ( p) = ( p +1)n−1 , |
|
|
|
|
|
|
(2.10) |
||||||||
|
|
|
1 |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a( p)a ( p) +b( p)b ( p) = b( p)( p +1)2n−m−1 . |
|
|
|
|
(2.11) |
|||||||||||
|
|
|
0 |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таким образом, принимая во внимание |
u2 |
|
|
|
|
|
|
|
|
|
|
|
y |
||||||
(2.7), (2.9)-(2.11), передаточная функ- |
|
|
|
|
|
W0 ( p) |
|
|
W ( p) |
|
|||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||
ция объекта |
управления |
охваченного |
u |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
отрицательной обратной связью по вы- |
|
|
|
|
|
|
|
|
|
|
W1 ( p) |
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
ходу будет иметь вид |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
Рис. 2.2. Структурная схема |
|
|
|||||||||||||||
|
|
1 |
|
|
|
|
|
|
|||||||||||
Wy (s) = |
|
|
, (2.12) |
|
|
|
|
|
|
системы управления |
|
|
|||||||
(s |
+1)n−m |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Структурная схема системы управления представлена на рис.2.2.
В заключении данного этапа строятся ЛАХ и ЛФХ передаточной функции Wy (s) .
2.3. Синтез следящей системы управления
С учетом того, что измерения проводятся с помехами, управляющее воздействие объектом (2.6) строится в форме:
|
|
|
|
|
|
u2 |
~ |
(2.13) |
|
|
|
|
|
|
|
=W2 ( p)ε, |
|||
|
W2 |
( p) = |
b2 |
( p) |
|
~ |
|
|
|
где |
|
|
и |
ε = y * (t) |
− y − f (t) = ε − f (t) . Структурная |
схема системы |
|||
a2 |
( p) |
||||||||
|
|
|
|
|
|
|
(2.2), (2.4), (2.5), (2.13) приведена на рис.2.3.
14
Задача этого этапа со- |
|
u2 |
|
W0 ( p) |
u |
W ( p) |
y |
y * |
|||||||||
стоит в выборе такой переда- |
|
|
|
||||||||||||||
точной |
|
функции |
регулятора |
|
|
|
u1 |
|
|
|
|
ε |
|||||
|
|
|
|
|
|
|
|
|
|
||||||||
W2 (s) , которая обеспечит бли- |
|
|
|
|
|
|
W1 ( p) |
|
|
||||||||
зость |
ошибки |
управления |
|
|
|
|
|
|
|
~ |
f |
||||||
ε(t) = y |
|
(t) − y(t) |
к |
нулю |
и |
|
|
|
|
|
|
W2 ( p) |
ε |
|
|||
|
|
|
|
|
|
|
|
|
|
||||||||
требуемые |
характеристики |
|
|
|
Рис. 2.3. Структурная схема сле- |
|
|||||||||||
замкнутой |
системы |
управле- |
|
|
|
дящей системы управления |
|
||||||||||
|
|
|
|
|
|
|
|
|
|||||||||
ния. Для решения этой задачи воспользуемся методом динамической компенсации. |
|||||||||||||||||
Используя уравнения (2.2), (2.4), (2.5) и (2.13) найдем выражение для ошибки |
|||||||||||||||||
управления |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
ε(t) = [1 −W |
( p)]y (t) +W |
( p) f (t) , |
|
|
|
(2.14) |
|||||||
|
|
|
|
|
|
|
|
з |
|
|
|
з |
|
|
|
|
|
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Wз( p) = |
|
Wp ( p) |
|
, |
Wp ( p) =Wy ( p)W2 ( p) . |
|
(2.15) |
|||||||
|
|
|
1 |
+Wp |
( p) |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Функция Wз(s) называется передаточной функцией замкнутой системы, а функция |
|||||||||||||||||
Wp (s) - передаточной функцией разомкнутой системы. |
|
|
|
||||||||||||||
Выбор передаточной функции регулятора определяет вид передаточной функции |
|||||||||||||||||
замкнутой системы. Простейший подход к выбору обратной связи заключается в |
|||||||||||||||||
том, чтобы предъявить требования к самой Wз(s) , к примеру, потребовать выполне- |
|||||||||||||||||
ние условия |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
W (s) |
|
|
|
|
|
|
|
|
|
|
W (s) =W (s) = |
p |
, |
(2.16) |
|||
|
|
|
|
|
|
|
1 +Wp (s) |
||||||
|
|
|
|
|
|
|
з |
|
з |
|
|
||
|
|
|
b |
(s) |
|
|
|
bp (s) |
|
|
|
|
|
где W |
з |
(s) = |
з |
|
и W |
p |
(s) = |
|
- соответственно, желаемые передаточные функ- |
||||
aз (s) |
ap (s) |
||||||||||||
|
|
|
|
|
|
|
|
||||||
ции замкнутой и разомкнутой системы, которые формируются на основе требований
ксистеме управления. Преобразуем требования, предъявленные в курсовой работе,
кзамкнутой системе в условия на передаточную функцию замкнутой системы.
15
Требование к устойчивости замкнутой системы (a) означает, что все корни харак-
теристического полинома |
|
∆(s) =[a(s)a0 (s)a1 (s) + b(s)b0 (s)b1 (s)]a2 (s) + a1 (s)b(s)b0 (s)b2 (s) |
(2.17) |
замкнутой системы (2.2), (2.4), (2.5), (2.13) должны иметь строго отрицательную вещественную часть.
Установившаяся ошибка, вызванная воздействием (1.1), равна |
|
||||||||
|
|
|
|
|
d |
|
|
|
|
ε |
∞ |
(t) = [1 −W |
(0)](a |
+ a t) + |
[1 −W (s)] |
a . |
(2.18) |
||
|
|||||||||
|
з |
0 |
1 |
ds |
з |
1 |
|
||
|
|
|
|
|
|
s=0 |
|
||
|
|
|
|
|
|
|
|||
Ошибка может быть ограниченной при a1 |
≠ 0 , только если выполнено условие ас- |
||||||||
татизма |
|
|
|
|
|
|
|
|
|
|
|
|
Wз (0) =1 |
|
|
(2.19) |
|||
или |
|
|
|
|
|
|
|
|
|
|
|
|
Wp−1(0) = 0 , |
|
|
(2.20) |
|||
но это возможно, только при условии, что передаточная функция разомкнутой системы имеет нулевой полюс. Если он простой, то справедливо представление
Wp (s) = |
kp |
~ |
~ |
|
|
|
|
Wp (s) , Wp (0) |
=1, |
(2.21) |
|||
s |
||||||
|
|
|
|
|
||
где kp - коэффициент усиления разомкнутого контура. Тогда используя выражения
(2.15) и (2.21), найдем |
|
|
|
|
|
||
|
d |
|
|
= k−1 |
|
||
|
[1 −W (s)] |
|
(2.22) |
||||
|
|
|
|||||
|
ds |
|
з |
s=0 |
p |
|
|
|
|
|
|
|
|||
и, следовательно, |
|
|
|
|
|
||
|
|
ε |
∞ |
(t) = k−1a . |
(2.23) |
||
|
|
|
|
p |
1 |
|
|
Таким образом, для выполнения требования (1.2) необходимо, чтобы коэффициент удовлетворял условию
k |
p |
≥ |
a1 |
= δ−1. |
(2.24) |
|
ε |
||||||
|
|
1 |
|
|||
|
|
1 |
|
|
||
Установившаяся ошибка, вызванная воздействием (1.3), является гармонической функцией с амплитудой
16
εy = |
|
1 −Wз ( jω) |
|
a y , |
(2.25) |
|
|
где ω- частота воздействия. Таким образом, для удовлетворения требования (1.4) необходимо, чтобы
|
1 −W |
( jω) |
|
≤ |
εy |
= δ |
y |
, ω ≤ ω . |
(2.26) |
|
|
||||||||
|
|
|
|||||||
|
з |
|
|
|
a y |
|
y |
|
|
|
|
|
|
|
|
|
|
|
Установившаяся ошибка вызванная воздействием (1.5), является гармонической функцией с амплитудой
|
ε f = |
|
Wз ( jω) |
|
a f . |
(2.27) |
||||
|
|
|||||||||
Таким образом, для удовлетворения требования (1.7) необходимо, чтобы |
|
|||||||||
|
Wз ( jω) |
|
≤ |
εf |
= δ f , ω ≥ ωf . |
(2.28) |
||||
|
|
|||||||||
|
|
|
||||||||
|
|
|
|
|
|
a f |
|
|||
Если передаточная функция объекта управления является минимально-фазовой и устойчивой, то для построения желаемой передаточной функции замкнутой системы, удовлетворяющей выше перечисленным условиям часто прибегают к методу основанному на соответствии между логарифмическими частотными характеристиками разомкнутой системы и ее статическими и динамическими свойствами в замкнутом состоянии.
Явные зависимости L(ω) , ϕ(ω) достаточно сложны. Поэтому часто ограничива-
ются построением асимптотических логарифмических частотных характеристик. Для примера рассмотрим произвольную минимально-фазовую передаточную функцию с вещественными нулями и полюсами
|
|
′ |
|
|
|
−ν |
∏(Tis +1) |
|
|
W (s) = ks |
i |
, |
(2.29) |
|
|
′′ |
|||
|
|
∏(Tis +1) |
|
|
i
где k - положительный коэффициент усиления, Ti′, Ti′′ - положительные постоян-
ные времени,. ν - натуральное число (порядок астатизма). Тогда соответствующие ЛАХ и ЛФХ даются формулами
′ |
2 |
′′ |
2 |
+1] |
, |
(2.30) |
L(ω) = 20lg k −20νlgω+ ∑20lg[(Tiω) |
|
+1] −∑20lg[(Ti ω) |
|
|||
i |
|
i |
|
|
|
|
17
ϕ(ω) = ν |
π |
′ |
′′ |
(2.31) |
2 |
+ ∑arctg(Tiω) −∑arctg(Ti ω) . |
|||
|
i |
i |
|
|
Асимптотическая ЛАХ есть кусочно-линейная функция, получаемая заменой в
(29) членов 20 lg[(Tω) 2 +1] на
0, ω ≤ T1 ,
20lgTω, ω ≥ 1 ,
T
где частота ω = T1 называется сопрягающей.
Асимптотическая ЛФХ есть кусочно-постоянная функция, получаемая заменой в
(30) членов arctg(Tω) на |
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
0, |
ω ≤ |
|
|
|
, |
|
||
T |
|
|||||||
|
|
|
|
|
|
|||
π |
|
|
|
1 |
|
|
||
|
, ω ≥ |
|
. |
|||||
|
|
T |
||||||
2 |
|
|
|
|||||
Нетрудно заметить, что для систем с передаточной функцией вида (2.29) по асимптотической ЛАХ можно восстановить асимптотическую ЛФХ и саму передаточную функцию. Данное свойство характерно для любой передаточной функции, не имеющей нулей и полюсов в правой полуплоскости.
Преобразуем требования к замкнутой системе в ограничения на свойства логарифмической амплитудной характеристики разомкнутой системы. Условие (2.21) на отработку линейно растущего задающего воздействия сводятся к ограничению на
поведение Wp ( jω) при низких частотах, близких ω = 0 : |
|
|||||||
Wp ( jω) ≈ |
kp |
|
(2.32) |
|||||
jω |
||||||||
|
|
|
|
|
|
|||
и, следовательно, |
|
|||||||
Lp (ω) = 20lg |
|
Wp ( jω) |
|
= 20lg kp − 20lg ω. |
(2.33) |
|||
|
|
|||||||
Иначе говоря, низкочастотная асимптота Lp (ω) должна иметь наклон -20 дБ/дек,
причем в силу (2.24) ее уровень определяется условием
18
|
20 lg k p ≥ −20 lg δ1 . |
(2.34) |
|||||||||||
Условия (2.26), (2.28) при достаточно малых δy , δ f можно заменить на |
|
||||||||||||
|
W |
p |
( jω) |
|
|
≥ δ−1 , ω ≤ ω , |
(2.35) |
||||||
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
y |
y |
|
|||
|
Wp ( jω) |
|
≤ δf |
, ω ≥ ωf |
(2.36) |
||||||||
|
|
||||||||||||
и, следовательно, |
|
|
|
|
|
|
|
|
|
|
|
|
|
Lp (ω) = 20lg |
|
Wp ( jω) |
|
≥ −20lg δy , ω ≤ ωy , |
(2.37) |
||||||||
|
|
||||||||||||
Lp (ω) = 20lg |
|
Wp ( jω) |
|
≤ 20lg δf , ω ≥ ωf . |
(2.38) |
||||||||
|
|
||||||||||||
Условие (2.37) задает ограничение на поведение логарифмической амплитудной частотной характеристики в области низких частот, а условие (2.38) - в высокочастотной области.
Поведение ЛАХ в области средних частот определяет запасы устойчивости по фазе и амплитуде и в значительной мере качество системы в переходном режиме (время переходного процесса и перерегулирование). Для обеспечения приемлемых запасов устойчивости наклон ЛАХ на частоте ωc такой, что Lp (ωc ) = 0 обычно выбира-
ется равным -20 дБ/дек, причем длительность этого участка должна быть не менее декады, что соответствует изменению частоты в 10 раз. Частота ωc называется час-
тотой среза. Значение частоты среза надо выбирать наиболее большим из всех возможных для того, чтобы увеличить быстродействие замкнутой системы.
Итак, все требования, которые были сформулированы в задании курсовой роботы, сведены к ограничениям на допустимое поведение ЛАХ разомкнутой системы (см. рис. 2.4, где заштрихованы границы зон, в которые не может заходить ЛАХ). Оста-
ется подобрать передаточную функцию разомкнутого контура Wp (s) , для которой эти ограничения выполнены. Если Wp (s) минимально-фазовая передаточная функ-
ция, то вначале строят асимптотическую ЛАХ, удовлетворяющую всем ограничени-
ем, а затем по ней находят саму Wp (s) .
Приведем в готовом виде сводку допустимых передаточных функций разомкнутой системы:
19
Lp (ω) |
|
|
|
20 lg k p |
-20 дБ/дек |
|
|
|
|
|
|
−20lg δy |
|
|
|
|
-20 дБ/дек |
|
−lg δ1 |
20 lg δf |
lg ωy lg ωc |
lg ωf |
lg ω |
|
|
|
Рис. 2.4. Ограничения на допустимое поведение ЛАХ
1) Пусть |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
≤ δ1 |
≤ |
δy |
. |
(2.39) |
|||
|
δ |
|
ω |
|
ω |
|
||||||
|
f |
f |
|
y |
|
|
||||||
|
|
|
|
|
|
|
||||||
Тогда передаточная функция |
|
|
|
|
|
|
|
|
|
|
|
|
Wp (s) = |
kp |
, kp = δf ωf |
(2.40) |
|||||||||
|
||||||||||||
|
|
|
|
|
|
s |
|
|
|
|
|
|
удовлетворяет всем ограничениям. При тех же условиях допустима, но требует меньшего усиления на частотах ω ≥ ωy и обеспечивает лучшее подавление помех передаточная функция вида
|
|
|
Wp (s) = |
|
|
|
|
kp (T2s +1) |
|
|
|
, |
|
(2.41) |
||||||||||
|
|
|
|
s(T1s +1)(T3s +1) |
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
δy |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
k |
p |
= δ−1 |
, T = |
|
|
|
|
|
|
, T = |
|
|
|
µδy |
|
|
, T ≈ 0.1T , |
(2.42) |
||||||
|
1 |
1 |
|
δ1 ωy2 |
2 |
|
|
ωy |
|
|
|
|
3 |
2 |
|
|||||||||
а параметр µ, характеризующий расположение частоты среза ωc |
= µ / T2 , может вы- |
|||||||||||||||||||||||
бираться в пределах 2÷4. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
2) Пусть |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
δ1 |
≤ |
|
δy |
|
, δ1 < |
|
|
|
|
1 |
|
|
|
. |
|
|
(2.43) |
||||
|
|
|
|
ω |
y |
δ |
f |
ω |
f |
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
20
