Поэтому для неубывающей функции l(v) при v ^ O
J |
l{v){v*-D)f~~S° ' d v ^ : l { V D ) Y |
Vd |
|
|
|
уъ |
_ |
v- |
----- |
I |
l ( v ) ( D - v 2)e |
20 |
d v ^ l ( V D ) y |
причем знак равенства в этих формулах одновременно не может иметь места, так как функция / ip) не является постоянной во всем диапазоне изменения v. Следовательно, правая часть формулы (12.29) строго положительна, т. е. [F"iQ) ]*=о > 0. Это служит доказательством того, что при х — 0 функция F(Q) достигает мини мума. Легко видеть, что х = 0 является в этом случае единствен ным решением уравнения F'(0) = 0. Поэтому других экстремумов
функция /•’(б) не имеет. |
|
|
|
0 (t, |
гг, . . |
zN) |
в формуле |
Для нахождения искомой функции |
(12. 24) достаточно положить х — 0. В результате получаем |
|
|
N |
|
|
|
N |
|
|
|
|
|
|
0 — |
2 |
v V r - |
|
2 |
|
|
|
= °- |
|
|
|
р, |
г=1 |
|
|
р, <7, / • = ! |
|
|
|
|
|
|
Отсюда находим |
|
|
|
|
N Г N |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
е = е [t, |
zl t ... , |
zN) — |
S |
2 |
|
M >p(0 |
zr + |
|
|
|
|
|
|
N |
|
r=i Lp=i |
|
|
|
|
|
|
|
|
+ |
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
P , |
q, r=1 |
|
|
|
|
|
|
|
|
Положив для |
краткости |
|
|
|
|
|
|
|
|
|
|
|
N |
|
|
|
|
|
N |
|
|
|
|
|
К Ф = |
|
2 |
«ргфр (0. Д (0 = |
2 |
|
K p rC rq H lq b (*). |
(12.30) |
|
Р = I |
|
|
|
Р . |
9 . |
г = 1 |
|
|
|
перепишем полученную формулу в виде |
|
|
|
|
|
|
B(t, |
2, ........ zN) = |
|
S Я,(0гг + Д(0. |
|
(12.31) |
|
|
|
|
|
|
|
r=I |
|
|
|
|
|
|
Из последней формулы следует, что функция 0 в данном случае линейна и не зависит от конкретного вида функции потерь. Однако она отличается от линейной функции, полученной в гл. 10, на личием дополнительного слагаемого A(t). Можно показать, что это слагаемое компенсирует систематическую ошибку, возникающую вследствие неравенства нулю математических ожиданий т ъ . . ., mN параметров сигнала Uх, . . UN.
Докажем, что при равных нулю математических ожиданиях т и . . ., mN случайных величин Ult . . ., UN функция (12. 31) сов падает с линейной функцией, полученной в гл. 10 для оптимальной
линейной системы. Прежде всего заметим, что в этом случае Д (0 = = 0 и формула (12. 31) имеет вид
9(*. .........zjV) = |
Xr(i)zr |
(12.32) |
|
r=l |
|
Остается доказать, что функции Xx(t), . . ., ^Л'(0> определяемые первой формулой (12.30), удовлетворяют в этом случае системе уравнений (10. 92). Так как кргпредставляют собой элементы матрицы, обратной по отношению к матрице В + С, то первая формула (12.30) дает решение системы линейных алгебраических уравнений
(Р= 1,--А^)
Заметим, что вторые начальные моменты ypq {р, q = 1, . . .,N) случайных величин Uх, . . ., UN представляют собой в данном слу чае элементы корреляционной матрицы случайного вектора U, обратной по отношению к матрице С в выражении (12.15) для нор мальной плотности вероятности. Поэтому, умножив уравнение (12.33) на угр и просуммировав по р, на основании известного соотношения между элементами обратных матриц получим
Е ( Е |
|
yrpbps + |
Е |
VrPcps ) К = |
Е уГР% (0 > |
s—1 \р=1 |
|
р=1 |
|
/ |
|
р=1 |
или |
|
|
|
|
|
|
|
Е ( Е уrpbPS+ |
1 |
к = |
£ |
УгрЬ (0- |
s=l |
\р=1 У |
|
|
р=1 |
|
(/•= 1 ,
Эти уравнения совпадают с уравнениями (10.92), что и требо валось доказать.
Полученные, результаты, дают основание во многих практических задачах ограничиться определением оптимальной линейной системы. А именно во всех случаях, когда функция потерь представляет собой неубывающую функцию абсолютной величины ошибки и есть основания считать совместное распределение входного и требуемого выходного сигналов нормальным, следует ограничиться определе нием оптимальной линейной системы по критерию минимума сред ней квадратической ошибки.
12.4. Оптимальная система при нелинейной зависимости входного сигнала от неизвестных параметров
Задача определения оптимальной системы в классе всех воз можных систем при линейной зависимости входного сигнала от случайных параметров Uu . . ., UN была рассмотрена в п. 12.2. При этом требуемый выходной сигнал мог зависеть от этих пара метров и нелинейно. Рассмотрим общую задачу, когда входной и
требуемый выходной сигналы являются произвольными, в общем случае нелинейными функциями параметров. Обозначив через U. вектор параметров с составляющими (/,, . . UN, полезный входной и требуемый выходной сигналы можно выразить следующими фор мулами:
5 (0 = ср (t, U), Гт (/) = ф(/, U), |
(12.34) |
где ср (t, U) и ф (/, U) — некоторые определенные функции указан ных аргументов. Помеху по-прежнему будем считать аддитивной и нормально распределенной. Тогда входной сигнал
* (0 = ф (/, U) + N(f). |
(12.35) |
Чтобы обобщить метод, изложенный в п. 12.2, на этот случай, дискретизируем векторный параметр U. Для этого разобьем об ласть его возможных значений на конечное число М малых эле ментов A«s и выберем в каждом элементе Aus точку us. Построим согласованные фильтры для сигналов
• ф(i, u.s) = ф(/, и!, ..., и%).
(s= 1____ М)
Весовая функция g (t,- т,- и) согласованного фильтра, соответ ствующего сигналу ф (t, и) при произвольном значении и, опреде ляется интегральным уравнением
t
т> и) dx = у (С, и). |
(12.36) |
Определив весовые функции согласованных фильтров, соответству ющих всем значениям us, выразим на основании формулы (12.35) выходной сигнал s-ro согласованного фильтра в фиксированный момент t в виде
t t
Zs = | g (t, t , us) X (t ) dx = |
J g (t, |
t , us) ф (t , |
(J)dx + |
to |
|
|
to |
|
|
+ Jt |
g (t, |
Г, |
11s) N |
(x) dx. |
(12.37) |
to |
|
|
|
|
|
(s = |
1 ........ M) |
■ |
|
Из этой формулы следует, что при любом значении и вектора U условное распределение выходных сигналов согласованных фильтров нормально, так как они представляют собой результаты линейных преобразований нормально распределенной случайной функции N(t). Поэтому для определения условного распределения выходных сиг налов согласованных фильтров при данном значении и вектора U достаточно найти их условные математические ожидания, дисперсии
и корреляционные моменты. Положив |
в выражении (12.37) U = и, |
о |
333 |
получим следующую формулу для условного математического ожи дания выходного сигнала s-ro согласованного фильтра:
ms (и) = М [Zs|и] = |
t |
|
J g (/, т, us) ср (т, и) dx. |
(12.38) |
(S= 1 , |
to |
|
• ■•, М) |
|
Дисперсии и корреляционные моменты выходных сигналов со гласованных фильтров выразятся формулой
kpq = Jt Jt |
g {t, T, |
UP) g (t, |
a, |
UP) Kn (t, 0 dx da. |
ta |
(p, |
q = 1, |
. . ., |
N) |
|
Но на основании формулы (12.36)
Jt Kn (t, a)g (t, a, up) da = 9(t, up). 10
Подставив это выражение в предыдущую формулу, получим t
kpq = \ g (t, x, up) ср (т, up) dx. |
(12.39) |
^0 |
|
(P. Я —1. •••. N)
Обозначив через Z вектор, составляющими которого являются выходные сигналы всех согласованных фильтров Z1, . . . . ZM, ус ловную плотность вероятности вектора Z при данном значении и вектора U можно выразить формулой
к (г\и) = |
— ■exp f■— -i-[z — т (и)]тК ~1 [г — m(u)]\, |
|
У (2л)" | к | |
l 1 |
> |
где m(u)— вектор с составляющими т^и), |
. . ., тм (и), а К — услов |
ная корреляционная матрица вектора Z, |
элементы которой |
опре |
деляются формулой |
(12.39). |
Как известно, |
для нахождения |
опти |
мальной системы достаточно знать условную плотность вероятности вектора U при данном значении z вектора Z. Пользуясь формулами теории вероятностей, для этой условной плотности вероятности по лучаем формулу
f (и) ехр { |
[г — т (и)]т К ~ 1 [г — т (и)]) |
f(u\z) = — --------------------------------------------------- |
1----, (12.40) |
j / (®) exp I— |
(z - m (®)]T /С-1 [z — tti (m)]| dv |
— 00 |
|
где f(u) — безусловная плотность вероятности вектора параметров сигнала U.
Найдем значения плотности вероятности f(u\z) в выбранных точках us (s = \, . . ., М). Имея в виду, что согласно (12.38) и (12.39)
tns (ur) = |
Jt g (t, |
т, |
us) ф (т, ur) dx = |
ksr |
|
|
to |
|
|
|
|
|
|
и обозначив через kjq элементы матрицы /С-1, |
можем написать |
[z — т (ur)]T K ~ l [z — m (иг)] = z^K - 'z — z JK~lm (иг) — |
|
— тп (ur)T K ~ xz + |
m (игУ K ~ xtn (u r) = |
|
м |
|
|
м |
|
|
|
|
= z TK~lZ — 2 Yi |
zpk~mq(ur) + S |
mp(ur) k~qmq(ur) = |
P,q=1 |
|
P .?=I |
|
|
|
|
|
M |
|
|
M |
|
|
|
= ZrK~lZ — 2 2 |
ZPKqkqr+ |
Yl |
kprkjqk{ |
|
|
P.q—1 |
P. <7=1 |
|
|
|
Ho |
Af |
|
|
|
|
|
|
|
kpqkqr — 6pf |
|
|
|
|
|
S |
|
|
|
|
|
</=I |
|
|
|
|
|
|
|
(p, r = 1 , . : . , M ) |
|
|
|
|
Следовательно, |
|
|
|
м |
м |
|
|
|
|
|
|
[z — m (ur)]T K~l [z — m {ur)\ = |
z ^ K ^ z — 2 Yi zP6Pf + £ V |
sp/- = |
|
|
|
|
p = |
i |
p = i |
|
|
= z TK~1z — 2zr -f |
|
|
|
|
Подставив это выражение в формулу |
(12.40), |
сократив множи |
тель ехр | ---- 2т/е- 1.г| и |
заменив |
интеграл |
суммой, получим |
сле |
дующее приближенное выражение для условной вероятности зна чения и г вектора U:
/ |
( u r) exp \ z r ----- ^ krr\ Ди ' |
|
/(«Ч«)Д»Г------Af----------L |
------1-------- • |
(12‘41) |
^ |
f (u s) ехр | zs ----- i - jfess\ Д а5 |
|
S=1 |
|
|
|
Эта формула определяет условную вероятность попадания слу чайного вектора U в элемент t±.ur с точностью до малых второго порядка относительно тах{Дяг).
Вычислим условное математическое ожидание функции потерь при данном вектрре Z. Для этого умножим значение I (Y , ф (/, иг))