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

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

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

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

Добавлен: 11.04.2024

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

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

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

где

t^lo,

с начальным

условием

в виде

матрицы

пере­

дачи

 

оптимального

фильтра,

сглаживающего

в

закреп­

ленной

точке to для

момента

t0+T

: С(t0

+ T) = В (U + T).

3)

Ошибка

сглаживания

 

с

постоянным

запаздыва­

нием {x(t\t + T),

t^zto}

является

гауссовский

марковским

процессом

второго порядка

с нулевым

средним

и

корре­

ляционной

матрицей,

удовлетворяющей

 

линейному

мат­

ричному

дифференциальному

 

уравнению

 

 

 

 

 

 

Р (t\t -\-T)

\F (t) +

А (01Р (ft

+

Т) +

 

 

 

 

 

+ P (t\t +

T)

[F (t) +

A (і)] ' ~

С (t -f- T) X

 

X К (t -f

T) R (t + T) K' (t - f

T) С

(t +

T) -G(t)Q(t)U>

(t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(7-73)

при

t ^ t 0

с начальным

условием

в

виде

корреляционной

матрицы

ошибки

сглаживания

в закрепленной

точке to

для

момента t0 + T:

 

 

 

 

 

 

 

 

 

 

 

 

 

P(t*\t0

+ T)=E[x(t0\t0

 

+ T)x'(t0\t0

 

+ T)].

 

Д о к а з а т е л ь с т в о . Доказательство теоремы

пре­

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

Как и в дискретном случае, алгоритм оптимального сглаживания с постоянным запаздыванием в непрерыв­ ном времени требует наличия двух корректирующих членов в сглаживающем фильтре и зависит как от ре­ зультатов оптимальной фильтрации данных измерения, так и от оптимального сглаживания этих данных в за­ крепленной точке. Эта зависимость представлена в табл. 7-1 и еще более подчеркнута в структурной схе­ ме на рис. 7-7.

На рисунке отмечено, что

для сглаживания

нужны

две шкалы

времени.

Шкала

времени а используется

в фильтре

Калмана

и оптимальном фильтре,

сглажи­

вающем в закрепленной точке. Шкала времени t исполь­ зуется для сглаживающего фильтра с постоянным за­ паздыванием. Фильтр с постоянным запаздыванием бездействует, пока a<to + T. Причина этого заключается в том, что он должен «ждать» в течение времени Т, пока время измерения не достигнет значения tQ + T и можно будет начать обработку данных. В течение этого вре­ мени фильтр, сглаживающий в закрепленной точке, обрабатывает данные для получения необходимых на-

319



x(t0\ )=0

В(в) с

X(t\t+T) +

G = tn + T

Fft)

A(t)

x(to\to)=0

x(e s)

z(e)

t—s-т

Рис. 7-7. Структурная

схема

оптимального

сглаживающего

фильтра с постоянным

запаздыванием (ключ К замкнут толь­

ко в момент a = to+T

для получения

начальных условий).

 

 

 

 

Т а б л и ц а 7-1

Требования

ко входным

данным

фильтра

с постоянным

запаздыванием

Источник дгнных

 

 

Требуемые

данные

Оптимальный фильтр

Оптимальный фин. ф, сглаживающий в за­ крепленной точке

 

