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

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

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

i'l

Но в Х \ (J Mj входят те и только те точки, которые либо не

/=!

принадлежат ни одному из Fj, либо принадлежат и (а',р (а')). Поэтому

I Я, (А() - Я,

(А"() I <

Я, (Х \

1J F,) +

Я,

(a'. Р (а'))) =

 

 

Я

 

 

 

 

/=1

 

Я

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= I -

S

Я, (Fj) +

Я, (U (a', р (а'))) <

1-

Ц

(Я, (А,) - ф п )

+

 

/=I

 

 

.

 

 

 

/=I

 

 

 

 

 

 

 

+ Р(а') = «/2 + р(а')<».

 

 

 

Точно так же получаются неравенства с Я2. Лемма доказана.

 

Если Ѳ (а) — некоторая функция на множестве

всех разбиений, то

мы будем писать

Ѳ (ajp^^g —» Ѳ, если для любого е > 0

можно ука­

зать 5 >

0,

что

для любого разбиения е > 0 , для

которого р(а)<8;

|Ѳ(а) — Ѳ| <е . Аналогичный смысл имеетсоотношенне Ѳ(а)--------

*■0.

Предложение

1.

Н Ра (Я,, а ) ^ ^

Я Ра(Я,).

 

 

 

Д о к а з а т е л ь с т в о .

Зафиксируем

некоторое е > 0.

Используя

определение

Я р (Я,),

найдем такое

разбиение

а0,

что

 

 

 

 

 

Я Рі (Р.) - <2 <

Н Ра j,а«).

 

(3-24)

Будем считать, что а0 — невырожденное разбиение (если а0 вы­

рождено, то можно перейти к невырожденному разбиению а0, полу-' ченному из а0 объединением всех вырожденных и одного из невы­ рожденных элементов разбиения а0 в один элемент). Очевидно, что

 

Нра(Я,, а0) =

ffp 2(Я„ а0).

 

 

По лемме 1 найдется такое у>0, что для любого разбиения Ь =

= {В,}, і — 1........

для которого выполнены

неравенства

 

 

|P i (A °,)-P 1(B ,)|< y; |Я2(А°0 —

 

 

—Ра(Ві) I < y;

(»=1....... П),

 

(3-25)

справедлива оценка

 

 

 

 

I Н р%(Я,, а0) -

Я ( Я 2, b) I <

е/2.

(3-26)

Из леммы 2 следует, что найдется такое q>0, что для каждого

разбиения а, для которого р(а) <<у, найдется такое разбиение Ь, бо­ лее крупное, чем а, н такое, что выполнены неравенства (3-25). Тогда, используя свойство 2 и оценки (3-26) и (3-24), получаем, что

Н Рі (Я,) ^ Н р%( л , а) > Н р%(Я,, Ь) 5= Я Ра (Я,, а») -

- е / 2 > Я Яі( Я , ) - « .

Предложение доказано.

123

С л е д с т в и е

1.

H p

(P t , а ) - -----H p lP ,) .

самом

деле,

р (а) -» 0, так как р (а) <

 

'clinma-vO

 

'

 

 

 

diam а).

 

 

 

 

 

 

С л е д с т в и е

2.

Если

а7*-—последовательность разбиений такая,

чтор(а’‘) -----*• 0 (или же diam а”----- *0),

и

при всех п

а71+1

мельче,

П-ИЗО

 

 

 

к-изо

 

г

 

 

 

 

чем а", то

(Я,, а71)---- >

Н р (Я,),

не убывая.

 

 

 

 

З а м е ч а н и е

1.

