и его решение, очевидно, |
|
|
|
|
B*(t)=e-'; |
t^O, |
|
|
поскольку |
ß*(0) =Р(010) = 1. Так как функции |
H(t) и R~l(t) |
тож |
дественно |
равны единице, из уравнения (8-94) |
следует, что |
коэф |
фициент передачи оптимального фильтра, сглаживающего в закреп
ленной точке A(t, t)=B*{t). |
Тогда из уравнения (8-93) получаем: |
x(Q\t) = |
e-*[z{t)—x(t\t)\. |
Уравнение (8-96) для дисперсии ошибки сглаживания в закреп ленной точке примет вид:
Р ( 0 | < ) = — е - " .
Так как Р (0 | 0) • ] , то
и для достаточно больших значений t ясно, что P(Q\t)—И/2. За метим, что в этом случае соответствующая дисперсия ошибки филь трации равна единице, и за счет сглаживания в закрепленной точке достигается выигрыш в точности, равный 50%. _
2. а > 0 при |
| а | > 1 . Здесь / > ( ф ) « 2 а и |
K(t)~2a, |
так что |
F(t)—K(t)H{t)=—а. |
Следовательно, |
В* = —аВ*, |
что дает |
|
|
B*(t)=2ae-at; |
f2?0, |
|
|
так как В*(0) =Р(010) =2а. Уравнение для дисперсии имеет вид Р(0|г)= — 4a z e~ 2 a t ; его решение, как легко убедиться,
P(0\t)=2ae~2a'.
Для достаточно больших t дисперсия P(0\t)—К). Следователь но, теоретически при сглаживании в закрепленной точке можно точ но определить х(0), в то время как соответствующая дисперсия ошибки фильтрации равна 2а при | а | > 1 .
3. а < 0 при |
[ о | < 1 . Система уравнений |
(8-97), |
очевидно, |
устой |
чива. Заметим, |
что F((\t)=a+ |
Ѵа2+\, |
и окончательный результат |
совпадает со случаем |
1, а именно, |
|
|
|
|
|
|
/>(0|о = -^"(1+*"">• |
|
|
4. о < 0 при | а | > 1 . В этом случае |
|
|
|
|
|
P (t |0 |
= |
а + |
VäF+~\ ^ а + |
I а | = 0, |
|
так как а отрицательно, |
и |
очевидно, |
улучшение |
точности |
оценки |
с помощью сглаживания в закрепленной точке |
является |
в этом |
случае излишним. |
|
|
|
|
|
|
|
З А Д А Ч И |
К ГЛ. 8 |
|
|
|
|
|
|
8-1. Доказать |
достаточность |
уравнения |
Винера — Хопфа, |
исследуя d2J/de2 |
при е = 0 . |
|
|
|
|
|
8-2. Записать уравнение |
Винера — Хопфа через |
соответствующие |
корреляционные |
функции. |
|
|
|
|
|
8-3. |
Как будет выглядеть |
уравнение Винера — Хопфа, |
если слу |
чайные |
процессы {x(ti), ti^to} |
и {z(x), |
to^x^t} |
имеют |
ненулевые |
математические ожидания? |
|
|
|
|
|
8-4. |
Как будет |
выглядеть |
уравнение Винера — Хопфа для за |
дач: |
|
|
|
|
|
|
|
а) |
оптимального |
предсказания? |
|
|
|
б) |
оптимального |
сглаживания |
на |
закрепленном интервале? |
в) оптимального |
сглаживания |
с постоянным |
запаздыванием? |
8-5. |
Записать и решить уравнение |
Винера — Хопфа для модели |
системы (7-1), (7-2) в задаче оптимального предсказания с посто янным упреждением x(t+T\t), где t^t0, а Т — положительная по стоянная («время упреждения»). Вывести также дифференциальное уравнение для корреляционной матрицы ошибки предсказания с по
стоянным упреждением |
x(t+T\t) |
=x(t+T)—x(t+T\t). |
8-6. В задаче оптимальной фильтрации показать, что весовая |
матрица оптимального фильтра имеет вид: |
|
A(t, x)=W(t, |
т)Д(т) |
для t, x^to, где Л'(т) |
определяется |
с помощью уравнения (8-53), |
a *¥(t, т) —переходная матрица состояния системы вида |
|
y=[F(t)-K(t)H(t)]y, |
где у — любой действительный я-вектор.
8-7. Пусть |
матрица передачи фильтра определена в предполо |
жении, что S 0 |
= 0 , т. е. с помощью соотношений |
|
|
K(t)=P(t\t)H'(t)R-i(l); |
|
|
P = F(t)P |
+ PF'(t)—PH'(t)R-4t)H(t)P |
+ |
G(t)Q(t)G'(t) |
для t^to, где |
P=P(t\t) |
и P(ta\ta)<=P(h) |
и |
пусть матрица S(t) |
в действительности не равна нулю. Какое дифференциальное урав
нение следует в этом случае решить, чтобы |
определить |
истинную |
корреляционную |
матрицу |
ошибки фильтрации? |
|
|
|
|
8-8. |
Если |
матрица R(t) |
неотрицательно |
определена для всех |
t^to, |
будет ли уравнение |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
дА (t |
і) |
|
|
|
|
|
|
|
|
F(t)A(t, |
х ) _ |
±~1~-A(t, |
t)H(t)A(t, |
x) |
= 0 |
|
при t0 < |
x < t |
эквивалентным |
уравнению |
|
|
|
|
|
|
t |
|
|
|
|
дА |
it |
і\ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
X |
|
|
F(t)A(t, |
|
x ) - |
± - L - A ( t , t)H{t)A{t, |
x) |
|
|
|
|
|
XE[z |
(x) z' |
(a)] c?x = |
0 |
|
|
|
|
при |
to^o<.t? |
Исследовать этот |
же вопрос |
в случае, |
когда |
матрица |
Q(t) |
положительно |
определена |
для t^t0, |
а |
матрица |
R(t) |
неотри |
цательно |
определена. |
|
|
|
|
|
|
|
|
|
|
8-9. |
Рассмотреть |
скалярную |
систему |
x = |
m(t), |
|
z(t)=x(t)+v(t) |
для |
t^O, |
где |
Q(0=9/4; |
R{t) = l; |
S(0 = 1/2, |
/>(*<>) =28, |
|
а) |
Определить P(t\t) и K(t) для всех |
t^O. |
|
|
б) |
При t—»-оо показать, что |
P(t\t)—И. |
|
|
в) Если коэффициент передачи фильтра определен |
без учета |
того, что S(0 = 1/2, т. е. в предположении |
S(0=0, чему |
равна фа |
ктическая дисперсия ошибки фильтрации при t—>-оо? |
|
|
8-10. |
Показать, что уравнение |
|
|
|
|
|
|
|
àA |
(t, т) |
|
|
|
|
возникающее |
при доказательстве |
теоремы |
об оптимальном сглажи |
вании |
в закрепленной |
точке, при ti^x^t |
эквивалентно |
уравнению |
С г дА (t, z) |
|
|
|
Е [г (т) z' (a)] dz = 0 |
\ |
I |
|
+A(t, |
t)H(t)A°(t, |
г) |
d t |
|
|
|
U |
|
|
Предположить, что матрица |
R(t) положительно |
опре |
при ti^a<t. |
делена для всех |
t^ti^to. |
|
|
|
|
8-11. Частица, |
которая может |
двигаться только в одном |
изме |
рении, покидает начало координат со случайной скоростью, распре
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
деленной по гауссовскому закону с нулевым математическим |
ожи |
данием и дисперсией |
"g = const |
0. |
|
|
|
|
|
|
|
|
|
|
Частицу |
ускоряет |
случайная сила, |
имеющая |
вид гауссовского |
белого |
шума {w(t), |
t^-.O} |
с |
нулевым |
средним и дисперсией з^—- |
= const>0. |
Положение |
частицы непрерывно |
измеряется |
в |
присутст |
вии аддитивного |
скалярного |
гауссовского |
белого |
шума |
{v(t), |
|
t^0} |
с нулевым |
средним |
и дисперсией |
o2 „ = const>0, |
независимого от |
{w(t)t |
t^0}. |
Для простоты |
предположим, |
что частица |
имеет |
еди |
ничную массу. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
а) |
Вывести, но не решать |
уравнения |
для оптимальных |
оценок |
положения и скорости частицы при |
t^0. |
|
|
|
|
|
|
|
|
|
б) |
Вывести, но не решать |
уравнения для оценки |
скорости, с ко |
торой |
частица покидает начало координат |
при t=0 |
с |
использова |
нием |
алгоритма |
оптимального |
сглаживания |
в закрепленной |
|
точке. |
8-12. Как выглядят уравнения оптимального сглаживания в за |
крепленной точке, если случайные процессы |
{w(t)t |
|
t^to} |
и |
|
{v(t). |
t^to} |
коррелированны и E\w(t)v'(x)]=S(t)à{t—т) |
|
для всех |
t, |
TJS^O? |
8-13. Для задачи |
оптимального |
сглаживания |
на закрепленном |
интервале to^t^tt |
|
показать, что: |
|
|
|
на закрепленном |
|
интер |
а) |
оптимальный |
фильтр, сглаживающий |
|
вале, описывается |
уравнением |
|
|
|
|
|
|
|
|
|
|
|
|
x(t |
К ) = F (0 x (t I U) + •-!.(/., i)\x(t |
| <,) - x |
(t |0] |
|
|
|
для |
г, приграничном |
условии x (tt |
\ tt) |
и при условии, что |
матрица передачи оптимального сглаживающего фильтра удовлетво
ряет соотношению A(t, <)='G(0Q(0G '(0^_ 1 (n0;
б) Корреляционная матрица соответствующей ошибки сглажи вания удовлетворяет матричному дифференциальному уравнению
^(i\ii)=lF(f) + 0(t)Q(t)Cr(t)P-i(t\t)]x.
Х Р ( ф і ) + £ ( ф і ) [P(() + G(t)Q(t)G'(t)P^(t\t)Y-r~G(t)Q(t)a'(t)
|
|
|
|
|
|
|
|
|
|
|
|
|
|
для |
Z o ^ / s s ^ i |
с граничным условием |
P(tt\tt). |
|
Предполагается, |
что |
процессы |
{w(t), |
t^to} |
и М О . t^to} |
взаимно |
независимы. |
|
|
|
8-14. |
Показать, |
можно ли модифицировать |
результаты |
задачи |
8-13 |
в случае |
коррелированных |
процессов |
{w(t), |
t^t0} |
и |
{v{t), |
t^k} |
при E[w(t), |
v'(%)}=S(t)b(t—x) |
для |
всех /, |
x^zto, |
и |
если |
можно, то каким образом. |
|
|
|
|
|
|
|
|
|
8-16. |
Вывести, |
но |
не решать |
уравнения для |
сглаженных |
на |
закрепленном интервале оценок положения и скорости частицы из
задачи 8-11. Рассмотреть систему на |
интервале |
времени |
О ^ г ^ Г , |
где Т постоянно. |
|
|
|
|
|
|
|
8-16. Из соответствующего |
интегрального |
уравнения |
Винера —• |
Хопфа получить уравнения оптимального сглаживания с |
постоян |
ным запаздыванием для г^Го, где запаздывание |
Г—const>0, в пред |
положении |
независимости |
|
случайных |
процессов |
{^(0, |
г ^ г 0 } и |
М О , t>to}. |
|
|
|
|
|
|
|
8-17. Для модели системы |
|
|
|
|
|
|
x = |
F (О X + G (t)w (0; |
|
|
|
|
z{t) |
= |
H |
(t)x(t) |
+ v(t), |
|
|
|
где t^to, |
предположить, |
что процесс |
М О . |
|
не гауссовский |
белый шум с нулевым средним, а коррелированный процесс, опи сываемый соотношением
|
|
|
|
|
i=A(t)v |
+ |
|
l(t), |
|
|
|
|
|
|
где t^to, |
A(t)—непрерывная |
|
матрица |
размера |
mxm, |
а |
{£(0. |
t^to}—m-мерный |
гауссовский |
белый |
шум с |
нулевым |
средним. |
Предположить, что |
вектор |
v(to) |
«е |
зависит |
|
от |
x(to), |
{£(0> |
i^-h), |
{w{t), |
i^to} |
|
и имеет положительно |
определенную |
корреля |
ционную матрицу E[v(to)v'(t0)]=V(to). |
|
|
Пусть |
для |
{|(0. |
|
справедливы соотношения |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
E[x(t0)l'(t)] |
|
= 0; |
|
|
|
|
|
|
|
|
|
|
|
E[W)¥W=W(t)ö(t-x); |
|
|
|
|
|
|
|
|
|
|
|
|
Ei[w(t)l'(x)]=0 |
|
|
|
|
|
|
|
при любых |
t, |
x^to, |
|
где W — матрица |
размера mx>n, |
положительно |
определенная и непрерывная. |
|
|
|
|
|
|
|
|
|
|
а) |
Можно |
ли |
в |
этом |
случае |
при |
определении |
оптимальной |
оценки |
x(t\t) |
для |
r^fo применить |
метод |
расширения вектора со |
стояния и теорему 8-1? Если можно, то |
составить |
соответствующие |
уравнения, а если нельзя, то объяснить почему. |
|
|
|
|
|
б) |
Ранее |
(см. задачи 4-13 и 6-19) |
была |
введена |
процедура |
формирования |
разностных |
измерений |
для |
дискретных |
измерений |
с коррелированной во времени ошибкой. Полученное разностное из мерение имеет ошибку в виде белого шума. Аналогично этой про цедуре в настоящей задаче для получения того же результата мож
но использовать |
дифференцирование. Определяя новое измерение |
как |
|
|
z* (0 — z(t) — A (0 г (О, |
показать, что это |
новое измерение можно записать в виде |
|
2*(0=я*(0*(0 + £(0. |
где |
|
|
|
|
|
|
|
|
И* (0 = H (t) + H (t) F(t)-A |
(t) Я (О, |
|
a {£(/), t^ta]—m-мерный |
|
гауссовский белый |
шум с нулевым |
сред |
ним, для которого |
|
|
|
|
|
|
Е [I (t) V (x)] = |
[H (t) G (t) Q (t) G' (t) M' (t) + W (t)] d (t - |
x); |
|
E [w (t) V (x)] = |
Q(t)G'(t)H> |
(t) S (t - x); |
|
|
|
E[x |
(<„) Ç' (01=0 |
|
|
при всех t, |
x^to. |
Заметим, |
что эти формулы |
требуют дифферен |
цируемое™ |
матрицы И(і) |
[Л. 8-2, 8-6—8-8J. |
|
|
|
в) Применить |
теорему |
8-1 к формулам |
п. «б» и получить |
урав |
нения оптимального фильтра. Составить структурную схему опти
мального |
фильтра, |
в котором |
не |
используется операция |
дифферен |
цирования |
z(t) |
[Л. 8-2, 8-6—8-8]. |
|
|
|
г) Показать, |
что для такого |
метода |
начальными |
условиями |
являются: |
|
|
|
|
|
|
|
x (to \to) = |
Р (to) H' (t„) [H |
(t,) |
P(t0) H' |
(t„) + V (t0)\- |
>z (t0) |
P(t0\ tt) = P (t,) - P (*,) H' (r0 ) [H (t0) P (t0) H' (t0) +
+V(t0))-<H(t0)P(tl>),
где матрица P(to)=E[x(t0)x'(t0)] |
предполагается |
известной |
[Л. 8-2, |
8-6—8-8]. |
|
|
|
|
д) Перечислить дополнительные условия, необходимые |
в п. «в», |
такие как дифференцируемость Н(і) и пр. Можно |
ли предположение |
о положительной |
определенности |
матрицы W(t) |
в п. «в» |
ослабить, |
потребовав, чтобы |
эта матрица |
была лишь неотрицательно опре |
делена? Объяснить |
ответ. |
|
|
|
Г л а в а д е в я т а я |
|
|
|
СТАТИСТИЧЕСКИ ОПТИМАЛЬНОЕ |
УПРАВЛЕНИЕ |
В ДИСКРЕТНЫХ ЛИНЕЙНЫХ СИСТЕМАХ |
|
Общая задача управления, оптимизирующего неко торый критерий качества поведения системы при воздей ствии возмущений и ошибок измерения, была кратко рассмотрена в гл. 1. В настоящей и следующей главах исследуется частный класс задач, входящих в подобное общее описание, в которых состояние системы представ ляет собой гауссовский марковский процесс, рассмотрен ный в § 4-2 и 4-3, а в модель динамики системы добав ляется аддитивный сигнал управления. За критерий качества принимается математическое ожидание квадра тичной формы от переменных состояния и управления на фиксированном интервале времени. Хотя этот крите рий качества и не обладает для задач управления той