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

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

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

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

йРч ПЗУ

Раз

її

'02

 

PC, ПЗУ,

a(><)=ailk)

~п(к)

РС2

ПЗУ,

ДРЧ

 

ИЛИ

 

PC,

ПЗУ,

 

РСИ

ПЗУ.

 

дш

x(k-t)

 

 

Рис. 1.1. Преобразователь равномерного распределения (а) и генера­ тор дискретной марковской последовательности (б).

Разрядность т датчика равномерных чисел и преоб­ разователей зависит от точности задания переходных ве­ роятностей в виде соотношения

Так, при т = 8 максимальная абсолютная погрешность

| Драр j m a x ^ « 0,004.

Соответственно относительная погреш­

ность есть | Дра р

| т а х / р а р .

50

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

Наряду с марковской аппроксимацией дискретных пе­ ременных возможно многомерное нормальное приближе­ ние. В частности, в монографии (38] приведено выраже.- иие для двумерного биномиального распределения, асим­ птотически нормального. Аналогично можно записать выражения для произвольных распределений последова­ тельностей дискретных переменных при многоуровневом квантовании в виде многомерного нормального распреде­ ления. Основная трудность такой аппроксимации состоит в определении матрицы ковариацнй дискретных пере­ менных.

Наряду с рассмотренными возможна иная аппрокси­ мация многомерного дискретного распределения для последовательности статистически связанных нулей и единиц:

 

 

 

 

р{и)

=

Р, («)f (").

(1.87)

где

 

p 1 ( » ) = n ^ t ( 1 - ^ ) ( I

 

 

 

 

 

 

 

di

— индикатор

состояния

(1.8),

т. е. di =

щ — 0,1;

 

 

 

-»•

 

 

 

 

 

 

 

 

 

f (її) =

1 + £

ГнУгУз

+

£ ГцьУгУіУк

+ • • • 5

(1-88)

 

 

 

i<j

 

 

i<j<k

 

 

 

 

 

У г =

(Ui

Pi)fVpi(l

—Pi)

 

 

нормоцентрированная

переменная;

 

 

 

 

ГЦ =

' " і і [УіУі]

• Ггjft =

m „ , [УгУіУк].

• • • .

 

парные, тройные

и т.. д. корреляции

соответственно.

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

(1.87)

тождественно

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

(1.20)

при г = 2.

 

 

 

 

 

 

 

С

помощью

(1.87) и ((1.88)

удобно рассчитывать

пе­

реходные вероятности при бинарном квантовании негаус­

совых

процессов.

 

 

 

 

 

 

 

1.7.2. Определение переходных вероятностей и началь­

ного вектора цепи Маркова. В общем случае квантования

4*

 

 

 

 

 

 

 

51

стационарного случайного процесса u(t) компоненты 'начального вектора р0 и вектора финальных вероятностей РФ определяются выражением

 

 

 

Р 0 . = Р а =

J t B , ( " ) r f " .

 

 

 

(1-89)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

->

 

 

где

р0 а —компонента начального

 

вектора

р0 =

0 1 , р 0 3

р 0 г) т ; Р0 компонента

финального

вектора

р ф = (р1 ,

р . , , . . . ,

р,-)т

односвязной

цепи;

 

ау, («) — одномерная

плотность

вероятности процесса

 

u(t).

 

 

 

 

 

 

 

 

 

 

 

 

Для двухсвязной

цепи:

 

 

 

 

 

 

 

 

 

 

 

 

 

Ро(«?)=/'(«р)=

 

J

J Ш Л " , . ua)dutdut,

 

(1.90)

где

р 0 ( а р ) —компонента

 

начального

вектора

р0 =

0 1 1 ,

Poi2> •••• Рогг)т ;

р( а р ) —компонента

 

финального вектора двух­

связной цепи Маркова; Даа , Да? —интервалы

квантова­

ния, а, р=:1, г ; £ Д« а =

u f f l a x

— u m I n ; ш, (и„

us ) — двумерная

 

 

а

 

 

 

 

 

 

 

 

 

u(t);

 

г , число уровней

плотность вероятности

процесса

 

 

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

определяемое

из

(В.б), г =

г, -f- 1.

 

 

Элементы

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

цепи

равны

 

 

 

 

 

 

j

1

 

1^ « J

2

(и,, и

) i/«,rfa

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

_ Z < 2

l

L

^

i

^

jШ,ви,(U)(и)rf«,

 

,

(1.90а)

 

 

 

Pi

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Л и »

 

 

 

 

 

 

 

 

где ур( а р ) —безусловная

 

вероятность

СОСТОЯНИЙ ыв

и мр

 

Соответственно для двухсвязной

цепи

 

 

 

 

 

 

|"

j"

£

tt>3 (и,, ц2. «з)

dUjdu^du,

 

 

Р(а?1)

Д Ц Д

Дир

 

А и т

.

 

 

 

 

 

,

(1.91)

^ • = - ^ 3 ! ) =

а аи1

"7

 

 

ш2

 

 

 

 

 

 

Р(<*?)

 

 

 

 

^

 

(и,, на) rfu,rfa2

 

 

 

 

 

 

 

 

Д « 0

 

Д и 3

 

 

 

 

 

 

 

 

 

где

м), («,, и2 , и3) — трехмерная

 

