Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Теория автоматического управления

..pdf
Скачиваний:
3
Добавлен:
12.11.2023
Размер:
16.96 Mб
Скачать

Рис. 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 (м) жидкости в баке, а в качестве входных — расходы и 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)].