ние можно оценить с такой же точностью, как и текущее, если начать оценку с момента /=0 при произвольно большой дисперсии
начальной ошибки а 2 |
0 > с Г ю О Ѵ |
|
|
|
|
|
|
Теперь известна |
вся необходимая информация для вычи:ления |
C(t + T) |
и P(t\t + T). Заметим, что |
|
|
|
|
|
A(t) |
= |
G (О QJt) G' (0 Р - > ( Ф ) = |
з2—l— |
th ^ |
= V- th v-t. |
(7-76) |
|
Подставляя |
это |
выражение |
в уравнение |
(7-72), |
получаем: |
|
|
С (t + |
Т) = |
fx [th |xr — th [л. (f + |
T)] С (t + |
Г) |
|
для |
г > |
0 , где C ( 7 ) = |
ß (7") = sch fx7\ |
|
|
|
|
|
Тогда |
|
|
|
|
|
|
|
|
dC |
|
|
|
|
|
- ç - |
= [th |л/ — th fx (/ -f |
7")] rf (fj-0- |
|
Интегрирование |
этого уравнения |
приводит к |
соотношению |
|
ch \i.t |
|
|
|
|
С (t + П -= Y ch ^(t+T) |
= |
Y c h |
^ S c h ^ ( ; + |
r ) ' |
где y — постоянная интегрирования. Для £=0, C(T) =В(Т) =sch \iT.
Следовательно,
откуда *y=l, так что |
Y ch 0 sch [x7 = sch [хГ, |
|
|
|
|
C(t+T)=àmT |
sch ii(t + T) |
|
|
|
(7-77) |
при t^O. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Для |
рассматриваемой |
задачи |
уравнение |
(7-73) |
приводится |
к виду |
|
|
|
|
|
|
|
|
|
|
|
|
P (t\ |
t + |
T) |
= |
2А {t) P (t |
11 + |
Г) |
- |
|
|
|
— R (t + |
7") С 2 (< + |
Г) K 2 (< + |
7) — G 2 |
(0 Q (0- |
|
Подставляя в это соотношение уравнения |
(7-74), |
(7-76) и (7-77), |
получаем |
уравнение |
|
|
|
|
|
|
|
|
|
|
|
P (t \t + 7) = (2[х th н . / ) Р ( ф + П - |
|
|
|
— a2 ch2 |х* SCh2 (J. (t + T) [X 2 |
Cth2 (X (* + |
7") - |
o2 |
, |
которое можно представить в виде |
|
|
|
|
|
|
|
Р (* И + |
7) — 2(х th [xrP (t\t+T) |
|
= |
|
|
|
= — а 2 |
[1 + |
ch2 [х^ csch2 |
|х (г + |
Г)]. |
|
(7-78) |
Уравнение (7-78) имеет стандартную форму обыкновенного ли нейного дифференциального уравнения первого порядка, решение которого можно записать в виде
где
h (0 = - ° 1 П + c h 2 ^ c s < * 2 M ' + 7)];
Ц — постоянная интегрирования, a с(0—интегрирующий множитель
вида
с (t) — exp J* (— 2,a th \xt) dt -•= exp |
2 th ptd (ixt) |
1
= exp (— 2 ln ch \xt) = -^Щц = sch2 \xt.
Следовательно, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
P (t |
11 + |
T) |
= |
— o^ch2 |
jxr ^ sch2 Ц.Г fl |
+ |
ch2 ^ |
csch2 (J. |
+ |
7)] df |
+ |
+ |
-r) ch2 |
|
= |
— |
|
ch2 |
ixt |
j" [sch2 |
jx/ - f csch2 |
fJ. |
(< |
+ |
7")] ^ (Н-О |
+ |
-f- |
vj ch2 |
[xr = |
— awov |
ch2 |
(J. |
< [th |
ixt |
— cth fx |
+ |
7")] |
+ |
v) ch2 |
|хЛ |
|
Но для / = |
0 согласно |
уравнению |
(7-75) P |
(0 | 7") = |
а щ о , cth |
[хТ. |
Следовательно, |
т) = |
0 и |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
для г^О. |
Р ( Ш + |
71) = |
з ш а в |
ch2 (х< [cth [x (г -4- 7) |
— th [xr] |
|
|
|
|
|
|
|
|
г, для которых [it ^2, |
|
|
|
|
|
|
Рассматривая |
значения |
|
можно |
получить |
полезное соотношение между |
P(t\t) |
и |
P(t\t |
+ |
T). |
|
Для |
рассматри |
ваемых значений t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
P(t\t)^ew°v; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
cth ix |
|
ch2 ixt |
д= |
— » |
|
|
|
|
|
|
|
|
|
|
|
|
(< -f |
T) ^ 1 + г*?-2^-2^; |
|
|
|
|
Тогда |
|
|
|
|
|
t h ^ s = 1 — |
2е~2^. |
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
P(t\t+T): |
|
|
|
|
[1 |
+ 2 е - 2 |
^ У - 2 |
^ |
- |
1 + |
|
le'2**} |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
r w , |
(1 + e |
-2|хГ |
|
|
|
|
|
|
|
|
Следовательно, |
для |
(х^5г2 |
имеем: |
|
|
|
|
|
|
|
|
|
|
|
|
|
Р ( * 11 + Т) = \ - ( 1 + e-W |
) Р (* |0. |
|
|
|
|
Из этого |
соотношения ясно, что сглаживание с постоянным за |
паздыванием |
позволяет |
добиться |
уменьшения |
дисперсии |
ошибки |
оценки до 50% по сравнению с возможной дисперсией ошибки оптимальной фильтрации.
Для |
данного |
примера соответствующие уравнения |
фильтрации |
и сглаживания имеют вид: |
|
|
|
|
Х ( а |
I а) = Л' (а) [г (а) |
— х (а | о ) ] , |
с 5* 0; |
|
|
X (0 I а) = |
В (а) К (о) [z (о) — х ( о | о ) ] ; |
0 < а < |
Т; |
x(t\t+T) |
= C(t |
+ T)K(t+T)[z |
(t+T)-~x(t+T\t+T)] |
+ |
|
+ |
A(t)[x(t\t+T)-x(t\t)], |
|
r > 0 , |
|
где x (0 10) = 0 и t = а — T.
x(o\0) = 0
В (а)
x(t\t + T)
C<t + T)
s = r
A(t)
£(o\o)=a
x(s I s)
БЗ T
Рис. 7-9. Оптимальный сглаживающий фильтр с постоянным запаздыванием из примера 7-6 (ключ К замкнут только в мо мент 0 = 7").
Структурная схема фильтра, соответствующая этим уравнениям, приведена на рис. 7-9.
Следует заметить, что здесь для К(а) следует применить какую-
либо аппроксимацию при малых 0, так как К(0) = оо в силу пред
положения Яд ^ > аюаѵ.
З А Д А Ч И К ГЛ. 7
7-1. Применить процедуру предельного перехода, использован ную в настоящей главе, к результатам задачи 5-11 и показать, что уравнения оптимальной фильтрации в этом случае принимают вид:
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x = F(t)x |
+ |
K (0 [г (0 - |
И (t)x\; |
|
|
|
|
К (t) =[P(t\ |
t) H' (t) + |
G(t)S |
(t)} R-• (0; |
|
|
P = |
[F (t) - |
G (t) S(l)R-> |
(t)] P |
+ |
P[F{t)-G |
(t) S (t) R- 1 (t)]' |
- |
- |
PH' (t) R-' (0 |
H(t)P+G |
|
(t) [Q (0 — S(t)R-* |
(t) S' {t)] G' (t) |
для системы (7-1), (7-2), |
если |
E [w (t) vf (x)] = 5 |
(t) S {t — x). |
В этих |
уравнениях |
t > tB, |
x = |
x (t \t), |
x(t0\ |
t0) = |
0, P =P (t \ t), |
P(t01 |
t0)= |
= P (t0) = |
E [x (*„) x' (/.)] |
(см. [Л. 7-1]). |
|
теорем |
7-1 |
и 7-2, |
|
7-2. Как следует |
изменить |
формулировку |
если |
x(to), |
{w(t), |
t^to) |
|
и |
{v(t), |
t^to] |
имеют |
ненулевые матема |
тические ожидания и присутствует известное управление u(t)7 За метим, что уравнение (7-1) для этой задачи принимает вид:
x=F(t)x+G(t)w(t)+C(t)u(i).
7-3. Если в уравнении (7-7) фигурирует произвольная матрица
передачи К* (t) |
вместо |
оптимальной матрицы, описываемой |
уравне |
нием (7-8), то |
какой |
вид будет иметь дифференциальное |
уравне |
ние для соответствующей корреляционной матрицы ошибки? Явля
ется ли в этом случае ошибка фильтрации гауссовским |
марковским |
процессом? |
|
|
|
|
|
|
|
|
|
|
|
7-5. В предположении, |
что при |
численном |
решении |
уравнения |
(7-9) на ЭВМ используется |
шаг Л>0, следует |
ли при вычислениях |
использовать |
матрицы: |
|
|
|
|
|
|
|
а) |
Q(t) |
и |
R(t)? |
|
|
|
|
|
|
|
б) |
Q(t)/h и |
R(t)/h? |
|
|
|
|
|
|
|
в) |
Q(t)h |
и |
|
R(t)h? |
|
|
|
|
|
|
|
Объяснить |
ответ. |
|
|
|
|
|
|
|
7-5. В |
системе |
связи |
передаваемое сообщение |
|
имеет вид |
x(t) = V sin Мог, |
где |
t^O, wo=const, |
а |
V — гауссовская |
случайная |
величина с нулевым средним и дисперсией а2 о. Принимаемый |
сигнал |
имеет |
вид z(t)—x(t)+v(t), |
где {v(t), |
t^O}—гауссовский |
|
белый |
шум с нулевым средним и постоянной |
дисперсией о 2 к , |
независимый |
от V. Вывести, но не решать уравнения, необходимые для реализа |
ции оптимального фильтра, вычисляющего оценку x(t\t), |
и составить |
структурную схему |
фильтра. |
|
|
|
|
|
|
|
7-6. Как будут выглядеть уравнения оптимального |
сглаживания |
на закрепленном |
интервале, |
если x(tB), |
{^(0. |
и |
{^(0, |
t^U} |
имеют ненулевые математические ожидания и присутствует извест
|
|
|
|
|
|
ное упраівление u(t) |
(см. задачу 7-2)? |
|
7-7. |
Вывести дифференциальное >равнение для взаимной корре |
ляционной матрицы |
Pii(i\ti) |
=E[x(t\ti)x'в |
случае оптималь |
ного сглаживания на закрепленном интервале. |
7-8. |
При некотором |
произвольном ^і>0 вывести уравнения опти |
мального |
сглаживания |
на |
закрепленном |
интервале [0, ti] для си |
стемы из примера 7-6 и составить структурную схему сглаживаю
щего фильтра. |
Показать, что P(t\tt)^P(t\i) |
для |
0 = £ ^ ^ і . |
|
7-9. Какому |
дифференциальному |
уравнению |
удовлетворяет |
ошибка сглаживания в закрепленной |
точке |
{x(ti\t), |
t^ti}? |
Пока |
зать, что этот процесс является гауссовским марковским второго порядка.
7-10. Показать, |
что |
утверждение теоремы |
7-4 |
можно также |
представить в виде |
|
|
|
|
|
|
x(t, \ t) = |
В* (0 |
H' (*)/?-• (t) [z (t) - |
H (t) x{t |
I /)]; |
|
B* |
= |
B* [F (t) — К (t) H (t)]'; |
|
|
|
|
|
ß * ( ' i ) = ^ ( M ' . ) ; |
|
|
|
P(tt \t) |
: |
= _ |
ß * (t) H' (t) R-l(t) |
H |
(t)B*'(t), |
где t ^ t i . Заметим, что эта формулировка позволяет избежать вы
числения Р _ 1 ( ф ) .
7-11. Рассмотреть систему из примера 7-1 при а = 0 и О д = о № в щ .
а) Для задачи оптимальной фильтрации показать, что в этом примере
P (t 11) = Яд и К (t) = .ч = а щ 5 , для всех t 5> 0.
б) Замечая, что оптимальный фильтр представляет линейную стационарную систему, определить ее передаточную функцию и по казать, как можно реализовать фильтр в виде простой RC-uenn.
в) Для оптимального сглаживания в закрепленной точке при
= |
0 показать, что |
В (t) |
= |
и Р (0 \t) = |
(о£/2) (1 + <?_21*') для |
^ О . |
Указать |
возможность |
уменьшения |
дисперсии начальной |
ошиб |
ки до |
50%. |
|
|
|
|
|
|
|
|
|
г) Для |
оптимального |
сглаживания с |
постоянным запаздыванием |
показать, |
что |
С (/ + |
Т) |
= |
и Р (t |
| * + |
T) = |
(og/2) (1 + |
е~2^) |
при |
^^0. |
Указать |
возможность уменьшения |
дисперсии |
ошибки |
оценки до |
50% |
при любом |
t^O. |
|
|
|
|
д) Составить структурную схему объединенного оптимального фильтра и оптимальных фильтров, сглаживающих в закрепленной точке и с постоянным запаздыванием.
7-12. Доказать теорему 7-5.
7-13. Показать, что при нулевом запаздывании (Г=0) утверж дение теоремы 7-5 сводится к утверждению теоремы 7-1.
Г л а в а в о с ь м а я
ОПТИМАЛЬНЫЕ ОЦЕНКИ В НЕПРЕРЫВНЫХ ЛИНЕЙНЫХ СИСТЕМАХ. АНАЛОГОВЫЕ
АЛГОРИТМЫ
В настоящей главе рассматривается иной подход к задаче оценки в непрерывных линейных системах по сравнению с методами предыдущей главы. Здесь задача с самого начала решается в непрерывном времени без обращения к процедуре предельного перехода. Полу ченные результаты совпадают с результатами, получен ными в гл. 7.
После постановки задачи будет получено матричное интегральное уравнение, представляющее собой необхо димое и достаточное условие оптимальности оценки. Это интегральное уравнение, называемое уравнением Винера — Хопфа, выводится здесь с использованием классических вариационных методов. В остальных пара графах главы уравнение Винера — Хопфа решается для задач оптимальной фильтрации и оптимального сглажи вания в закрепленной точке. Решение остальных задач предоставляется читателю в качестве упражнения, по скольку эти результаты известны из гл. 7.
8-1. ИНТЕГРАЛЬНОЕ УРАВНЕНИЕ В И Н Е Р А — Х О П Ф А
Вначале рассмотрим задачу оптимальной оценки, несколько отличную от задачи § 7-1. Предположим, что состояние является непрерывным гауссовский процессом