Файл: Медич Д. Статистически оптимальные линейные оценки и управление.pdf

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

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

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

Добавлен: 11.04.2024

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

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

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

Пример 2-4. Применим изложенную методику для получения ди­ скретной модели стационарной системы второго порядка из приме­ ра 2-1. Для простоты предположим, что дискретные моменты отсчета отстоят друг от друга на постоянные промежутки времени Т. Тогда k+l соответствует моменту времени (£+1)7", a k— моменту kT.

Так как

О 1

О

то Fn=0 для всех п^2. Следовательно, переходная матрица со­ стояния

ф(*,

т) = ф(г т) = e f ( ' - T ) =

/ + f

1

— *)

=

 

1 t —

 

О

1

 

 

 

 

 

 

 

 

 

 

 

Поэтому

 

 

 

 

 

 

1

Т

 

 

 

 

 

 

 

ф - (Л+

1.

k)

=

 

 

 

 

 

 

 

О

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда, согласно уравнению (2-33),

 

 

 

 

 

 

 

 

 

 

 

 

1

+

1)7"—т

 

0

dz

=

 

 

 

 

 

 

 

1

 

 

1

 

 

 

kT

 

 

 

 

 

 

 

 

(k + l)T

 

 

 

 

 

 

 

 

 

 

{k+l)T-z

 

 

 

1/2 7"г

 

 

 

 

 

 

 

 

 

 

 

kT

 

 

1

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Аналогично

 

(k+\)T

 

 

 

 

 

 

 

 

 

 

 

 

 

1 (k+l)

T z

 

0

 

 

 

 

 

 

 

 

dz

=

 

 

 

 

 

0

 

1

 

1

 

 

 

kT

 

 

 

ь

 

 

 

 

 

 

 

1/2 ЬТг

 

 

 

 

=

bT (k +

1,

k)

=

 

 

 

 

 

bT

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Следовательно,

динамика

системы описывается

уравнением

x(k+

1

Т

 

 

 

1/2 Т2

 

 

 

1/267™2

1):

 

 

 

 

Т

 

» ( * )

+

 

67-

a(Ä).

 

011

 

 

 

 

 

 

 

 

Соответствующее

уравнение

модели* измерения

 

 

 

 

Z(k +

 

1

0

 

x(k+l)+v(k+l).

 

 

 

 

 

1):

1

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

Непрерывная

 

модель

как предельный

 

случай

дискретной

 

модели

 

 

 

 

 

 

 

 

 

При получении некоторых результатов теории оцен­ ки и управления в непрерывных линейных системах ока зывается удобным дискретизовать уравнения системы,

53


получить алгоритмы оценки и управления для дискрет­ ных систем и затем рассматривать предельное поведе­ ние этих алгоритмов при стремлении к нулю интервала между дискретными моментами времени.

Методика получения дискретной модели из непрерыв­ ной приведена выше. Теперь составим методику обрат­ ного перехода. Для этого обозначим дискретные момен­

ты времени k и k+l

через t

и t+At

соответственно,

где

Д / > 0 .

 

 

 

 

 

t)

 

 

 

 

 

 

 

 

 

Затем,

разлагая

Ф(г + Д£,

в

ряд

Тейлора

по

степе­

ням At, получаем:

 

 

 

 

 

 

 

 

 

 

 

 

 

ФЦ+М,

t) =<b(t,

t) + Ô(t,

t)At

+

 

0(At2),

 

 

где О (At2)

—матрица

размера

пХп,

все

элементы кото­

рой имеют порядок

(At2)

и выше.

 

 

 

 

 

 

 

 

Поскольку Ф(і,

і ) = /

и <b(t, t)=F(t)<D(t,

t),

имеем:

 

<î>(t+At, t)=I

+ F(t)At

+ 0(M2).

 

 

 

(2-34)

Из четвертого равенства в (2-33) ясно, что

 

 

 

 

 

t+At

 

 

 

 

 

 

 

 

 

 

 

 

 

Г(* + Д*. t)=

 

j*

Ф(* +

Д/, z)G(T)dz

=

G(t)At

+

0(Af).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2-35)

Аналогично

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t+м

 

 

 

 

 

 

 

 

 

 

 

 

 

W(t-\-At,

f)=

J Ф(* +

ДЛ х)С(т)Л =

