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

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

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

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

Добавлен: 11.04.2024

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

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

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

ние можно оценить с такой же точностью, как и текущее, если начать оценку с момента /=0 при произвольно большой дисперсии

начальной ошибки а 2

0 > с Г ю О Ѵ

 

 

 

 

 

 

Теперь известна

вся необходимая информация для вычи:ления

C(t + T)

и P(t\t + T). Заметим, что

 

 

 

 

 

A(t)

=

G QJt) G' (0 Р - > ( Ф ) =

з2l

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) = (2th н . / ) Р ( ф + П -

 

 

 

— a2 ch2 |х* SCh2 (J. (t + T) [X 2

Cth2 (X (* +

7") -

o2

,

которое можно представить в виде

 

 

 

 

 

 

 

Р (* И +

7) — 2th [xrP (t\t+T)

 

=

 

 

 

= — а 2

[1 +

ch2 [х^ csch2

|х (г +

Г)].

 

(7-78)

Уравнение (7-78) имеет стандартную форму обыкновенного ли­ нейного дифференциального уравнения первого порядка, решение которого можно записать в виде

где

h (0 = - ° 1 П + c h 2 ^ c s < * 2 M ' + 7)];

21*

323


Ц — постоянная интегрирования, 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.

324


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).

325


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.

326


в) Для оптимального сглаживания в закрепленной точке при

=

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. Предположим, что состояние является непрерывным гауссовский процессом

327