Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Лихарев В.А. Цифровые методы и устройства в радиолокации

.pdf
Скачиваний:
47
Добавлен:
25.10.2023
Размер:
14.41 Mб
Скачать

имеющая предельное ^-распределение с г 1) [п1) степенями свободы. Как обычно гипотеза На принимает­ ся, если с доверительной вероятностью р %2 (1р) ^

^%Z<\I4P).

И, наконец, перейдем к вопросу проверки гипотезы о порядке связности цепа Маркова. При этом ограничим­ ся однородными цепями. Необходимо проверить две ги­

потезы:

 

 

# о — цепь Маркова

одиосвязная, v = 1,

# i — цепь Маркова

двухсвязная, v = 2.

Наряду с

оценками переходных вероятностей ра? одно­

связной цепи

(1.13), нам потребуются оценки переходных

вероятностей

ра^ двухсвязной цепи

 

 

 

 

 

Л з т = ^ з т

Ithsr

 

 

(I - 1 1 4 )

В

то же

время

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а = 1

 

/

а = 1

/=1

 

 

 

 

 

Для

проверки

гипотезы

Н0

вычисляются статистики

 

 

 

 

 

г

г

 

 

 

 

 

 

 

 

 

 

 

 

« = 2

Я

^ І

Р

^

-

Р М

*

-

 

(1Л16>

 

 

 

 

 

a = l

7=1

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7=1

 

 

 

 

 

 

 

 

 

Если

гипотеза

На

верна,

имеет

предельное

/^-рас­

пределение

с 1)а

степенями

свободы. При

проверке

совместной

гипотезы

Н0:

р а

? = р

для

всех

а,

