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

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

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

ötkyÄä ■

=

r °rf j

f

 

+ 1)

ь

 

^ux" — 2 Kna

J

*>2^ + l

 

X

 

X cos (Щ )

 

Tlrl

2 (•

 

 

 

 

 

 

 

 

 

 

( tt*

\

 

 

TlTl

 

 

Xexpl —4^r) cos Mf/rfco +

/ —

 

 

00

 

 

 

 

 

 

 

X«;* '

to2T2p +

1

exp

G)2

\

 

(6-204)

J COS Ci)// rfco.

 

 

 

Второй из приведенных в формуле (6-204) интегралов тождест­ венно равен нулю. Тогда

С учетом малости величины у и в предположении, что аТ2р<С 1,

можно произвести разложение некоторых функций в формуле (6-205) и далее пренебречь членами с у в третьей степени и выше.

Преобразованные таким образом величины имеют следующий вид:

ехр (—у2а) ~ 1—у2а;

ch(y/Tp) ~ l+ y 2/2T2p.

Разность двух интегралов вероятности можно преобразовать, применяя разложение Тейлора

М г т ’рУ ^Г - ^

“ ) - ®1( щ 7 ^ +

+ у У ^ г ) ъ 2

).

■При больших значениях 'аргумента 4>і(х) может быть пред­

ставлена в виде

293

В результате преобразований формула (6-205) принимает вид

Rux„ (у) ^ 2 а2Л ja (1 - Зар2) + ~ у = ^ ~ ехр +

+ Г2 ехР 2Тѵу

а ^

’г 2Тр\Г

(6-206)

Взаимная корреляционная функция первой и второй производ­

ной процесса на интервале у

равна:

 

 

RX’x" (У) :

â3P* (у)

 

ду3

 

При автокорреляционной

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

вида р.х(т) =

=ехр(—ат2) и преобразованиях, аналогичных принятым при выводе формулы (6-206), имеем:

Rx'x"(u) « 4 s“

а2у (3—2ар2)ехр (—ар2) =

 

 

 

 

 

= 12а| а2у.

 

(6-207)

С учетом (6-206)

и (6-207)

взаимная

корреляция

процессов

Уя.г (і) равна:

 

 

 

 

 

 

 

 

*«„

 

(!')= 2s' r o k 4

x

 

ХеХР [зГрУ

а ( І + 2ГРК a

) . }

 

12°л-“2</-

2 ^ 'р 2

 

ехр\4Ъ ) \

(6-208)

 

 

 

 

 

 

Введем для краткости обозначения

 

 

Y==2a27’5|a + ^r

exp

2 T v

У.

+ ;

 

 

 

 

V a \

Т 2Т

 

 

 

К —

Г2а^а2;

 

I

(6-209)

 

За -

 

Уі

 

 

 

6 = - 2 °1А

 

т ІѴ ~ • ехр А ,Т \

 

 

с учетом которых имеем:

 

 

 

 

 

 

К»

 

к” (и) —'( + ?{/ + Ь/1-

(6-210)

 

Jn..y

 

 

 

 

 

Закон распределения Y запишется в виде

 

 

 

 

 

 

kV a+ by + cy2

( 6- 211)

 

W

( У )

I _ j_

, П у ----- П у 2

 

294

где

а = о*

o*»+Ys; ö= 2Y?;

c=

?* + 2y6:

и

"I- ах'

 

^д.у

m = — 2у; /г = — 2^-

•'я.у

 

 

 

 

 

 

 

Норміфующии

множитель

 

 

 

 

 

 

 

 

со

)Лг +

by +

с(/г

 

 

 

 

 

/г =

Jсо

 

 

( 6-212)

 

* +

ту + /?;/-

у

 

Для нахождения интеграла

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

его

в виде

 

 

 

 

 

 

 

 

 

 

 

 

+ by + су2 dy,

 

 

 

 

 

+ Р) (У +

<7)

 

 

 

где

тУт2— 4пі

.

 

т+ 2 — -іпі

 

 

Р

 

 

 

