Файл: Д’Анжело, Г. Линейные системы с переменными параметрами. Анализ и синтез.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 08.11.2024

Просмотров: 201

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Эти уравнения в матричной форме имеют вид

 

 

 

x(*) =

A ( 0 x ( * ) + B(*)u(/),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.26)

где А(/), B(t),

С (t) и D(/) соответственно

матрицы типа

п-У^п,

а х ( 0 , и (0 и

y(t)

— n-, т- и

/г X т > Г X tt и

Г

Х т '

 

 

 

г-векторы.

 

 

 

 

 

 

Системы с несколькими

входами и выходами, в данном случае с

т входами и г выходами, будем называть многомерными.

Рас­

сматривая произведение

B(t)u(t)

как один

/г-вектор

f(t),

запи­

сываем уравнения (2. 26) в виде

 

 

 

 

 

 

 

 

 

 

 

x(n---A[7)x(/) + f(0,

 

 

 

 

 

(2. 27)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y(/)=.-C(/)x(0 + D(/)u(0.

 

 

 

 

 

 

Во многих современных работах по теории

автоматического

управления вектор х(/) называется вектором

 

состояния,

а урав­

нения (2. 26) и (2. 27) уравнениями

состояния.

 

 

 

 

Часто бывает,

что система n-го порядка

с

одним

входом и

одним

выходом,

т. е. скалярная

система,

описывается

одним

обыкновенным

 

линейным

дифференциальным

уравнением

га-го порядка

 

 

 

 

 

 

 

 

 

 

 

 

 

«*(t)ywV)

+ ai(t)yln-14t)

+

- + a«V)y(t)

=

 

 

 

 

 

 

= PeWr<»>(0 + ..- +

P „ W - W ,

 

 

 

 

(2-28)

где у(t)

—скалярный выход, г(t)

—скалярный

вход, a

 

производные по времени г'-го порядка. При проведении

анализа

правую

часть

уравнения

(2. 28) можно рассматривать

как одну

возмущающую

функцию u(t).

Тогда

скалярное

дифференциаль­

ное уравнение принимает вид

 

 

 

 

 

 

 

 

 

 

^{t)yw{t)

 

+ a1{t)y^{t)

+

...-Yan{t)y{t)^u{t).

 

 

 

(2.29)

Вводя линейный дифференциальный оператор

 

 

 

 

 

/: ^

а„ (0

+

ох (/)

 

+ • • • + « »

W ,

 

 

(2.30)

 

 

 

dt

 

dt"-1

 

 

 

 

 

 

можно уравнение

(2. 29)

переписать

в более

простой

форме

L\y{t)\=u{t). (2.31)

Желательно показать, почему уравнение (2. 29) является ча­ стным случаем матричных уравнений состояния (2.27). Это можно сделать, вводя п переменных согласно соотношениям

Xl{f) = yV-V{t), /==1,2, ... , л.

(2.32)

г , . 25


Используя уравнения (2.29) и (2.32), получаем тогда сле­ дующую систему из п линейных дифференциальных уравнений первого порядка:

x1{t)

=

x2(t)

 

x2(t)

=

x3(t)

(2.33а)

*я-1 (') = •*,, (О

x A t ) = - ^ x l {

«о

где

t ) - . . . ~ ^ x

n { t n u { t )

«о

«о

0 (/) = *!(/).

 

(2.336)

Записывая уравнения (2.33) в матричной форме, получим уравнение состояния *

 

 

 

 

х(/) = А(0х(*) + Ь ( 0 » ( 0 .

 

 

 

(2.34)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

г/(0 =

с(/)х(/),

 

 

 

 

 

* При такой записи необходимо

получать функцию u{t),

представляющую

собой сумму производных действительного входного сигнала

 

 

 

 

 

 

и (О = Ро (О /-( я ) (0 + • • • + К г (О-

 

 

 

(2. ЗЗА)

С точки зрения физической

реализации это означает, что в дополнение к п ин­

теграторам, необходимым для моделирования

системы n-го порядка,

требуется

п дифференциаторов для получения u(t).

Кроме

того, при данной

форме за­

писи

решение

задачи

управления

сопряжено

с

неудобством

рассмотрения

управляющей

функции

u(t)

как функции,

связанной с действительным вход­

ным

сигналом

r(t) соотношением

(2. ЗЗА). Этих трудностей

можно

избежать,

рассматривая другую,

более удобную, форму

уравнений

состояния

 

 

 

 

 

 

x(t)

=

 

k(t)x(t)+br(t)

 

 

 

(2. ЗЗБ)

 

 

 

 

у (О = с (0 х (О + d (0 г

(f).

 

 

 

 

 

 

 

 

 

 

 

Можно показать (см. задачу 2. 16), что при такой форме

записи

А(^) и с(^)

определяются выражениями

(2. 35), в то время как

 

 

 

 

 

 

 

 

 

 

 

 

 

• Ьх

(t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h(t)

 

 

 

 

 

 

 

 

 

 

 

 

ь'(/)

= |

 

 

 

 

 

 

(2. 33В)

 

 

 

 

 

 

 

 

bn (t)

_

 

 

 

 

 

Здесь

bQ (0 ==

«о (О '

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i - \

 

t-j

 

 

 

 

 

 

 

 

*/ (О =

 

h -

2

£

 

« г

ч

-

 

 

 

 

(2. ЗЗГ)

 

«о

 

 

;-0

ft-0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.33Д)

26


где

О

1

О

О

О

1

A(t)-

О

О

О

ап

Q«-i(0

<*л-2(0

«о

«о

«о

 

 

о

 

 

о

4 0 =

с ( 0 = [ 1 0 . . . 0 0 ] ,

* 2 w

х ( ф

о

о

(2.35а)

1 «1 «о

(2.356)

(2.35в)

(2.36)

Уравнение

состояния,

содержащее

матрицу

A(t)

 

вида

(2.35а), называется каноническим уравнением в фазовых

 

коор­

динатах.

Уравнение (2.28)

эквивалентно

уравнениям

(2. ЗЗБ)

и (2.34), которые, в свою очередь, являются частным

случаем

уравнений

состояния (2.27).

Поэтому любой результат,

дока­

занный

для уравнений состояния (2.27),

справедлив

для

урав­

нений

(2. ЗЗБ)

и (2.34), а следовательно,

и для их

скалярного

эквивалента (2.28).

2. 2. 2. Существование и единственность решения однородной системы [9]. Для неоднородной системы высокого порядка, опи­ сываемой уравнением

(2.37)

27


можно доказать теорему о существовании и единственности, аналогичную теореме 2. 2 для системы первого порядка *. Вы­ кладки значительно упрощаются, если эту теорему доказать сна­ чала для однородной системы

x(t) =

A(t)x(t),

х(/0 ) = х0 .

(2.38)

Поэтому рассмотрение общей теоремы о существовании и единственности для неоднородной системы начинается с доказа­ тельства этой теоремы и изучения свойств решения в случае од­

нородной системы.

 

 

 

 

 

 

 

 

Теорема 2.3. Пусть А(^) матрица типа

пХп

с элементами

aij{i),

i, / = 1 , 2, ... , п в

виде

непрерывных

на интервале [а, Ь]

функций. Тогда на интервале [а, Ь] существует

единственный

непрерывно дифференцируемый д-вектор

х(^),

удовлетворяю­

щий

уравнению

(2.38), если

только t0

принадлежит

интерва­

лу [а, Ь].

 

 

 

 

 

 

 

 

 

Доказательство. Предполагая, что вектор x(t)

удовлетворяет

условиям теоремы 2. 3, перепишем

уравнение (2. 38) в виде

 

 

 

 

 

t

 

 

 

 

 

 

 

 

х [t) =

хо +

J" A (t) х (t) dr.

 

 

(2. 39)

 

 

 

 

to

 

 

 

 

 

 

* Было бы очень желательно результаты, полученные для скалярного

диф­

ференциального уравнения первого порядка

 

 

 

 

 

 

 

х (0 =

A (t)x

(t),

Л(0) = ЛГ0

 

 

(2.38А)

распространить на

матричное

дифференциальное

уравнение

первого

порядка

(2.38).

Поскольку

теорема 2.2

устанавливает, что решение

скалярного

урав­

нения первого порядка (2. 38А) определяется формулой

 

 

 

 

 

 

х (t)

= ехр

Л (7])йГт)

 

 

 

(2. 38Б)

можно бы надеяться, что решение матричного уравнения первого порядка

(2.Э8) определяется аналогичной формулой

 

х (с) — ехр

А (т))йГт) х0-

(2. 38В)

Ксожалению, это справедливо лишь для очень ограниченного класса систем,,

аименно [8]:

а) А(^)—постоянная матрица; б) более общее условие

 

А (О = 2

Asa* (О.

(2. 38Г)

где а1(1)Фй](Ц

при 1Ф] и А ; А, - =А 3 А{ при всех i и /;

 

в) эквивалентное условие B(t)B(t)

B(t)B(t),

(2. 38Д)

где B(t)=\

A(t)dt.

 

 

 

28


Поскольку

A(t)

и x(t) непрерывны, j

A(x)x(x)dx—диф-

ференцируемая

функция и, следовательно,

вектор x(t),

опреде­

ляемый выражением

(2.39), дифференцируем. Полагая

в фор­

муле (2. 39), что t—t0, видим, что, как это и требуется, х(/0 ) = х 0 . Таким образом, вектор x(t), определяемый выражением (2.39), может быть решением уравнения (2.38), удовлетворяющим на­

чальным условиям.

 

 

 

 

 

 

 

Для удобства

записи

введем следующие обозначения*:

| А ( / ) | = т а х У

\aiJ{t)\,

 

tt[a,b],

 

(2.40)

 

 

1

U

 

 

 

 

 

 

||А(/)|| =

Л и . 6 . | А (0| ,

 

(2.41)

 

 

 

t 6[а, Ь]

 

 

 

 

 

|х (ОИ max | JC, (0|,

t$\a, Ь],

 

(2.42)

 

 

 

i

 

 

 

 

 

/. и. Ь.\

х

а >

0.

 

 

te[a,b]

 

 

 

 

(2.43)

Заметим, что

|А(^)|

и | х ( / ) | с к а л я р н ы е

функции

t, в то

время как ||А(£)||

и \\x(t)

[| — скалярные

константы. Из

выраже­

ний (2. 40) и (2. 42) можно видеть, что

 

 

 

 

| А ( / ) х ( 0 1 < | А ( 0 |

|х(*)|,

 

(2.44)

а из (2.41) и (2.43) следует, что

 

 

 

 

 

| | А ( / ) х (0 | | < | | А ( / ) | |

||х(01|.

 

(2.45)

Уравнение (2. 42) обосновывает интегральное неравенство

 

I х 'т)

dx

< j

\x{x)\dx.

 

(2.46)

 

 

 

 

to

 

 

 

 

Определим интегральный оператор Г следующим образом:

 

Гх(/) == |

 

A{x)x(x)dx.

(2.47)

Пусть U — некоторая

точка

в интервале [а, Ь], подчиняющая­

ся условию ti>to.

Тогда

 

 

 

 

| Г х & )

( А(т)хОx)dx

< j |А (t) х (t)| dx

<

 

 

 

 

<0

 

 

<

f |A(t)|

|x(t)|rfr.

(2. 48)

* l. и. b. — наименьшая

верхняя

грань (прим. редактора).

 

29