Файл: Калинчук, Б. А. Анализаторы инфразвуковых случайных процессов.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 19.10.2024
Просмотров: 88
Скачиваний: 0
где Т 0 — длительность исследуемой реализации, тх — переменная интегрирования.
Анализ показывает, что составляющая погрешности измерения, связанная с ограниченностью интервала интегрирования Т 0 — т, для произвольных законов распределения отличается от известных результатов для нормального распределения.
Пусть стационарный случайный процесс имеет закон распреде
ления ординат |
W (а*), |
отличный от нормального. Тогда |
|||
|
2 |
|
Rxx (т) |
7*„—X |
|
а2 (т) = у \ѵ Т 0 |
|
■^“ 11 |
IR.XX(Ті) Rxx (т 1 + |
||
т |
|
R xx (0 ) |
|||
|
|
|
|
Г„—т |
|
+ |
т)] гітх + |
T — % [Rxx (Tl) + Rxx (T1 “Ь |
|||
+ т ) К х х Ы — т ) ] * ! , |
|
|
(1-5) |
||
где yw = ^ — |
3, |
—эксцесс, |
Ox = l - D x x , |
|
|
CT., |
|
-fc o |
(x —inx)i W(x)dx, tnx |
M (A (0). |
|
(-4= |
J |
||||
|
|
—00 |
|
|
|
Заметим, что из выражения (1-5) как частный случай (у^- = 0) может быть получена формула (1-4).
Таким образом, погрешность от несоответствия принятой ги потезы о гауссовом виде распределения истинному закону имеет вид:
|
|
Rxx (т) |
Г0—X |
К М] «7 (л) |
2Уп/ |
Г |
|
Т о - Т |
R (0) |
J |
|
|
|
|
О |
X Rxx (Ti -I- т)] dxv ( 1-6)
При т = 0 формулы (1-4) и (1-5) дают оценки статистической дисперсии.
Другой метод определения минимального допустимого с точки зрения заданной погрешности интервала интегрирования Тп = = Т — ттах заключается в аппроксимации на бесконечном интер вале стационарного случайного процесса х (t) линейной комбина цией синусоидальных сигналов
A'v (t) = av sin (оѵЧФѵ) |
(1-7) |
со случайной амплитудой и фазой [97, 6, 104, 1051. Из-за отсутст вия взаимной корреляции между сигналами вида (1-7) с различ ными V корреляционная функция исходного процесса может быть записана в виде:
Rxx^) = y ,R xx( т), |
(1-8) |
V |
|
9
где |
|
|
|
|
|
Я.ѵ.ѵ(т) = 1іт-^- |
J a2vsin (w^t + cpv) sin [cov (/ + |
г) -f фѵ] dt. (1-9) |
|||
T -+CO A l |
— T |
|
|
|
|
Предельный переход в (1-9) обеспечивает получение |
|||||
|
Rxx (т) = |
а2 |
|
|
(1-10) |
|
cos соѵг. |
|
|||
В практических |
случаях (Т =f=со) |
имеем: |
|
||
а2 |
|
|
|
|
|
Rxx (т) = — cos соѵт— |
2Фу) + |
sin ( ( ѴГ - |
( ѵт — 2фѵ) |
||
°ѵ sin (2мѵг + мѵт + |
|||||
4 |
|
|
2 |
0 |
0 |
|
2озѵ |
|
|
||
Т |
|
|
|
|
|
Погрешность измерения корреляционной функции ѵ-й гармо ники
Sv |
|
(£>vT |
( 1- 12) |
2 |
cos (соѵТ -|- 2фѵ). |
||
2%T |
|
||
|
sin 2 |
|
Выражение (1-12) представляет собой косинусоиду со случайной начальной фазой 2срѵ и амплитудой, изменяющейся по закону, sin 2соѵ77(2сйѵГ). Погрешность бѵ убывает с увеличением Т. При заданной допустимой погрешности измерения 6ѴД0П определим ми нимальный допустимый интервал интегрирования Т — ттах для низшей гармоники:
[ Т - т„ |
(1-13) |
|
V доп |
где 7"vmax = 2я/соѵтіп — период |
низшей гармоники. |
2лб„ |
Выражение (1-13) позволяет оценить необходимую (с точки зре ния допустимой погрешности) длительность реализации случай ного процесса. Максимальный интервал корреляции тшах опреде ляется здесь выражением (1-3). Следует отметить, что максималь ный интервал корреляции зависит от ширины спектра случайного сигнала. В настоящее время в публикациях, посвященных вопро сам корреляционного анализа, нельзя обнаружить единого подхода к определению максимального времени задержки сигнала т,„ по
известным |
граничным частотам |
Fmin—Fmax. Большинство |
работ |
базируется |
на ставшем уже классическим определении Ланге |
[331: |
|
|
~ — < |
т ш « Т • |
( 1 - 1 4 ) |
|
г m in |
|
|
Близкие или тождественные формулы для определения т,„ при водят В. В. Солодовников [97], А. Ф. Котюк [6] и др. Полемизи
10
руя по этому вопросу, Г. А. Балл [111 справедливо отмечает, что трудно согласиться с рекомендацией выбирать максимальную ве личину задержки, исходя только из нижней спектральной состав ляющей сигнала. Действительно, максимальная задержка тт (ап паратурный параметр) очевидно должна выбираться по крайней мере не меньшей максимального интервала корреляции случайного сигнала:
Тш<А'І’тах> |
(1-15) |
при этом тшах определяется либо в виде (1-3), либо из соотношения:
|
I Р.ѵ.ѵ (т) |
< 8 п р н |т |> т д.в, |
(1-16) |
тогда |
; здесь е наперед заданная величина: е |
1. |
|
С другой |
стороны, приI |
использовании коррелометра для оты |
скания скрытой в шумах периодичности, измерение корреляцион ной функции на интервале, равном лишь одному периоду перио дической составляющей, явно недостаточно. Поэтому для опреде ления тт и тпіах будем исходить из соотношений (1-3), (1-15) и (1-16). Тогда для случайного стационарного процесса, образованного про
пусканием |
белого шума |
через идеальный фильтр нижних частот |
|
с полосой |
прозрачности |
AF, максимальный |
интервал корреляции |
т ш а х п0 уровню е определится из (1-16) в виде |
[11 ]: |
||
|
|
_1_ |
(1-І7) |
|
|
Б п а х |
|
|
|
е я Д F |
|
Вопросам классификации методов измерений и способов построе ния приборов корреляционного анализа посвящены многочислен ные исследования [11, 106, 89,91,32, 33 и др. 1. Не задаваясь целью разработки новой классификационной схемы, разделим все корре ляционные анализаторы на приборы, использующие пропорцио нальные умножители (мультипликационные коррелометры), и при боры, структурная схема которых строится без использования устройств пропорционального умножения (знакосочетательные, ре лейные, квазимультипликационные и др. коррелометры).
Обобщенная структурная схема мультипликационного корре лометра с пропорциональным умножителем МУ показана на рис. 1-1. Структурные схемы реальных коррелометров, как пра вило, не содержат показанных на рис. 1-1 переключателей П1 и П2. Здесь они приведены для иллюстрации возможности как пред посылки устройства задержки (УЗ) входным преобразователям ка налов ВПК*, ВПК, (положение 2 переключателей Ш , П2), так и включения его непосредственно перед множительным устройством коррелометра (положение / переключателей П 1, П2). Последний вариант наиболее часто используется в приборах дискретного типа, так как дискретная и кодовая форма представления информации облегчают в ряде случаев техническую реализацию устройства за держки.
11
Исполнение входных преобразователей каналов определяется конкретным типом коррелометра и формой представления подлежа щей обработке исходной информации (аналоговые или дискретные сигналы, кодовые последовательности в виде электрических импуль сов, информация на перфокартах, кинопленке, в виде графиков ит. п.). Исполнение множительного (МУ) и накопителыю-усредняю- щего устройств (НУУ) определяется формой представления инфор мации на их входах. В качестве МУ коррелометров используются различные время-импульсные устройства, каскады с квазиквадра тичными характеристиками, цифровые регистры, преобразователи Холла и др., а в качестве НУУ — аналоговые интеграторы, интег рирующие устройства дискретного действия, счетчики импульсов и т. п.
Рис. 1-1. Структурная схема мультипликационного корре лометра
Снятие результатов измерений осуществляется с помощью уст ройств вывода результатов (УВР), выполняемых на стрелочных приборах, самописцах, цифровых индикаторах, на цифропечатаю щих и магнитных регистраторах и т. п. В особых случаях исполь зования результатов корреляционных измерений в качестве проме жуточной информации УВР в структурной схеме коррелометра может отсутствовать. Синхронизация работы всех узлов коррело метра осуществляется блоком синхронизации (БС).
Знакосочетательные коррелометры (называемые также просто знаковыми или полярными) строятся без применения пропорцио нальных умножителей и выполняют измерение нормированной функции корреляции (коэффициента корреляции) р (т) случайных сигналов на основании соотношения:
Q+ + (т) = 0,25 1 + 2 arcsin р (т) |
(1-18) |
где q++ (т) — вероятность совпадения положительных |
полярно |
стей исследуемого сигнала в функции от аргумента т. Соотношение (1-18) справедливо лишь для гауссовых сигналов.
Измерение функции корреляции сигналов с произвольными за конами распределения с помощью знакосочетательной аппаратуры связано с необходимостью отыскания и решения уравнения связи между q++ (т) и р (т) и соответствующей калибровки прибора.
12
Обобщенную структурную схему знакосочетательного коррело метра можно видеть на рис. 1-2. Принцип работы схемы заключается
в следующем. Центрированные |
входные случайные сигналы х {t) |
|
и у (/) поступают через |
входные |
преобразователи каналов ВПКЛ., |
ВПК;/ на ограничители |
ОгЛ., Ог;/ (сигнал у (t) — дополнительно |
через устройство задержки УЗ). Блоки Ог формируют последова тельности прямоугольных импульсов стандартной амплитуды и случайной длительности гиі. Длительности импульсов в зависимо сти от положенного в основу работы прибора уравнения связи ме жду р (т) и q или q,__ равны временным интервалам дискретных
выбросов сигналов определенной полярности. Как правило, в блоки
Рис. 1-2. Структурная схема знакосочетательного кор релометра
Ог входят усилительные элементы, обеспечивающие увеличение первой производной сигнала в точках пересечения оси времени.
Рост позволяет уменьшить по абсолютной величине раз-
•ѵ=0 ность между длительностью г-го импульса т„(- на выходе ограничи
теля и соответствующим ему временным интервалом между двумя ближайшими моментами пересечения сигналом х (t) оси времени с противоположными знаками производных.
Сигналы с выхода ограничителей поступают на схему совпаде ния (СС), выполняющую логическую операцию конъюнкции. Вы ходной сигнал блока СС представляет собой импульсную временную функцию, мгновенные значения которой пропорциональны величине
ОО
sgn .у (/) sgn у (/ — т). Накопителы-ю-усредняющее устройство НУУ обрабатывает за интервал (Т—т) импульсные последовательности с выхода СС и образует оценку нормированной функции корреля ции, регистрируемую устройством вывода информации УВР.
Релейные коррелометры представляют собой предельную моди фикацию коррелометра Стильтеса, в одном из каналов которого используется только знаковая информация о сигнале. Алгоритм работы релейных коррелометров основан на измерении релейной функции корреляции:
Я рел'О гН ^г f X { t ) s g n y { t — x ) d t . |
(1-19) |
1 о
13