Как следует [из [Л.

1-54],

если Н р (Я ,)< с о

то мера Я,

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

äP,/dP 2 относительно меры А,

и

 

 

 

Нр, (^>і)=

*°2 ~Jp^ (х) (^х)-

 

 

 

 

 

 

 

 

х

 

 

 

 

 

 

Если, кроме того, Л, и Рг имеют плотности рх и р2 относительно

некоторой меры р, то

 

 

 

 

 

 

 

 

 

 

^ ( ' D1) = j

І0%рР'Хк) Р'

 

 

 

 

 

 

 

 

 

•V

 

 

 

 

 

 

З а м е н а н и е

2.

Возникает естественный вопрос, всегда ли су­

ществуют

разбиения

а с произвольно малыми р(а) и диаметром

diam а. Что касается р(а), то при сделанных предположениях отно­

сительно множества X это действительно так. Для того чтобы это

было верно и для diam а нужно потребовать, чтобы (помимо сде­

ланных предположений) X было компактно.

 

 

 

 

Отсюда ясно, что соотношения типа Н D (Я,,

а )-------►

Н „ (Я.)

содержательны только

 

 

 

 

 

diam а-»0

 

 

в том случае, когда множество X компактно.

3.

Теперь

мы рассмотрим случай, когда

Х = ХіХХ2, приче

X] и Х2 — полные

сепарабельные метрические пространства

с рас­

стояниями г 1, г2, и X снабжено расстоянием г = Ѵ г2\ + г22. Пусть на

некотором вероятностном пространстве заданы случайная величина —1со значениями в Xt и случайная величина Е2 со значениями в Х2, тем самым определена случайная величина Е =(Е |, Н2) со зна­

чениями в Х = ХіХХг.

Пусть Я-

,

 

Я-, ,

Я ,-

-

, — распределения

(т.

е. меры)

случай-

-1

 

 

-2

1-1, - 21

 

в X, , X2, X.

В X можно ввес­

пых величин Е,,

 

Е2, Е соответственно

ти также меру

Р -

,

которая

являлась бы распределением Е,

если бы Е, и

Е2

1—2

 

 

(иначе говоря,

Яс

— прямое

были независимы

произведение мер Яс

и Р-.

 

 

 

 

 

1/\ —2

Я. =

). Если положить Я, = Я ,-

-

= Р ^- у ^ - у

 

 

-1

 

-2

Н р^{Р^)

переходит

1-1. -2І

 

то определение

в определение ин­

формации случайной величины

Е,

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

случайной

величины

Е2, обозначаемой через / (Е,,

Е2)

(см. [Л. 1-54])

или

w .

Пусть а 1 =

{А] },

і =

1.......пи а 2 =

{А? },

 

-1

-2

/г2 —

у == 1 , . . . ,

разбиения пространств Х ,,Х 2.

Построим

разбиение

а’Х а 2

пространс­

тва X, относя две точки (х,, х2) и (х[, х 12) к одному и тому же эле­

менту разбиения а 1 X а2, если Xj и xj находились в одном элементе разбиения а 1, а х2 и х \ — в одном элементе разбиения а2. Разбне-

124

»ne a1 X а2

состоит

из

/z, X пг

множеств

вида

А' х Ау , где А! —

элемент

а 1,

А? — элемент а2.

 

 

 

 

 

 

 

 

 

 

Л е м м а

3.

Пусть

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рг-

(а1) =

inf

max (е, p^

(u (a',

 

е))),

 

 

 

 

 

 

 

 

 

e > o

 

 

Q 1

 

 

 

 

 

 

 

 

 

 

 

 

(а2) =

inf

шах (e, p-, (u(a2,e))).

 

 

Тогда

 

 

- 2

 

E > 0

 

 

w2

 

 

 

 

 

 

 

 

 

p(a’ X a 2X p -

(a ')-fp -

 

 

 

 

 

 

 

 

 

 

 

(a2).

 

 

Кроме того,

 

 

 

 

‘-i

 

‘-г

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

diam а <

V (cliam а 1)2 +

(diam а2)2.

 

 

Д о к а з а т е л ь с т в о .

Из неравенства

 

 

 

 

 

diam (А \

Х А/ ) > V

 

 

следует,

что либо

diam (А;')> е ,,

либо diam (Ay ) >

е2. Поэтому

 

 

 

 

 

 

 

 

 

 

 

 

 

 

inf max (е, Я™

 

. (и (а, е))) =

 

 

 

 

 

 

е>0

 

 

V—1>“ 2f

 

 

 

 

 

 

= inf

max)]/

в?-I-

Р , ^

ъ

,(u<a,]/

e? +

Eo ) ) X inf max(e,-j-

 

в ,, Ea> 0

 

 

VT

 

-

 

l - i . - г . )