-=--------2п----------- (І = -

 

2п

 

 

Тогда

 

 

 

 

 

 

пк

 

 

$ х= к У с In [2 Ус ('і +

 

 

 

 

 

 

Ьу + су2) + 2cy + b\-\-y. acin2— 4пі

X

X [(« +

bq + cq=) In

2'?- + by + bq + 2

Va

(a -f by + cy2)

 

 

 

y + q

 

 

 

 

2a + by + Ö/7 -f- Va (a +

by + cy2)

(6-213)

(abp + cp2) In

 

 

'/ + />

 

.

 

 

 

 

 

 

 

 

Решение уравнения вида к =

дает искомое условие норми­

ровки. Дисперсия •инструментальной 'погрешности, усредненная по всем У, при любом виде закона распределения равна:

2

z \ = 2 I z~

- р* ІТ)

 

ач

Т ) ааи + Т2

Т2

(6-214)

где принято ог2о = сг2:с, что справедливо при определенных условиях, оговоренных в § 6-5. .

В случае, когда w(y) имеет вид (6-203),

2

 

+ Ьу + су2 dy — — =

аи

 

 

к (2су + Ь) У а2 + Ьу + су2

 

(4ас — Ь)2 к

4пс

1_

со

8пс У с

In [2

Ус(а+Ьу + су2)+ 2су + Ь}

П

 

 

(6-215)

295

Окончательное выражение для а2у имеет вид:

о У тт + 4 -)

(6-216)

24

Следует отметить, что полученный результат хорошо согласу­ ется с физической стороной явления. Действительно, дисперсия по­ грешности определения момента времени наступления экстремума

тем значительнее, чем больше дисперсия погрешности дифферен-

о

ниатора яил. у и второй производной о2*”. По мере сглаживания

процесса, т. е. уменьшения коэффициента а, величина а2ѵ растет.

Возвращаясь к формуле (6-124), замечаем, что с ростом дис­ персии погрешности'определения временного положения экстремума

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

Для нахождения экстремальных значений а2 (т) продифференци-

ѵг

руем ее по и и приравняем эту производную нулю:

 

да- (т)

+ aXay — “7'ад1=

 

 

 

 

 

дх

 

 

 

 

Решая это уравнение относительно т, находим,

что от имеет эк-

стремум в точке

 

 

 

 

 

 

 

 

у%

 

 

 

 

 

 

 

 

 

 

 

 

%= Т ■

+

Ѵ

 

 

 

(6-217)

 

 

 

 

 

 

аор

 

 

 

 

В этой точке имеет место минимум, так как

 

 

 

 

 

д2а\ (т)

Aal

 

 

 

 

 

 

 

 

------ W [1 + ^ ];> о .

 

 

 

 

 

да,2

 

 

 

 

 

 

 

 

Этот минимум равен:

(т)'= 2<ХДК

(2 + а°1)

 

 

 

 

 

min а2

 

 

(6-218)

 

 

 

 

 

 

1+*'1

 

 

 

Поскольку ход кривой"а2

(г) на интервале (О,

Т)

представляет

собой параболу второго порядка относительно т и кривая

имеет ми­

нимум внутри

интервала,

следует

ожидать

появления

наибольшего

значения на интервале либо в

точке т = 0,

либо

і = Т .

Нетрудно

получить, что

 

 

 

 

 

 

.2 п

 

 

 

 

 

 

 

 

 

 

(6-219)

 

 

 

 

аа^ + (і

■ е - аТг)1Т г

 

 

 

1

