![](/user_photo/_userpic.png)
книги из ГПНТБ / Кавалеров, Г. И. Введение в информационную теорию измерений
.pdfПерейдем далее к рассмотрению погрешностей интер поляции случайного процесса по дискретным во времени измерениям с учетом квантования в момент измерения. Используя методику, принятую в гл. 4, при ступенчатой интерпсляции имеем:
1 в (ік -)- т) = X (tu -|- т) — |
(th) — Y'вх(tu) — |
(tu); |
(5-34) |
|||||||
a l ^ h |
"О — а’ ( t k + х) + а1 Ѵ г ) + % |
( t h ) + |
|
|||||||
( t h ) |
- - |
% R x ( t h “ Ь т > til) |
— 2 /? |
( t i i ~ \ - t , |
t h ) |
— |
||||
j k |
|
|
|
|
|
*y oX |
|
|
|
|
~ 2 R vy(tk + i , t h) + 2R |
(tu,tk) + |
|
|
|||||||
|
‘ |
|
|
|
' -'ox |
|
|
|
||
4“ 2/? |
(tu, |
tu)-\-2R |
(tu, tu), |
|
(5-35) |
|||||
где О^тг^гГ. |
|
' ‘ 'K |
|
|
' / в Х*/ |
к |
|
|
|
|
и |
Увх(0— стационарные |
и стационарно |
||||||||
Если X(t) |
||||||||||
связанные случайные процессы и г/вх(0 = 0, то |
|
|
||||||||
|
|
1 |
(tu -f- т) I = I Y K(іи) |; |
|
|
(5-36) |
||||
|
+ |
’) = |
24 1 1 — Р* <-)] + |
я'ца + |
’», + |
|
||||
+ 2 К » Ол„ (0) - |
DA ( - |
’)] + |
2 К ,К. (0) - |
2R„К ( - |
’Я + |
|||||
|
|
|
+ 2R |
(0). |
|
' |
• |
(5-37) |
||
|
|
|
|
-/вх-/ к |
|
|
|
|
Из входящих в формулу (5-37) величин новыми являются только взанмнокорреляцпонные функции по грешности дискретности и x(t) и yBx(t) соответственно:
Rsu. СО= |
R* Ь) + |
Rxe.. (**) - |
R**СО; |
|
|||
|
|
|
|
|
|
(0 ); |
|
V |
( |
° |
) ä=, . ( |
° |
) +ä w |
. (Р) = |
(5-38) |
= |
N - |
1 |
(І+І)? |
|
|
||
s |
гур (зу) |
f |
uw (и I zy) du. |
|
|||
|
1=1 |
|
“7 |
|
|
1 |
Можно воспользоваться также выражением для взаимнокорреляцпонной функции сигнала и шумов квантования, полученным в |[Л. 1-3]. Отметим, что разли чие постановок задали в [Л. 1-3] и нашем случае состоит
13—301 |
193 |
в учете нами погрешностей і/вх. В наших обозначениях без учета //вх имеем:
D /т\ |
л_ |
00 |
_L ( |
1\j+ i |
('J' 2я//у, т) |
|
|
|
|
||||||
К х У к і } ~ |
^ |
Ъ |
|
І { |
} L |
д » |
J u = o ’ |
|
|
/=—оо |
|
|
|
|
где 02.-V— двумерная характеристическая функция кван туемого сигнала с произвольным вещественным параме тром и. Используя это выражение, в предположении о нормальности X(t) и Увх(/) и их некоррелированности имеем:
R |
|
(т) =о 2 [а2 рх (т) + а2 |
ри M IX |
|
|||
|
(•Ѵ+(/пх) !/кw |
I * ' * ' |
' ‘ |
„ |
|
||
|
"вх "вх |
|
|||||
|
00 |
|
|
2л=/= |
|
+ 0wdX) |
|
|
Х $ ] ( ~ |
l)j exp |
|
(5-39) |
|||
|
|
Ф |
. |
||||
|
/=1 |
|
|
|
|
||
|
|
|
|
|
|
|
|
Можно показать, что модули членов ряда, стоящего |
|||||||
под знаком суммы в формуле |
(5-39), |
удовлетворительно |
аппроксимируются гиперболической зависимостью, так как ox2^>q2-
Кроме |
того, |
опираясь |
на результаты, |
полученные в |
||||
[Л. 1-3], |
нетрудно выявить малость R |
(-с) |
по |
сравнению |
||||
с |
|
|
Тогда |
|
|
|
|
|
О |
|
і ) г « С і (т) -j- а~ -j- а“ |
— 2 (з" + |
о" |
) X |
|||
(th + |
||||||||
|
|
|
exp |
2"2/2(«X + ° £ ) |
|
(5-40) |
||
|
|
|
|
|
|
|
||
|
|
/=1 |
|
|
|
|
|
|
где Ск(т )— структурная |
функция |
X(f) |
(см. |
формулу |
||||
(1-13)]. |
|
по |
формуле |
(5-40) показывают, что при |
||||
Расчеты |
||||||||
большом |
числе областей |
квантования (Д1> |
100) послед |
ним членом формулы можно заведомо пренебречь. Подобным же образом может быть определена по
грешность при других видах интерполяции.
Отметим две особенности анализа дисперсии погреш ности квантования в динамике. Во-первых, при этом анализе не играет роли, каков порядок обобщенного марковского процесса, полученного после дискретиза ции. Во-вторых, дисперсия погрешности восстановления
194
дискретизированного и проквантованного процесса суще ственно зависит от применяемого метода интерполяции. В этом отношении расчет дисперсий погрешности суще ственно отличается от информационного анализа, являясь более простой, но более частной задачей
Даже при стационарности Х(1) погрешность интер поляции представляет собой нестационарный процесс.
В ряде случаев можно указать оптимальные соотно шения между интервалом времени между измерениями Т и шагом квантования по уровню q. Например, если измерения следуют непосредственно друг за другом че рез время, необходимое для считывания показаний, то имеется конкретная связь Т с q. Для развертывающих
АЦП и время-импульсных цифровых вольтметров |
|
Т — TqN, |
(5-41) |
где Т,, определяется быстродействием нуль-органа. Если применяется ступенчатая интерполяция, то при условии монотонного спада автокорреляционной функции мини мум дисперсии сг/в.макс достигается при Т— Т0, удовле творяющем условиям:
д?х ( Т ) |
|
1 |
|
|
д Т |
Е- |
124 |
|
|
о |
|
|||
|
|
|||
|
ч |
|
||
|
чІ- II |
|
|
|
d=P, ( Г ) |
|
|
L-T~q |
|
д Г - 1 |
|
< |
|
|
|
|г= Г 0 |
4 а х |
) |
|
|
|
|
|
В частности, при рЦт) — ехр(— а-г) имеем:
Та& ехр ( аГд) — JLTqj !2а^ .
(5-42)
(5-43)
Для преобразователей с поразрядным уравновешива нием при применении двоичного кода
T=Tqlog N, |
(5-44) |
откуда минимум максимальном дисперсии погрешности восстановления обеспечивается Т.= Т0, для которого
dpХ (Т) |
г=г0 |
L2 In 2 |
(5-45) |
|
дТ |
6о£ Гд-22Го/7« |
|||
|
||||
|
|
В рассмотренные примерах фактически показывается целесообразная связь быстродействия аппаратуры и чув-
13* |
195 |
ствительности. Интересно отметить, что статистические характеристики YBX(t) не влияют на эти соотношения при ступенчатой интерполяции.
5-3. СПОРАДИЧЕСКАЯ ДИСКРЕТИЗАЦИЯ ПРИ КВАНТОВАНИИ ПО УРОВНЮ
Большинство существующих измерительных информационных систем построено на основе цикличес кого опроса, рассмотренного выше. В некоторых систе мах применяют измерение по вызову. Вместе с тем для решения ряда задач перспективным представляется спо радический принцип дискретизации.
Под спорадическим опросом [Л. 5-8] понимается изме рение в случайные моменты времени, определяемые хо дом конкретной реализации случайного процесса.
Примером воплощения такого принципа является рас сматриваемый ниже метод дискретизации при кванто вании по уровню, реализуемый в АЦП следящего дей ствия.
Пусть, как и раньше, истинное значение измеряемой
величины |
X(t) |
есть стационарный |
случайный процесс, |
|
причем |
0 |
Требуется квантовать |
величину X |
|
с шагом |
q = L/N = const. На первом |
этапе |
рассмотрения |
погрешности измерения в момент снятия отсчета будем считать пренебрежимо малыми. В § 5-1 было показано [см. (5-13)], что при определенных условиях среднее чис ло пересечений в единицу времени процессом Х(і) гра ниц областей квантования по уровню равно:
Дя == ] / 2%а_ѵ, -J-.
В спорадической системе каждое такое пересечение является как бы заявкой на обслуживание [Л. 5-10]. Дру гими словами, в системе не возникает никакого сообще ния, пока X находится внутри кванта (г—l)q<x<iq, где і=1, 2, ..., N. В момент времени th>определяемый из ра венства x(til)=iq, возникает сообщение. Время обслужи вания одного сообщения, т. е. время считывания показа ний (и, если требуется, математической обработки и ре гистрации), обозначим через тп. Время пребывания процесса X{t) в границах ((і—\)q, iq] — случайная ве личина, обозначаемая как T(iq).
196
При нормальности процесса X (/) среднее значение
Т (iq) = Г Ѵ
2я
exp
(0)
, ) q -
(/<7 —
T*iq == |
? - |
T% = |
|
Л')2 |
,n ,/g2 |
(2 /- |
\ ) - X , \ _ j |
|
exp |
|
+ |
|
|
2^ |
+ Ф, {{ІЯ- X) ox '] - Ф, [(/<7 - q - X) зв '] j , (5-46)
где T*a_-l)q и T*iq — длительности выбросов соответст венно за уровень (i— \)q и iq; Ф, — форма интеграла вероятности (см. приложение 1). При больших значениях N последняя формула преобразуется к виду
Т (iq) |
Ѵ2к ( iq — X |
exp ( i q - x y |
+ |
|
К |
2<- |
|
|
|
|
|
|
exp |
(fr-ff)« |
(5-47) |
|
V2k |
2®?. |
|
Точность приближенного равенства в формуле (5-47) весьма велика и зависит от N. В частности, при іѴ= 100 и правильном выборе соотношения между ох й L по грешность не превышает 1,5%. Средняя по всем уров ням длительность интервалов времени между ближай шими пересечениями границ квантов
ТСѵ = ^ Т ( і д)р(гі) = ) Г \ |
(5-48) |
і=і |
|
где |
|
р fa) = Ф, \(iq —X) а“ 1] — Ф, [(iq — q — LY) ох]. |
(5-49) |
Если T(iq)>Тп, то устройство к моменту возникнове ния новой заявки успевает освободиться от обслужива ния предыдущей. Если Т (iq) < т1Ь то в момент получения второй заявки система занята. В этом последнем случае возможны следующие варианты: либо система способна запомнить заявку, либо память отсутствует. В первом случае заявка, возникшая в то время, когда система была занята, спустя некоторое время і3 дождется момен та освобождения системы и вызовет новое сообщение; во втором случае сообщение будет потеряно. Возможно воз-
197
никновение белое чем одной заявки за время занятости системы то. Тогда для системы с запоминанием возмож на память по первому из возникших сообщении, по по следнему или по нескольким.
Успешность функционирования системы может быть оценена по критериям, принятым в теории массового обслуживания, например среднему времени задержки и вероятности потери сообщения (см., например, (Л. 5-4, 5-5]). Однако в данном случае эти критерии хотя и полу чили некоторое применение (см., например, [Л. 5-10, 5-11]), но отражают менее существенные особенности ИИ С по сравнению с развиваемыми ниже информацион ными оценками пли дисперсией погрешности восстанов ления процесса.
Следует оговорить, что сам по себе алгоритм спора дической дискретизации предполагает лишь процедуру определения моментов времени измерения {//,}. В эти мо менты может измеряться как полное значение х (//,), так и приращения. В последнем случае, если требуется опре
делить полное значение х (4), |
то имеет место накопле |
|
ние |
погрешностей (например, |
из-за сбоев). Кроме того, |
для |
восстановления процесса |
А'(/д.) по дискретным во |
времени измерениям могут применяться различные ме тоды интерполяции, рассмотренные в гл. 4.
Рассмотрим прежде всего систему без потерь сооб щений. Теоретически в любой системе с конечным и от личным от нуля временем обслуживания т,і, если искус ственно не наложить на нее условие иі(7\) = 0 при Ті,< <т„, всегда имеются отличные от нуля математическое ожидание времени задержки и вероятность потери сооб щений. При прочих равных условиях величина их тем меньше, чем меньше N и и чем медленнее изменяется процесс. Поэтому при наблюдении медленных процессов (например, давление в магистральном газопроводе) спо радическая система может работь практически без по терь сообщения.
Другими примерами возможных систем без потерь являются измерение температур тел с большой тепловой инерцией и наблюдение за осадкой грунтов в лаборато риях строительных учреждений.
Если для восстановления исходного непрерывного процесса после дискретизации вышеописанным образом применяется ступенчатая интерполяция, то в соответст вии с [Л. 2-9] автокорреляционная функция ня пыхопе
198
равна:
А/ |
I |
iq |
PS W |
|
,('с) = |
(«) |
(5-50) |
||
Ф1 |
|
/г! |
и
,/=1
где ф(”)і — /г-я производная от интеграла вероятностей. Усредненная по всем интервалам квантования дис персия погрешности интерполяции приблизительно,равна q2l 12, а автокорреляционная функция погрешности
интерполяции в соответствии с [Л. 1-3]
R |
2и2 U іѵ ехр |
J |
4/1=71= 4 [ 1 — Рх (t)J |
1 |
<72 |
||
|
/1=1 |
|
|
причем формулы (5-50) и (5-51) выведены в предполо жении о нормальности процесса X(t).
Отметим, что во всех приведенных в данной главе формулах возможна замена х(і) на 2вх(0 =х(і) +у ПХ(0 . что позволяет увязать результаты данного параграфа с предыдущим материалом.
Перейдем далее к информационному описанию систе мы без потерь. Отличительной особенностью информа ционного анализа является возможность описания пото ка на выходе системы без учета методов интерполяции. В этом смысле информационный анализ описывает как бы предельные возможности восстановления. В данном случае, как во всякой спорадической системе, информа цию песет сам момент появления сообщения и характер сообщения. В данном случае характер сообщения — это знак сообщения, т. е. знак приращения. Первую состав ляющую количества информации назовем фазовой, вто рую— амплитудной. Для системы без потерь сообщений амплитудная составляющая / а не превышает 1 дв. ед.
Рассмотрим количество информации в фазе поступле ния измерения. Естественно, что эта величина /ф отлич на от пуля только в том случае, если интервалы времени между измерениями Т случайны. Предположим, что по ток заявок является простейшим, т. е. ординарен, стацио нарен и не имеет последействия. Запишем прежде всего дифференциальную энтропию величины Т в виде
СО
// (Г) = - J да (Г) log W (Г) dT, |
(5-52) |
где w ( T ) — плотность |
вероятностей |
Т. Так как |
для |
пуассоновского потока |
с параметром |
Хх плотность |
ве |
роятностей |
|
|
|
w(T)=XXexp (—ХхТ), |
|
||
Н (Т) =log(eA,v). |
(5-53) |
В связи с тем, что любой измеритель времени имеет квантующее устройство, то фактически может быть заре гистрирована не непрерывная величина Г, а результат ее аналого-цифрового преобразования (время дискрети зации) Г д . Обозначим величину кванта через То- Анало говой ошибкой измерителя времени при измерении ма лых интервалов, как правило, можно пренебречь. Тогда вся ошибка измерения фазы поступления, сводится к по грешности дискретности Ут- Среднее значение условной энтропии равно:
со |
Т0 |
|
|
|
Н (Т I Гд) == — £ Р (Ги) J W (ут) log W (ут) dyT, |
(5-54) |
|||
А=1 |
О |
|
|
|
где |
кТ |
|
|
|
Р (Ги) = |
w(T)dT. |
|
|
|
'j |
|
(5-55) |
||
|
(*—I) То |
|
|
|
В частноегп, для пуассоновского потока в силу от |
||||
сутствия последействия |
|
|
|
|
w(yT) =Хх ехр [—І х{То—ут)1 |
(5-5G) |
|||
00 |
|
|
|
|
Учитывая, что ^ , р ( Т к) — ], |
и меняя |
порядок суммн- |
||
А=1 |
|
|
|
|
рования и интегрирования, имеем: |
|
|
||
То |
|
|
|
|
Н (Т I Гд) = J Я, exp I - Я, (Г, - |
ут)] X |
|
||
О |
|
|
|
|
X log {Я* ехр[—Я,(Г0—г/г)]} dyT=[ 1 - exp ( - l xT0)] l o g — |
|
|
''■je |
- Е Т ехР < - ;,-*7'»>- |
<6-67> |
При этом предполагалось, что величина уд- достаточ но мала, чтобы вероятностью появления нового измере ния за это время можно было пренебречь, т. е. 1—ро<С
< 1.
200
Количество информации в фазе поступления измере
ния при пуассоновском потоке |
|
Іф = Н(Т)—Н(Т\ Тя) =ехр (—кхТо) X |
|
X Iogfexp (1 +ЯЖГ0)1_ІЛ-]. |
(5-58) |
Нетрудно видеть, что при увеличении Хх фазовая |
|
информация уменьшается, а |
|
lim / ф — log- (е/Іх) = Н (7), |
(5-59) |
Го->0 |
|
что и следовало ожидать.
Таким образом, скорость передачи информации в си
стеме без потерь сообщений |
|
|
Е —АД/а + Ліі) |
+ехр(—АдТо) X |
|
X logfexp (1 +КхТо) АгД]. |
(5-60) |
|
В реальных спорадических системах всегда |
сущест |
вует отличная от нуля вероятность потери заявки, как это будет показано ниже. Однако полученные в данном
параграфе |
|
результаты |
|
|
оказываются |
|
полезными |
|
|
в ходе дальнейших выво |
|
|||
дов. |
|
|
|
|
|
Исследование систем |
|
||
с потерями начнем с уче |
|
|||
та |
конечного |
времени |
|
|
обслуживания |
заявки. |
|
||
Пусть время, |
необходи |
|
||
мое для измерения в спо |
|
|||
радической системе, хиф' |
|
|||
ф.О. Если в системе нет |
|
|||
накопителя (памяти) за |
|
|||
явок, то все заявки, воз |
Рис. 5-6. Временной график рабо |
|||
никающие |
за |
время тп |
ты спорадической одиоканальной |
|
с |
момента |
предыдущего |
системы без памяти. |
измерения, будут потеря ны. Следующее измерение будет ’ вызвано первой же
заявкой, возникшей после освобождения от обслужива ния (рис. 5-6). При этом характер каждого измерения меняется — счет единичных приращений невозможен. Возможно либо измерение приращений, что связано с накоплением погрешности, либо измерение полных значений. Допустим, что реализуется второй вариант.
201
Количество информации в одном сообщении по-пре>к- нему будем описывать через амплитудную и фазовую компоненты. Нетрудно видеть, что
/а = Е Р „ К )//д ,(« + 1 ), |
(5-61) |
H=U |
|
где рп(ті,)— вероятность появления п пересечений |
гра |
ниц областей квантования по уровню процессом А'(/) за время обслуживания одной заявки; HN{n-\-\) опреде ляется формулой (5-23).
Компонента /ф определяется формулой (5-53), если поток пересечений простейший. При соблюдении усло вия пуассоновского потока имеем также: •
Рп Ы = Р % Г еХР -(- А*хи)- |
(5-62) |
Для нахождения скорости получения информации определим также интенсивность потока обслуживаемых заявок Запишем:
Ѵ= Ы 1 —Рп), |
(5-63) |
где рп — вероятность пеобслуживапия заявки, т. е. |
появ |
ления такого пересечения, которое не послужит основа нием для измерения.
Вероятность рп можно определить как |
|
рц= 1 1 (1 + Н ) - \ |
(5-64) |
где Н — математическое ожидание числа потерянных за явок за время тп.
Для пуассоновского потока |
|
||
Ы = £ 11Рп Ы = Яя-Еи. |
(5-65) |
||
/ 1=1 |
|
|
|
Подставляя (5-65) в (5-64) и (5-64) в (5-63), имеем: |
|||
X’ |
К |
(5-66) |
|
1+ Ххіа |
|||
|
|
202