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

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

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

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

Добавлен: 16.10.2024

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

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

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

функции потерь I (Y,

Yr) при и = иг на / (ur\Z) Аиг и просумми­

руем по г. В результате получим

 

AI

M[l{Y,

YT)\Z]=--

 

 

 

 

^ / ( Y ,

ч|) (/, a ') ) f

(«О exp U r ------ЛГг | Амг

 

= — ------м------------------------------------------•

(12-42)

2

/ (« 0 exp

f z r — - i - Arrj Диг

 

Определив Y как функцию выходных сигналов согласованных фильтров Zr минимизацией числителя в последней формуле, найдем оптимальное нелинейное преобразование выходных сигналов согла­ сованных фильтров. Таким образом, заменив непрерывное распре­ деление случайного вектора U дискретным с достаточно мелким разбиением области его возможных значений, как и в п. 12.2, получим оптимальную систему, состоящую из ряда согласованных фильтров и оптимального безынерционного нелинейного преобразователя их выходных сигналов.

Для нахождения точного решения задачи необходимо перейти к пределу при max J Aurj —>0. При этом число согласованных фильт­ ров будет возрастать и мы получим бесконечное множество (непре­ рывный спектр) согласованных фильтров.

Переходя к пределу в формуле (12.42) и принимая во внимание, что условное математическое ожидание функции потерь при данном векторе Z переходит при этом в условное математическое ожидание

функции потерь при данном входном сигнале Х(т),

получаем

 

 

 

00

 

 

 

м [/(Г, У,)|Х] =

х J l(Y, .f (/ u))f{ti)

х

 

 

X exp j J g (t,

T,

u)X(x)dx---- 4_P(tt)|rf«,

(12.43)

где

 

 

 

 

 

х = J

f (a) exp IJ g (/,

т,

и) X (x) dr — - L p (u)\du

, (12.44)

—03

\10

 

J

-

 

 

P (и) = Jt

g (t, T, и) cp (t , u) dx.

 

(12.45)

Для определения выходного сигнала оптимальной системы ин­ теграл в выражении (12.43) следует рассматривать как функцию переменной Y и найти то значение Y, которое обеспечивает минимум этого интеграла. Величина % не зависит от Y, и поэтому ее можно считать постоянной при минимизации M[l(Y, УТ)|Х ].

336


Формулы (12.41) и (12.37) при переходе к пределу дают следу­ ющее выражение для условной плотности вероятности случайного вектора U при данной реализации х (т) входного сигнала X (т):

/(я|д;) = х/(и)ехр jjg (f,

т, u)x(x)dx---- ^-[5(и)|. (12.46)

Wo

>

Итак, при нелинейной зависимости входного сигнала от пара­ метров оптимальная система состоит из бесконечного множества согласованных фильтров, соответствующих всем возможным зна­ чениям случайного вектора параметров U, и функционального преобразователя, осуществляющего безынерционное нелинейное пре­ образование выходных сигналов согласованных фильтров путем минимизации интеграла в (12.43).

Можно показать, что при линейной зависимости входного сиг­ нала от параметров вида

ср {I,

U) = ср (/, и ъ . . ., UN) = £

и г(рг (/)

 

Г=1

 

интеграл в формуле (12.43) принимает вид (12.46).

На основании

полученных результатов, можно заключить, что

определение оптимальной системы сводится к выполнению следу­ ющих операций:

решение интегрального уравнения (12.36); вычисление функции {5(и) по формуле (12.45);

определение выходного сигнала оптимальной системы, точнее нахождение выражения ее выходного сигнала через входной сигнал

(оператора

оптимальной системы) путем минимизации интеграла

в формуле

(12.43).

Изложенный метод дает оптимальную систему не только в сред­ нем, но и для каждой конкретной реализации входного сигнала Х(т), так как основан на минимизации условного математического ожидания функции потерь при данном входном сигнале.

12.5. Определение оптимальной системы по критерию минимума средней квадратической ошибки

При определении оптимальной системы по критерию минимума средней квадратической ошибки функция потерь выражается фор­ мулой

l(Y, Yr) = ( Y - Y T)\

Подставив это выражение в формулу (12.43), получим

СО

M [ l ( Y , Y r) \Х] =

х J [7 - ф ( г , u)ff(u) х

 

 

— со

 

Xexp j J g(t, т,

u)X(x)dx ---- ^-|3(#)jdtt.

(12.47)

22 В. С. Пугачев

337


Приравнивая нулю производную этого интеграла по Y, полу­ чаем уравнение для определения выходного сигнала оптимальной системы У* =

СО

2% J [Y* — ф (t, и)] / (и) X

СО

X exp (иJ. g (t, т, и) X (т) clx— -у р (»)

Решая это уравнение, находим

СО

П* = К* (f) = х j ф (*, и) /(») X

СО

X exp I J g (t,

т, a) X(x)dx---- Р (и) | du

(12.48)

и»

1

 

Таким образом, при использовании критерия минимума средней квадратической ошибки метод, изложенный в предыдущем параграфе, дает явное выражение для оператора оптимальной системы.

Правая часть равенства (12.48) на основании формулы (12.46) для условной плотности вероятности случайного вектора U при данном входном сигнале Х(т) представляет собой условное матема­ тическое ожидание функции ф(/, U), т. е. требуемого выходного сигнала YT(t) при данном входном сигнале Х(т). Следовательно,

формула (12.48) может быть

переписана в следующем виде:

 

У*(0

= М[Ут (0|Х].

(12.49)

Таким образом, оптимальной системой среди всех возможных систем по критерию минимума средней квадратической ошибки является система, которая дает на выходе условное математическое ожидание требуемого выходного сигнала при данном входном сиг­ нале, т. е. апостериорное математическое ожидание требуемого выходного сигнала. В каждом конкретном случае это условное математическое ожидание вычисляют по формуле (12.48).

Если требуемый выходной сигнал YT(i) является линейной функцией параметров

П (0 = !>(*,

^ М О ,

 

то формула (12. 49) принимает вид

Г— 1

 

 

 

Y*(t)= S M[Ur\X]^r{l).

(12.50)

Л=1

 

 

Но величины

 

 

U* = M[Ur\X]

(12.51)

(r= 1, . . . .

А/)

 

согласно той же формуле (12.49) представляют собой оптимальные

338


оценки параметров U ъ . . ., UN по критерию минимума средней квадратической ошибки. Следовательно, для нахождения оценки сигнала, представляющего собой линейную функцию параметров U .......... UN, по критерию минимума средней квадратической ошибки

достаточно найти

оценки U\, . . ., U*N параметров

Uv . . ., UN

и заменить этими

оценками параметры U ........... UN в выражении

требуемого выходного сигнала:

 

 

r ( 0 = S t / > | V(0-

(12.52)

 

Г=1

 

В заключение рассмотрим частный случай, когда не только требуемый выходной сигнал, но и входной сигнал X(t) линейно зависят от нормально распределенных параметров U ъ . . ., UN. В этом случае

Ф(Л и ) = Ъ и л Л 1 ) -

г=1

Заметив, что величины р^, . . ., |%, введенные в п. 12.3, пред­ ставляют собой реализации случайных величин U\, . . ., U*N, опре­

деляемых формулой (12. 51),

при данной реализации х(т) входного

сигнала Х(т), из формулы (12.

22)

получаем следующую систему линей­

ных алгебраических

уравнений

для определения

U\.......... U‘N\

£ ( » „ - К ,)

у; = ■ ? ,« +

(12.63)

< 7 = 1

4

 

< 7 = 1

 

(Г= 1 ____ N)

Таким образом, при линейной зависимости входного и требуемого выходного сигналов от нормально распределенных параметров и критерия минимума средней квадратической ошибки оператор оп­ тимальной системы определяют чисто алгебраически.

12.6. Оптимальные системы обнаружения сигналов

Изложенный в п. 12.4 общий метод определения оптимальной системы в классе всех возможных систем можно, в частности, при­ менить для синтеза оптимальных систем обнаружения и распозна­ вания сигналов. Если перенумеровать сигналы разных классов, подлежащих распознаванию, и принять номер сигнала за дополни­ тельный параметр входного сигнала и за требуемый выходной сигнал, то задача распознавания сведется к задаче, рассмотренной в п. 12. 4, причем один из параметров сигнала будет иметь только конечное число возможных значений. Определив соответствующим образом функцию потерь, можно применить метод, изложенный в п. 12. 4 для определения оптимальной системы обнаружения.

Рассмотрим задачу обнаружения сигнала. В этом случае рас­ познаванию подлежат сигналы двух классов: сигнал, представляю-

22*

339


щий собой чистый шум, и сигнал, представляющий собой сумму полезного сигнала и помехи. Если решать задачу обнаружения по критерию минимума вероятности ошибки, то, как было показано в п. 10.2, функция потерь равна нулю, когда выходной сигнал сис­ темы равен требуемому выходному сигналу, и единице, когда вход­ ной сигнал системы не равен требуемому выходному сигналу:

(0

при У =

Ут,

(12.54)

ЦУ, Гт)~ 1 1

при Y ф У г.

Если определить функцию ср (/, и) так,

чтобы она

обращалась

в тождественный нуль при и = 0, то можно считать, что случай отсутствия полезного сигнала соответствует нулевому значению параметра и , а случай наличия сигнала соответствует любому отлич­ ному от нуля значению параметра и. При этом требуемым выходным сигналом является номер класса полезного сигнала. Условимся приписывать номер 0 сигналу, представляющему собой чистый шум, и номер 1 сигналу, содержащему полезный сигнал. Тогда требуе­

мый выходной сигнал определится

формулой

_

[ 0

при

U = 0,

т

\ 1

при

(12.55)

U =j=0.

Пусть р — вероятность того, что полезный сигнал есть в при­ нимаемом сигнале, q — 1 р — вероятность отсутствия полез­ ного сигнала, h (а) — условная плотность вероятности случайного вектора U при наличии полезного сигнала. Тогда плотность вероят­ ности случайного вектора U

f(u) — q8(a)-\-ph(u).

(12.56)

Имея в виду что ср (t, 0) = 0 и, следовательно, g (t, т,

0) = 0,

и принимая во внимание, что условное математическое ожидание функции потерь вида (12.54) представляет собой вероятность ошибки

системы,

перепишем формулу

(12.43)

в

виде:

 

 

р0Ш(Х,

Y) = M [/(У,

У7)\Х] =

 

 

 

со

f

t

 

ql (Y,

0)+ р /(У ,

1) J

/г(и) exp

Jg (f, т, и)Х{т)йх —

 

 

 

— со

 

П о

 

 

 

2

Р (и) | da .

 

 

 

 

 

 

 

Вводя

обозначение

 

 

 

 

 

 

со

 

 

 

 

Q (X) =

j h (и) exp

J g(t, X, u)X{x)dx — -±-$ (u)\da,

 

— со

U

 

 

(12.57)

 

 

 

 

 

340