(3, у =

=

1, г,

следует

вычислить

статистику

 

 

 

 

 

х 2 = Е x j

= 2

2

2 ( Р а н - К У І К *

 

(1Л!7>

 

 

 

P = l

 

a = l

(3=1

7=1

 

 

 

 

 

 

 

имеющую

предельное

^-распределение

с

г (г—I)2

степе­

нями свободы.

 

 

 

 

 

 

 

 

 

 

 

 

60

Аналогичным образом образуются статистики для проверки гипотезы о том, что цепь Маркова имеет связ­ ность v > l при альтернативе v = l . Соотношения (1.116) и (1.117) легко обобщаются на неоднородные цепи Мар­

кова. В

общем случае при проверке гипотез:

# 0 ( v = i>) и

Hi(v—u),

где

и < « , статистика %2 имеет

предельное

Х2 -распределение с (/'"—rv )<(r—1) степенями свободы.

Пример. Пусть производится бинарное квантование стационар­

ных гауссовых

процессов с коэффициентами' корреляции R(i) —

=е»р{—.а|х|}, іУ?(т)=ехр{ — ат 2 } . Переходные верятности цепи Мар­ кова определяются по формулам (1.97), (1.98), (1.100). Проверка

гипотез о связности была

произведена

методом Монте-Карло

с использованием критерия х2 -

Результаты

моделирования приведены

на рис. 1.2. Следует отметить, что оценка связности v, полученная

Г "

г Г

• і

і

Рис. 1.2. Оценка связности цепи Маркова.

выборочным методом, хорошо согласуется с формулой (1.83). Эта оценка зависит от доверительной вероятности н вида коэффициента корреляции процесса R(x). Приведенные результаты будут справед­

ливы для стационарных случайных процессов, получаемых-'из нор­ мального нелинейным безынерционным преобразованием: равномер­ ного, релеевского, логарифмически нормального. У всех названных процессов коэффициент корреляции близок к коэффициенту корреля­ ции исходного нормального процесса [152].

1.7.4. Информационные критерии в задачах аппрокси­ мации распределений. Вопросы аппроксимации распреде­ лений последовательностей дискретных переменных на выходе аналого-цифровых преобразователей конечными

цепями Маркова были рассмотрены

в п. 1.7.3 в

рамках

статистических критериев согласия

(критерий

Одна­

ко такое решение задачи может

оказаться

не до­

 

 

бі

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

каковы потери

за

счет дискретизации

и квантования?

как

выбрать

оптимальные параметры

квантизатора?

чему

равны потерн

при аппроксимации

точных распре­

делений приближенными? Ответы на эти вопросы могут быть получены на основе подхода, связанного с исполь­

зованием

информационных

мер в задачах

статистиче­

ских решений

[18—21].

 

 

 

 

Общая

постановка

задачи

и ее решение таковы. Зада-

 

 

 

 

 

->

->

 

на многомерная плотность вероятности w {и,

Я)

процесса

u[t,X(t)\

на

входе

аналого-цифрового

пр:образователя.

Выборочное пространство U,

u(t)aU,

непрерывно. Зада-

 

 

 

—•

 

 

 

 

на решающая

функция у (и),

минимизирующая

средний

—>

риск г решения при заданной функции штрафов с (Я, у).

Структура решающей функции у(ы) определяется распре-

дел?нием w (и,

Я)

и функцией

штрафов

с (Я, у).

Средний

риск

решения

равен

[6]

 

 

 

 

 

 

 

 

 

 

 

 

->

->

-*•-+-»•-+

=

 

 

 

 

 

 

 

(Я, у (и)) w (и,

Я) dudX

 

 

 

 

 

и Ф

 

 

 

 

 

 

 

 

 

 

=

f

j

с (Я,у (и)) w (и | Я) w (X)dudX,

 

(1.118)

 

 

 

й Ф

 

 

 

 

 

 

 

 

 

где

—>

 

 

->

 

 

 

ПраВДОПОДОбИЯ Парамет­

Я ( і ) с Ф ,

ш (ы | Я) функция

 

ру

- у

 

 

 

 

 

 

 

 

 

 

-V

ров

Я, ш (Я) априорная

плотность

вероятности

для Я.

Произведем теперь редукцию выборочного пространства

U =±>U с помощью

операций аналого-цифрового преобразо­

вания

(дискретизации и квантования). Теперь

мы имеем

 

 

 

 

 

 

 

 

 

 

 

->

->

 

дискретное многомерное

распределение

Р(и\Х),

которое

будем

считать

дискретной

цепью

Маркова

с г

состоя­

ниями

и связностью v. Средний риск,

связанный с распре-

делением Р(и\Х),

 

равен

 

 

 

 

 

 

 

 

 

7=-- J

j

с (Я, у (u))P(u| Я)о»(A) d u d l .

(1.119)

62

В

работах

[18,

19] показано, что оценка

изменения

байесовского риска при редукции выборочного

простран­

ства дается

соотношением

 

 

 

 

 

 

 

Vr(c2)2I{P:w)

 

 

< 7 — г <

Vr{c"-)2I{w:P),

(1.120)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7(с2 )

== j

j V ( l , Y (и))Р

(u\X)w(X)d!ldX,

(1.121)

 

 

 

 

и ф

 

 

 

 

 

 

 

 

 

 

'(w:Р)

=

f

f wСи,

X) In

(s(й, \\)

dlidl.

(1.122)

 

 

 

 

 

 

J

J

 

Р

 

 

 

 

 

 

 

 

 

и Ф

 

 

 

 

 

 

 

Соотношение

 

(1.122)

имеет

смысл

„информационного рас-

 

 

 

 

 

 

 

 

 

 

->• -+

~

 

стояния"

между

распределениями

w(u,X)

и Р(и, Я).1 ' „Ин­

формационное расстояние"

I(w:P)

необходимо

выразить

через

информационную меру Шеннона [20]

 

 

 

 

 

 

 

 

I(w:P)^J(w)

 

— J{P),

 

(1.123)

где

 

 

 

 

 

 

 

 

 

 

—У - У

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7 (ш) =

(

( w (!) ш (и11)1п

 

та(Х)

rfudt;

(1.124)

 

 

 

 

 

v

Ф

 

 

 

У

'

 

 

 

 

 

/ (Я) =

У

f w (І) Р (и | Я) InЯ

(

d l

(1.125)

 

 

 

 

 

^ J J '

 

 

Р ('ц( В))

 

 

 

 

 

 

 

 

V Ф

 

 

 

 

 

 

 

В свою очередь, J{w)=-H

 

(Я) — Я

(Я j «),

 

 

/(/?) = Я (Я) — Я

(Я | uj, Я (Я) = — j

w (Я) In ти (Я) d l

априорная

 

энтропия,

 

 

 

 

 

 

 

*> Интеграл

I(w:P)

существует,

если

вероятностные

меры

->

 

 

. ->

 

 

 

 

 

 

 

 

 

 

au(u, >,) и

Р(н ,

Я,)

абсолютно непрерывны.

Это требование

выпол-

 

 

 

 

 

 

 

 

 

 

->

 

 

 

няется

для

дискретного

распределения Р(и,

X),

кусочно-постоянного

на интервалах

 

Ди,-,

i = l , г, являющегося

достаточной

статистикой

 

 

 

 

 

 

 

- V

 

 

 

 

 

 

 

цля распределения Р(и, X).

-

63

 

 

Н (я I и) = J f w{и, Я) In w (Я I и) dudX

иФ

апостериорная энтропия. Значит,

->

I(w:P)^H

(Х\и) — Н{Х\и),

(1.126)

т. е. больше или равно разности апостериорных энтропии. Пусть колебание на входе приемника имеет вид

u(t) =

bS{t,

X(t),a{t))-\-n{t)

+ p(t),

(1.127)

где її (t) — внутренние

шумы приемника; р

(t) — внешние

—>

->

 

 

 

помехи; 5 (t, X (t), a (t)) — сигнал, заданная функция времени

-> ->

и параметров; Я(/) вектор измеряемых параметров; a (t)—

вектор неизменяемых параметров; 0 — дискретный слу­ чайный параметр, принимающий значения: 0 — отсутст­ вие сигнала, 1—наличие сигнала от цели. Задача про­ верки двух гипотез сводится к оценке параметра 0. При

1 1 этом функция штрафов обычно задается в виде

(

1 при 6

=

0 =

90 ,

Y

= 1 .

 

 

 

c ( 6 , y ( u ) ) = J

I

при 0 = 1 =

6,,

Y = 0,

(1.128)

[

0

при 0

=

Т =

1

или 8 = Y =

°-

Тогда оценка изменения среднего риска (1.120) принимает вид [21]

— ]Г2г1{Р:ш)

< r — r<V2rI

(w:P) .

(1.129)

Спомощью соотношений (1.120) — (1.125) и (1.129)

можно оценить потери в среднем риске при выбранной

~>. ->

нами аппроксимации распределения Р (и | Я). Если эти потери окажутся неприемлемыми, следует усложнить аппроксимацию, т. е. увеличить число состояний и связ­ ность цепи Маркова, путем увеличения разрядности пре­ образования и уменьшения периода дискретизации.

При проверке двух гипотез

(обнаружение сигнала)

средний риск равен

 

r=pM+qF,

(1.130)

64

где М = \—D — вероятность

пропуска цели; D — вероят­

ность правильного обнаружения; F — вероятность

лож­

ной тревоги; р — априорная

вероятность

наличия

цели;

q — априорная

вероятность

отсутствия цели (p-\-q = l).

-При измерении скалярного параметра X, квадратич­

ной функции штрафов

 

 

 

 

с(Я,Я) =

( Я —If

 

( і . і з і у

и оптимальной

решающей функции Y 0 ? t ( " ) '

соответствую­

щей оценке по центру тяжести апостериорного распреде­ ления

 

X*

( " ) = V Й

=

 

J

« » № ) *da,

 

(1-і 32)

 

 

 

 

 

 

 

 

 

—со

 

 

 

 

 

 

 

 

 

 

средний риск

равен

 

дисперсии

а2

[Я] оценки

параметра

Я.

Из уравнения

(1.120)

 

следует,

что

минимизация

потерь

в

среднем риске

эквивалентна минимизации

„информационно-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

->

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

->

->

 

~ ->

го расстояния"

между распределениями

w(u\X) и Р(и\Х),

что, в свою очередь, эквивалентно

максимизации

инфор­

мационной меры Шеннона (1.123)—(1.125).

 

 

 

 

 

Максимизируя

(1.125),

можно

 

решить

 

задачу

об

опти­

мальной расстановке

 

порогов квантования

аналого-цифро­

вого

преобразователя.

В

самом

деле,

пусть

область

определения

переменной и,

равная

и

— и

= £ / , ,

раз­

бита

на г подынтервалов

 

с помощью r, =

r — 1

порогов

квантования

0

 

01

а2

 

«

0 Г і

)

т

- Необходимо найти

век­

ы =

 

 

 

 

тор «0 , максимизирующий

 

функцию

(1.125),

т. е. решить

систему уравнений

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/(Р(н0 ))== max/(/>(«)).

 

 

 

(1.133)

 

 

 

 

 

 

 

 

 

->

 

 

 

 

 

 

 

 

 

 

и

Такая задача может быть решена методом неопределен­ ных множителей Лагранжа. Условие связи параметров таково:

г,

S 4 » = ^ . Auoi=U0i—«„;_,.

(1.134)

Составим функцию Лагранжа

Ф & Ф ) = / ( Р Й , ) ) + * Д и . < - u } j .

(1.135)

5—1410

65

Точка максимума определится

из системы /', - f - 1

уравне­

ний

с г, -|- 1 неизвестными Дн0 1 , Д«0 2

Ды

 

 

 

 

V

Д«0 1 _ £ У, = 0

 

 

 

 

 

 

йФ(а„ф)

 

0.

 

 

(1.136)

 

 

 

диа

 

 

 

 

 

 

 

 

 

 

 

 

Когда

неизвестны

априорные

вероятности

гипотез

/У и И0, равные соответственно р

и q,

обычно

применяет­

сяІ

критерий Неймана — Пирсона,

а использование соот­

ношений

(1.133) — (1.136) при

оптимизации

параметроз

аналого-цифрового преобразования невозможно. В этом случае потерн информации при проверке гипотез # i И Я 0

за счет перехода от распределения гс'„(м) к распределе-

нию Р(и) оцениваются величинами

Y, = /_>(! :2;я)//^(1 : 2; л) при Я 0 ,

 

 

 

 

 

 

(1.137)

 

Т8 =

/_Д2: 1;я)//_Д2:

при Я „

 

где

 

 

 

 

 

 

 

 

^ ( 1 : 2 ; « . ) = Г ш . Л « | 8 „ ) 1 п ^ Ш 1 ^ ;

 

 

 

J

 

ш„(и|в,)

 

 

 

 

 

 

 

 

(1.138)

 

M 2 : l ; / z ) = U n g | 6 , ) l n а , а ( ! 1 " ' )

^

 

 

"

J

 

Ш я

(И | 60 )

 

 

— информационные числа Кульбака—Лейблера [20].

 

Величины /_»(1: 2; я) и / _ (2 : 1;«)

имеют

аналогичную

 

и

 

и

 

 

 

 

форму

записи.

 

 

 

 

_^

 

При

независимых

выборках,

т.

е. при

wn(u\Q0)

=

п

 

«

 

 

 

 

= П w, (Ui 1 в„)

и ш„ (а 1 в,) = Ц

®i ("*

I б>) п о

т е Р и №

ф о р "

i = l

 

 

i = I

 

 

 

 

66

мации

(1.137)

имеют вид:

 

 

 

 

 

 

 

I / £ 8 i l

= Y 1 = i ( u ) / / i ( u )

=

7 t t ( l :2)//Jl

:2) при Я 0 ,

(1.139)

 

 

 

 

 

 

 

 

 

 

U

 

 

 

 

1 / £ ' г 1

1 = у „ =

/1(и)/л(«) =

/„ (2: \)/IJ2:

1) при Я , ,

(1.140)

где

 

/ u ( l : 2 ) =

J » l (

 

 

 

 

 

 

 

 

 

« | f l . ) l n 5 | | ^ r f « ;

 

 

(1.141)

 

 

/ u ( 2 : l ) = j » 1 ( u | e i ) l n = i | l ^ d « .

 

 

 

а ЕІЛ,Е'—асимптотическая

 

 

относительная

эффектив­

ность

[см. (2.169)];

/ _ (1:2) и / _ (2:1)

имеют

аналогич-

ную с

 

 

 

 

и

записи

с

заменой

ш,(и|б0 ) на

(1.141) форму

P(u|G0 ),

В У

, (

И

| 0 , )

на

P(u|9,)

и интегралов

на

суммы.

При

использовании

 

критерия

Вальда

соотношения

(1.139) , (1.140) остаются в силе; под величинами /г (и) и п (и) теперь следует понимать средние объемы выборок.

Объемы выборок п (и) и п(и), определяемые из (1.139), обеспечивают одни и те же вероятности ложной тревоги F и пропуска M^>F при использовании распределений

Р(и) и

Wi(u); при

M^F

справедливо

соотношение

(1.140) .

Соотношения

(1.139)

и (1.140)

получены для

больших

выборок, т. е. при П

*0О.

 

Теперь процедура

оптимизации параметров аналого-

цифрового преобразования сводится к минимизации объ­

ема

выборки /і (її) при заданной

вероятности

ложной

тревоги F, т. е. при /r = const

и М~^>Р к

максимизации

функции

 

 

 

 

/

(1 :2|и0 ) = т а х У І Р ( ц | ц ,

6„)ln р

& 1

^ .

(1.142)

и

Определение максимума (1.142) можно выполнить ме­ тодом неопределенных множителей Лагранжа, как при решении уравнения (1.133) при условиях связи (1.134).

5*

67

В общем случае при зависимых выборках, оптимизи­ руя параметры 'преобразования, в качестве информацион­ ной меры можно использовать дивергенцию Кульбака

[ 2 2 ]

£ > ( 1 , 2 ) = / ( 1 : 2 ) + / ( 2 : 1 ) ,

( 1 . 1 4 3 )

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

->

-> —>

D(\, 2) = V [РСи\80).-Р(Ц|

0,)] l n ^ I M . .

( 1 . 1 4 4 )

При гауссовых статистиках распределений, соответствую­ щих гипотезам Н0 и Ни дивергенция равна отношению энергии сигнала к спектральной плотности шума [20].

—У

Поэтому выбор вектора и0 можно производить, максимизируя D(\,2\u) по и методом неопределенных множите­ лей Лапранжа:

 

Z ) ( l , 2 | u 0

) = m a x £ ) ( l , 2 | t t )

 

 

( 1 . 1 4 5 )

 

 

 

 

—V

 

 

 

 

 

 

 

 

и

 

 

 

 

с учетом

( 1 . 1 3 4 ) .

 

 

 

 

-+

 

 

 

 

 

 

 

 

 

При

оценке

векторного

параметра

Я = (Я,,

Я/)т и

квадратичной функции

штрафов

 

 

 

 

->

з

'

 

 

 

 

 

 

с(Я,

Я) = 2

Rii(*i

li){K

h)'

 

 

 

;,/=i

 

 

 

 

 

где [ R i j ]

— некоторая

симметричная неособенная

матрица,

выражение для

среднего

риска имеет

вид [ 5 5 , 1 7 1 ]

г >

det1 / 2 tf ехр

(X)

ехр | -

- f

[J(w(Z,

Я))]}-

 

 

 

 

 

 

 

 

 

 

( 1 . 1 4 6 )

Выражение ( 1 . 1 4 6 ) тем точнее, чем выше точность измерения параметров X, т. е. когда совместная плотность

w (X, X) является многомерным нормальным распределением.

Аналогично записывается выражение риска для ре­ дуцированного выборочного пространства. Относитель­ ное изменение риска за счет редукции выборочного про68

странства,

равное относительной

эффективности (3.38),

с учетом

(1.146)

принимает вид

 

 

«

ехр | -

- f [•/ И

- •/(/>)] }• - ехр / А

(Я |") -

 

-

Я (11 u)] J «

ехр { - - f

/ (ш : Р) } •

(1.147)

Выражение (1.147) удобнее при расчетах, чем формула (1.120).

Итак, мы рассмотрели задачи оценки увеличения байесовского риска за счет дискретизации и квантова­ ния, когда дискретное распределение Р{и) точно соот­ ветствует истинному распределению дискретных пере-

менных и, а также проблему оптимизации параметров квантизатора. Не менее важной является задача аппро­ ксимации многосвязных цепей Маркова цепями невысо­ кой связности, т. е. вопросы замены точного распреде-

ления Р(и) приближенным

Р'{и). В работе [77] найдено,

что

относительное

увеличение

среднего

риска,

вызван-

 

 

 

 

->•

 

 

 

 

ное

заменой

Р(и)

на

Р'(и),

ограничено

сверху

значе­

ниями

 

 

 

 

 

 

 

 

V2I(Р':

Р) У

21 (Р':

Р)

+V2!(Р

: Р').

(1.147а)

С помощью (1.147а) можно решать задачи аппроксима­ ции многосвязных и бесконечно усложняющихся цепей

Маркова Р(и) цепями невысокой связности Р' (и). При

этом, осуществив замену

P(u)z$P'(u),

оценим относи­

тельное увеличение риска

по (1.147а).

Если потери ока­

жутся приемлемыми, то

решающая функция может син­

тезироваться далее для

распределения Р'(и).

Соседние файлы в папке книги из ГПНТБ