
книги из ГПНТБ / Лихарев В.А. Цифровые методы и устройства в радиолокации
.pdfооразователь соответствует определенно!! строке пере ходной матрицы. Выбор строки на данном такте зависит от предшествующего состояния, запоминаемого блоком задержки, и осуществляется с помощью дешифратора и регистров совпадения.
йРч ПЗУ
Раз
її |
'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