r

 

 

 

 

z

e,.ea>0

+ e2. P-s (U (a’, e,))+

 

(u (а2, ег)))

< in f шахе,, Я ,

(u(a>,e,))) +

 

“ l

 

 

 

 

w’2

 

 

 

e,p ea> 0

 

 

 

 

 

 

 

+

 

max (e2,

 

(u (a2, e2)))) = P-,

(а’) +

Яя

(а2).

 

 

 

 

 

 

ы 2

 

 

 

‘- l

 

 

 

 

w2

 

 

Диалогично доказывается,

что

 

 

 

 

 

 

 

 

 

 

inf

max (e, P „

v r

. (u (a, e))) < p_,

 

(a1) + p,

(a2).

 

 

в>о

 

 

Ѵ-1

/Ч-2У

 

 

 

- 1

 

 

 

■-г

 

Второе из доказываемых неравенств следует из того, что

 

 

diam (А! X А?) <

 

(diam (Aj ))2 +

(diam (И?))2 .

 

Лемма доказана.

 

 

 

 

 

 

 

 

 

 

 

 

 

Будем писать вместо

 

 

 

 

 

 

 

 

 

 

Н

£.j

ulj

(2

н у а 1X а2) величину X

1*^2/

(а1X а2).

 

 

» 1

^2/

 

 

 

 

 

\

 

 

 

 

П р е д л о ж е н и е 2.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/ / Я

^ \ (а1X а2)------- —

/(3„

 

Е2).

 

 

 

 

 

 

1“ Х' -2^

 

 

 

Р„

(а1)->0

 

 

 

 

 

Р* (а2)—о

‘—‘2

Д о к а з а т е л ь с т в о . Нужно применить предложение 1 илемму 3.

125

1),

З а м е ч а н и е

3. Как указано

в [Л. 1-54]

(см. также замечание

если / (Hj, Sj) <

со,

то

 

 

 

 

Г Г

log

dP 13

Н )

(-’Сц х 2)

(dXi X dx2).

I

(йі> аі) = \ I

cfPZ

^

x,x> S,XS2

Есл^і на X, и X2 заданы такие меры А, и А2, что Р (3 ,,3 2) имеет

плотность р (:ct, х 2) относительно \ , Х \ 2, то тогда Р ^ и Р-, пмс-

'“І —42

ют плотности Рі (х,).

и р 2 (х2) относительно Аі и А2 и справедлива

формула

 

 

I (“1. ^г) = ^ |*

log /?, (у,) р„ (.С,) Р

Х-^ ^-іХ^-2 ( d x ^ d x 2).

х,'ха

Вприложениях обычно Аі, А2— меры Лебега на прямой (или на некотором ее отрезке) или в некоторой области многомерного про­ странства, а 3 j и Н2 являются скалярными пли векторными вели­ чинами. В качестве г в этих случаях можно взять обычное рас­ стояние.

Вусловиях предложения 2 справедливы также утверждения, аналогичные следствиям 1, 2 и замечанию 2.

4.

Начиная отсюда, мы будем

считать, что Хі

= Х2 = Хп,

а величи

ны 3, и

3 2 будем

интерпретировать

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

как истинное

значение

измеряемой

величины и смесь его с погрешностью на

вхо­

де; относительно характера их взаимодействия никаких спе­ циальных предположений не делается. Предположим, что 3, и Е2 подвергаются дискретизации. Математически дискретизацию можно описать заданием на Х0 некоторой функции G с конечным множе­ ством различных значений (Яі, . .., g „). Эти значения не обязаны