плотность

вероятности

процесса u(t); а,

р, у " = 1, г.

 

 

 

 

 

 

 

 

 

 

52

Вычислить вероятности (1.87) — (1.91) в аналитиче­ ской форме, как правило, невозможно. Поэтому целесо­ образно использовать приближенные методы. Отметим некоторые из них.

1.Численное интегрирование на Э Ц В М с использова­ нием стандартной процедуры вычисления определенных интегралов методом Монте-Карло.

2.Численное интегрирование на Э Ц В М с использова­ нием стандартных процедур вычисления кратных инте­ гралов.

3.Использование многомерного ряда Эджворта [1'2] для разложения многомерных плотностей и почленного интегрирования ряда.

4. Моделирование на Э Ц В М последовательностей с заданными многомерными плотностями и определение

соответствующих вероятностей

методом

Монте-Карло

[см.

(1.13)].

 

 

5.

Использование теоремы

о среднем

яри большом

числе уровней квантования. Пусть число уровней кван­

тования велико, так что Ды< (0,5-7-0,7)0",

где а — средне-

квадратическое

значение процесса

u(,t).

Тогда, применив

к (1.89) — (1.91)

теорему о среднем [13], получим

 

 

/ ? 0 в ~ Д и в г М " * ) .

 

 

(1.92)

Р о ^ , ~ Л и « Д и р ш Я ( и * , . и * а ) .

 

 

(1-93)

Л * ~ Л " е " ' f f o f f i =

 

 

 

(1-94)

где и*и и*2, и*3 — середины соответствующих

интервалов

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

 

 

 

 

 

Рассмотрим

теперь другой предельный

случай — би­

нарное квантование. Пусть wn{u)—n-мерная

 

нормальная

плотность с параметрами N[0, В(х)],

где В(х) =

azR(х)—

корреляционная

функция процесса

u(t).

В

этом

случае

при квантовании по нулевому уровню для

переходных

вероятностей односвязной и двухсвязной цепи сущест­ вуют аналитические соотношения.

Для односвязной цепи переходная матрица имеет вид

/>=ГЛв

Л »1.

(1.96)

L Pit

Pa J

 

53

При этом в силу симметрии распределения [1]

 

 

 

А . =

Л о = І

4 ~ a r c c

o s /

? ( 7 , « ) .

(1.97)

где

і? (Г д ) — коэффициент

корреляции

процесса

и (/)•

В случае двухсвязной цепи переходная матрица имеет

вид

(1.99),

а

переходные

вероятности ра?

, а,

р, у = 0,1

вычисляют

из

соотношений

[1]:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

arccos R (27"д )

 

 

 

 

 

 

 

 

 

 

 

2 (л—arccos У? (7"д ))'

 

 

 

 

 

 

 

 

 

 

arccos

/? (2 Уд)

 

(1.98)

 

 

 

Рои —

Ptoo

2 arccos R(Ta)

'

 

 

 

 

 

 

Р п о = 1

 

— A n .

A o i =

1

PIOO =

J D

O . O = 1

Рои.

 

 

 

 

 

Pooi:

= :

^

Рооо>

 

 

 

 

 

 

 

 

 

 

 

Лоо

Лої

О

 

О

 

 

 

 

 

 

 

 

 

 

О

0

р010

 

рои

 

(1.99)

 

 

 

 

 

 

 

/>юо

/>ю.

 

 

 

 

 

 

 

 

 

 

 

0

 

0

 

 

 

 

 

 

 

 

 

 

0

0

Рив

Ріп

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В

свою

очередь,

 

абсолютные вероятности равны

 

 

 

 

 

 

 

 

со

 

 

 

 

 

 

 

 

 

 

 

p ( l

)

=

| a>j («)rf« =

0,5,

 

 

 

 

 

 

 

 

 

 

о

 

 

 

 

 

 

 

 

оо со

 

 

 

 

 

 

 

 

 

 

 

 

Р(„) = j

|'ЧУ2 ("І> u^du^diu—-^-^

 

 

arccosЯ (Г„) ^,

 

о

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.100)

 

 

 

P(.n) =

j

j

 

 

("i .

 

 

u3)duldu2du3^

 

 

 

 

 

 

0

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 arccos R(TS)

 

arccos /?

( 2 Г д )

 

5*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

)•

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При бинарном квантований гауссова процесса не rid нулевому уровню переходные вероятности можно рассчи­ тывать с использованием аппроксимаций (1.87), (1.88).

Рассмотрим еще одни пример бинарного квантования. Пусть wL(u) и ш 2 ( " і , « 2 ) —соответственно одномерная и двумерная плотности вероятностей косинуса фазы, соот­ ветствующие случайному процессу, получаемому на вы­ ходе системы идеальный ограничитель — фазовый детек­ тор при .подаче на ее вход нормального случайного про­ цесса:

ш, (u,) = l / 1 l y r \—и],

(1.101)

 

w2 (и„

и.,,

t ) =

1 +

 

 

 

 

* 2 Y\~u\]/

\-а\

 

 

+

 

8*-ч£Аг(,)Гг(иЛТг(іи)

(1.102)

Здесь

|и,|<1;

2 |<1; a = cos^;

 

 

 

 

/і=0

 

 

Tr(u) =

cos(r arccos и)—полиномы

Чебышева;

V{z) —

гамма-функция.

 

 

 

 

Ряд

(1.102)

есть разложение

двумерной плотности

вероятности косинуса фазы стационарного нормального процесса по ортогональным полиномам Чебышева.

 

Элементы начального вектора ро и переходной мат­

рицы

Р

определяются из (1.89) и (і . 90а) . В

соответствии

с (1.89) и (1.101) имеем

 

 

 

 

р

0 а =

1 -

F,

(uJ - l +

F t (ua+l) =

Fl а+1)

- Л а),

где

 

Fl

а)