К (t +

T) [z (t +

T)-

 

-H{t

+

T)x(t +

T\t +

T)\

для

всех t ^

t0

 

 

 

I

x(t\t)

для всех

t 3 Ï

tB

x(t0\t0 + Т) в

качестве

•••ачг.пі ного усло­

вия,

требуется

 

сглаживали f закрепленной

точке

ta в течение времени

 

320


чальных условий x(t0\t0 + T). Ясно, что фильтр, сглажи­ вающий в закрепленной точке, не должен работать при

a>to+T, а

фильтр Калмана должен работать и при

этом условии.

 

Матрица

A (t) = G(t) Q (t) G'(t) P-^ (t\t),

являющаяся

матрицей передачи фильтра, сглаживающего на закреп­ ленном интервале, присутствует в уравнениях (7-71 ) — (7-73). Здесь, однако, ее следует вычислять в прямом времени. Обращения матрицы P(t\t) можно избежать, непосредственно решая уравнение (7-52) с начальным

условием M (t0)

= Р - 1 ( ^ 0 | г 0 ) . Для

этого,

очевидно, тре­

буется, чтобы

матрица

P(t0\t0)

=P(t0)

была

несингу­

лярной.

 

 

 

 

 

Уравнение

(7-72) для

матрицы передачи

фильтра

с постоянным запаздыванием представляет собой систе­ му п2 обыкновенных линейных дифференциальных урав­ нений. Начальными условиями для этих уравнений яв­ ляются п2 элементов матрицы B(t0 + T), которая может быть получена только после предварительного решения уравнения (7-58) на интервале t0^t<^t0 + T.

Аналогичная ситуация возникает при расчете на­ чальных условий для уравнения (7-73).

Пример 7-6. Рассмотрим сглаживание с постоянным запазды­ ванием для простой системы связи, в которой передаваемое сооб­ щение можно представить как выходной сигнал интегратора с вход­ ным сигналом в виде белого шума с нулевым средним. Уравнение

передаваемого

сообщения имеет вид:

 

 

 

 

 

 

 

 

x=w(t)

 

 

 

 

для t^O,

где

x (0)—гауссовская

случайная

величина

с

нулевым

средним

и дисперсией

а^,

a {w(t),

t~^0}—гауссовский

белый шум

с нулевым средним и дисперсией

, независимый от х ( 0 ) .

 

Предположим, что рассматривается следующая модель прини­

маемого

сигнала:

 

z(t)=x(t)+v(t)

 

 

 

 

 

 

 

 

 

 

 

 

для t^O,

где

{v(t),

t^O}—гауссовский

белый

шум

с

нулевым

средним

и дисперсией

,

независимый от

х (0) и

{w(t),

t^Q}.

Также

предположим,

что

Од^>

атзѵ.

 

 

 

 

Заметим, что эта модель является частным случаем модели из примера 7-1 при а—0. Из примера 7-1 имеем:

Рі = »«j't,; Рг = — Vw°v>

21—85

321


Тогда дисперсия ошибки оптимальной фильтрации составит:

P(t\t)=

а ш а г YZT^Wt -= а ^0 -" ^ _ с - и - /

ИЛИ

для ^ О . Кроме того, имеем:

 

 

1

 

 

 

 

 

 

Д ' ( / ) = — Р(< 10 =

l J - cthfA/ .

(7-74)

 

 

°*

 

 

 

 

Зависимость

/'(ОО/0 "*0

о т

изображена

графически

на

рис. 7-8, откуда

можно видеть,

что для \it^2

P{t\t)/aKav~\.

 

Pit

I t)

 

 

 

 

 

 

 

 

Рис. 7-8. Дисперсия ошибки

 

 

 

оптимальной

фильтрации

из

 

 

 

примера 7-6 в

виде функции

 

 

 

времени.

 

 

Теперь получим используемые в дальнейшем выражения для коэффициента передачи и дисперсии ошибки оптимального сглажи­ вания в закрепленной точке. Из уравнения (7-58) имеем:

É ' . = -

°1

ѵ Ѵ В t h W = -

В (Mh f*0

 

 

для t^O при начальном условии ß(0) = l .

 

 

 

 

Решение этого уравнения

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

В

W = сТТй = s c h

 

 

 

 

получается с помощью метода разделения

переменных.

 

 

При / і = 0 уравнение

(7-59)

здесь

принимает вид:

 

 

Р (0 I 0 = — R В2

(0 Я 2

(0 =

°l sch2 Mi ((Д.2 cth2

[л0

=

 

 

= —o^csch2 ^

 

 

 

 

для <^=0. Отсюда следует, что

 

 

 

 

 

 

Я (0 I 0 = -

«£ - ^ - j " esch2

f*W ((о-О =

cth

 

(7-75)

Попутно заметим,

что P(0\t)=P(t\t)

для всех

t^O.

Это озна­

чает, что при сглаживании в закрепленной точке начальное

состоя-

322