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

книги из ГПНТБ / Кавалеров, Г. И. Введение в информационную теорию измерений

.pdf
Скачиваний:
11
Добавлен:
22.10.2023
Размер:
12.18 Mб
Скачать

гдё

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+1th и ве­ роятность pn(t)- Указанное распределение вероятностей в общем виде определить в форме, пригодной для даль­ нейших расчетов, не удалось. Однако при очень большом числе областей квантования (У > 5 000) можно показать, что

а»(7\)>

• ехр 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

Зависимость от 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

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