=

р {1 <

cos <р < » а } =

1

]r arccos и

— интегральная функция распределения. Тогда

р = — [arccos ы _,. — arccos и 1 =

= ~ arccos [ u e + , u e У1 — u* j / 1 — и2а+1].

55

При бинарном квантовании иа = 0,

к , = 1,

 

р, =

arccos 0 =

0,5.

 

 

ЇЇ

 

 

С учетом соотношения [\]

 

 

I

 

 

 

 

и

 

 

 

 

где Ur-i (") — полином

Чебышева

2-го рода; г = 1, 2,

можно записать,

что

 

 

 

Г

Г М

^ _

г ^

.

J

J * ] Л - и ?

« і Л - а ?

^ ^ [ - г / ї ^ ^ - . к ) - ^ К Г

-и ' + .

] х

X [

7

^ і / г .

, ( « р ) -

l /

l - ^ f

v . ) ]

Тогда

двумерная

(безусловная)

вероятность

состояний

иа и ы?

на выходе квантизатора будет

равна

 

 

 

 

00

 

 

 

- V i - u i + l с / г . , к ^ ш / ! - " ; ^ - ,

( « р )

-

- l / i - u ^ . C / r - . f V i H -

 

(1.104)

При бинарном квантовании

« a = «g = 0 ,

=

ы р + ) =

1,

тогда получим

 

 

 

 

р ( 1 1 ) = 0 , 2 5 +

Л г (,) - i - l U , ( 0 ) ] *

=

 

 

=0 , 2 5 + 8 ^ Л г ( х ) з т 2 ( г - | - ) .

56

Ограничившись двумя первыми членами ряда и подставив выражения для Д (і), Л3 (т), получим

А , = РЫ)1Р, = 0 , 5 + ^ - [ l + 4 - ^ 3 W +

Если

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

1гф0,

тогда

 

 

/><„> =

/>? + 8 ( 1 - А ' ) Л 1 + 1 6 А а ( 1 - А » ) Л а +

...,

где

=

arccos /г; Л, (-с) и

Л, (т)

определяются в

соот­

ветствии с (1.103).

Зная условную вероятность рц с учетом симметрии закона распределения w2(uu и2) относительно порога «о = 0, получим

Pu = Poo, Pio = P o i = l —Ри.

(1.106)

1.7.3. Статистические выводы о цепях Маркова. Опре-

делив начальный вектор ро, пеіреходную матрицу Р и связность v цепи Маркова по описанной выше методике, необходимо убедиться в правомочности таких аппрокси­ маций. Проверка связности цепи основана на сравнении переходных вероятностей. Так, если цепь односвязная, то должны выполняться соотношения:

Ром) = Ро* Рог Ро№і) = Ро« 'Ро? Рог-•

