GOS for Iphone / mobile / ТАУ / Книги / ТАУ учебник
.pdf
где τ – длительность импульса на выходе модулятора; Tn – период следования импульсов на выходе модулятора, величина которого зависит от сигнала x(t) на входе модулятора.
2.2. Постановка задачи синтеза и общая схема ее решения
Задача синтеза систем автоматического управления, содержащих модуляторы различного вида, рассматривается в следующей постановке. Предполагается, что известна структура синтезируемой САУ и параметры объекта управления. Параметры регулятора (оператора управления), структура которого задана в самом общем виде, определяются из условия приближенного обеспечения заданных показателей качества работы САУ в переходном режиме (времени переходного процесса – Tп.п; перерегулирования – σ ; колебательности – µ). При этом, безусловно, должна обеспечиваться устойчивость и грубость системы по варьируемым параметрам.
Ввиду того, что число искомых параметров может быть любым, оператор управления структурно может быть задан со значительной избыточностью. В этом случае после определения значений искомых параметров в результате применения методов теории чувствительности, определяющих координаты системы, чувствительные к варьируемым параметрам, часть этих параметров может быть принята равной нулю (бесконечности), что приводит к упрощению оператора управления и выявлению тем самым наиболее целесообразной его структуры.
Как правило, задача синтеза решается при технических ограничениях, которые накладываются на значения варьируемых параметров:
ck− ≤ ck ≤ ck+ k= , … m |
(2.11) |
где ck+ – максимально допустимые значения варьируемых параметров; ck– – минимально допустимые значения варьируемых параметров.
Ограничения на грубость системы по варьируемым параметрам имеют следующий вид:
∆ = |
δ, ck ≤ ∆ 0 |
(2.12) |
|
ck |
|
где ∆ 0 – заданное значение грубости системы; δ ck – вариации параметров, в пределах которых обеспечивается устойчивость системы.
31
Для определенности задачу синтеза рассмотрим при внешнем скачкообразном входном воздействии f(t) = H1(t) и нулевых начальных условиях для момента времени t = –0, т. е. до приложения к системе воздействия амплитудой H
x |
|
= |
x 0, |
= |
x 0, |
= |
(n−1) |
(2.13) |
−0 |
… x0, , = |
|||||||
|
|
−0 |
|
−0 |
|
−0 |
|
Так как при синтезированных параметрах система должна быть устойчива, то
x(∞ )= |
|
|
… |
, x |
(n−1) |
0. |
(2.14) |
H , x∞( =) |
0, ∞x(= ) 0, |
|
= |
|
Выбираем систему из m непрерывно дифференцируемых линейнонезависимых координатных функций
ϕ 1 (t )ϕ, 2 (t ),…ϕ , q (t ),…ϕ , m (t ). |
(2.15) |
В соответствии с требуемыми показателями качества работы синтезируемой системы управления в переходном режиме зададимся желаемым программным движением
x0 (t) = Ω0 (t) + ∑l |
aiΩi (t), i = 1,2,…,l, |
(2.16) |
i=1 |
|
|
где Ω 0(t) = ω 0(t)1(t) – функция, удовлетворяющая заданным граничным (начальным (2.13) и конечным (2.14)) условиям; Ω i(t) = ω i(t)1(t) – функции, удовлетворяющие однородным граничным условиям; ai – известные коэффициенты.
Анализ структурных схем импульсных систем автоматического управления, содержащих один нелинейный элемент, которые представлены на рис. 2.2, показывает следующее:
– уравнение движения системы, структура которой показана на рис. 2.2, а, записанное относительно координаты входа нелинейного элемента, имеет вид
N (ck , D) x (t) + M (ck , D)F x (t), x (t) = N (ck , D) f (t);
– динамика САУ, приведенной на рис. 2.2, б, относительно входа нелинейного элемента описывается уравнением
N (ck , D ) x (t ) + M (ck , D )F x (t ), x (t ) = N (ck , D ) f (t );
32
а) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x*(t) |
|
|
|
|
M (ck , p) |
|
z(t) |
||||
|
f(t) |
|
|
|
|
ИЭ |
F(x,px) |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
N (ck , p) |
|
|
|
||
|
|
|
– |
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
б) |
f(t) |
f*(t) |
x*(t) |
|
M (c , p) z(t) |
|
ИЭ |
|
F(x,px) |
k |
|
|
|
|
N (ck , p) |
||
|
|
|
– |
|
|
|
|
|
|
|
|
|
|
|
|
z*(t) |
ИЭ |
|
|
|
|
|
в) |
f(t) |
ИЭ |
M (ck , p) x(t) |
F(x,px) |
L (ck , p) |
z(t) |
|
|
|
N (ck |
, p) |
|
K (ck , p) |
|
|
|
|
|
|
|
|||
|
|
– |
|
|
|
|
|
г) |
f(t) |
M (ck , p) |
L (ck , p) |
x(t) |
|
||||
|
|
ИЭ |
K (ck , p) |
|
|
|
N (ck , p) |
|
|
|
|
– |
– |
|
F(x,px) 
д) |
f(t) |
x(t) |
|
M (ck , p) |
z(t) |
|
|
ИЭ |
|||||
|
|
F(x,px) |
N (ck |
, p) |
|
|
|
|
- |
|
|
||
|
|
|
|
|
|
|
Рис. 2.2
– для структурной схемы рис. 2.2, в уравнение движения будет следующим:
N(ck , D)K (ck , D) x (t) + L(ck , D)M (ck , D)F x (t), x (t) =
=L(ck , D)M (ck , D) f (t);
33
– в том случае, если структура системы управления имеет вид, показанный на рис. 2.2, г, динамика САУ относительно входа нелинейного звена описывается следующим образом:
N (c , D)K (c , D)x(t) + L(c , D)M (c , D)x (t) + |
|
|
||||||
k |
|
k |
|
k |
|
k |
|
|
+L(c , D)N (c , D)F |
x(t), x(t) = L(c , D)M (c , D) f |
|
(t); |
|||||
k |
k |
|
|
|
k |
k |
|
|
|
|
|
||||||
– наконец для САУ, показанной на рис. 2.2, д, уравнение движения будет
N (ck , D) x (t) + M (ck , D)F x (t), x (t) = N (ck , D) f (t).
Обобщение результатов анализа частных случаев структур систем управления, приведенных на рис. 2.2, показывает, что система автоматического управления, содержащая модулятор и нелинейный элемент, в общем виде описывается следующим дифференциальным уравнением:
Q (c , D) x (t) |
+ Q (c , D) x |
(t) + |
|
|
|||||||
+R (ck , D) y (t) + R |
|
k |
|
|
|
k |
|
|
|
|
(t) , |
|
|
|
|
|
|
|
|
|
|
||
|
(ck , D) y |
|
(t) = S (ck , D) f (t) + S |
(ck , D) f |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
y (t) = F |
x (t), x (t) |
, |
y (t) = F |
|
x |
(t), x |
(t) , |
(2.17) |
|||
|
|
|
|
|
|
|
|
|
|
|
|
где x(t), x*(t) – исследуемая координата на входе и выходе модулятора, соответственно, относительно которой записано уравнение движения синтезируемой САУ; f(t), f*(t) – внешнее входное воздействие
на входе и выходе модулятора соответственно; |
y (t ) = F x (t ), x (t ) , |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
– нелинейные функции; |
|
|
|
|
|
|
|||||||
y (t ) = F x |
(t ), x |
(t ) |
|
|
|
|
|
|
|||||||||
|
Q (c , D ) = |
n |
|
(c |
)Di ; Q (c , D ) |
|
|
n |
a |
(c )Di |
|
|
|
||||
|
∑ |
a |
|
= |
|
; |
|
|
|||||||||
|
|
k |
|
i |
k |
k |
|
∑ |
|
|
i |
|
k |
|
|
|
|
|
|
|
|
i=0 |
|
|
|
|
i=0 |
|
|
|
|
|
|
|
|
|
R (c |
|
, D ) = |
u |
|
(c |
)Di ; R (c , D ) |
|
u |
b |
|
(c |
)Di ; |
|
|
|
|
|
|
∑ |
b |
= |
|
|
|
|
|
||||||||
|
k |
|
i |
k |
k |
∑ |
|
i |
|
k |
|
|
|
|
|||
|
|
|
|
i=0 |
|
|
|
|
i=0 |
|
|
|
|
|
|
|
|
|
S (c |
|
, D ) = |
v |
|
(c |
)Di ; S (c , D ) |
|
v |
e |
|
(c |
)Di |
|
|
|
|
|
|
∑ |
e |
= |
|
|
|
|
|
||||||||
|
k |
|
i |
k |
k |
∑ |
|
i |
|
k |
|
|
|
|
|||
|
|
|
|
i=0 |
|
|
|
|
i=0 |
|
|
|
|
|
|
|
|
34
– полиномы оператора обобщенного дифференцирования D с вещественными постоянными коэффициентами степеней n, n*, u, u*, v, v* соответственно.
Очевидно, что при описании динамики импульсных САУ с одним нелинейным элементом в частных случаях часть слагаемых уравнения (2.17) может отсутствовать. Также необходимо отметить, что запись уравнения движения относительно координаты входа нелинейного звена, которая впервые была предложена И. А. Орурком, дает несомненные преимущества при реализации метода синтеза систем на основе обобщенного метода Галеркина. Это связано с упрощением процедуры определения соотношений вида «вход – выход» интегралов Галеркина, что будет подробно показано в следующих главах книги.
В уравнении (2.17) применяется универсальная координата времени. Это дает возможность использовать дискретно-непрерывные модели систем, определяющих их описание на каждом из интервалов дискретности, и позволяет без перехода к разностным уравнениям, которые требуют получения аналитических решений нелинейных нестационарных дифференциальных уравнений, решать задачу синтеза обобщенным методом Галеркина с единых математических позиций для САУ широкого класса.
Поставим желаемое программное движение (2.16) в уравнение движения системы (2.17) и образуем невязку
ψ(ck ,t) = Q (ck , D) x0 (t) + Q (ck , D) x0* (t) +
+R (ck , D)F x0 (t), D{x0 (t)} +
+R (ck , D)F x0 (t), D{x0 (t)} −
−S (c , D) f (t) − S |
(c , D) f |
(t). |
(2.18) |
k |
k |
|
Если предположить, что система с синтезированными параметрами заведомо устойчива, то значения искомых параметров определяются из условия ортогональности невязки (2.18) координатным функциям (2.13)
∞ |
|
∫ ψ (ck ,t )ϕ q (t )dt= 0, k, q= 1,2,…,m, |
(2.19) |
0 |
|
что приводит к следующей системе алгебраических уравнений:
35
∞ |
Q(c ,D)x0 |
(t)ϕ |
q |
(t)dt +∞ |
∫ |
Q (c ,D)x0 |
(t)ϕ |
q |
(t)dt + |
∫ |
k |
|
|
k |
|
|
|||
0 |
|
|
|
|
0 |
|
|
|
|
+∞∫ R(ck ,D)F x0 (t),D{x0 (t)} ϕ q (t)dt+
0
|
+ |
∞ |
|
|
(c ,D)F |
|
|
0 |
(t),D x |
0 |
(t) |
|
ϕ |
|
(t)dt− |
||||||||
|
∫ |
R |
x |
|
|||||||||||||||||||
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
k |
|
|
|
|
|
|
|
{ |
|
|
} |
|
q |
|
|
|
|||
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
−∞ |
S (c ,D)f (t)ϕ |
q |
(t)dt −∞ |
∫ |
S |
(c ,D) f |
(t)ϕ |
q |
(t)dt = 0, |
||||||||||||||
∫ |
k |
|
|
|
|
|
|
|
|
|
|
k |
|
|
|
|
|
|
|||||
0 |
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
k,q = 1,2,…,m. |
|
|
|
|
|
|
(2.20) |
||||||||
Решая систему из m алгебраических уравнений (2.20), определяем значения варьируемых параметров оператора управления. Так как задача синтеза решается при ограничениях на значения искомых параметров, наложенных исходя из возможности их технической реализации, ограничениях на устойчивость и грубость САУ с синтезированными параметрами, а также в силу того, что имеет место нелинейная зависимость между варьируемыми параметрами, то строго равенство (2.19) выполняться не будет. Поэтому задача синтеза параметров обобщенным методом Галеркина в вычислительном плане представляет собой задачу нелинейного программирования с целевой функцией, построенной на основе уравнений (2.20) и имеющей вид
m |
∞ |
|
|
2 |
|
|
|
J = ∑ |
|
|
ψ (ck ,t )ϕ |
|
, minc |
J→ 0, |
|
|
∫ |
q (t )dt |
(2.21) |
||||
q=1 |
|
|
|
|
k |
|
|
|
|
0 |
|
|
|
|
|
оптимум которой определяется при ограничениях, отмеченных выше, путем использования известных методов поиска экстремума функционала [89, 90].
2.3. Построение математической модели желаемого программного движения произвольно высокого порядка
Проблема обеспечения заданных показателей качества функционирования систем управления при решении задачи параметрического синтеза вари-
36
ационными методами, в том числе и обобщенным методом Галеркина, непосредственно связана с необходимостью решения следующих вопросов:
–задание желаемого программного движения произвольно высокого порядка, в общем случае соответствующего порядку синтезируемой системы;
–установление взаимосвязи между параметрами программного движения и показателями качества работы синтезируемой системы в переходном и установившемся режимах.
При синтезе нелинейных САУ обобщенным методом Галеркина в первом приближении целесообразно задавать желаемое программное движение для гармонически линеаризованной системы, поскольку задать в нелинейной системе управления заведомо реализуемый процесс достаточно сложно. В этом случае динамика системы будет описываться уравнением
Q (ck , p ) x (t ) = S (ck , p ) f (t ), |
(2.22) |
где p = d/dt – оператор дифференцирования.
Характеристическое уравнение однородного дифференциального уравнения (2.22) в случае известных значений параметров ck будет
λn + an−1λn−1 +…+ a1λ + a0 = 0, |
(2.23) |
где λ – корни уравнения; an−1 – коэффициенты уравнения.
Уравнение (2.23) имеет n корней, часть которых может быть вещественными, часть – комплексно-сопряженными. Для устойчивой системы, как известно, вещественные корни и вещественные части комплексно-сопря- женных корней должны быть отрицательными, т. е. быть расположенными в левой части комплексной плоскости. Тогда можно записать следующее:
λ1,2 = −α1 ± jβ1; λ3,4 = −α2 ± jβ2;
λk −1,k = −αk ± jβk ;
22
λk +1 = −αk +1;
2
|
|
λi = −αi ; |
|
λ j = λ j+1 = = λ j+r = −αm , |
(2.24) |
|
37 |
где k – четные числа.
Корни, которые имеют наименьшее (наибольшее) абсолютное значение вещественной части, в дальнейшем будем называть младшими (старшими).
Решением уравнения вида (2.22) будет следующее выражение [190]:
|
k |
|
cosβ i t + Ci |
|
|
|
|
|
||||
x0 (t) = ∑ |
e−αit Ci |
sinβ i |
t |
+ |
|
|||||||
|
i=2 |
|
|
|
|
|
|
+1 |
|
|
|
|
|
2 |
2 |
2 |
2 |
|
|
|
|||||
+ ∑j |
Cle−αlt + e−αmt (C j+1 + C j+2t + + C j+r+1tr ), |
(2.25) |
||||||||||
l=k+1 |
|
|
|
|
|
|
|
|
|
|
|
|
где i – ряд четных чисел от двух до k; l – ряд натуральных чисел; αi,...,αm
– вещественные части корней характеристического уравнения (коэффициенты затухания составляющих); β1,..., βk/2 – частоты собственных колебаний системы управления; Ci/2,..., Cj+r+1 – постоянные, определяемые начальными условиями.
Таким образом, в уравнении (2.25) первая сумма от i = 2 до i = k определяется комплексно-сопряженными корнями λ 1,2,..., λ k–1,k; вторая сумма от l = k+1 до l = j определяется вещественными корнями λ k+1,..., λ j; последняя сумма определяется вещественными кратными корнями
λj+1,..., λ j+r (рис. 2.3).
Вряде работ [191–197] предлагается оценивать качество переходно-
го процесса интегральными критериями качества вида
t |
|
|
x (t ) |
|
|
|
|
t |
|
|
|
I = ∫ t |
n |
|
|
n |
dt, |
либо |
I = ∫ t |
n |
|
n |
|
|
|
||||||||||
|
|
|
|
|
|
x (t ) |
dt, |
||||
0 |
|
|
|
|
|
|
|
0 |
|
|
|
где n у различных авторов принимает значение 0, 1, 2.
Однако обратный переход от показателей интегральной оценки к параметрам желаемого программного движения достаточно сложен, а подчас невозможен.
Как показывает практика проектирования систем автоматического управления различных по своему назначению и функциональным возможностям, затруднительно дать общие рекомендации по заданию желаемого программного движения при решении задачи синтеза оператора управления.
В то же время, учитывая специфику работы проектируемой системы, можно задавать квазиоптимальные переходные процессы, т. е. про-
38
|
|
|
|
|
|
|
|
|
λ |
|
|
|
+j |
|
|
|
|
|
|
|
|
|
|
4 |
|
−α2 |
|||
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
+jβ2 |
|
|
|
|
|
λ 2 |
|||
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
+jβ |
|
−α1 |
|||||
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
||||||||
λ m |
λ i |
λ k+1 |
|
|
|
|||||||||
|
|
|
|
|
|
|||||||||
|
|
0 |
||||||||||||
|
|
|
||||||||||||
|
|
|
||||||||||||
|
|
|
||||||||||||
|
|
|
||||||||||||
|
|
|
||||||||||||
|
|
|
|
|
−jβ2 |
|
|
−jβ1 |
||||||
|
|
|
|
|
|
|
|
λ |
||||||
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
1 |
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
λ 2 |
|
λ k/2+1 |
|
||
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|||
λi
λm
−j
Рис. 2.3
цессы, оптимальные по одному из показателей качества работы САУ в переходном режиме: величине перерегулирования, длительности переходного процесса, быстродействию [190].
Известно, что, чем меньше величина перерегулирования σ , тем качество переходного процесса лучше ( при выполнении требований к другим показателям качества). Перерегулирование зависит от корней уравнения (2.23) и начальных условий. Причем в случае внешнего скачкообразного воздействия f = H1(t) перерегулирование в системе будет лишь в случае комплексно-сопряженных корней при достаточно больших значениях собственных частот β.
Как показано в [190, 195], при вещественных корнях, комплексносопряженных корнях с малыми значениями β, а также при смешанных корнях, когда младшими являются вещественные корни, переходный процесс протекает монотонно.
Для многих САУ длительность переходного процесса Tп.п (промежуток времени, в течение которого регулируемая величина с заданной точностью достигает установившегося значения x является одним из важнейших показателей их оптимальности. При этом минимальная длительность переходного процесса при f(t) = H1(t) достигается при малых
39
мнимых частях комплексно-сопряженных корней, либо при смешанных корнях, при которых σ ≈ (1–4)% [195].
При увеличении мнимых частей (собственных частот колебаний β) комплексно-сопряженных корней возрастают как Tп.п, так и σ . Поэтому с точки зрения скорейшего затухания желаемого программного движения необходимо, чтобы вещественные части всех корней (коэффициенты затухания составляющих) были, возможно, большими. Однако наилучшие результаты получаются при равенстве вещественных частей всех корней уравнения (2.23), поскольку их сумма численно равна первому коэффициенту a0 характеристического уравнения.
В случае только вещественных корней САУ становится более устойчивой, однако при этом существенно возрастает длительность переходного процесса, причем тем больше, чем значительнее различие в коэффициентах затухания α отдельных составляющих процесса.
Под быстродействием САУ будем понимать скорость нарастания переходного процесса, что характеризует способность системы управления за минимальное время достигнуть первого согласованного с задающим устройством положения. В работах [190, 192] быстродействие характеризуется максимальной скоростью xmax и средним ускорением xср изменения регулируемой величины x во времени. Как отмечается в [190], максимальные значения характеристик быстродействия xmax и xср могут быть получены при комплексно-сопряженных корнях уравнения (2.23) с большими значениями β и малыми α, т. е. при больших перерегулированиях.
Малым быстродействием будет обладать САУ, характеристическое уравнение которой дает лишь вещественные составляющие процесса (2.23) при выполнении условия
α1 << α2 << α3 << << αn .
В ряде работ [192, 193, 195, 196] предлагаются некоторые стандартные (типовые) уравнения, дающие оптимальные переходные процессы. Приводятся не только значения коэффициентов стандартных характеристических уравнений [143, 197] и численные значения корней [190, 195], но также и значения показателей качества переходных процессов, соответствующих рассматриваемому распределению корней. К сожалению, в известных работах не рассматриваются уравнения выше 8-го порядка, что несколько снижает их ценность для практического применения в случае синтеза САУ более высокого порядка.
40