быть числами и могут пониматься как абстрактные символы. Дискре­ тизация G отвечает разбиение ас пространства Аг0, состоящее из п

элементов, причем Л,- содержит те и только те точки хеА'о, для

которых G (x)— gi. Положим А'с =

( Яь. . ,

Я«)

51 определим в Х 0

расстояние формулой r ( g g j )

= б ,2, где

б,-2 — символ Кронекера,

равный нулю при і=И=/ и единице

при і= /.

 

 

Условимся также писать р,-,

(G)

вместо

p-j

(aG), diam G вместо

1—'i

 

“і

 

diam aG и т. д.

Вопрос, которым мы будем заниматься, таков: как ведет себя информация о Еі (или о дискретизированном сигнале G(Hi), содер­ жащаяся в дискретизированной смеси G(E2), когда элементы раз­ биения а0 становятся достаточно малыми?

П р е д л о ж е н и е 3.

/ ( 3 I, G( S2) ) — Т с)/ ( 3 і,Н2).

ѣ

До к а з а т е л ь с т в о . Заметим, что

/(3,, G (Н2)) = H pE^G (E2) ( Я(3„ G(32)))-

Если мы обозначим разбиение^ пространства Х а на отдельные

точки через ас, то, очевидно, diam ас=р(ас) =0 и для произвольного

12R

разбиения

а1 пространства Х\ согласно лемме 3 выполнено неравен­

ство