( 1 Л 0 8 )

Если же эти равенства выполняются приближенно, то можно говорить лишь о 'приближенном соответствии ди­ скретного распределения цепи Маркова со связностью v = l . Степень близости этих вероятностей и справедли­ вость аппроксимации в общем виде определить затруд­ нительно. Выходом из положения является решение ти­ повой задачи обнаружения или измерения на основе вы­ бранной аппроксимации и сравнение результатов реше­ ния с полученными методом Монте-Карло на ЭЦВМ .

Проверка тех или иных статистических^ гипотез мо­ жет осуществляться как с помощью критерия отноше­ ния правдоподобия, так и по критерию %2.

57

Чаще всего''Подобные задачи решаются выборочным методом і[14]. Как известно [14, 15], оценками максималь­ ного правдоподобия для абсолютных и переходных ве­ роятностей являются частости (1,13):

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

ное

значение

р°а? , а, р = 1 , / \

Альтернативная гипотеза

Я , :

ра? не соответствует р°а? .

 

 

Отношение

правдоподобия будет равно [20]

 

 

^

где PQ и .Р, — распределения, соответствующие гипотезам

Я 0 и Я р Я — оценка в?ктора параметров методом макси­ мального правдоподобия.

Как показано в монографии [16], статистика —21пЛ имеет асимптотически центральное ^-распределение с соответствующим числом степеней свободы, если выбо­ рочное и гипотетическое распределения совпадают. При

несовпадении гипотез статистика —21пА имеет асимпто­ тически нецентральное ^-распределение и оба критерия асимптотически эффективны. На практике чаще исполь­ зуется критерий Xі.

 

Вернемся к

задача

проверки

гипотезы

о том,

что Paf=

=

р°а? , а, р =

1, г. В

этом

случае

для проверки

гипотезы

Я„

образуются

статистики

 

 

 

 

 

 

•£=УК

( ^ ~

^ )

2 .

а,р =

Т77,

(1.109)

 

 

р =1

р°-ї

 

 

 

 

58

имеющие асимптотически -/--распределение с г — І степе­

нями свободы. Поскольку переменные

Яа (/?а 3 —Р°з ) 2

Д л я

различных о. = 1, г

асимптотически

независимы,

то

для

проверки

гипотезы о том, что все переходные вероятности

ра% имеют

заданные

значения

/?°3 , образуется статистика

 

 

•**"

О

 

 

 

с г (г1) степенями

свободы.

 

 

 

 

Критическая область «непринятия гипотезы Н0»

обыч­

но выбирается двухсторонней, соответствующей довери­ тельным вероятностям р и (1—р) . Гипотеза отвергается,

если %2>%2{р)

и %2 <%2 (1—Р)> и

принимается,

если

Х 2 ( 1 Р)^12^%2{р)-

Значения х 2 (р )

и %2 (1 —р)

опреде­

ляются «з таблиц или номограммы ^-распределения [17].

Доверительная вероятность выбирается обычно

равной

р = 0,9н-0,99. Методика практического

применения

крите­

рия х2 описана в работах (68, 156, 157].

 

 

 

Рассмотрим следующую задачу: является ли данная

цепь Маркова однородной, т. е. Pa9(k)

= p a ? ;

k—1,2,....

п, —нулевая гипотеза Н0, или цепь неоднородна, т. е переходные вероятности зависят от времени pap{k) — аль­ тернативная гипотеза Ну. В этом случае оценки переход­ ных вероятностей получаются методом максимального правдоподобия по выборкам требуемого объема:

Я,р= VA «==S А«Р (*)

/ І t KtW=t

 

h?W

j'y-xjk),

 

k=l

I 1=1

ft=.l

k=\

I

ft=0

 

 

 

 

 

 

 

(1.111)

P ^ ) ^ ( k

) l X a { k

- \

)

= X^(k)

I І

KM)..

(1.112)

 

 

 

 

 

I 1=1

 

Для проверки

гипотезы

Я 0

образуется

статистика

х 2 = 2

E . E U * - ^ « p ( * ) - u * 1 7 A , p .

(1-113)

a = I

ft=l p=I

 

 

 

 

 

 

59

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