12.2. Системы для измерения законов распределения вероятностей случайных процессов

Одномерный интегральный закон, или интегральная функция распределения вероятностей, определяется вероятностями нахож­дения исследуемого процесса ниже уровня х, который может из­меняться от—оо до +оо: F(x)=p[—оо<£Гл:(^) sgrjc] (рис. 12.2). По определению F(—оо)=0 и /⁷(+с»)=1.

^хтах

^xmin

а) Рис. 12.2. Функция распределения вероятностей:

а — к определению функции распределения вероятностей; б — вид функции распределения вероятностей

Ятах х+йх X	х
			^"р	V	fflCZZJl	1—
хт'ш
	t

а)

kf

Xmin

б)

^хтах

Рис. 12.3. Плотность распределения вероятностей:

а — к определению плотности распределения вероятностей; б — вид плотности распределе­ния вероятностей

Двумерный интегральный закон распределения вероятностей стационарного случайного процесса x(t) определяется вероятностя­ми нахождения x(t) и x(t-\-x) ниже уровней х\ и х₂:

F{x_u x₂, т)=р[л:(0<л:,; x{t+x) sCx₂].

Одно- и двумерные дифференциальные законы распределения вероятностей случайного процесса x(t)—плотности распределения вероятностей (рис. 12.3)—определяются следующими выраже­ниями:

f(x)=dF(x)/dx; f(x_u х₂, т) =

= d²F(X\, Х₂, Т)1йХ\йхг

или

f (x) dx=p [x^x(t) ^x-\-dx]; f(x_u x₂, x)dx_u dx₂=p[x_l.<^x(t)<^ ^Xi+dxu x₂^.x(t-\-t)z^x₂zkdx₂].

Определение вероятностей p[x(t)^x] и p[x^x(t) ^x-r-dx] можно производить путем суммирования отрезков времени, в те­чение которых выполняются указанные неравенства, и отнесения полученной суммы к времени наблюдения:

для функции распределения

F{x) = lim^-,

для плотности вероятности

f(x, Дх:) = lim

г->оо Т Ах

Если исследуемый процесс представлен решетчатой функцией (в виде дискрет, размешенных через интервал дискретизации Д^), то выражения для интегральной функции и плотности распре­деления примут следующий вид:

F(x, N, Д/) = Шп

эо П/ 1

Л/-»=о ^N

f(x, N, At, Дяг) = Шп

N->oo ^N & Дл:->0

В этих выражениях At—интервал дискретизации; п—количест­во дискрет, уровень которых меньше х для функции распределения или находится в интервале х-\-Ах для плотности распределения; N==T/At—количество дискрет в исследуемой реализации—объем выборки.

В большинстве случаев Д^ равно или больше интервала кор­реляции исследуемого процесса.

Для определения функции и плотности распределения должен быть получен ряд значений F(x) и f(x, Ax) в пределах динами­ческого диапазона исследуемого процесса. Для этого динамичес­кий диапазон разбивается на несколько интервалов квантования и для каждого интервала квантования определяются значения F и f.

При экспериментальном определении закона распределения вероятностей неизбежно возникают методические погрешности обусловленные конечной длительностью наблюдения (Т<с°о) или выборки (Л/<оо) реализаций и конечным значением интервала квантования по уровню АхфО. Именно ввиду наличия в первую очередь методических погрешностей в результате измерительного эксперимента получаются не точные, а приближенные выраже­ния—оценки законов распределения вероятностей:

Т '

St,- 1

F*(x, T)= ^&1

Т Дх Г(х, N, Д*) = "

f*(x, Ax, Т) =

N '

Определение интервала квантования по параметру х произво­дится известными способами. В зависимости от степени аппрокси­мирующего многочлена в теории вероятностей различают пред­ставление плотности распределения гистограммой (приближение на интервалах квантования степенным многочленом нулевой сте­пени) или полигоном (приближение многочленом первой степени).

Неравенства Чебышева [12.5] позволяют при заданной по­грешности оценки законов распределения приближенно определить время наблюдения при некоррелированной выборке, когда интер­вал квантования At равен или больше интервала корреляции:

p[max\F*(x, T)-F(x)\^b₁]^ ^Ы

47V '

р [max | f* (х, их, T)Ax-f (х) Д* | > в,]

47V

Более точное определение времени анализа Т или объема вы­борки N (при некоррелированной выборке) может быть выполне­но, если задана относительная средняя квадратическая погреш-ость измерения а в фиксированных точках F*(x) и f*{x, Ax) по формулам [12.2, 12.5]:

_N _ l-F*(x) l__

^F~~ o*_FF*(x) ~ o2_FF*(x) '

—-f*(x. Ax) #, = > ^l

^f a^_ff(x, Ax) o2_ff*(x, Дх)Лх "

Приближенные значения N_F и Nf получены при условиях F*W«I uAxf*(x, Д*)«1.

Если между выборочными значениями имеется корреляционная связь (Д?<то, где то—интервал корреляции), то объем выборки нужно увеличить.

Поскольку F*(x, T) и f*(x, Ax, T) связаны между собой из­вестными зависимостями, то в статистических анализаторах, види­мо, можно ограничиться только измерением F*{x, T) или f*(x, Ax, T). Аппаратурные решения для измерения F*{x, T) проще, чем для f*(x, Ax, T). Однако большинство ИС делается для измерения плотности распределения. Это объясняется тем, что при преобразовании плотности распределения в функцию рас­пределения погрешности получаются существенно меньшими, чем при обратном преобразовании [12.6].

Рассмотрим основные структуры и алгоритмы ИС, измеряю­щих дискреты функций или плотностей распределения вероят­ностей. Такие ИС (будем называть их анализаторами вероятнос­тей) могут быть одноканальными и многоканальными.

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

для представления законов распределения, следует последователь­но изменять х или месторасположение интервалов квантования-по уровню Да: и производить измерение величин F*(Xi) и f*(xt, Дх)- Естественно, что при последовательном измерении всех дискретных величин F* и f* на проведение анализа затрачивает­ся большое время.

x(t)

х (П

x(t)^x_z

Рис. 12.4. Анализатор с ШИМ и аналоговым интегрированием (х, и х₂—уровни: анализа)

			1		—*•	СТ
X} Ж.		\		Л- / /JT-	—	СТ
. з»		\
	X

r(t)

x(ty

Ш)

_П

Рис. 12.5. Одноканальная система для измерения функции распределе­ния вероятностей

Рис. 12.6. Система с телевизионным датчиком ТД для измерения распре­деления вероятностей

Структурные схемы одноканальных анализаторов представле­ны на рис. 12.4—12.6. Наиболее часто в одноканальных анализа­торах, оперирующих с непрерывными сигналами и основанных на измерении относительного времени пребывания процесса x{t), ис­пользуется образование ШИМ-сигналов и их интегрирование в те­чение времени Т.

Логическая схема алгоритма анализатора вероятностей, изо­браженного на рис. 12.4 и на'строенного на получение дискрет плотности распределения вероятностей, охватывает:

начальное состояние: Ф(В) q>[x\; х?, Кл{$)\ I[*(0JX; образование ШИМ-сигналов: X{I[*(0<^i/"o; A*(*i)I II I[*(0 =C <;х₂/ио; At{x₂)]}X;

формирование дискреты плотности распределения: У,1{[М(х_г)—

М(х,)]; ы₀/П*.. ux]\l[S:r(x_v Д*)]Ф(Е),

~М(х_г)]; u₀} IjElMfo)

Разрешающая способность таких анализаторов определяется дрейфом уровня амплитудных селекторов и формой импульсов. Она обычно ограничивается примерно 1/20—1/100 частью макси­мального значения измеряемого напряжения.

Если исследуемый процесс представлен в квантованном по вре­мени виде, используются амплитудный дискриминатор с изменяе­мым уровнем сравнения х (рис. 12.5) и суммирующие счетчики, на выходе которых получаются значения п и N.

СЛСА такого анализатора вероятностей:

начальное состояние: ф(х; N=0; п=0; /=0;)Х;

формирование х,- и х: X Ф(/ : =/+1) 1(xj || х)Х

счет суммарного количества дискрет /' и I (с уровнем Xj^.x): X{I(S:D/) О [I(CR :*,,*) а(х^х) I(S:Di]}X;

выдача информации: x<a(/=M) [l(Zi)\\I(Di)] Ф(Е).

Здесь /—текущее суммарное количество всех дискрет проанали­зированной части исследуемой реализации, а I—текущее количест­во дискрет с амплитудой, меньшей или равной уровню анализа х.

На рис. 12.5 показан анализатор, предназначенный для изме­рения F* (х). Предоставляем читателям возможность составить схе­му и алгоритм анализатора такого типа, позволяющего произвести .измерение f * {х, Ах).

По такому принципу работают одноканальные анализаторы .импульсов АИ-3—АИ-8, обладающие высоким быстродействием ■{до 10⁵ имп./с), повышенной чувствительностью (максимальное значение амплитуды исследуемых процессов до 3—10 В).

Простотой и высоким быстродействием отличаются анализато­ры с преобразованием исследуемого временного процесса x(t) в лространственный сигнал х(1). Это преобразование может быть ■осуществлено с помощью электронного осциллографа при подаче ла отклоняющие пластины электронно-лучевой трубки исследуе­мого процесса и несинхронного развертывающего напряжения. Возможно также использование оптических систем с перемещаю­щимися диафрагмами и фотоэлектронным воспринимающим и ин­тегрирующим устройством либо промежуточное преобразование с ломощью телевизионного датчика ТД (рис. 12.6) с последующим анализом полученного потенциального изображения.

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

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

В аналоговых анализаторах используются дискриминаторы, выделяющие каналы—интервалы значений измеряемых величин,

аналоговые накопители (индивидуальные для каждого канала) и устройства вывода (рис. 12.7). В качестве амплитудных дискрими­наторов могут использоваться пороговые устройства (типа тригге­ров Шмитта), фотоприемники с диафрагмами, выделяющими ка­налы, и т. п. Обычно при использовании дискриминаторов число-каналов не превышает 30—50. Для накопления данных в этом слу­чае часто используются конденсаторные устройства.

x(t,l)

Рис. 12.7. Многоканальная аналого­вая система для измерения распре­деления вероятностей

		1	-*■		—»■	ст
x(t,l)
		п	—*-		—*-	ст

Рис. 12.8. Многоканальная цифровая система для измерения плотности распределения вероятностей

Рис. 12.9. Структурная схема ИС типа АИ

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

В многоканальных анализаторных ИС типа АИ, предназначен­ных для статистического анализа импульсов [12.5J, используются в качестве интегрирующих устройств магнитные устройства памя­ти (рис. 12.9). Результат измерения в виде кодоимпульсного сиг­нала I (Zj) в таких ИС поступает в регистр RG^ и определяет ад­рес ячеек памяти Adz,-, в которых накапливается количество зна­чений исследуемой величины, соответствующих каждому результа­ту измерения. С помощью RG\ к числу, хранящемуся в ячейках памяти данного интервала АЦП, после поступления очередного результата измерения добавляется единица. Промежуточные ЗУ (на магнитных лентах, перфолентах, перфокартах и т. п.) исполь­зуются для увеличения объема выборки, уменьшения потерь ин­формации и т. д.

Весьма удобной для построения многоканальных анализаторов является мультиплицированная структура (рис. 12.10), пригодная для оценки параметров законов распределения как случайных ве­личин, так и случайных функций. Здесь число ступеней выходного

напряжения и_к ЦАП (число ступеней «развертки») равно числу интервалов группировки f*(x) или F*(x); при этом на каждом такте и_к проводятся определение числа сработавших устройств сравнения и запись этого числа в соответствующие сумматоры, количество которых равно числу ступеней и_к. Если определение числа сработавших устройств сравнения осуществляется путем опроса, то в качестве сумматоров можно использовать обычные счетчики. В зависимости от того, сбрасываются или нет каналь-

х,

—»-

\

' Т

( '

—»■

1

—*-ь

—>

с^ стг

L—»-

.

2С2

>—»■

\

* т г

СП ~

i

I

<

-*■/

г-*■

*

ЗначениЯ у f*(x) unu F*(x)

Рис. 12.10. Мультиплицированная система для измерения распреде­ления вероятностей

ные триггеры в конце каждой ступени, будет измеряться f*(x) или F* (х). Анализаторные ИС позволяют получать законы распре­деления амплитуд А, временных интервалов В, номеров источни­ков информации Д в различных сочетаниях (обычно не превышаю­щих трехмерного): А, В, АА, ААА, АВ, ААВ, ВВ, АД, ВД.

Для получения ^-мерных законов распределения ИС содержат соответствующее количество АЦП и промежуточное ЗУ (на рис. 12.9—RAMi). При этом основное запоминающее устройство ЯАМ₂ разделяется на k частей.

В Советском Союзе разработаны и выпускаются анализаторы импульсов типа АИ. Отечественные промышленные анализаторы АИ-128-1М (А), АИ-256 (А, В, Д), АИ-1024 (АА, В, АД), АИ-2048 (А, В, АВ), АИ-4096 (АА, ВВ, АД, ВД) имеют число каналов, равное 'соответственно 128, 256, 1024, 2048 и 4096, число уровней квантования по амплитуде от 100 до 512 (в АИ-4096-3-А6—до 4096), ширину временных каналов от Ю^-9 до 64-Ю^-6 с, число каналов координат датчиков от 8 до 100, емкость каналов от 10⁴ до 10¹⁸ бит, цикл сортировки от 10 до 56 мкс.

Анализаторы АИ-4096-3 и А-2048-1 выполняют дополнительные операции обработки сигналов: сглаживание, нормирование, раз­ложение на элементарные компоненты и т. п. Для этого анализа­торы имеют арифметические устройства и устройства управления, позволяющие производить помимо основных арифметических one-

раций операции сравнения, модификации адреса, условного пере­хода.

Анализаторы типа АИ выпускаются в блочном исполнении. К отдельным блокам, входящим в состав современных анализато­ров АИ, относятся АЦП, ЗУ, АЛУ, управляюшие, выходные устройства, устройства коммутации и передачи информации. Принцип их действия и технические характеристики изложены в [12.5].

Ъ&Ю]

Функцио­нальный преобразо­ватель

Xlfj

Рис. 12.11. Система для измерения распределения вероятнос­тей с аппроксимацией

Для контроля правильности функционирования АИ и опреде­ления их характеристик используются генераторы образцовых по­токов импульсов. Разработаны генераторы точных амплитудных значений, распределенных, например, по нормальному закону, ге­нераторы белого шума и т. п. Предусмотрен контроль входных и выходных устройств АИ. Разработаны тестовые программы, слу­жащие для проверки работоспособности узлов обработки инфор­мации АИ. В [12.5] приводятся данные измерительных комплек­сов и центров, в которых с помощью коммутаторов и согласующих устройств объединяются для совместной работы несколько АИ и ЭВМ.

Следует отметить, что при работе АИ возникают специфические погрешности, которые следует учитывать. Для этого используется аппарат теории массового обслуживания.

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

00

Если учесть, что M[x(t)] = f к (t) f {х) dx, коэффициенты раз-

СО

ложения Фурье функции C_k — Г <р« [х Ш ^dx< ^a f (^х) =

—сю со

=2 ^Ck9k\t]. то можно принять, что C_kr=M{f_k[x(t)]} [12.2].

k=0

Сигналы на выходе фильтра ф/_г(^) ортогональны и зависят от времени, а на выходе функционального преобразователя ЧМ*(01—^от исследуемого процесса.

Для графического представления необходимо, чтобы время бы­ло связано с уровнем сигнала. Интервал между импульсами 6(0 при определении Си должен быть больше времени усреднения.

Порядок приближения т может быть определен из анализа средней квадратической погрешности оценки:

со т

—со /г=1

При получении f* (x) описанным методом отсутствует необхо­димость квантования процесса по времени, могут быть сокращены емкость используемой памяти и время измерения.