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

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

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

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

Добавлен: 11.04.2024

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

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

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

состояния

x(k), минимизирующую

критерий

качества ви­

да J[x(k\j)]

= E{L[x(k\j)]},

где L —допустимая

функция

потерь.

 

 

 

 

 

 

Заметим, что в

силу

свойств

условного

математиче­

ского ожидания

 

 

 

 

 

E7{L [x (k I /)]} =

Et.[E7{L

[x (k I /) J| г (1),

. .,

 

где 2* обозначает /га-мерный вектор вида

ІИ>)||

IIz С') I

Так как внешнее математическое ожидание в пра­

вой части равенства не

зависит от выбора

x(k\j),

то

минимизация Е {L [x (k \ /)]} эквивалентна

минимизации

Е {L[x(k\j)]\z(l),

2(2),

2(/)}. Поэтому уравнение

(5-1)

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

 

 

J[x(k\i)]

= Ex{L[x(k\j)]\z(\),

z(j)}.

(5-2)

Следовательно, значение критерия качества, как можно было интуитивно ожидать, обусловлено факти­ ческими значениями измерений. Уравнение (5-2) можно рассматривать как определение критерия качества.

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

Пример 5-1. Предположим, что одна из составляющих передавае­ мого сообщения представляет собой случайную величину x ( k ) , где£ =

= const, неизвестную получателю.

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

гается воздействию мультипликативного шума v t ( i ) в канале

переда­

чи и аддитивного

шума и2 (0

в приемнике, так что принимаемый сиг­

нал в / дискретных точках

может быть представлен в виде

г ( і ) =

— x(k)Vi(i)+v2(i),

i'=l, 2. ..., /, где &=const. Вся система изображе­

на в виде структурной схемы на рис. 5-3. Задача заключается

в опре­

делении значения x(k) по известным измерениям 2(1), z(2), .. ., z ( j ) .

Сформулируем

ее как задачу

оптимальной

оценки,

утверждая,

что любая оценка

x (k) вида x {k\j) =

<р [г ((')• і =

1, 2

/] являет»

180


u,(i)

u,(i)

 

Принимаемый,

x(k)

сигнал

Передатчик

г(і)

 

i=l,2,...,j

Рис. 5-3. Структурная схема системы связи.

'ся оптимальной, если у—скалярная функция / переменных, миними­ зирующая критерий качества вида / [х (k\j)] = £ [\х (k) — х

В этом примере функция потерь представляет собой просто модуль ошибки оценки.

Основы теории оценок

Представим условную функцию распределения ве­ роятностей x(k) при условии z ( l ) , ..., z(y') в виде

Р[х(Ь)<1\г*(І)]

=

Г[1\г*(І)]

z*(j)—/m-мерный

для всех n-мерных векторов

\, где

вектор вида

 

 

 

 

 

 

 

 

z(/)

 

 

 

Решение задачи оценки при выполнении

определен­

ных условий, накладываемых на F[l\z*(j)],

дается сле­

дующей теоремой.

 

 

 

 

 

Теорема

5-1. Если функция

L есть допустимая

функ­

ция потерь

и

функция

F[l\z*(j)]

 

симметрична

относи­

тельно

своего

 

математического

ожидания

I,

т. е.

F{%—1|г*(/)]=1

 

F ff— l\z*(i)]

 

для любого

g, и

выпукла

для всех l ^ i

т.е. F{1^+(\-К)Щг*

 

(j)]^lF[^\z*_{j)]

+

+ ( 1 — К)Filz\z*(j)]

для

любых

п-векторов I 1 , 12<СІ при

0<czX<c\,

то оптимальная оценка

имеет вид

 

 

 

 

 

 

x(k\j)

=

E[x{k)\z*(j)}.

 

(5-3)

Этот фундаментальный и мощный результат принад­ лежит Шерману [Л. 5-1]. Другими словами, в теореме утверждается, что при выполнении требований симмет­ ричности и выпуклости условной функции распределения оптимальная оценка представляет собой условное мате­ матическое ожидание. Очевидно, оптимальная оценка

ISi


определяется

сразу,

если

известна

условная

функция

распределения

F{x(k)

\z*(j)]

или, что то же самое, услов­

ная плотность

распределения f[x(k)

\z*(j)], при

условии,

что она существует. Однако более важно заметить, что

E[x(k)\z*(j)]

— оптимальная

оценка

при

любых комби­

нациях допустимых функций

потерь

для

симметричных

и выпуклых условных функций распределения.

Пример функции распределения, симметричной отно­

сительно

своего математического ожидания и выпуклой

1,0

Выпуклость

Симметрия

І

Рис. 5-4. Иллюстрация симметричной и выпуклой условной функции распределения вероятностей случайной величины.

для всех значений аргумента, меньших или равных ее математическому ожиданию, показан на рис. 5-4.

Очевидным следствием теоремы 5-1 является то, что при независимых случайных процессах {x(k), kŒl} и

{z(i), / = 1 , ..., /} оптимальная оценка равна x(k\j)

=x(k)

при любых k и /.

 

Доказательство теоремы основано на следующей лем­ ме, известной из теории вероятностей [Л. 5-2, 5-3].

Лемма

5-1.

Пусть

L — допустимая

функция потерь,

а у — n-мерный

случайный

вектор с нулевым

математи­

ческим

ожиданием и плотностью

распределения,

симмет­

ричной

относительно

этого математического

ожидания

(нуля)

и

выпуклой

для

всех

у^О.

Тогда

E[L(y)]^

^E{L(y—а)]

для любого п-вектора а.

 

 

 

Доказательство теоремы 5-1 с помощью леммы пред­

ставляется читателю в качестве упражнения.

 

В

связи

с теоремой 5-1 заметим,

что

оптимальная

оценка также имеет вид E[x(k)\z*(/)],

если

требование

выпуклости

функции

/ 7 [ | | г * ( / ) ]

заменить

условием вы­

пуклости функции потерь L {Л. 5-4].

 

 

 

182


Теперь рассмотрим следствие теоремы 5-1, которое относится к частному случаю более общей задачи, по­

ставленной выше.

 

 

 

 

 

 

Следствие 5-1. Если

функция

L

является

допустимой

функцией

потерь, a

ix(k),

k<=I} и

{z(i); t = l , ...,

/} — га-

уссовские

случайные

процессы с

дискретным

 

временем,

то

 

 

 

 

 

 

 

 

x{k\jf)==E[x{k)\z*{j)\.

 

 

 

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

Следствие

сводится

к

теореме

5-1, если удастся показать, что гауссовская

условная

функция

распределения

F[x(k)\z*(j)]

симметрична от­

носительно своего математического ожидания и выпукла

для всех |<;' ; | .

 

 

 

 

 

 

Опуская для простоты

аргументы,

имеем:

 

F(Ç| 2 *) = j

j aexp [

- - i

- ( x - ï ) ' Q - > ( x - f ) dx,

 

 

 

 

 

 

— О

 

 

 

 

 

(5-4)

г д е

 

 

 

 

 

 

t=X+P

.P"1

lz* — ?•);

 

 

1

xz* z*z*v

 

 

 

 

P

 

JP~..P.

;

(5-5)

 

 

 

xz* z*z*

z*x'

 

 

 

 

 

Ѵ{2ъГ

I

Q I

 

 

в предположении, что матрица Рг-г* размера jmxjtn и матрица Q размера пХп несингулярны.

Сначала докажем симметричность. Заменяя g в урав­ нении (5-4) на I\ получаем:

 

 

a ехр

(x-t.yQ->X

-ОО

оо

 

(5-6)

Аналогично

Х(х-Щах.

 

 

 

P ( È - S | 2 * ) = j . . . j

a e x p [ - 4 - ( x - ^ Q - ' X

—00

—00

 

 

X(x-l)]dx,

183