Файл: Пугачев, В. С. Основы статистической теории автоматических систем.pdf

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

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

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

Добавлен: 16.10.2024

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

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

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

Поэтому для неубывающей функции 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 для оптимальной

331


линейной системы. Прежде всего заметим, что в этом случае Д (0 = = 0 и формула (12. 31) имеет вид

9(*. .........zjV) =

Xr(i)zr

(12.32)

 

r=l

 

Остается доказать, что функции Xx(t), . . ., ^Л'(0> определяемые первой формулой (12.30), удовлетворяют в этом случае системе уравнений (10. 92). Так как кргпредставляют собой элементы матрицы, обратной по отношению к матрице В + С, то первая формула (12.30) дает решение системы линейных алгебраических уравнений

Е

+

О 2-33)

s=1

 

 

(Р= 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. При этом требуемый выходной сигнал мог зависеть от этих пара­ метров и нелинейно. Рассмотрим общую задачу, когда входной и

332


требуемый выходной сигналы являются произвольными, в общем случае нелинейными функциями параметров. Обозначив через 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.

334


Найдем значения плотности вероятности 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 , ф (/, иг))

335