книги / Теория автоматического управления
..pdfРис. 2.17 Модель динамического объекта (2.163) в переменных состояния
“ 0 |
1 |
0 |
О |
0 |
0 |
1 |
О |
А =
— ап 7 _ а0
"0 ~
0
В
1
_ |
_ пх 1 |
—ал-1. |
|
■ai |
а0 7 |
|
ny.lt |
|
Ь/П |
|
|
Ьщ-1 |
|
Ст= |
D = О |
(2.164) |
|
Ьг |
|
_ |
_ (m+ l) X 1 |
|
Соответствующая матрицам (2.164) алгоритмическая схема по» казана на рис. 2.17. При этом переменные состояния объекта свя заны соотношением (2.158), т. е. являются фазовыми.
Если в передаточной функции (2.163) полиномы числителя и зна менателя имеют одинаковый порядок, т. е. т = л, то матрица D = b,Ja0 Ф 0. Но при этом матрицы А, В, С имеют вид, отличный от (2.164), а переменные состояния уже не удается выразить как фазовые.
Покажем, что объекту с передаточной функцией (2.163) в про странстве состояний действительно соответствуют матрицы (2.164). Преобразуем дифференциальное уравнение
(Оорп+ а1рп~*+ |
+ ап) хв (t) = (Ь0рт+ |
|
+ b1ffn~i + |
+ Ьт)У (0 (m e n ) |
(2.165) |
к нормальной форме Коши (2.156) или (2.157). Для этого введем
вспомогательную |
переменную |
х г (t) и запишем уравнение |
(2.165) |
в виде пропорции: |
|
|
|
*в ШЬор™+ Ьхрт~' + |
+ bm) = y (/)/(аоР" + |
|
|
+ a1pn~l + |
+ an) = x1(t). |
(2.166) |
Очевидно, что соотношению (2.165) соответствуют два уравнения:
(b0prn + b1pm~l + |
+ bm)x1(t) = xB(t), |
(2.167) |
||
(a0pn+ a1pn~i + |
+ a„)x1(t) = y(t). |
(2.168) |
||
Введем обозначения |
|
|
|
|
x,{t) = x2{t), x2(t) = *з(0. |
xn_1(t)= xn(t), |
(2.169) |
||
тогда вместо уравнения (2.168) получим |
|
|||
х п (0 = |
ао 1[a ix n (t) + агх„_1 (/) + |
фй'пХi (01 |
°о V (0- |
|
|
|
|
|
(2.170) |
Объединяя последнее уравнение с выражениями (2.158), получим систему уравнений в форме Коши
xi (0 = хы (0. |
1= 1; 2; |
п — 1, |
|
|
х п (0 = |
а 0 |
[a.iXn (t) -f- а,-2.хп_1 (t) -f- |
(2.171) |
++ a nx1(t)]+cbly(t).
Применяя соотношения (2.169) к уравнению (2.167), найдем выра жение для выходной переменной
хв(0 — ЬдХт+1 (t) -(- biXm (t) -f- |
-(- bmx-y (t). |
(2.172) |
Нетрудно заметить, что именно при матрицах вида (2.164) по лученные уравнения (2.171) и (2.172) эквивалентны уравнениям состояния (2.156) и выхода (2.159).
82
От модели объекта, записанной по способу ПС, можно перейти к описанию по способу ВВ. Например, матричная передаточная функция объекта между вектором управления у (/) и вектором вы хода х в (t), согласно уравнениям (2.156) и (2.159), равна (при 2 (0 = 0 и g (t) = 0)
W 0(Р) = Хв (р)/у (Р) = С (pr— A)-1 B + D. |
(2.173) |
Изменение вектора состояния х (t) во времени при известном векторе управления у (t) определяется следующим решением не однородного уравнения (2.156):
л: (0 = Ф (0 JC (0) + $ Ф (t— Ъ) By (Ф) сЮ, |
(2.174) |
|
о |
|
|
где х (0) — начальное положение |
вектора х (t) в момент t = 0; |
|
Ф (/) = [Ф,-/ {t)\ny.n — переходная |
(фундаментальная) |
матрица, |
определяющая так называемую фундаментальную систему решений однородного уравнения х (t) = А х (t):
(0 |
" Ф ц(0 |
фщ (0 |
" -* i(0 ) |
- |
|
= |
|
|
(2.175) |
-Xn(t) _ |
_ Ф„1(0 |
Фnn(t) |
_ _ * я (0) |
_ |
Каждый элемент Ф£|- (/) фундаментальной матрицы описывает процесс перехода i-й переменной состояния, когда /-я переменная состояния имеет начальное значение х,- (0) = 1.
Первое слагаемое в выражении (2.174) переходного процесса определяет свободную составляющую х св (0, зависящую только от начальных условий х (0) и собственных свойств объекта, т. е. от матрицы А, а второе слагаемое — вынужденную составляющую
*вын (t), |
зависящую |
как |
от собственных свойств объекта |
(матриц |
А и В), так и от входных воздействий (вектора у (/)). |
|
|||
Через |
переходную |
матрицу Ф (t) можно выразить матричную |
||
весовую |
функцию w (t) |
для выходного вектора х в (при |
D = 0): |
|
w(t) = CO>(t)B. |
|
|
(2.176) |
Свободную и вынужденную составляющие вектора х (t) можно определить с помощью преобразования Лапласа. Продемонстрируем этот способ на примере определения свободной составляющей. Однородное дифференциальное уравнение, описывающее свобод ное движение, имеет вид
x(t) = Ax(t), |
(2.177) |
соответствующее ему преобразованное уравнение — |
|
рх{р) — х(0) = Ах{р) или х (р) = (р/—Л)-1 х (0), |
(2.178) |
отсюда
х (t) = х св (0 = .S?'1 {х (р)} = i?"1 {[р/— Л г 1 л: (0)} = eAl х (0).
(2 Л79)
Из сравнения выражения (2Л79) с первым слагаемым в формуле (2Л74), характеризующим лгсв(0. следует, что фундаментальная матрица всегда имеет вид
Ф (/) = еА‘. |
(2.180) |
Элементы матрицы Ф (t) можно определить разложением мат ричной экспоненциальной функции (2Л80) в степенной ряд
|
оо |
|
|
|
ф(<) = е'4' = ^ Г ^ - = = / + |
Л / + ^ р - + |
(2Л81) |
||
|
t=0 |
|
|
|
Собственные |
динамические |
свойства |
объекта |
можно оценить |
и без решения |
уравнения (2.177) — по |
собственным значениям X |
матрицы А , которые определяются по ее характеристической мат рице
Л —XI = |
ап |
X] |
а 1п |
(2.182) |
а 2Ъ |
й22 Х\ |
а2п |
||
- |
&пЪ |
а п2ч |
О'пп |
^ - |
Эти собственные значения отыскиваются как корни характери
стического уравнения |
матрицы |
\ A - X I | - 0 , |
(2.183) |
левая часть которого представляет собой многочлен степени п относительно X
(ап X) (а22 А,) |
(агт-—X) —0 |
(2.184) |
или |
|
|
(с№ + сгХп- ' + |
+ сп) (— 1)" = 0. |
(2.185) |
Собственные значения Xt матрицы А совпадают с полюсами pi передаточной функции (2.173) объекта. Это утверждение следует из известного свойства определителя квадратной матрицы
\А —XI\ = (— 1)" IXI— А |. |
(2.186) |
Пример. Составим уравнения состояния и выхода смесительного бака (рис. 2.18) емкостью V (м3), в который одновременно подаются два регули руемых потока жидкости Qx и Q2 (м3/с), содержащие один и тот же компонент с постоянными концентрациями с10 и с20 (кг/м3). В качестве выходных пе ременных этого объекта управления будем рассматривать расход Q, концен трацию с компонента в выходном потоке Q и уровень h (м) жидкости в баке, а в качестве входных — расходы Q± и Q2.
Рис. 2.18. Смесительный бак как объект управле ния
Уравнения материального баланса имеют вид: по массе жидкости
d V (/)/d t = Qx (i) + Q2 (t) — Q (/) |
|
|
(2.187) |
|||||||
|
по массе компонента |
|
|
|
|
|
||||
|
d [с (/) V (/)]/d t = |
|
c,,,!?! (/) + |
c2 ,Q2 (t) — c(t)Q (/). |
|
(2.188) |
||||
|
Расход |
жидкости |
из бака согласно законам |
гидравлики |
|
|||||
|
Q ( t ) = k |
л / Щ |
= |
k V V (/)/S |
|
|
|
(2.189) |
||
где |
/г — коэффициент |
истечения; |
5 — площадь |
поверхности жидкости, м2. |
||||||
q (/) |
С учетом (2.189) |
последнее слагаемое в (2.188), соответствующее расходу |
||||||||
компонента из |
резервуара, |
равно |
|
|
|
|||||
|
q (t) = с (/) Q (/) - |
|
с (/) k У V(t)/S |
|
|
(2.190) |
||||
|
Для некоторого фиксированного установившегося режима, характери |
|||||||||
зуемого значениями |
У0, Qo> со» |
h0 и соотношениями |
|
|
||||||
|
Рм + Рг> = Р». |
|
Q„ = k ^ V J s ’ |
|
|
(2.191) |
||||
можно нелинейную |
функцию (2.189) двух переменных линеаризовать: |
|||||||||
|
Q = f (С, V) « |
р„ -1- с» 1 ^ |
^ |
~^Г ДУ + |
* Д / - J - |
Дс- |
(2.192) |
|||
|
Аналогично линеаризуются |
и остальные слагаемые |
уравнений |
(2.187) |
и (2.188), представляющие собой произведение двух переменных или нели нейную зависимость от V.
Введем следующие обозначения для отклонений всех переменных от их
значений в рассматриваемом режиме: |
|
AQi W = Pi (/) - Qv, = У1 (0; Д<?2 (0 = |
W - Р*. = у* (0;' |
А V (0 = |
V (/) - V„ = |
it, (/); |
Дс (/) - Со = х2 (/) = *в2 (/); |
(2.193) |
АЛ(/) = |
А(0 — A.. = |
*BSW; |
AP(0 = PW — Р« = *в1(0- |
|
Тогда, учитывая соотношения (2.191), можно исходные уравнения (2.187) и (2.188) записать в отклонениях от установившегося режима:
W |
У\ (0 + Уг (*)----- |
|
|
|
ZVп |
|
|
V QX2 (О “Ь СКХ1 (0 ^ С10У1 ( 0 “1” ^2 |
У2 (О |
(2.194) |
|
|
gpOo X! (f)-QoX2(t).
2К0
Обозначая V0/Q0 = Т и исключая из второго уравнения х г (/), можно получить искомые уравнения состояния:
xi (t) = |
------- —1 Х\ (/) + |
У\ (0 + Уг (*); |
|
2 Т |
(2.195) |
|
|
|
Х2 (0 = ----- |
--- х2 (0 ---- |
--------- У\ (О Н “ “--- “ ^2 (0. |
Выходные переменные Q (t), с (/) и h (t) однозначно выражаются через переменные состояния (в отклонениях):
*В1(/) = до(0= 0 (0- Qo» |
(0; |
|||
*в2 (0 = |
Ас (0 = |
с (() — с0 = |
х2 (0; |
(2.196) |
*вз (0 = |
ДА (/) = |
ДУ (t)IS = |
- i - хх (0. |
|
|
|
|
о |
|
Очевидно, что для рассматриваемого объекта матрицы А, В я С равны:
|
Г -1Й Г; |
0 I |
|
|
-1/2 Т; |
0- |
|
В = \ 1; |
Ч ; |
С = 0; |
1 |
||
|
|_0; |
— 1/7Т |
||||
|
1_^21» |
^22J |
.1/S; |
0. |
||
|
|
|
|
|
||
|
|
|
|
|
|
(2.197) |
где |
b21 = (c10—c0)/V0; b22 = |
(c20—c0)/V0. |
|
|
|
|
2.10. Элементарные операции аналогового |
|
|
||||
и |
цифрового моделирования |
|
|
|
Для экспериментального исследования динамических свойств систем управления широко применяется машинное моделирование, осуществляемое на различных видах ЭВМ. Используются два ос новных метода машинного моделирования: метод аналогового моде лирования, реализуемый на АВМ, и метод цифрового моделирова ния, реализуемый на ЦВМ.
При аналоговом моделировании решение динамической задачи воспроизводится в виде электрического переходного процесса, адек ватного свойствам исследуемой системы управления. Результат
Рис. 2.19. Операционные элементы аналоговых моделей:
а — потенциометр; в — операционный усилитель; д — операционный усилитель интегри рующий; б, г, е — соответствующие условные обозначения элементов
решения наблюдают на экране осциллографа или регистрируют с помощью записывающих приборов.
Цифровое моделирование заключается в выполнении определен ной последовательности вычислительных процедур, соответствую щих численному решению дифференциальных уравнений системы. Результат решения получается на экране дисплея или выводится на печатающее устройство.
При аналоговом и цифровом моделировании систем управления наиболее широко применяется методика структурного моделиро вания, согласно которой в машинной модели воспроизводится фи зическая структура исследуемой системы, расчлененная на эле ментарные математические операции (интегрирование, умножение на постоянный коэффициент, суммирование и др.). Эти элементар ные операции реализуются при аналоговом моделировании на опе рационных усилителях, а при цифровом — соответствующими про граммными блоками.
Как было показано в 2.9, любая сложная динамическая система п-то порядка, описываемая обыкновенным дифференциальным уравнением (2.165) или передаточной функцией (2.163), может быть представлена в виде структурной схемы, состоящей только из ин теграторов, сумматоров и блоков умножения сигналов на постоян ный коэффициент (см. рис. 2.17).
Простейшим элементом аналоговой модели является потенцио метр (рис. 2.19, а), умножающий входное напряжение их на по стоянный коэффициент а:
иу= аиХ1 |
(2.198) |
где а = г2/гх = 0 |
1. |
Основным элементом при аналоговом моделировании является операционный (решающий) усилитель — электронный усилитель
постоянного тока с большим коэффициентом усиления (k =
= 104 10е), охваченный отрицательной обратной связью. Ха рактер математических операций, выполняемых таким усилителем, зависит от операторных сопротивлений ZBX(р) и Z0. с (р), включен ных на входе и в обратную связь усилителя. Передаточная функция операционного усилителя, связывающая входное и выходное на
пряжения, |
равна |
|
| W (р) = |
Uy (р)/их (р) = - Z0. с (p)IZsx (Р), |
(2.199) |
где знак «минус» указывает, что одновременно с математическим преобразованием входного сигнала усилитель меняет его поляр ность на противоположную.
Так, если на входе и в обратную связь включены активные со противления Ri и R 2 (рис. 2.19, в), то операционный усилитель выполняет функции масштабного блока — умножает входной сиг нал их на постоянный коэффициент
kBX= — г2/гъ |
|
|
|
(2.200) |
||
который можно изменять в широких пределах. |
|
масштабный |
||||
В частном случае, |
когда |
г х = г 2 и kBX = — 1, |
||||
усилитель |
используется как |
инвертор — только |
для |
изменения |
||
полярности |
входного |
сигнала. |
|
|
|
|
Если на входе стоит резистор R ly а обратная связь образована |
||||||
емкостью |
С |
с операторным |
сопротивлением |
Z0. с (р) = 11рС |
(рис. 2.19, д), то операционный усилитель является интегратором с передаточной функцией
W (р) = г 0яс (P )/Z BX (р) = - ( 1 /рС)/гх = - 11ггСр. |
(2.201) |
||
Коэффициент |
kBX = 1/ г х а р а к т е р и з у е т |
скорость |
интегриро |
вания входного |
сигнала. |
|
|
Масштабный |
и интегрирующий усилители |
могут одновременно |
выполнять суммирование нескольких входных сигналов. При та ком совмещении основной функции усилителя и функции суммиро вания операторное сопротивление обратной связи для всех сумми руемых сигналов всегда одинаково, а входные сопротивления для отдельных сигналов могут быть различными. Поэтому коэффици енты kBX для отдельных масштабируемых или интегрируемых сла гаемых также могут быть разными.
Рассмотрим основные операционные элементы цифровых мо делей. Цифровое моделирование непрерывных систем основано на
дискретизации |
обыкновенных дифференциальных уравнений вида |
||||||
(2.165) — приближенной замене их |
разностными |
уравнениями, |
|||||
связывающими |
между собой |
дискретные |
значения |
л: (/t) = х |
(i&t) |
||
и У (О = У (*Д0 входного |
и выходного |
сигналов |
х (/) и |
у (t). |
|||
При |
этом непрерывная |
независимая |
переменная 0 <С / < 00 |
||||
заменяется дискретной переменной tt = |
0; 1А/; 2Д/; |
; Ш ; |
; |
||||
оо. |
Переход к дискретным значениям сигнала х (t) |
и дискретному |
88
а, |
в |
|
|
At |
/ У£ |
|
|
o-^Y <°~1 |
|
|
|
A t |
|
|
|
xll) |
SL |
|
|
|
|
|
|
|
е |
|
|
|
Уб П |
At |
Hi |
|
|
|
Рис. 2.20. Операционные элементы цифровых моделей
времени может быть осуществлен с помощью идеального ключа (рис. 2.20, а), который через равные интервалы времени At на ко роткое мгновение замыкается и выдает дискретные значения х,- = = х (iAt) и Ус = у (iAt) = kxc (рис. 2.20, б).
Вразностных уравнениях вместо обычных дифференциалов dу
иdx, используемых в обыкновенных дифференциальных уравнениях,
применяют разности (отсюда название уравнений!):
d У (0 « |
Ау (ti) = y ( t t) — y {tui)=yt— yi-i, |
(2.202) |
d х (t) « |
Ax (t() = x (tt) — x (tc_i) = хс— хс_ъ |
(2.203) |
а вместо производных по времени — отношение конечных прира щений (2.202) и (2.203) к шагу дискретизации по времени At. Поэ тому элементарным операциям непрерывного дифференцирования
y(t) = dx(t)/dt |
(2.204) |
и интегрирования |
входного сигнала |
dy(t) = x { t ) d t |
(2.205) |
в дискретной форме будут соответствовать разностные уравнения
Ус= (*<— |
Д£, |
(2.206) |
Ус = Ус~1~\~ XcAt. |
(2.207) |
|
На схемах цифрового моделирования уравнения дискретного |
||
дифференцирования |
(2.206) и интегрирования (2.207) изображают |
в виде полных схем вычислений (рис. 2.20, в, д) или в условном виде (рис. 2.20, г, ё).
Как видим, для цифрового моделирования операции дифферен цирования или интегрирования требуется выполнение лишь ариф метических операций суммирования, умножения на постоянный коэффициент и запоминания предыдущего значения сигнала на один шаг At, которые легко программируются. Процедура запоминания
89
или задержки дискретных значений реализуется при программиро вании вычислений с помощью оператора присваивания
УГ.=Уиъ |
(2-208) |
согласно которому в ячейке памяти происходит пошаговая замена предыдущего значения на последующее.
Формула (2.207) и схема на рис. 2.20, д соответствуют алго ритму приближенного численного интегрирования непрерывного
сигнала |
х (t), |
представленного своими дискретными |
значениями |
х (ti) = |
xi9 по |
наиболее простому методу — методу |
Эйлера (ва |
риант обратной разности или прямоугольников с упреждением). Поскольку в этом алгоритме каждое последующее значение выра жается через предыдущее, он называется рекуррентным.
Приближенное решение исходных дифференциальных уравнений с помощью их цифровых моделей дает удовлетворительные ре зультаты только в тех случаях, когда шаг интегрирования А/ достаточно мал по сравнению с постоянными времени исследуе мой системы и темпом изменения сигналов х (t) и у (t) во вре мени.
Например, для дифференциального уравнения первого порядка
с постоянной времени Т = а0/ах шаг интегрирования |
рекомен |
дуется принимать равным |
|
At < (0,1 ч-0,2) Т |
(2.209) |
Контрольные задания и вопросы
1. В чем заключается свойство однонаправленности передачи воздейст вий отдельным элементом (например, электрическим четырехполюсником)?
2.Назовите внешние и внутренние переменные, а также параметры ма тематической модели электромагнита (см. рис. 2.2, б).
3.Какие типовые воздействия используются при изучении динамики элементов и систем?
4.Поясните понятия переходного и установившегося режима на гра фике процесса,у (/), возникшего на выходе статического элемента после сту пенчатого воздействия на его вход.
5.Как будет выглядеть график амплитудного спектра отдельного гар
монического |
сигнала с x im = 10 В и сох = 2 |
рад/с? |
6. Как |
изменяется график спектральной |
плотности | X (/со) | любого |
конечного импульса при увеличении его длительности? |
||
7. Проинтегрируйте квадрат экспоненциального импульса (2.37) по |
времени и его спектральную плотность мощности (2.41) по частоте и убеди
тесь |
в |
этом частном |
случае в справедливости равенства Парсеваля |
(2.36). |
|||
п = |
8. Постройте линеаризованную статическую характеристику |
двигателя |
|||||
f (и), соответствующую |
k = |
0,05 (об/с)/В. Укажите на ней |
точку но |
||||
минального режима — и0 = |
200 |
В, п0 = 10 об/с. |
|
|
|||
|
9. Как классифицируются элементы систем управления по виду стати |
||||||
ческих |
характеристик? |
|
множительного элемента у = /j (xlt |
х 2) — |
|||
|
10. |
Линеаризуйте |
уравнение |
||||
= 5x^2 |
в точке х 10 = |
1 В, хао = |
2 В. Запишите линеаризованное уравне |
||||
ние для |
абсолютных значений [см. формулу (2.46)1 и для отклонений [см. |
формулу (2.47)].