(П '=:2а2

1 + (1 — е ~ аТг)

 

 

(6- 220)

т- «• а1 (П >«2

(0).

 

 

 

 

 

 

 

 

г

 

 

 

 

 

 

 

 

 

Отсюда

 

 

sup ° l

M = 2 e *

2в *(1 + « < $ .

 

*5

Ъ

( 6 ‘ 2 2 1 )

В предыдущем параграфе было показано, что методическая погрешность была .минимальна при т=0 я х=Т , т. е. именно там,

где максимальна инструментальная погрешность. Принимая во вни­ мание, что последняя является .величиной второго порядка малости, можно рекомендовать учитывать ее при расчете суммарной погреш­ ности только на концах интервала интерполяции.

Рис. 6-7. К возникновению инструментальной погрешности из-за неточности фиксации мо­ мента экстремума.

Рассмотрим также инструментальную погрешность, обусловлен­ ную конечной точностью измерения и фиксации моментов времени измерений. На рис. 6-7 эти погрешности обозначены у *і и у*2. Зна­

чения этих погрешностей обусловлены типом применяемого измери­ теля времени. Естественно, что анализ погрешностей часовых меха­ низмов выходит за рамки данной работы. Поэтому ограничимся исследованием .влияния этой погрешности на погрешность интер­ поляции. С учетом этой погрешности интерполирующая прямая име­ ет вид:

z 2 (*)= *„ (0) - T^ ^ y*s [X. (Т) ха (0)].

(6-222)

где у* 1и у*2— погрешности в определении моментов времени, к ко­

торым следует отнести экстремумы. Учитывая, что часовым меха­ низмам свойственна накапливающаяся погрешность, можно принять:

\у *і— У*Л < Т .

 

Математическое ожидание Уг (т) равно нулю, а

дисперсия

2

2у*в\

 

С®! У*) —

"гг П — ?х (Г)]-

(6-223)

297

Усредняя по //*, имеем:

 

(6-224)

а при р* (х) = ехр (— ах2)

(6-225)

о

т. е. диаперсия дополнительной погрешности не зависит от т. Не­ трудно видеть, что величина дисперсии этой погрешности еще мень­ ше, чем в 'предыдущем случае. По-видимому, на практике ею можно пренебрегать.

ГЛАВА СЕДЬМАЯ ПРЕДЕЛЬНЫЕ ВОЗМОЖНОСТИ ИЗМЕРЕНИЙ

7-1. О НЕИЗБЕЖНЫХ ПОТЕРЯХ ИНФОРМАЦИИ

ПРИ ИЗМЕРЕНИЯХ

Как указывалось выше, формирование измеритель­ ного сигнала п его преобразования в измерительном канале неизбежно связаны с потерями информации. Однако эти потерн не являются единственной причиной неполноты информации об истинном значении измеряе­ мой физической величины. Сам первый этап измерения— взаимодействие чувствительного элемента датчика с изу­ чаемым физическим телом — уже неизбежно связан с неполнотой отражения, с потерями информации. Здесь сказываются общие закономерности природы.

Действительно, чтобы изучить любое физическое тело, необходимо познать его свойства. Другого пути нет. Свойства же проявляются во взаимодействии с другими физическими телами, в нашем случае—-с датчиками из­ мерительных устройств. При этом мы неизбежно нару­ шаем замкнутость исследуемой термодинамической системы, т. е. искажаем какие-то параметры. Эти изме­ нения могут быть сравнительно малыми, но они неиз­ бежно имеют место. Так, внося термометр в котел с во­ дой, мы очень несущественно изменяем интересующую

нас величину — температуру воды.

Но это

изменение

будет значительным, если тот же

термометр

вносится

в каплю воды.

Покажем это на примере. Обозначим удельную теп­ лоемкость воды Си, ее массу тв и начальную температу­ ру воды хв■Аналогично для термометра введем обозна-

298

чения ст, Шт, л'т (в данном случае для упрощения не учитывается, -что термометр может состоять из частей с различной теплоемкостью, учет этого обстоятельства не меняет общей картины, но усложняет формулы). На

основании

законов

физики

конечная

температура Ѳ

воды и термометра выразится как

 

 

 

 

Q £ц//гпхп-|- сттт$т

 

IV

 

 

cTmT+ cBmB

 

'

Разность между Ѳ и начальной температурой воды,

определяющая погрешность измерения,

 

 

 

У== О

в

сттт+ свтв

 

( 7 - 2 )

 

 

 

 

 

 

а дисперсия

этой величины может быть

записана в форме

 

 

(1 +я=) ,

где а--

свтв

( 7 - 3 )

 

 

сгтТ

 

 

 

 

 

 

Безусловная дифференциальная энтропия температу­ ры воды в предположении о нормальности распределе­ ния вероятностей значения измеряемой величины в дан­ ном диапазоне равна:

Н (Х Я) = log

хп

Предположим, что начальная температура воды име­ ет плотность вероятностей, также описываемую с по­ мощью гауссовской кривой. Тогда условная энтропия после измерения

tf(X B|0) =

log

У 2т.е

стгпт

1

(7-4)

Cntnв

I >

 

 

 

** I

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

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

х в имеем в

виде

 

 

 

 

 

в

<sXJTcrmT

 

 

( 7 - 5 )

 

 

 

причем масса воды предполагается постоянной.

На практике при измерении температуры в большин­ стве случаев нельзя говорить о постоянстве массы. Для , учета этого обстоятельства введем обозначения Ь=а~і. Предположим теперь, чтр. величина b распределена по

299,

нормальному закону. Тогда условная энтропия после из­ мерения

 

Н (-^в 1Ѳ)= log {]/2ite X

 

X

(Ѳ -

*т)2 Н-<т К +

(£)2] },

(7-6)

а частное количество

информации,

