книги из ГПНТБ / Кавалеров, Г. И. Введение в информационную теорию измерений
.pdfгдё
L — X — |
Gxo V p (iq — Z) |
|
А = |
ь п г |
|
У 1— p2 |
9 |
|
|
|
— X — axa~' p(iq— Z)
Ь п Г
'BX
°* V 1 — p2
Подставляя эту зависимость в (5-і) и (5-2) и учиты вая, что p(Zi) определяется соотношением (3-2), имеем формулу для вычисления количества информации в инте ресующем нас случае. При N— >-оэ имеем:
(5-8)
что совпадает с результатами, полученными в гл. 3. Та ким образом, даже для марковского процесса, получен ного путем дискретизации и квантования непрерывного нормального процесса при конечном числе областей квантования по уровню N, выражения для количества и скорости передачи информации принимают исключитель но громоздкий вид. Вид формул еще более усложняется, если измерения образуют обобщенный марковский про цесс. В целом можно сделать вывод, что метод вычис ления количества информации через условные плотности вероятности приводит при определении информационных характеристик с учетом квантования по уровню и во вре мени к непригодным для практического применения ре зультатам, причем основную трудность представляет вы числение условной энтропии
H[X(tk)\Z(tk- T ) , Z(th- 2 T ) ...].
В связи с этим в [Л. 5-1] был предложен другой спо соб вычисления верхней оценки информационных харак теристик при коррелированных измерениях. В основе этого метода лежит замена непрерывного процесса х(і) на входе АЦП случайной последовательностью, образо ванной моментами пересечения x ( t ) с границами облас тей квантования по уровню (рис. 5-1). Эту последова тельность назовем входящим потоком, а каждое пересе чение— заявкой на обслуживание.
Отвлекаясь от содержания каждой заявки, можно охарактеризовать поток заявок на обслуживание метода ми, принятыми в теории массового обслуживания. Сфор-
183
Мулируем некоторые очевидные положения, Не остапяйливаясь на их строгом доказательстве:
1) при стационарности X(t) поток заявок стациона рен;
2) при сингулярности Х(і) поток заявок ординарен; 3) последействие потока зависит от корреляционной
функции р*(тг);
4)интенсивность потока равна конечному числу Хх, определяемому из записанных ниже формул (5-11)— (5-13);
5)в силу теоремы В. С. Королюка {Л. 5-5] из условий стационарности и ординарности потока вытекает, что па раметр этого потока также равен Ä*;
6) вероятность п заявок за время t при стационар ности потока равна:
t |
|
|
A>W= **J['P«-i(f) — <Р» (*)]<#. я = 1 . |
2, |
3 ...; |
о |
|
|
t |
|
|
Po(0 = 1 —* j ?о (г) dz, |
|
|
о |
|
|
где фп( 0 — функция Пальма — Хинчина |
{Л. |
5-5], кото |
рая представляет собой вероятность поступления за вре мя t ровно п заявок при условии, что в начальный мо мент этого промежутка времени заявка поступила.
Для дальнейших вычислений существенное значение имеют такие характеристики потока пересечений (зая-
184
вок), как распределение интервалов Tk= th+1—th и ве роятность pn(t)- Указанное распределение вероятностей в общем виде определить в форме, пригодной для даль нейших расчетов, не удалось. Однако при очень большом числе областей квантования (У > 5 000) можно показать, что
а»(7\)> |
2д |
• ехр 2<ОТ"«(0) |
(5-9) |
; У2*?"х (0 ) |
При числе областей квантования б ^Л /^300, как по казано в [Л. 5-10, 5-11] на обработке большого экспери ментального материала, удовлетворительная аппрокси мация достигается при применении модели пуассоновско го потока.
Используя формулу (6-19) из [Л. 5-21], можно полу чить математическое ожидание числа пересечений гра ниц квантов за время (t0, /о+7’) в виде
Я* &>. t0+ Т) = |
j1 |
|
11 x' (t) I w2 [iq, x' (f)] dx' dt. |
(5-10) |
|
І= 1 1 о |
—oo |
|
|
Если X(t) — стационарный процесс, то среднее число |
||||
выбросов в единицу времени |
|
|||
|
|
N |
со |
|
|
Ях= |
2 |
j x'w2{iq, x')dx', |
(5-11) |
і=1 —oo
где Wo — двумерная плотность вероятностей самого процесса и его первой производной X' в совпадающие мо менты времени. При гауссовском характере процесса X(t) параметр
Я* = r ~ f « (0) 2] exp [ ~ 0 . 5 ( - f ^ y ] . (5-12)
і=I
Зависимость 1Х от N=l/q представлена на рис. 5-2. При большом числе областей квантования (N1>16)
с приемлемой точностью можно представить последнее выражение в виде
1/ |
Р^ ( 0)= |
У ^ . |
(5-13) |
|
Я У |
- 2тс |
I |
||
|
где — щкада прибора.
!§§
Процедуру вычисления верхней оценки информацион ных характеристик рассмотрим для случая, когда дис кретные отсчеты образуют марковский процесс.
Напомним, что случайный процесс называется обоб щенным м а р к о в с к и м процессом порядка т, если для конечного множества моментов времени t \ < . h ,
.,., <_tn- i < t n, п ^ т он полностью определяется совмест ной плотностью вероятности в т моментов времени, так что условная плотность вероятности
W (Хп, t?i 1Хт\—1, tn—l, ■• |
X-otz', Xi, tij — |
= ay(Xn, tn\%n—1, tn—i, ■• |
Xn—m, ^n—m)- |
Марковский процесс первого порядка называют прос
то марковским процессом (см. стр. |
35). |
|
|
|
|||||||
Если при предыдущем измерении получено значение |
|||||||||||
Zi {tu—Т), |
то |
в момент |
времени |
первого |
пересечения |
||||||
|
|
|
|
|
|
уровней, следующего пос |
|||||
|
|
|
|
|
|
ле th—Т, возможны зна |
|||||
|
|
|
|
|
|
чения Zi+i и Zi-1. Энтро |
|||||
|
|
|
|
|
|
пия |
в |
этом |
временном |
||
|
|
|
|
|
|
сечении зависит от соот |
|||||
|
|
|
|
|
|
ношения p(Zi-i) и p(Zi+1), |
|||||
|
|
|
|
|
|
причем |
наибольшее |
зна |
|||
|
|
|
|
|
|
чение |
энтропии |
реали |
|||
|
|
|
|
|
|
зуется |
при |
p(Zi-1) = |
|||
|
|
|
|
|
|
= p(z,-+1) =0,5. |
Хотя |
на |
|||
|
|
|
|
|
|
самом деле эти вероят |
|||||
|
|
|
|
|
|
ности |
почти |
при |
всех і |
||
|
|
|
|
|
|
не равны между собой, |
|||||
|
|
|
|
|
|
для |
|
искомой |
оценки |
||
Рис. 5-2. |
|
Зависимость |
парамет |
сверху принимается усло |
|||||||
ра %х от числа областей квантова |
вие |
равенства. Тогда |
ве |
||||||||
ния при нормальном законе изме |
роятность |
отклонения |
|||||||||
нения |
измеряемой |
величины |
процесса на |
пг |
квантов |
||||||
(Ь/ах=6). |
|
|
|
|
|
||||||
|
|
|
|
|
|
от предыдущего |
отсчета |
||||
при наличии п пересечений процессом X(t) |
уровней |
||||||||||
квантования |
выразится |
как |
|
|
|
|
|
|
|||
. . . |
(2~пС(п т)/2 |
при /i-4-m |
четном и т < я ; |
|
|||||||
р (т I п) = і |
п |
|
|
1 |
|
|
|
|
|
||
|
|
I |
0 |
|
при п-\-ш нечетном и т'у-п. |
||||||
|
|
|
|
|
|
|
|
|
|
(5-14) |
Возможные варианты переходов показаны на рис. 5-3. Оценка энтропии в сечении, соответствующем п перехо-
186
дам с уровня на уровень, имеет вид:
п
Я* (а) = — 2 ^ р(т I и) log р (т | п) —р (01п) log р (0 | а).
т—1
(5-15)
Оценка энтропии в момент времени tu
СО
Я* (X) = - р 0(Т) log- р0 (Т) - 2 р*(Г) Я* (п), (5-16)
Л=І
где Р п ( Т ) — вероятность п пересечений за время Т. Следующее допущение связано с видом входящего по
тока. В § 5-3 будет показано, что при определенных условиях возможно описание его как пуассоновского.
Рис. 5-3. Сетка переходов и вероятности перехода на сосед ние уровни при различном числе шагов п.
Принятие этого допущения, вообще говоря, |
приводит |
||
к дальнейшему завышению оценки энтропии. |
Принимая |
||
Рп(Т) = ^ р е х р ( - |
1ХТ), |
|
(5-17) |
где .Я* — математическое ожидание |
числа |
пересечений |
|
в единицу времени (интенсивность потока) |
и подставляя |
187
(5-15) и (5-І7) в (5-І6), после ряда преобразований при условии \ ХТ< 2 имеем:
|
|
|
Я*(Х) |
г . ^ Г е х р (— Хх Т) |
|
(5-18) |
|||||
|
|
|
|
1 - |
кт/ъ |
|
|
||||
|
|
|
|
|
|
|
|
|
|
||
Сравнительные расчеты по формулам |
(5-18), с одной |
||||||||||
стороны, |
и формулам (5-2) и |
(5-7), с другой, показы |
|||||||||
вают, |
что |
удовлетворительная |
оценка |
достигается |
при |
||||||
|
|
|
|
|
|
|
сравнительно малых |
N и |
|||
|
|
|
|
|
|
|
Г (условие /\,Г<0,15). За |
||||
|
|
|
|
|
|
|
висимость |
количества |
ин |
||
|
|
|
|
|
|
|
формации |
от |
интервала |
||
|
|
|
|
|
|
|
времени |
Т показана |
на |
||
|
|
|
|
|
|
|
рис. 5-4. Для расшире |
||||
|
|
|
|
|
|
|
ния зоны действия оце |
||||
|
|
|
|
|
|
|
ночного соотношения |
вве |
|||
|
|
|
|
|
|
|
дем некоторые уточнения. |
||||
Рис. 5-4. Зависимость количества |
|
Учтем |
прежде |
всего |
эф |
||||||
информации |
на |
одно |
измерение |
|
фект |
«отражения» |
от |
||||
от интервала |
опроса |
(<?/0^ = 50, |
|
краев |
шкалы. |
Так |
как |
||||
Хх=0,1). |
|
|
|
|
|
принято условие O ^ x ^ L |
|||||
нечно, |
то |
после |
|
|
|
и число квантов N ко |
|||||
|
определенного |
числа |
переходов |
і—m процесс X(t) может достигнуть пределов, уста новленных шкалой. При этом действие формулы (5-14) прекращается. Обозначим вероятность попадания в точ
ку і—m за п |
шагов при N областях квантования |
как |
||
Pi(N, i, |
іг, пг), |
а вероятность |
попадания в точку |
і+ т |
при тех |
же условиях как pz{N, |
i, п, m ) . Тогда |
|
п + т
p,(N, і, п, т) = 2 |
п |
2 |
|
П—П -(лг-о |
п + т |
— N |
|
+ С 2 |
+ С« |
|
|
|
п + т |
- ( i-Д ') |
|
|
+ с ~ ^ |
||
|
|
+ . J |
‘п
|
|
|
|
|
|
|
п + т |
p2(N, |
і, п, т) — 2~п с п |
+ |
С / |
||||
|
п —т |
|
|
|
п + т |
|
я—т |
+ |
С ~ |
‘ + |
‘ |
С |
*1 |
+ |
с * |
' |
*» |
|
|
п+т
2+
—N
+
(5-19)
+
- (лг-0
+
188
Оценка условной энтропии запишется в виде
Я дг(л) = 2 р ( ^ ) Я дг(і , я ), |
(5-20) |
/=і |
|
где
і
Яд, (і, п) = — £ А (Я, г, я, т) logр, (Я, г, я, яг) —
Ш=1
УѴ—і |
|
|
|
і, я, т) — |
|
— £ Ра (Я, |
г, |
я, |
/и) log р2(Я, |
|
|
m=I |
|
|
0)logр (Я, г, |
|
|
— р(Я, |
г, |
/г, |
я, 0). |
(5-21) |
Индекс г в формуле (5-21) означает зависимость от начального уровня 2j. Учет вероятностей д(х,) представ ляет собой еще одно уточнение оценки энтропии. Отме тим, что
Рі (Я, іи я, т) = д 2(Я, k, я, т) |
при іі = Я—і2. |
(5-22) |
|||||||||
Подставляя (5-22) и |
(5-21) |
в |
(5-20), имеем: |
|
|||||||
N |
Г |
і |
|
|
|
|
|
|
|
|
|
HN(я) = — £ р (Хі) |
|
2 £ |
Pi (Я, |
г, |
я, |
яг) log а (Я, |
г, я, яг)+ |
||||
(=1 |
_ |
Ш=1 |
|
|
|
|
|
|
|
|
|
+ р (Я, г, я, |
0) log р (Я, |
|
г, я, |
0) |
, |
|
(5-23) |
||||
где |
|
|
(і+і) я |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Р ( Х і ) = |
^ |
W |
( Х ) СІХ. |
|
|
|
|
||||
|
|
|
ІЯ |
|
|
|
|
|
|
|
|
Результаты расчета условной энтропии ЯЛ(я) |
по этой |
||||||||||
формуле на ЭВМ приведены в табл. 5-1. |
|
пересечений |
|||||||||
Пример зависимости |
Н^{п) |
от |
числа |
||||||||
приведен на рис. |
5-5. С увеличением |
числа |
шагов п |
средняя энтропия на отсчет увеличивается, так как рас пределение вероятностей при этом становится более рав
номерным. |
Подставляя |
полученное |
значение |
HN{n) |
|||
в формулу |
(5-16) |
вместо |
Я* (я) |
при |
условии |
(5-17), |
|
имеем: |
|
|
|
|
|
|
|
|
|
|
СО |
|
Яд, (П) |
|
|
Я(Х) = ЯхГ ехр ( - |
ХхТ) 1 , 4 + 2 |
{ХхТ)П |
(5-24) |
||||
ПІ |
п —1
189
Т а б л и ца 5- i
|
|
|
H'N (/1). |
дв. ед. |
|
|
|
N = 10 |
/V = 0 |
N = 30 |
N = 50 |
N = 100 |
N = oo |
1 |
0,9000 |
0,9500 |
0,9667 |
0,9800 |
0,9900 |
1,0000 |
3 |
1,6051 |
1,7084 |
1,7311 |
1,7706 |
1,7913 |
1,8113 |
5 |
1.9179 |
2,0580 |
2,1056 |
2,1420 |
2,1693 |
2,1984 |
7 |
2,1078 |
2,2772 |
2,3344 |
2,3734 |
2,4118 |
2,4468 |
9 |
2,2342 |
2,4322 |
2,4950 |
2,5508 |
2,5924 |
2,6302 |
11 |
2,3244 |
2,5506 |
2,6258 |
2,6858' |
2,7316 |
2,7762 |
13 |
2,3892 |
2,6454 |
2,7290 |
2,7958 |
2,8470 |
2,8966 |
15 |
2,4384 |
2,7246 |
2,8164 |
2,8898 |
2,9448 |
3,0006 |
17 |
2,4752 |
2,7908 |
2,8906 |
2,9704 |
3,0326 |
3,0900 |
19 |
2,5040 |
2,8490 |
2,9560 |
3,0416 |
3,1086 |
3,1702 |
21 |
2,5262 |
2,9006 |
3,0150 |
3,1066 |
3,1742 |
3,2432 |
23 |
2,5426 |
2,9458 |
3,0668 |
3,1632 |
3,2336 |
3,3088 |
25 |
2,5548 |
2,9860 |
3,1134 |
3,2152 |
3,2894 |
3,3680 |
27 |
2,5638 |
3,0224 |
3,1562 |
3,2632 |
3,3184 |
3,4236 |
29 |
2,5702 |
3,0550 |
3,1950 |
3,3074 |
3,3926 |
3,4750 |
31 |
2,5752 |
3,0858 |
3,2318 |
3,3486 |
3,4348 |
3,5236 |
33 |
2,5780 |
3,1136 |
3,2658 |
3,3874 |
3,4760 |
3,5692 |
35 |
2,5796 |
3,1396 |
3,2974 |
3,4234 |
3,5202 |
3,6126 |
37 |
2,5802 |
3,1634 |
3,3274 |
3,4580 |
3,5576 |
3,6534 |
39 |
2,5796 |
3,1852 |
3,3550 |
3,4898 |
3,5934 |
3,6916 |
45 |
2,5718 |
3,2392 |
3,4256 |
3,5724 |
3,6854 |
3,7928 |
55 |
2,5524 |
3,3080 |
3,5394 |
3.6S82 |
3,8134 |
3,9078 |
|
|
|
Т аблица |
5-2 |
T, c |
V r> |
H (Л), дв. ед. |
E, дв. ед. |
|
|
|
|
|
|
|
0,05 |
0,1139 |
0,1139 |
|
|
0,1 |
0,2194 |
'0,2194 |
|
1 |
0,4 |
0,7062 |
0,7062 |
|
1.0 |
1,2163 |
1,2163 |
|
|
|
1,2 |
1,4279 |
1,4279 |
|
|
1,6 |
1,4427 |
1,4427 |
|
|
2,0 |
1,5331 |
1,5331 |
|
|
0,05 |
0,1139 |
0,01139 |
|
|
0; 1 |
0,2194 |
0,02194 |
|
10 |
0,4 |
0,7062 |
0,07062 |
■ |
1,0 |
1,2162 |
0,12163 |
||
|
1,2 |
1,4279 |
0,14279 |
|
|
1,6 |
1,4427 |
0,14427 |
|
|
2,0 |
1,5331 |
0,15331 |
|
190
Расчет |
|
по |
формуле |
|
|
|||
(5-24) показывает, что |
|
|
||||||
данная оценка дает удов |
|
|
||||||
летворительное |
прибли |
|
|
|||||
жение при /ЦГ<2. При |
|
|
||||||
мер результатов |
расчета |
|
|
|||||
приведен в табл. 5-2. |
что |
|
|
|||||
Нетрудно |
видеть, |
|
|
|||||
учет краевого эффекта |
не |
|
|
|||||
меняет качественной |
кар |
|
|
|||||
тины |
зависимости |
# ( X) |
|
|
||||
от Т. При этом |
скорость |
|
|
|||||
прохождения информации |
|
|
||||||
по измерительному тракту |
|
|
||||||
Ех^ . г |
= |
4 -{Я[Х (**)І- |
Рис. 5-5. Зависимость условной |
|||||
ПЫХ |
1 |
|
|
|
|
|||
- H [ X ( t k) \ Z ( t u)]}, |
|
энтропии от числа пересечений. |
||||||
где H[X(tk)] |
|
|
(5-25) |
|
|
|||
приблизительно определяется |
формулой |
|||||||
(5-25), |
а H[X(th)\Z(tk)] — формулами § 3-3. |
В заключе |
||||||
ние отметим, |
что при Т->-0 последовательность измерен |
|||||||
ных значений стремится к марковской. |
|
|||||||
5-2. ДИСПЕРСИЯ ПОГРЕШНОСТИ |
|
|||||||
Перейдем далее |
к средиеквадратическнм оценкам |
погрешностей при совместном квантовании по уровню и дискретизации во времени. Значение этих оценок особен но возрастает в связи с тем, что среднеквадратпческий критерий по существу не требует знания порядка мар ковского процесса, в чем мы убеждались в гл. 4.
Пусть исходный процесс x(t) после промежуточных преобразований аналоговыми блоками и дискретизации по временным сечениям поступает на квантующий пре образователь. Так же, как и в предыдущем параграфе этой главы, будем предполагать, что на идеальный пре образователь поступает сигнал x(t) + y Bx(t), где уВх(і) — суперпозиция погрешностей, накопленных на предыду щих блоках, и погрешностей самого АЦП, приведенных ко входу; в данном рассмотрении предполагается, что интервал дискретизации Т = const и шаг квантования q= L/N. Моменты времени, в которые поступают сигналы на вход АЦП, обозначим tu t2, ■.., А, • • •, причем
191
Используя обозначения § 3-2, для погрешности дис кретности в момент считывания показаний можно запи сать:
|
y\<(h) =u(th)—Zi(tk), если /<7<ы(М<(і+1)<7, |
(5‘26) |
||||||||
где |
u(tk)=x(tk)+yBx(tk)- |
|
|
|
|
|
||||
Перепишем это выражение в виде |
|
|
||||||||
|
|
Ук(іи) = и (Ui) = 7 ent[u(lk)lq], |
(5-27) |
|||||||
т. е. |
примем для определенности Zi = iq. |
|
||||||||
Математическое ожидание и дисперсия этой погреш |
||||||||||
ности соответственно равны: |
|
|
|
|
||||||
|
|
|
У к ( 4 ) = |
Ä ' (tu) + |
F „х (tu) — |
|
||||
|
|
|
- |
q ent {[X(4) + |
r BX(4)lM; |
(5-28) |
||||
|
ук |
= |
x |
+ |
4 |
(k) + |
2 R xg |
(h, tu) + |
|
|
|
|
|
Vax |
|
Лі/вх |
|
||||
|
+ |
\ z (t)i |
— 2R xZ(tu> |
tu) — 27?"bx (tu, tu), |
(5-29) |
|||||
При стационарности X(t) |
и Y Bli(t) |
имеем: |
|
|||||||
|
|
|
|
|
|
|
|
N—\ |
|
(5-30) |
|
|
У к (tu)-- tx -f- <Рвх — Q 'Eitp (zi), |
||||||||
где |
|
|
|
|
|
(<+i) я |
i= 0 |
|
|
|
|
|
|
p(Zi)= |
|
|
|
||||
|
|
|
|
J |
w(u)du\ |
|
||||
|
|
|
|
|
|
iq |
|
|
|
|
|
jV—1 г |
V |
^ |
— 3x3+ |
° L + |
|
(0)+ |
|
||
|
|
iV—1 |
|
|
|
|
|
|
||
+ < гs£ |
i - |
X ip Ш |
|
27?xz (0) — 27?»nx* (0). |
(5-31) |
|||||
|
i= 0 |
|
i=0 |
|
|
|
|
|
|
Полученные зависимости легко раскрываются при за дании конкретных законов распределения X и Увх. На пример, при соблюдении условий (3-11)
У к (tu) = 0,5q, |
(5-32) |
а если добавить условие N>2L/a, то
О2 (tu) Ъ <7712. |
(5-33) |
Отметим, что формулу (5-33) иногда применяют без учета условий, при которых она справедлива,
192