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

Teoriya avtomatichnogo keruvannya

.pdf
Скачиваний:
75
Добавлен:
02.03.2016
Размер:
9.97 Mб
Скачать

Г л а в а 3

МАТЕМАТИЧНЕ ОПИСАННЯ ЛІНІЙНИХ НЕПЕРЕРВНИХ СИСТЕМ

 

АВТОМАТИЧНОГО КЕРУВАННЯ

 

 

1

-(<> +1) 52 + 5

+1)

-1)

У 1974 р.

Г. Розенброк заклав основи

методу автоматизованого

проектування САК, який з часом було вдосконалено і розвинуто іншими вченими. Цей метод дістав назву методу змінних стану.

Метод змінних стану дає можливість визначити залежність вихідних (керованих) параметрів багатовимірної системи х,, х 2 , х г (в узагальненому вигляді х) від деяких змінних стану у{,...9уп (в узагальненому вигляді у) з урахуванням дії керуючих впливів щ ,..., ит (в узагальненому вигляді и).

Отже, в узагальненому спрощеному вигляді рівняння стану можна записати так:

У = Ау+ В и ;

( 3 . 1 0 1 )

х = Су,

 

де А, В, С — відповідно матриці стану, керування, виходу (вихідних координат).

Розглянемо принципові загальні положення даного напряму теорії автоматичного керування.

Для описання динамічних процесів у системі автоматичного керування необхідно в загальному випадку встановити характер зміни в часовому просторі керованої величини х системи з параметрами Ь{,...,Ьп від збурень / , , . . . , / „ (до складу збурень належать

задавальні і керуючі функції и(і)).

Поширення дістав спосіб описання систем, що базується на диференціальних рівняннях і передаточних функціях, які застосовують

водновимірних системах.

Убагатовимірній системі з п змінними хх,..., хп при дії в загальному випадку збурень / , , . . . , /п для описання системи необхідно мати систему п рівнянь вигляду

СІи1

+СІ12(р)Х2

+ ...

+ (/1п(р)хп=/{\

сі21(р)хї

+сі22(р)х2

+...

+

сі2п(р)хп =/2;

СІ

" п2 (р)х2

+ ...

+ СІпп

де с!ік (/;) — поліном від оператора р = — .

160

3.14.Багатовимірні системи та метод змінних стану

Уразі застосування перетворення Лапласа при ненульових початкових умовах ці поліноми є функцією деякої комплексної величини

5 і мали б записуватись у вигляді с1ік(з). При нульових початкових умовах, як зазначалося раніше, поліноми можна записувати у вигля-

ді СІік (р).

Перехід від рівнянь ланок конкретної системи і визначення системи рівнянь та їх коефіцієнтів вигляду сІік(р) для описання конкретної САК досить детально викладено раніше на конкретному прикладі застосування теореми Крамера.

У разі багатовимірних систем їх описання і дослідження значно ускладнюються. Тому для дослідження таких систем розроблено так званий метод змінних стану, який може застосовуватись також і для

дослідження

одновимірних систем.

 

 

 

Метод базується на тому, що для розв'я-

щ(г)

зання диференціального рівняння п-го по-

и2(ї)

рядку, яке характеризує динамічні процеси в

 

деякій

системі

автоматичного керування

 

(рис. 3.40), необхідно знати п початкових зна-

Рис. 3.40

чень регульованої величини х — її початково-

го

значення

х(і0) = х(0) та п- І похідних при

 

і =

і0,

а саме:

х(і0)

= х(0); х(і0) = х(0);...;

1

(/0) = хп~1(0).

Ці початкові умови визначать подальший рух системи.

Для характеристики внутрішнього стану системи вводять деякі абстрактні характеристики системи у{ (і), у2(і),..., уп(і)з початковими умовами ух (0), ^(О),..., >>„(0), за допомогою яких з урахуванням вхідної дії и(/) можна однозначно визначити динамічні процеси в системі.

Змінні ух (і),..., уп(і) називають змінними стану системи.

У загальному вигляді деяка змінна стану ук (і) записується так:

УА0 = ЛІУі (0)... у„(0); щ (і)... ит

Стан системи в довільний момент часу і визначають за допомогою вектора стану у(/), що є функцією змінних стануук (і)(к = 1 ,...,п):

У ( ї ) = ІУЛО,У2«),.-,УМт.

(3.102)

Множину всіх значень вектора стану у (і) на деякому інтервалі

зміни часу і називають простором стану або фазовим простором. Й о -

го можна розглядати як деякий ^-вимірний абстрактний простір, що відповідає порядку диференціального рівняння системи.

6 Т е о р і я а в т о м а т и ч н о г о к е р у в а н н я

161

Глава З

МАТЕМАТИЧНЕ ОПИСАННЯ ЛІНІЙНИХ НЕПЕРЕРВНИХ СИСТЕМ

 

АВТОМАТИЧНОГО КЕРУВАННЯ

Кінець вектора стану при русі в просторі стану має множину по-

л о ж е н ь , я к а називається траєкторією вектора стану або фазовою

траєкторією. Кожна фіксована точка цього положення називається

зображуючою точкою.

Фазовий простір, який відповідає рівнянню динаміки першого порядку, називають одновимірним. Його можна зобразити як пряму на деякій площині, а фазову траєкторію — як сукупність точок (відрізок) на цій прямій.

При п = 2 вектор стану є функцією двох змінних ух (і) і у2). Фазовий простір при цьому називають двовимірним. Він є площиною з координатами ух - у2, а фазова траєкторія — крива на цій площині.

Тривимірний простір відповідає системам при п= 3.

Фазові траєкторії і простори стану одно-, дво- і тривимірних систем зображено на рис. 3.41, а—в.

Рис. 3.41

При п - 4 матимемо чотиривимірний простір і т. д.

При застосуванні методу змінних стану виділяють поняття просто-

ру керувань, простору виходів і простору збурень. При цьому під прос-

тором керувань розуміють множину можливих значень их (і), и2(і),...

...,ит(і), під простором

виходів

множину хх (/), х2 (і),..., хг(і), а

під

простором збурень

множину /, (і), / 2 ( і ) , / { ( і ) . Ці

множи-

ни

є координатами

векторів:

керувань

в

^-вимірному

просторі

и(0= [Щ (і), и2 (/),..., ит(і)У\

виходів

в

 

/-вимірному

просторі

х(0 = х (/), х2(/),

(/)]1

і

збурень

в

/-вимірному

просторі

т=[мо,мо,...,/,(()]'.

162

3.14. Багатовимірні системи та метод змінних стану

Фізична суть змінних стану (які, як було вказано раніше, в загальному випадку можуть бути абстрактними величинами) залежить ні і вибору базису, при зміні якого можна змінити і змінні стану.

ІІрп дослідженні конкретних систем за змінні стану приймають •и'які фізичні величини, які характеризують поведінку даної системи.

Рівняння стану. Якщо вважати, що похідна сіуі (і)/сіі = г, залежить тільки від поточного стану системи, то систему автоматичного керування можна описати системою диференціальних рівнянь першого порядку в нормальній формі Коші, згідно з якими похідні змінних стану залежать від значення змінних, часу і керу-

пані. //,(/).

Згідно з рис. 3.40 можна записати:

 

 

Уі(*)=

ФіЬ>і (*),у2(*),...,уп(*У, Щ (0,

М 0 , . . . ,

итЦУ,*У,

 

 

 

(3.103)

уп(о =

ФЛ^і(0, у2(*)>уп(*У, щ(0,

М 0

, / ] .

Знаючи змінні стану, знаходимо кожний з вихідних сигналів хДг)

м (функцію у ДО

змінних стану і керувань (входів):

 

( 0 =

У

І

І У

І ( 0 , ц

(

0 ,

м 2 ( 0 , « „ ( 0 ; ' ] ;

* 2 ( > ) =

У 2 Ь ' і

(*),У2(*),...,Уп(*У,

Щ

( 0 ,

и 2

( ' ) , . » ,

ит{і)\і]\

 

 

 

 

 

 

 

 

 

(3.104)

хп«) = УПІУІ (0,

у2(*)>—> У Ж Щ (0,

М

0 ,

/ ] .

Якщо рівняння (3.103), (3.104) вважати лінійними (лінеаризовапимп), а систему детермінованою, то ці рівняння однозначно визначають поведінку системи. Ці рівняння називають рівняннями стану І истеми. При зазначених вище умовах змінні стану і вихідні величини не залежать від моменту появи вхідного сигналу щ(і).

І Іри цьому рівняння (3.103), (3.104) можна записати у вигляді

>1 (0 =

аИу1 (і) +

а12у2 (/) + ... + ау„(ґ) +

+

Ьп

щ (і) +

Ьї2и2(і) +... + ЬХтит(і)\

У2 (0

=

а21уї (і)

+

а22у2 (/) + ...+ а2пуп(і) +

163

Глава З

МАТЕМАТИЧНЕ ОПИСАННЯ ЛІНІЙНИХ НЕПЕРЕРВНИХ СИСТЕМ

 

 

АВТОМАТИЧНОГО КЕРУВАННЯ

 

 

 

 

+ Ь2[ щ (ґ) + Ь22 и2(і) +

... + Ьит(ґ);

 

 

 

 

 

 

(3.105)

 

У,М) =

Ух (0 + ап2у2 (і) + ... +

аппуп(і) +

 

 

+ Ьп{ ил (0 + Ьп2и2(ї) +... + Ьтп

 

 

 

(0 =

с,, у, (/) + спу2 (()

+... +

с, „у„(() +

 

 

+ сіпи1 (() + <іпи2(1) +

... + с1ііт(();

 

 

*2 (0 = с21 У і (0 + с22^2 (/) + ...+ с2„у„(0 +

 

 

+

щ (0 + Л22и2(1) + ... + с1ит(();

(3.106)

хг(() = сГІу1 (() + сг2у2(0 + ... + сгпу„(() +

+(0 + </г2 «2 (0 + • • • + <*тин(і),

де аікікік,(1ік сталі коефіцієнти, знайдені при визначенні і перетворенні відповідних функцій ф, , у, в рівняннях (3.103), (3.104).

Для складних об'єктів запис рівнянь вигляду (3.105), (3.106) та їх перетворення при дослідженнях системи утруднені, тому значне поширення дістала раціональніша компактна векторно-матрична форма рівнянь стану. Згідно з цією формою рівняння стану можна записати у вигляді

У = Ау(П+ В и(/);

х = Су(ґ)+ £и(ґ),

(3.107)

де

(3.108)

є матриця стану розмірністю пхп, яка визначає вільні і вимушені рухи системи;

В =

(3.109)

Ь„,Ь „-,...Ь,

164

3 . 1 4 . Багатовимірні системи та метод змінних стану

\ит\рицн керування (матриця входу) розмірності пхт, яка визначає

п»а< ми іанежпість входу системи і змінних стану;

 

С\ і Сп.. .Сі

 

 

"1112*

'\п

 

С

С 2 \ С 2 2 " ' С 2 п

(ЗЛЮ)

 

 

сгХ\сх г2"

- і ш три ця вихідних координат (матриця виходу) розмірності г х п, яка ми шачає характер взаємозв'язку вихідних величин зі змінними ста-

ну,

~сіпсі{2...сі{п'

 

СІ1{СІ22...СІ2п

(3.111)

 

сіг{сіг2...сіп

 

> маїриця розмірністю г х /7, що характеризує прямии зв язок вихіднії координат х/ (і) з керуванням; визначає безпосередній вплив керн мані ія на різні складові вихідних координат; для багатьох систем

/>

0

< і рук гурпу схему багатовимірної системи, що відповідає наведе-

ним

рівнянням стану, подано на рис. 3.42. Подвійні лінії відобража-

І-НІ.

маїрично-векторний характер математичної моделі системи.

Рис. 3.42

Знаходження рівнянь стану і аналогових моделей сис-

тем Розглянемо одновимірну систему, динаміка якої описана дифе- і" шпачинім рівнянням А/-ГО порядку, зведеним до нормованого ви-

165

Глава З

МАТЕМАТИЧНЕ ОПИСАННЯ ЛІНІЙНИХ НЕПЕРЕРВНИХ СИСТЕМ

 

АВТОМАТИЧНОГО КЕРУВАННЯ

гляду, при якому коефіцієнт за найвищого порядку похідної від регульованої величини х дорівнює одиниці (о0 = 1); и(()— вхідна величина. Отже, вихідне рівняння має вигляд

сГ'х(і)

+

а

сІ"~іх(і)

+... + ап-і

4хЦ)

+ а„х(0 = Ь0и(і).

( 3 . 1 1 2 )

Сіґ'

сії"

 

сії

 

 

 

 

 

 

 

Якщо початкові умови відомі, то, як вказувалося раніше, змінні

стану можна записати так:

 

 

 

 

 

 

 

 

 

У\

(Г)=х((),

 

 

 

 

 

 

у1(()

=

х(() = у1

(І),

 

 

 

 

 

У3(<)=Х(()=У2(1),

 

 

(3.113)

 

 

 

уп(ґ)

 

=

х("-1}(Г)

= у„^(ґ).

 

З рівняння (3.112) знайдемо значення похідної найвищого по-

рядку:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

^сГ^хіТ) + Ьои(г)_

( 3 ] ] 4 )

сіґ

 

 

 

 

 

сії

 

сії

 

На основі рівнянь (3.113) і (3.1 14) можна дістати систему диференціальних рівнянь першого порядку:

у, (0 = .М0;

 

Ш = У3(0;

 

 

(3.115)

УпіО = ~апУ\ (0 - а„-іУг(0

а\уЛО + Ь0и(*У,

х(0 = у1 (0-

 

Звідси знаходимо систему рівнянь стану (3.115) в матричній формі:

 

 

0

1

 

0

0

.

0

"

 

 

"0 ~

~и(Г)

ш

 

0

0

 

1

0

.

0

 

у2(п

т1

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0

 

0

0

...

1

 

У..-Л0

 

0

0

уп(0

 

га"

 

1

2

 

 

 

 

У,«)

 

А .

0

166

3.14.

Багатовимірні

системи

та

 

 

 

 

 

метод змінних

стану

 

 

 

 

 

 

х(г) =[ 1

0

0... 0]

[у,

(і)у2

(ґ)...

^(ґ)]т.

У пек і орній формі рівняння стану матимуть вигляд

 

 

у(і)

=

Ау(і)

+

Ви((У,

 

 

 

 

х(/) = Су(0,

 

 

.і иідновідні матриці

 

 

 

 

 

 

 

 

0

1

 

0

 

 

0

0

 

0

0

 

1

 

о

0

 

0

 

 

 

 

 

 

1

 

-А.. -а..

 

-а„

 

я-З

-а.

Матриця керування В розмірністю п х 1 має вигляд ' 0 0

В

0

К

( 3 . 1 1 6 )

(3.117)

(3.118)

( 3 . 1 1 9 )

лматриця вихідних координат С розмірністю їх п

С= [1 0 0... 0];

( 3 . 1 2 0 )

( 3 . 1 2 1 )

Схему аналогового моделювання одновимірної системи автомаиі'шого керування, що відповідає рівнянням (3.115), наведено на риє. 3.43. У модель введено початкові умови змінних стану »•„(()),..., ух (0), а також значення коефіцієнтів ах,...,ап. Змінюючи ичідну величину и(і) за допомогою даної аналогової моделі, можна аістати залежність х(і).

Загальна методика розв'язання рівнянь стану. Однорічні рівняння. Однорідне рівняння стану (яке не враховує дію вхідних величин и(/)) є векторно-матричним рівнянням вигляду

у(і) = Ау(і).

(3.122)

1 6 7

Глава З

МАТЕМАТИЧНЕ ОПИСАННЯ ЛІНІЙНИХ НЕПЕРЕРВНИХ СИСТЕМ

 

АВТОМАТИЧНОГО КЕРУВАННЯ

 

 

 

 

 

 

(/2(0)

У№

 

 

 

 

 

 

 

Уп-М

УїМ

УМ х(0

 

 

 

 

 

 

 

 

І

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 3.43

Воно визначає вільні коливання системи. Позначимо початкові умови: у, (0), >>2(0),... >>„(()). їхній вектор

Відповідне скалярне диференціальне однорідне рівняння

С І ^ = ау(П або

^ = асі!,

сіґ

у(Т)

Зінтегрувавши обидві частини останнього виразу, дістанемо

,<,»

у^

і

або

 

 

1п Ш

= а і ,

звідки

УФ)

 

Розклавши експоненціальну(Т)=е"'функці( 0)ю. е"' в степеневий ряд, одержимо

у(і) = е"'у( 0 ) =

X"- ^ 1- ( а'г )А*,к у ( 0 ) .

( 3 . 1 2 3 )

*

= о/(

 

Для матричної функції е згідно з формулами матричного числення маємо

еА' = Г + А1 + (АГ)2

, (АГ?

 

21

З!

к = о К

168

3.14. Багатовимірні системи та метод змінних стану

і- І, /

відповідно динамічна й одинична матриці розмірності п х п.

Згідно зі скалярним рівнянням ( 3 . 1 2 3 )

 

 

у(ґ)

= еА'у(0) =

^±Ак(к'у(0).

(3.124)

 

 

 

к = 0 & '

 

^диференціювавши обидві частини цього рівняння, дістанемо

 

у ( 0 =

~ !)•

І~(АПку(0),

(3.125)

 

к = 1

к = 1 !

 

ме к

1 = к{.

 

 

 

1 Іідставляючи в ліву частину виразу ( 3 . 1 2 5 ) відповідне значення з ( І. І 23), а в праву — з ( 3 . 1 2 4 ) , бачимо, що вони тотожні, а це дає змогу сіверджувати, що ( 3 . 1 2 4 ) — загальний розв'язок векторно-матрично- и> рівняння.

( ф у н к ц і ю еАг

н а з и в а ю т ь перехідною функцією стану а б о

( ч а с т о )

(фундаментальною

матрицею.

 

Отже, розв'язок однорідного векторно-матричного рівняння

( 4 22) матиме вигляд

 

 

у(Г) = еА'у(0).

 

Для лінійних стаціонарних систем розв'язок однорідного рівнян-

ії і через перехідну матрицю стану можна записати у вигляді

 

 

у(() = еА'у(0).

(3.126)

Основні властивості перехідної матриці стану наведено без відпоиі.'цшх доказів.

1. Вимірність перехідної матриці стану дорівнює вимірності динамічної квадратної матриці А. Тому вона також є квадратною матрице ю розмірності п х п, а її складові — функціями часу.

2. При і = 0

=1.

Перехідна матриця стану еАг є неособливою матрицею, а це можливо при умові, що сіег еАі Ф 0.

4. Для перехідної матриці стану еАіеВї = е(А+ В)г лише тоді, коли

І// ПА.

Я. \ е А і ]"'

=е-Аі.

(). еА' =еАГ А'А.

169

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]