Файл: Калинчук, Б. А. Анализаторы инфразвуковых случайных процессов.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). Следует отметить, что максималь­ ный интервал корреляции зависит от ширины спектра случайного сигнала. В настоящее время в публикациях, посвященных вопро­ сам корреляционного анализа, нельзя обнаружить единого подхода к определению максимального времени задержки сигнала т,„ по

известным

граничным частотам

FminFmax. Большинство

работ

базируется

на ставшем уже классическим определении Ланге

[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