С(0Д< +

О(ДО-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2-36)

Заменяя в уравнении (2-30) k

и k + \ на

t и t+At

со­

ответственно, имеем:

 

 

 

 

 

 

 

 

 

 

 

 

 

*(/+Д/)=Ф(*+1Д*,

t)x(t)

+

 

 

 

 

 

+ T(t+At,

t)w(t)+W(t

 

+ At,

 

t)u(t).

 

 

 

Подставляя

в это

соотношение

уравнения

(2-34) —

(2-36) и группируя члены, получаем:

 

 

 

 

 

 

 

 

x(t+At)

={I+T(t)At

 

+ 0(At2)]x(t)

 

+

 

 

 

 

+{G

(t) At + O (At2)]w

(t) +[C(t)At

 

+

 

 

 

 

+0(At2)]u(t)=x(t)+F(t)x(t)At+

 

 

 

 

 

 

 

 

+ G(t)w(t)At

+ C(t)u(t)At

 

+ 0

(At2).

 

 

 

54


Перенося x(t) в левую часть уравнения, Деля обе ча­

сти на M и переходя к пределу

при At—>-0, для t^to по­

лучаем уравнение

 

 

 

x =

F(t)x+G(t)w{t)+C(t)u(t),

совпадающее с уравнением (2-1). Уравнение

z(t + At) = H(t + At)x(t-{-At)

+v{t+At)

при At—>-0 принимает

вид:

 

 

z(t)

=H(t)x(t)

+v(t).

 

2-4. Н А Б Л Ю Д А Е М О С Т Ь И УПРАВЛЯЕМОСТЬ

Теперь перейдем

к рассмотрению

двух фундамен­

тальных понятий теории линейных систем, тесно связан­ ных с основными идеями оценки и управления. Эти по­

нятия,

называемые

наблюдаемостью

и

управляемостью,

введены Калманом ,[Л. 2-7, 2-8], работы

которого яви­

лись

основой для

других

исследований

(см. также

[Л. 2-9—2-15]).

 

 

 

 

Формулировка и изучение

понятий

наблюдаемости и

управляемости основаны, соответственно, на двух сле­ дующих вопросах, возникающих из очевидных физиче­ ских соображений:

1. При каких условиях можно восстановить поведе­ ние вектора состояния х динамической системы на ко­ нечном интервале времени, зная поведение вектора изме­ рения z на том же интервале?

2. При каких условиях можно перевести динамиче­ скую систему из заданного начального в требуемое ко­ нечное состояние за ограниченное время с использова­ нием кусочно-непрерывного управления и?

Перед тем как перейти к подробному обсуждению поставленных вопросов для двух классов линейных си­ стем, введенных в § 2-2 и 2-3, с помощью нескольких простых схем качественно проиллюстрируем понятия на­ блюдаемости и управляемости на системе общего вида.

Рассмотрим динамическую систему S с вектором со­ стояния X, вектором управления и и вектором измере­ ния z. Предположим, что возмущения и ошибки измере­ ния отсутствуют и система может быть как дискретной, так и непрерывной.

Вначале допустим, что структурная схема системы имеет вид рис. 2-9, где у — вектор, компонентами кото-

55


рого являются первые k компонент вектора х (хі,.. ,,Хь) или их часть. Из структуры системы ясно, что значения Xh+t, ..., xk нельзя определить, исследуя z, поскольку эти переменные не влияют на хи .. ., Хь. и не включены в z. Следовательно, система является ненаблюдаемой. В то

Рис. 2-9. Схема ненаблюдаемой

Рис. 2-10. Схема наблюдаемой

управляемой системы.

неуправляемой системы.

же время, если и воздействует

на

все элементы

х, систе­

ма является управляемой.

 

 

 

Аналогичным образом система

на рис. 2-10

будет на­

блюдаемой, но неуправляемой, поскольку и воздействует только на переменные хи ..., xk.

Ясно, что системы можно разделить на четыре сле­ дующие категории: наблюдаемые управляемые, наблю­ даемые неуправляемые, ненаблюдаемые управляемые и ненаблюдаемые неуправляемые.

2-5. Н А Б Л Ю Д А Е М О С Т Ь В НЕПРЕРЫВНЫХ И ДИСКРЕТНЫХ ЛИНЕЙНЫХ СИСТЕМАХ

Исследуем более подробно понятие наблюдаемости для двух классов линейных систем, введенных выше (см. [Л. 2-7—2-15]). Рассмотрим этот вопрос в предположении, что возмущения системы и ошибки измерения равны нулю, т. е. при идеальных условиях (в этом же предпо­ ложении в § 2-6 будет рассматриваться понятие управ­ ляемости) .

Непрерывные линейные системы

Вначале рассмотрим непрерывную линейную си­

стему

 

£ = F(t)x+C(t)u(t);

(2-37)

z(t)=H(t)x(t)

(2-38)

56

для i^U, где все члены определены в § 2-2. Предполо­ жим, что «(/) является известной функцией времени для

всех t^t 0 ,

но x(to) неизвестно.

Пусть

требуется определить x(t), исследуя z(t) на

некотором конечном интервале времени \t0, t^\. Очевидно,

если

матрица

H (t)

имеет

размер

пХп

и

несингулярна

для всех t^to,

то

x(t)=H-4t)z(t),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и вопрос о наблюдаемости решается

тривиально.

 

Этот вопрос

также

легко

решается,

если матрица

H(t)

имеет

размер

пХп

и

несингулярна

только

для

одного значения

t^40,

 

скажем, для іл.

Чтобы убедиться

в этом, вначале

заметим, что

 

 

 

 

 

 

 

 

 

 

x(ta)

=

H-^ta)z(ta).

 

 

 

