Файл: Медич Д. Статистически оптимальные линейные оценки и управление.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