доставляемое зна-

чением 0,

 

 

 

 

' Ѳч—►Х

= log-

а\ (Ѳ— хт)2 + а2^ fa“ (Б)2]

(7-7)

 

] /"

 

Нетрудно видеть, что при а2в= 0

формулы

(7-6) и

(7-7) превращаются в (7-4) и (7-5).

 

 

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

ры термометра

величина

также уменьшается.

Аналогично

можно

оценить количество инфор­

мации при измерении напряжения Е а некоторого источ­ ника постоянного тока, имеющего внутреннее сопротив­ ление Двн, с помощью вольтметра с входным сопротив­ лением Rbx- Присоединяя вольтметр к источнику, на зажимах вольтметра имеем:

U

(7-8)

 

R ax “Ь

Если измеряемое напряжение — случайная величина, то среднеквадратическое значение его при постоянстве R t>x и Двн равно:

 

О

A f lX “г * 'В Н

(7‘9)

Погрешность, обусловленная конечностью входного

сопротивления, может быть записана в виде

 

У

EgRim

= сЕа,

(7-10)

R a x + R aa

где

R a x + Ra

300

Дисперсия этой погрешности при случайном характе­ ре величины с равна:

а1ЁС,Е2 ~ К ес “Ь (С)2 вЕа + (Еа)23“ +

2 [СЯ 2 + 2CEaRcE -{-EaRc,E ] -f- (С)“{ЕаУ-

(7-11)

 

cCq

а

а

 

В случае нормального закона распределения С и Еа

выражение

(7-й 1) может

быть

преобразовано к

виду

< = К

+ (С )2] 4 а + К

+ (С)2] (.ЕаУ+ (С)2(ЕаУ+

 

 

+ (RcE + 4 CEa)R cE.

(7-12)

 

а

 

а

 

При этом безусловная энтропия

Н (Еа) =

log (]/г2тсЕаЕ ).

(7-13)

 

о

 

Условная энтропия при нормальности Еа

и d — c~l

равна:

 

 

Н (ЕаIV) =

log (V2-ке Uod),

(7-14)

а частное значение количества информации, доставля­ емое значением U,

ІЕ ~ ѵ = ^ аЕІѴ°л).

(7-15)

аа

Напомним, что d = \ + R BXIRBn.- Последняя формула получает ясное истолкование. Количество информации тем больше, чем больше дисперсия (изменчивость) изме­ ряемой величины и чем меньше разброс возможных со­ отношений между Двх и Rbb. Наличие величины U в зна­ менателе формулы (7-15) имеет тот же смысл, что и вве­ дение понятия относительных погрешностей.

Другим источником потерь на начальном этапе изме­ рений является тот факт, что измерения всегда по необ­ ходимости интегральны. Любой прибор измеряет не мгновенное значение, а усредненное за какой-то, пусть короткий, интервал-времени. Это приводит к необ­ ратимым потерям информации. Действительно, пусть измеряется некоторый стационарный случайный процесс X(t) на протяжении времени Гц. Если бы -мы могли получить полную информацию о процессе, то имели бы без учета погрешностей измерения [Л. 1-48]

/ . : = £ - ял-*),

(7-16)

LT0x

 

301

где

to.*— интервал корреляции X(t);

Н0(Х) — безуслов­

ная

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

отсчет. (Здесь для

простоты расчета считается,

что Ти кратно то.*.) В дейст­

вительности мы имеем:

 

х ср = ~

I* Л- (/) dt,

 

о

содержащее, как будет показано ниже, меньшее количе­ ство информации.

Представляя

 

 

 

k

где

 

 

<=і

ft =

 

ent [TJxaX),

 

 

имеем дисперсию в виде

 

 

9

9

 

 

(7-17)

 

 

 

т. е. а2 <

о2. Отсюда

 

следует, что энтропия на одно

^ср

х

 

 

измерение при усреднении меньше, чем Яо(Х), и анало­

гичным образом количество информации

меньше, чем Іи

Оценим это уменьшение на примере нормального

процесса, для которого

 

 

%ох

log (/2т«? ах),

(7-18)

 

 

а при усреднении

 

 

/ = - log

-,)■

(7-19)

 

Методические потери информации (аналогично мето­

дической погрешности)

 

 

Д / , =

 

(7-20)

Естественно, что эти потери тем больше, чем быстрее изменяется процесс X(t), т. е. чем меньше то*, а также, чем больше интервал усреднения Тщ и дисперсия процед-

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