p(a1 X aGX p Si(a1).

 

 

 

Зафиксируем е < 0 . Из предложения 2 следует, что найдется та­

кое (7> 0,

что для всякого разбиения а 1 пространства

такого, что

р_ (а') < Х 11 Для любой дискретизации G, для которой p^ (G)<Cq,

справедливо неравенство

Hr.

{Р(-р

пн

а ' X aG) I (-^1, іа) I

'2.X V (3‘-s*)

 

Но, как нетрудно видеть,

 

 

Н р -

(Р,

а' хао) =

 

,, с (Н2)’

а' X ас) < Я я3і хО (32)(Я(

 

= / (3 ,,G (S 2)).

Поэтому

/(H „6 '(S 2))5=/(31,E 2) - , .

£-

°(зг)))-

С другой стороны, как было показано в начале данного парагра-

фа,

/ ( 3 i,G (32) ) ^ / ( 3 2, 3 2).

Предложение доказано. П р е д л о ж е н и е 4.

/(G(3,)),

G(32)). PT (G )-*0; р - (G)

/ (—и “г)*

 

 

-1

—2

 

 

Д о к а з а т е л ь с т в о .

Предложение следует

из предложения

2 и того факта, что

 

 

 

 

 

/(0(3,);

0(32)) =

/-/р_

1Äw2

 

3gX ao).

 

 

^

 

 

З а м е ч а н и е

4. В условиях

предложений

3

и 4 справедливы

утверждения о монотонности сходимости по измельчающимся по­ следовательностям разбиений (см. следствие 2), о возможности за­ мены р(а) на сііаш а (см. следствие 2) и формула для #(Е і, S2) (см. замечание 3).

Установив таким образом общие закономерности квантования в статике, перейдем далее к формулам для инженерных расчетов. Воспользуемся системой обозна­ чений, принятой в предыдущем параграфе и в § 1-5.

Для оценки количества информации воспользуемся формулой (1-52), причем z в данном случае есть диск­ ретная величина, а к — непрерывная.

127

Так же, как и в § 3-2, будем различать

два метода

оценки I z ^ x и І г^ Хц в зависимости от

соотношения

единиц измерения х и рв.\- с квантом д. Входящие -в фор­ мулу (1-52) условные энтропии выразятся следующими формулами:

при первом подходе (т. е. с учетом 'погрешности диск­ ретности)

н (Z I Х )=

- (w (х)

£ р (г,: | х) log р & | х) dx\

(3-27)

 

О

1=0

 

 

N - X

Б + 1

(3-28)

Н {X I Z) = — £ р (z,-)

I“ w (XI Zi) log w (x | zf) dx;

 

i= 0

c . i

 

при втором

подходе

 

 

Я (Z 1Х )=

—Л£ р (хі) Л£ Р(?j 1Xt) log р (г, I Х іУ ,

(3-29)

 

;=о

/=о

 

H(X\Z) = —Л£ р (Zi) £ р (Xi I Zi) log р (Xj 1Zi), (3-30)

1=0 1=0

где p(Zj\x) определяется формулой (3-10); p(zj\xi) и p(Xi) — формулой (3-15); p(zd — формулой (3-2).

Безусловные энтропии можно записать в форме

H (Z ) '= — £ p (^)lo g p (?i);

(3-31)

і=0

при первом подходе

L

 

Н ( Х ) = — ^w(x)logw(x)dx;

(3-32)

6

 

при втором подходе

 

Н ( Х ) = — £ Р (Хі) log Р (Хі).

(3-33)

і=0

 

Для вычислений можно воспользоваться как левой частью формулы (1-52), так и правой. В [Л. 3-17—3-20, 3-22] использовалось выражение H(Z)H(Z\X), в то время как в опубликованных одновременно работах [Л. 3-13—3-15] использовалось выражение Н(Х)

128

H(X\Z). Оба подхода, естественно, дают совпадающие результаты.

Можно показать, что при выполнении условий (3-11)

UJ

= (log N - (aN[2L) log е при q> a\

д_34

x*~~iZ (log (Ljа) — Lj2aN ln 2 'при q<a,

где L — ширина интервала, в котором равномерно рас­ пределена измеряемая величина X, а — то же для г/Вх.

Ход кривой / х^ 7 в зависимости от N при различ­

ных а показан на рис. 3-3. Количество информации по мере возрастания числа областей квантования монотон­ но увеличивается, тіо скорость этого роста падает. При N'>-оо количество информации стремится к log {Lfa),

Рис. 3-3.

Зависимость количества информации

от числа

областей квантования.

т. е. к постоянному пределу. Выясним, что представляет собой этот предел. Отметим, что дифференциальная без­ условная энтропия X с учетом равномерности распре­ деления равна H(X) = \ogL; аналогично Н (Y^) = \og а. Таким образом,

Нт /*«__*= Hm Ix^ z =log(L/a) =

N —>oo

N->oo

/v

 

=

lo g L - lo g a =

tf ( X ) - tf ( 7 BX).

(3-35)

Отсюда вся формула (3-34) получает простое и есте­

ственное [объяснение: "количество информации

I x<_^.z

равно безусловной энтропии истинного значения измеря­ емой величины Н(Х) минус безусловная энтропия погрешности H(YBX) и минус потери информации на квантование L/{2aN ln 2).

9-301

129

Рассмотренный нами случай .квантования относится к способу аналого-цифрового преобразования с фикси­ рованной начальной фазой. Другими словами, мы пред­ полагали границы областей квантования заранее зафик­

сированными и не зависящими от

х.

Это имеет место

в обычных стрелочных приборах,

в

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

с пространственными кодирующими масками и при по­ разрядном уравновешивании. В преобразователях с про­ межуточным преобразованием в длительность импульса (развертывающего типа) фаза фиксирована только при синхронизации генератора развертки и генератора так­ товых импульсов (образцовой частоты). Если начальная фаза не фиксирована, то вид формулы изменяется, но общий характер зависимости / x*~*z от ^ сохраняется.

Так, для равномерного распределения начальной фазы при соблюдении условий (3-11) в [Л. 3-25] было получе­ но следующее выражение:

 

f 4- ^ - 1оК(У 2 + 1) -

/

(3-36)

9 ln 2 ( 9 - 3 ß + ^ ) ] при K ß < 2;

log (3 -j- 2 I/ 2) при ßs=2,

где $=a/q.

130

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

Перейдем к рассмотрению второго подхода к оценке информационных свойств, связанного с измерением х и Уъх в единицах, равных интервалу квантования q. (Ранее отмечалось, что это относится к анализу цифро­ вых измерительных приборов.)

Считаем, что квантование выполнено безошибочно,

если величина

я + i/nx

лежит

в том же

интервале кван­

тования, что и X .

В

предположении

о

справедливости

условий (3-11)

после ряда

преобразований имеем:

 

 

 

 

 

 

 

 

 

a N

 

 

 

 

 

 

 

 

 

 

4L

 

 

 

I

a N

,

 

a N

при

2L

 

 

 

 

4 L

log-g^-

а

 

 

 

 

 

 

 

 

 

 

 

 

log N

(г — 1) -Кт- log

1

 

 

 

 

 

 

 

 

 

 

a N

 

 

 

a N

1+ * - т - ) і ° е [ т й г ( і - +

 

+ 4

- т У

]

4

к -

L

log

L

 

4

a N

 

~8~

a N

 

 

 

 

при N

2L

 

 

 

I

 

 

 

а

 

 

 

 

 

 

 

 

 

 

 

(3-37)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где k и г имеют тот же смысл,

что и в формуле (3-21).

В частности, три /е = 0 из (3-37)

имеем:

 

Л

:= І0&'- Г — - Ж ПРН N > Ж

<3-38>

Формула (3-38) может ібыть

 

получена

и независимо

от (3-37). Действительно, если

a=rq, где г — любое

целое четное число, то заведомо

имеется

г—1 область

квантования, перекрываемая законом распределения погрешности УВх, т. е. зоной шириной а, и, кроме того, две области, условная вероятность попадания в которые равна q/2a. Соответственно условная вероятность по­ падания в г—1 область равна q/a. Так же, как и в предыдущем случае, все значения условных энтропий вида H(Z\xi) равны между собой и, следовательно, рав-

9*

131

ны среднему значению H(Z\X) в силу упомянутой спе­ цифичности законов распределения ‘Вероятностей X и

Твх. Тогда

= l°S‘N :l) - f 1о&' -Г +

что соответствует формуле (3-38).

Сопоставляя последнюю формулу с (3-34), видим, что характер зависимости количества информации от основ­ ных параметров L, а, N в обоих вы­

 

 

ражениях совпадает. Отличие име­

 

 

ется лишь в коэффициенте перед

 

 

вторым членом.

 

 

разность

 

 

 

На

рис.

3-4 показана

 

 

Г,

 

-

IXN*

Разность

эта

 

 

всегда положительна. Максималь­

 

 

ное

отличие

результатов

имеет

ме­

 

 

сто

при a/q = 0,86 и равно 0,37

дво­

 

 

ичных единиц. Это означает, что при

 

 

/,/<2=100

максимальная

разность

Рис. 3-4.

Разность

имеет место при N = 86.

При возрас­

значений

количества

тании

числа

областей

квантования

информации, получае­

(N-i-oo)

разность уменьшается и

мых при

различных

методах

вычислений.

 

lim (/,

 

x n *—*z n '

 

 

 

 

іѴ - > о о

 

 

 

Реальная область N

 

 

 

 

 

 

(3-39)

для большинства

измерительных

приборов соответствует сравнительно небольшим значе­

ниям разности Ix^ z N~ I x

^.z

. Положительные значения

разности соответствуют

тому

обстоятельству, что со

стрелочного прибора в принципе можно получить больше информации, чем с сопоставленного цифрового. В [Л.3-20]

показано,

что

разность

I x^ z

— /*

всегда неотри­

цательна

при

любых

законах

распределения X и УВх.

Можно показать, что при других комбинациях наиболее употребительных законов распределения .вероятностей для X и Увх характер зависимости I x^ z от ^ сохраняет­ ся. Так, при равномерном законе распределения X и тре-- УГОЛЬНОМ Увх

fi/M2-} -/-1

при K y < Q \

%JВХ(г/вх)— I— yl~2

при 0 < г /< /

132

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