Файл: Медич Д. Статистически оптимальные линейные оценки и управление.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 11.04.2024
Просмотров: 245
Скачиваний: 1
Но
У —7 |
II PxyPxz ||> |
|
Z 2 |
||
|
( I I I I J I K - ^ ' ( 3 ^ ' 0 r =
Рѵѵ |
0 |
- 1 |
р - 1 |
0 |
|
^ УУ |
|
||
0 |
р |
|
0 |
р - і |
|
|
|
|
' ZZ |
где использовано равенство Pj/Z = 0, поскольку случайные векторы у и z независимы. Следовательно,
|
|
|
|
|
|
|
|
«/ — У |
|
|
|
|
|
|
|
О |
Р.,—1 |
г — |
г |
|
= * + л * * 3 " 1 |
(у - у ) + ( z - s ) = |
|
||||||
|
|
|
= £ С ф ) + £ ( . * | г ) - * . |
|
|
||||
4. Если |
г/ и z не обязательно |
независимы, то |
|
||||||
|
|
£(х|г/, |
z ) = £ ( * | & z); |
|
(3-50) |
||||
|
|
Е(х\у, |
z)=E(x\y) |
|
+E(x\z)—x, |
|
(3-51) |
||
где |
z=z—E(z\y). |
|
|
|
|
|
|
||
Уравнение (3-51) сразу следует из свойств |
2 и 3, до |
||||||||
казанных |
выше. |
|
|
|
|
|
|
|
|
Чтобы |
доказать (3-50), убедимся, что в силу уравне |
||||||||
ния |
(3-41) |
|
|
|
|
|
|
|
|
|
Е(х\у. |
г) = |
х\+\\РіуРх |
РУУ |
РѴг -1 |
У — У |
|||
|
|
|
|
|
|
Ргу |
Р гг |
Z — 2 |
|
|
Из уравнений |
(3-36) — (3-38) |
получим: |
|
|
||||
|
|
|
|
Руу |
Руг |
|
|
|
|
|
|
|
|
Ргу |
Ргг |
|
|
|
|
|
|
РУУ + РУУ |
РугСРгУРш |
-Р^\Р«гС |
|
|
|||
|
|
|
|
|
,—1 |
|
|
|
|
|
|
|
~СРгуР\*УГуу |
|
|
|
|
||
где |
|
|
|
|
|
|
|
|
|
|
|
С |
(PzZ |
Pzy^uu Pyz) |
* ' |
|
|
116
Следовательно,
|
|
Е {х\у, |
z) |
= |
X + {РхуР-1 |
|
+ |
Piyp->PyzCPzyP^ |
|
|
|
- |
|
|||||||
|
- |
P*zCPZyP-1) |
|
(у |
- |
у) |
! + |
(Р«С - |
P*yP-'PyzC) |
|
(z |
- |
г) = |
|||||||
|
|
= |
* + I V " 1 |
+ |
iP^Pyz |
|
- |
P*z) |
CPzyP^\ |
|
X |
|
||||||||
|
|
|
|
|
Х(У-У) |
|
+ |
(P*z - |
P^Pyz)C(z—z). |
|
|
|
(3-52) |
|||||||
|
С |
другой |
стороны, |
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
E(x\y, |
г) = |
Е(х\у) |
+ Е(у\г')-х |
|
= |
|
|
|
|||||||
|
|
|
= |
* + |
P*»P~l |
(У-У) |
|
+ |
Р „PZl |
[г - |
|
E{z)\. |
|
|||||||
Ho |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
£(z) = |
£ [ z - £ ( z | t , ) ] = |
0; |
|
|
|
|
|||||||
|
|
|
|
|
P |
„ = |
E{(x~x)[z-E{z)}>} |
|
|
= |
|
|
|
|
||||||
|
|
|
|
|
|
xz |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= |
E{(x-x)[z-z-PzyP~l |
|
|
|
|
(y - |
y)]'} |
= |
|
|
||||||
|
|
|
|
|
|
|
|
— |
P*z |
|
PxyPyy |
Pyz- |
|
|
|
|
|
|
||
|
Кроме |
того, |
|
|
|
|
PzyP~x |
|
|
|
|
|
|
|
||||||
|
|
P^ |
z z |
= |
E{[(z- |
|
z) - |
{y - |
y)} [(z |
- |
z |
) |
- |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
- |
PryP'l |
{У - |
У)}'} |
|
= |
Pzz |
- |
PzyP^Pyz- |
|
|
|
|||||
|
Следовательно, |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
E |
(x\y. |
z) = |
X + |
|
P^P-1 |
(y - |
|
jr) + |
|
|
|
||||
+ |
(/>« |
- |
P*yP~lpyz) |
(Pzz |
- |
V " |
' |
Py r ) " 1 |
[Z — |
Я |
( ф ) ] |
= |
||||||||
|
|
|
= |
* + |
|
|
(У - |
|
|
+ |
(P« |
- P*yP~lPyz) |
|
X |
|
|||||
|
|
|
|
|
X C [ ( z - 2 ) - P n , P - , 0 / - ^ ) ] = |
|
|
|
|
|||||||||||
= |
* + |
|
l * V £ |
+ ( Л ^ Ч * |
- |
P*z) CPzyP-1] |
(У-У) |
+ |
||||||||||||
v |
|
|
|
|
+ ( ^ z - - P x y P - 1 P y 2 ) C ( z - 2 ) . |
|
|
(3-53) |
Из уравнений (3-52) и (3-53) сразу следует уравнение (3-50).
117
Заметим, что если бы все случайные векторы имели нулевые математические ожидания, то многие преобра зования существенно упростились. Поскольку математи ческое ожидание случайного вектора изображается сред ней точкой евклидова пространства, около которой рас пределены точки, изображающие значения случайного вектора, часто бывает удобно перенести начало коорди нат с помощью формулы % — х—x и проводить все опера ции со случайным вектором g, имеющим нулевое мате матическое ожидание. Ненулевое математическое ожида ние можно ввести вновь в любое удобное время. В дальнейшем обычно будет использоваться именно та кой подход.
З А Д А Ч И |
К ГЛ. 3 |
|
|
|
|
|
|
|
|
|
|
|
||||
3-1. Вывести уравнение (3-3). |
|
|
|
|
|
|
|
|
||||||||
3-2. |
Как |
выглядит функция |
f(x\y) |
в |
примере |
3-2 |
для |
случаев |
||||||||
ï = + l и у |
=—П |
|
|
|
|
|
|
|
|
|
|
|
|
|||
3-3. |
Изобразить функцию f(xi, |
xz) |
из примера 3-3 |
для Хі = х2 и |
||||||||||||
отметить на рисунке х. |
|
|
|
|
|
|
|
|
|
|
|
|||||
3-4. Доказать четыре основных свойства условного математиче |
||||||||||||||||
ского ожидания, |
приведенные на стр. 98. |
|
|
|
|
|
|
|||||||||
3-5. |
Вывести |
уравнения |
(3-36) — (3-38). |
|
|
|
|
|
|
|||||||
3-6. |
Пусть |
x — гауссовский |
случайный п-вектор с |
математиче |
||||||||||||
ским ожиданием х и корреляционной матрицей |
Рхх. |
Показать, что |
||||||||||||||
вектор у=Ах+Ь, |
|
где А—матрица |
размера тХп, |
а Ь — постоянный |
||||||||||||
гсг-вектор, является гауссовский |
случайным |
вектором |
и определить |
|||||||||||||
его математическое ожидание и корреляционную |
матрицу. |
|
||||||||||||||
3-7.. |
С помощью гауссовской плотности распределения |
(3-33), |
||||||||||||||
используя |
определение Е(х) |
и Е[(х—х)(х—г)1], |
показать, что Е(х) = |
|||||||||||||
= х и Е[(х—х) |
(х—х)'] = |
Р. |
|
|
|
|
|
|
|
|
|
|
||||
3-8. |
Пусть |
x — гауссовская |
случайная |
величина |
с нулевым мате |
|||||||||||
матическим |
ожиданием и |
дисперсией |
Р = |
о"2 >0. |
Показать, |
что А-й |
||||||||||
момент |
x, |
определяемый |
соотношением |
р-ь = £ ( * * ) , |
где k—\, 2, |
|||||||||||
имеет вид |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 X 3 X |
•• • X |
{k — 1) о* для четных k; |
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
для |
нечетных |
k. |
|
||
Указание. |
|
Использовать |
характеристическую |
функцию случайной |
||||||||||||
величины и вывести соотношение |
|
|
|
|
|
|
|
|
~ds* «=0
где }=Ѵ—1. Этот результат в n-мерном случае называется теоремой моментов [Л. 3-4].
3-9. Если x — гауссовский случайный п-вектор, то следует ли отсюда, что ^-вектор, включающий любые k компонент вектора х,
m
при k=\, |
2, ..., h—i также является гауссовским случайным |
векто |
|||||||||
ром? |
|
|
|
|
|
|
|
|
|
|
|
3-10. Показать с |
помощью |
функции |
плотности |
распределения |
|||||||
(3-35), что два некоррелированных |
гауссоЕСких |
случайных |
вектора |
||||||||
независимы в предположении, что матрицы РХх |
и Руѵ |
положительно |
|||||||||
определены. |
|
|
|
|
|
|
|
|
|
||
Остальные задачи имеют своей целью |
иллюстриро |
||||||||||
вать |
возможность |
применения |
некоторых |
результатов |
|||||||
главы для получения оптимальных оценок |
случайного |
||||||||||
вектора. |
|
|
|
|
|
|
|
|
|
||
3-11. |
Пусть X и у— |
случайные векторы с нулевыми математиче |
|||||||||
скими ожиданиями п- и m-мерные соответственно. |
|
|
|||||||||
Предположим, |
что |
вектор у |
известен, |
и |
обозначим через jT= |
||||||
= h (у) |
|
оценку х, |
где |
|
h — определенная |
л-мерная |
вектор-функция |
||||
вектора у; т. е. значение |
Л(-) известно при любом заданном |
значе |
|||||||||
нии у. |
Поскольку ошибка |
оценки, |
определяемая |
здесь |
как х |
=х—х, |
не будет в общем случае нулевой, требуется определить х таким об разом, чтобы минимизировать некоторую меру ошибки х. Одной из возможных мер ошибки является величина
называемая среднеквадратической ошибкой. Разумеется, возможны и другие меры ошибки, обычно называемые функциями потерь, функ циями штрафов, показателями или критериями качества и т. д. По казать, что значение х, минимизирующее среднеквадратическую
ошибку L, определяется выражением х = Е(х]у). Указание: Использовать тот факт, что
и |
положить |
|
£* [g (*)] = |
£« {£*[g |
(*)ЫЬ |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
g(x) |
~?'& |
= |
(х—х)' |
{х — х). |
|
|
|
|
|||
|
Этот результат показывает, что для квадратичного критерия ка |
|||||||||||||
чества наилучшей или оптимальной |
оценкой |
является |
условное мате |
|||||||||||
матическое |
ожидание. |
|
|
|
|
|
|
|
|
|
|
|
||
|
3-12. Как выглядит выражение для оптимальной оценки в задаче |
|||||||||||||
3-11, если X и у |
независимы? |
|
|
|
|
|
|
|
|
|
|
|||
|
3-13. Пусть |
X — «-мерный, |
а |
у — /«-мерный |
случайные |
векторы |
||||||||
с |
нулевыми |
математическими |
ожиданиями |
и известной |
совместной |
|||||||||
плотностью |
распределения f(x, |
у). |
Предположим, |
что |
значение у |
|||||||||
известно и нужно определить |
оценку х |
вида |
х=А°у, |
где |
А0—матрица |
|||||||||
размера пхпг. Предположим |
также, что А0 |
требуется |
отыскать среди |
|||||||||||
множества всевозможных матриц А размера пХпг |
таким |
образом, |
||||||||||||
чтобы минимизировать |
критерий качества |
|
|
|
|
|
|
|||||||
|
|
|
|
L = |
|
E[(x-Ay)'(x-Ay)). |
|
|
|
|
119