ö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>і(х) может быть пред
ставлена в виде
В результате преобразований формула (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 |
|
где
а = о* |
o*»+Ys; ö= 2Y?; |
c= |
?* + 2y6: |
(°и |
"I- ах' |
|
^д.у |
m = — 2у; /г = — 2^- |
•'я.у |
|
|
|
|
|
|
|
Норміфующии |
множитель |
|
|
|
|
|
|
|
|
со |
)Лг + |
by + |
с(/г |
|
|
|
|
|
/г = |
—Jсо |
|
|
( 6-212) |
|
* + |
ту + /?;/- |
у |
|
Для нахождения интеграла |
в последней формуле представим |
его |
в виде |
|
|
|
|
|
|
|
|
|
|
|
|
Vа + by + су2 dy, |
|
|
|
|
|
(У + Р) (У + |
<7) |
|
|
|
где |
т— Ут2— 4пі |
. |
|
т+ Vт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) |
(a— bp + 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)
Окончательное выражение для а2у имеет вид:
о У тт + 4 -)
(6-216)
24
Следует отметить, что полученный результат хорошо согласу ется с физической стороной явления. Действительно, дисперсия по грешности определения момента времени наступления экстремума
тем значительнее, чем больше дисперсия погрешности дифферен-
о
ниатора яил. у и второй производной о2*”. По мере сглаживания
процесса, т. е. уменьшения коэффициента а, величина а2ѵ растет.
Возвращаясь к формуле (6-124), замечаем, что с ростом дис персии погрешности'определения временного положения экстремума
дисперсия дополнительной погрешности также возрастает, что соответствует физической картине.
Для нахождения экстремальных значений а2 (т) продифференци-
ѵг
руем ее по и и приравняем эту производную нулю:
|
да- (т) |
4а |
Iх+ 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). |
|
|
|
|
|
|
|
|
г |
Jг |
|
|
|
|
|
|
|
|
|
Отсюда |
|
|
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) |
Усредняя по //*, имеем:
|
(6-224) |
а при р* (х) = ехр (— ах2) |
(6-225) |
о |
т. е. диаперсия дополнительной погрешности не зависит от т. Не трудно видеть, что величина дисперсии этой погрешности еще мень ше, чем в 'предыдущем случае. По-видимому, на практике ею можно пренебрегать.
ГЛАВА СЕДЬМАЯ ПРЕДЕЛЬНЫЕ ВОЗМОЖНОСТИ ИЗМЕРЕНИЙ
7-1. О НЕИЗБЕЖНЫХ ПОТЕРЯХ ИНФОРМАЦИИ
ПРИ ИЗМЕРЕНИЯХ
Как указывалось выше, формирование измеритель ного сигнала п его преобразования в измерительном канале неизбежно связаны с потерями информации. Однако эти потерн не являются единственной причиной неполноты информации об истинном значении измеряе мой физической величины. Сам первый этап измерения— взаимодействие чувствительного элемента датчика с изу чаемым физическим телом — уже неизбежно связан с неполнотой отражения, с потерями информации. Здесь сказываются общие закономерности природы.
Действительно, чтобы изучить любое физическое тело, необходимо познать его свойства. Другого пути нет. Свойства же проявляются во взаимодействии с другими физическими телами, в нашем случае—-с датчиками из мерительных устройств. При этом мы неизбежно нару шаем замкнутость исследуемой термодинамической системы, т. е. искажаем какие-то параметры. Эти изме нения могут быть сравнительно малыми, но они неиз бежно имеют место. Так, внося термометр в котел с во дой, мы очень несущественно изменяем интересующую
нас величину — температуру воды. |
Но это |
изменение |
будет значительным, если тот же |
термометр |
вносится |
в каплю воды.
Покажем это на примере. Обозначим удельную теп лоемкость воды Си, ее массу тв и начальную температу ру воды хв■Аналогично для термометра введем обозна-
чения ст, Шт, л'т (в данном случае для упрощения не учитывается, -что термометр может состоять из частей с различной теплоемкостью, учет этого обстоятельства не меняет общей картины, но усложняет формулы). На
основании |
законов |
физики |
конечная |
температура Ѳ |
воды и термометра выразится как |
|
|
|
|
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 распределена по
нормальному закону. Тогда условная энтропия после из мерения
|
Н (-^в 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- Присоединяя вольтметр к источнику, на зажимах вольтметра имеем:
Если измеряемое напряжение — случайная величина, то среднеквадратическое значение его при постоянстве R t>x и Двн равно:
|
О |
A f lX “г * 'В Н |
(7‘9) |
Погрешность, обусловленная конечностью входного |
сопротивления, может быть записана в виде |
|
У |
EgRim |
= сЕа, |
(7-10) |
R a x + R aa
где
R a x + Ra
Дисперсия этой погрешности при случайном характе ре величины с равна:
а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 |
|
где |
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), т. е. чем меньше то*, а также, чем больше интервал усреднения Тщ и дисперсия процед-