Файл: Медич Д. Статистически оптимальные линейные оценки и управление.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 11.04.2024
Просмотров: 235
Скачиваний: 1
х(0) и |
начального |
времени, |
для |
которой |
x(N)=0, |
где |
N |
|
конечно. |
Если это справедливо |
для |
любого х(0) |
и |
для |
|||
любого |
начального |
времени, |
то система называется |
пол |
||||
ностью |
управляемой. |
|
|
|
|
|
|
|
Докажем следующую теорему. |
|
|
|
|
||||
Теорема 2-4. Дискретная |
линейная |
система (2-67) |
яв |
|||||
ляется |
полностью |
управляемой |
тогда и |
только |
тогда, |
когда |
симметрическая |
матрица |
размера |
пХп |
|
|||
|
|
N |
|
|
|
|
|
|
Wd |
(О, ^ ) = |
}]Ф(О, |
О 47(і, |
і - 1)Щі, |
І-1)Ф' (О, 0 (2-68) |
|||
положительно |
определена |
для |
некоторого |
конечного |
||||
/Ѵ>0, |
где |
|
|
|
|
|
|
|
Ф(0, О = Ф ( 0 , 1)ф(1, 2) |
. . . Ф ( і — 1 , |
і); |
і = 1 , |
2,...,N. |
Д о к а з а т е л ь с т в о .
Достаточность. Используя уравнение (2-67), можно показать, что
|
|
|
|
|
|
N |
|
|
|
|
|
|
|
х(М) = Ф (N, 0) je (0) + |
£ |
Ф (N, t) 47 (/, |
/ _ |
1) и (і - 1 ) , |
(2-69) |
||||||||
где |
Ф(УѴ, i)=0(N, |
N— 1) . . . Ф ( / + 1 , |
i), |
i=N— |
1,..., 1, |
0. |
|||||||
|
Теперь |
пусть |
|
|
|
|
|
|
|
|
|
||
u(î-l) |
= |
-W |
(i, i - |
1) Ф' (0, i) W~l |
(О, N) x (0) |
(2-70) |
|||||||
для |
|
i=\,...,N. |
Подставляя уравнение (2-70) |
в уравне |
|||||||||
ние |
|
(2-69), |
получаем: |
|
|
|
|
|
|
|
|
||
|
|
|
|
|
x (N) = |
Ф (N, 0) je (0) |
- |
|
|
|
|
||
- S |
Ф;(^Ѵ, |
О T (І, І - 1) T ' |
(i, f - |
1) Ф' |
(0, î) W-x |
(О, N) |
x (0). |
||||||
|
Но Ф(ІѴ, |
і ) = Ф ( Л / , |
0)Ф(0, |
i) для |
любого |
начального |
|||||||
времени. Следовательно, |
|
N |
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x (N) = |
Ф {N, |
0) x (0) - |
Ф (УѴ, 0) X Ф (0, i) W (i, i |
- |
|
||||||
- 1) 47' (f, i - |
1) Ф' (0, i) [UT-1 (О, ЛО 'х (0)] = |
Ф (/V, 0) x (0) |
- |
||||||||||
|
|
|
- |
Ф {N, 0) r d |
(О, N) ИГ"1 (О, УѴ) je (Ü), |
|
|
|
68
откуда ясно, что x(N)=ö. Поэтому система является полностью управляемой, если матрица Wd(0, N) положи тельно определена для HeKOTqporo конечного JV>0.
Необходимость. Доказательство необходимости остав лено читателю в качестве упражнения, поскольку за ис ключением несущественных подробностей оно повто ряет доказательство для непрерывных линейных систем.
Для стационарных систем справедлив следующий ре
зультат. |
|
|
|
|
|
|
|
|
|
|
Следствие 2-4. Стационарная |
|
дискретная |
линейная |
|||||||
система x(k+1) =Флг(&)+Чг и(А), |
k—Ô, |
1, ... |
является |
|||||||
полностью |
управляемой |
тогда |
и |
только |
тогда, |
когда |
||||
матрица размера |
пХпг |
|
|
|
|
|
|
|
||
имеет ранг |
п. |
WW ФУ... |
<№-»Щ |
|
|
|
|
|||
|
|
|
|
|
|
|
|
|||
Пример 2-7. Для стационарной дискретной линейной системы |
||||||||||
второго |
порядка из примера |
2-4 |
|
|
|
|
|
|
||
|
|
Ф |
1 |
т |
[ |
1/267 |
|
|
|
|
|
|
= |
1 ; ч? |
= 1 |
ьт |
|
|
|
|
|
откуда |
ясно, что |
0 |
|
|
|
|
||||
|
1/267"« |
3/26 Л |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|||
|
|
|
|
67 |
|
67 |
|
|
|
|
Очевидно, что эта матрица имеет ранг, равный 2 |
для |
7>0 |
и |
|||||||
6=И=0. |
|
|
|
|
|
|
|
|
|
|
Вполне очевидно подобие между теоремами |
2-1 и 2-2 |
|||||||||
о наблюдаемости |
с одной стороны |
и теоремами |
2-3 |
и |
2-4 |
об управляемости с другой. Например, теорему 2-3 фор мально можно получить из теоремы 2-1, проделав в по следней следующие замены:
Наблюдаемость -* Управляемость |
|
||||
|
Ф(і, |
г 0 ) - +Ф'('„ . |
0 |
|
|
|
н |
(t) -» с (о |
|
|
|
Аналогичная ситуация наблюдается также для тео- |
|||||
ірем 2-2 и 2-4 и соответствующих следствий. |
|
||||
Это свойство |
было |
впервые |
замечено Калманом |
||
[Л. 2-7], который назвал |
его дуальностью. |
Итак, |
наблю |
||
даемость и управляемость являются дуальными |
свойст |
||||
вами линейных |
систем. |
Исследование |
многочисленных |
приложений понятия дуальности выходит за рамки на стоящей книги и в дальнейшем этот вопрос не будет
69
больше рассматриваться. Однако следует заметить, что многое из того, что было сказано ранее по поводу на блюдаемости, применимо теперь и к управляемости.
Уравнение (2-61) в доказательстве теоремы 2-3 и уравнение (2-70) в доказательстве теоремы 2-4 представ ляют собой алгоритмы управления с целью перехода от заданного начального состояния к началу координат за конечное время. Отсюда очевидна связь между понятия ми управления и управляемости. Действительно, можно показать [Л. 2-2], что управление вида (2-61) минимизи рует «работу» управления
г,
J и' (0 и (0 dt to
при переходе от x(t0) к началу координат на закреплен ном интервале времени [to, f j . Иными словами, здесь получено решение частной задачи оптимального управ ления. То же справедливо для уравнения (2-70) в ди скретном случае, где работа управления имеет вид:
N
2 u'(l- l)u(l- 1).
'=1
После изложенного введения в теорию наблюдаемо сти и управляемости читатель уже должен составить представление о характере задач оценки и управления. Дальнейшая работа будет связана с получением и при менением методов решения задач о наблюдаемости и управляемости в случае, когда присутствуют возмуще ния и ошибки измерения некоторого частного вида.
2-7. НЕЛИНЕЙНЫЕ СИСТЕМЫ
Во многих случаях, имеющих практическую важ ность, физические системы нельзя описать обыкновенны ми линейными дифференциальными уравнениями или линейными разностными уравнениями и для этой цели следует использовать системы нелинейных уравнений.
В таких задачах часто удобно линеаризовать уравне ния системы относительного некоторого заданного на бора номинальных условий и получить алгоритмы оцен ки и управления относительно этих номинальных усло вий.
70
Рассмотрим систему, которую можно описать соот ношениями
|
|
|
x = f[x, w(t), |
u(t), |
t]; |
(2-71) |
|
|
|
z(t)=h[x(t), |
u{t), |
t], |
(2-72) |
где |
t^t0, |
a |
векторы x, w, и |
и v имеют тот же |
смысл, |
|
что |
и ранее. |
Предполагается, что / — n-мерная |
вектор- |
функция указанных переменных, непрерывная и непре
рывно дифференцируемая по всем компонентам |
векто |
||||||||||
ров x, w и и; |
h — m-мерная |
вектор-функция |
указанных |
||||||||
переменных, |
непрерывная |
и непрерывно |
дифференци |
||||||||
руемая по всем компонентам векторов |
х |
и ѵ. |
|
||||||||
Для данного x(t0) |
и известных |
кусочно-непрерывных |
|||||||||
функций |
w(t)—w{t) |
и u(t)—ïï(t) |
из теории |
обыкновен |
|||||||
ных |
дифференциальных |
уравнений известно, |
что |
урав |
|||||||
нение (2-71) можно решить |
и получить |
некоторое x(t) = |
|||||||||
= x(t). Тогда |
для известной |
кусочно-непрерывной |
функ |
||||||||
ции |
v(t)=v{t) |
можно |
получить |
z(t)=z(t), |
|
подставляя |
|||||
x(t) |
wv(t) |
в уравнение |
(2-72). |
|
|
|
|
|
|||
Будем |
считать x(t), |
w(t), |
ïï(t), |
v{t) |
и z(t), |
где |
t^to, |
номинальными значениями, относительно которых тре буется линеаризовать систему (2-71), (2-72). Чтобы про
вести эту линеаризацию, |
обозначим: |
|
|
||
x(t) =x(t)+Êuc(t); |
w{t) |
=w(t) |
+Aw(t); |
||
u(t)=ïï(t)+&u(t); |
|
v(t) |
=v(t) |
+Av(t); |
|
z{t)=z(t)+&z{t) |
|
|
|
||
и затем разложим обе части |
исходных |
систем уравнений |
|||
в ряд Тейлора. |
|
|
|
|
|
Вначале рассмотрим і-е уравнение из системы урав |
|||||
нений (2-71) |
|
|
|
|
|
Xi=fi\x, |
w(t), |
u(t), |
t) |
|
|
для некоторого і=\,...,п. |
Разлагая |
его в ряд Тейлора, |
|||
получаем: |
|
|
|
|
|
|
-f- члены второго порядка, |
(2-73) |
|
где |
индекс 0 означает, что указанные |
частные |
производ |
ные |
следует вычислять при x=x{t), |
w = w(t) |
и u = ü~(t). |
71
Полагая, что члены второго и более высоких поряд
ков в уравнении (2-73) пренебрежимо малы, и |
исполь |
|||||||||||
зуя векторно-матричные обозначения, имеем: |
|
|
||||||||||
Ах,- |
дхі |
' 'дхп |
|
|
' ' ' dwp |
|
Aw (t) |
+ |
|
|||
|
|
|
|
|
|
|
|
|||||
|
|
|
+ |
|
|
|
Au (t), |
|
|
(2-74) |
||
для любого і= 1,..., п, где |
Дда, |
|
|
Д«і |
|
|
||||||
|
|
|
|
|
|
|
*- |
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Ах-. |
|
; |
Aw — |
; |
Au — |
|
|
||||
|
|
Дх„ |
|
|
àwp |
|
|
|
|
|
||
Систему уравнений (2-74) можно записать |
в более |
|||||||||||
удобном |
виде: |
|
|
|
|
|
|
|
|
|
|
|
|
Ax = F(t)Ax+G(t,)Aw(t) |
+C(t)Au(t), |
|
(2-75) |
||||||||
для t^t0j |
|
где |
F (t) |
— Wfij(t) II — матрица |
|
размера |
пХп |
|||||
{U){t) = |
(dfi/dXj)o; |
і, |
/ = 1 , . . . , я ) ; |
G (t) |
=\\gi} |
(t) || - м а т р и |
||||||
ца рамера |
nXp |
(gij(t) |
= (âfi/dwj)0; |
i = l , . . . , n ; |
/ = 1 , . . . |
|||||||
...,p); |
C(t) |
=\\Cij(t) Il—матрица |
размера |
nXr |
(сц(і) |
= |
||||||
= (dfi/dujjo; |
i= 1,..., |
«; |
/ = 1 , . . . , |
r). |
|
|
|
|
|
|||
Очевидно, уравнение (2-75) |
имеет |
ту |
же форму, |
что |
||||||||
и уравнение (2-1). |
|
|
|
|
|
|
|
|
|
|||
Аналогичным |
образом получим: |
|
|
|
|
|
||||||
|
|
Az{t)=H{t)Ax{t)+A{t)Av{t), |
|
|
|
|
|
|||||
где H (t) = \\hij(t) |
II — матрица |
размера |
m X / i |
(hij(t) |
= |
|||||||
= {дкі(дх})0; |
i = l , . . . , |
m; |
/ = ! , . . . , |
n) ; |
|
Л (г) = |
Ііа^ (^) || — |
|||||
матрица |
размеіра тХт |
|
(ciij(t) |
= (dhifdüj)0; |
i—l,...,m, |
|||||||
/ = 1 , . . . , т ) . |
|
|
|
|
|
|
|
|
|
|
|
|
Обозначая Аѵ*(t) |
=А(t)Аѵ(t), |
получим уравнение |
||||||||||
|
|
Az(t)=H(t)Ax(t)+Av*(t), |
|
|
|
(2-76) |
||||||
имеющее такую же же форму, что и уравнение |
(2-2). |
|||||||||||
Следует |
иметь в |
виду, |
что линеаризованная |
система |
уравнений применима только для достаточно малых от клонений от принятых номинальных условий.
Пример 2-8. Для иллюстрации описанной здесь процедуры рас смотрим движение спутника массы m около сферической планеты
72