Согласно уравнению

(2-19)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

X (0 =

Ф [t,

ta) X (t0)

+ J Ф (t,

x) С (т) и (х) dz.

(2-39)

Для

t =

ta

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

('.)«('„)

=

<">('„. Qx{t0)+

 

\Ф{іа,

x)C(x)a(x)rfx,

откуда

следует,

что

 

 

 

 

 

 

 

 

= Ф & . ОН-Чіа)гѴа)-

 

| Ф ( / 0 ,

-г)С(т)и (x)Л. (2-40)

 

 

 

 

 

и

 

 

 

 

 

 

Подстановка

значения

x(t0)

из

 

уравнения (2-40)

в уравнение (2-39) позволяет

определить x(t)

для

всех

t^t 0 ,

поскольку

и (г) известно

при

x^U.

сингулярна

Однако если

матрица H(t)

размера

пХп

для

всех t^tQ

ИЛИ если

H (4)

имеет

размер

тХп

при

тф'п,

то не совсем

ясно, как можно

определить х(4) из

z(x),

U^Lx^ti

при

некотором

конечном

Это и

есть

вопрос, который

требуется

здесь

исследовать.

 

 

 

 

 

 

 

 

 

 

 

 

 

$7


 

Заметим,

что в силу

уравнения

(2-39)

состояние

x(t)

известно для всех

t^t 0 ,

если

начальное

условие

 

x(t0)

можно определить из измерений. Имея

это в виду,

вве­

дем

определение

наблюдаемости

следующим

образом.

 

Непрерывная

линейная

система

(2-37), (2-38)

 

назы­

вается

наблюдаемой,

если x(to)

можно

определить,

 

зная

z(t),

h^.t^ti,

для

некоторого

конечного

ti. Если

это

справедливо

для

любого

to, то система

называется

 

пол­

ностью

наблюдаемой.

 

 

 

 

 

 

 

 

 

Необходимое

и достаточное

условие

полной наблю­

даемости системы дает следующая теорема.

 

 

 

Теорема 2-1. Непрерывная

линейная

 

система

(2-37),

(2-38) является полностью наблюдаемой

 

тогда

и

только

тогда,

когда

симметрическая

матрица

размера

пХп

 

 

Мс

О =

(V (t. Q H' (t) H (t) Ф {t, g dt

(2-41)

положительно

определена

для

некоторого

конечного

ti>t0.

 

 

 

 

 

 

 

 

 

 

 

 

 

Д о к а з а т е л ь с т в о . Так как предполагается, что u(t) известно для всех t^t 0 , легко определить воздей­ ствие управления u{t) на состояние x(t):

t

| Ф ( / , т)С(т)и(і)Л .

Следовательно, достаточно рассмотреть систему вида

x=F{t)x;

(2-42)

z(t)=H(t)x(t),

(2-43)

где t^t0.

Вначале докажем

достаточность. Подставляя реше­

ние уравнения (2-42)

 

 

 

х(і)=Ф(і,

t0)x(t0)

 

в уравнение (2-43), получаем:

 

 

z(t)=H(t)0(t,to)x(to).

(2-44)

58