Устройство для определения характеристик случайного процесса

 

Изобретение относится к специализированным средствам вычислительной техники и может быть использовано для определения выборочных квантилей. Цель изобретения - расширение функциональных возможностей за счет обеспечения возможности вычисления оценок дисперсий выборочных квантилей и оценок исходной плотности вероятности в точках квантилей. Устройство содержит регистр 1 буферной памяти, генератор 2 импульсов, переключатели 3, регистр 4, два коммутатора 5 и 9, счетчик 6 импульсов, счетчик 7 числа инверсий, схему 8 сравнения, группу элементов И 10, две группы блоков вычитания 11 и 15, группу масштабных усилителей 12, сумматоры 13, счетчик 14 циклов, группу блоков вычисления обратной величины 16, группу квадраторов 17. 1 ил.

СОЮЗ СОВЕТСКИХ

СОЦИАЛИСТИЧЕСКИХ

РЕСПУБЛИК ии 4 С 06 F 15/36

ОПИСАНИЕ ИЗОБРЕТ

М ABTOPCHOMV СВИДЕТЕЛЬСТВУ

ГОСУДАРСТВЕННЫЙ НОМИТЕТ

ПО ИЗОБРЕТЕНИЯМ И ОТКРЫТИЯМ

ПРИ ГКНТ СССР (61) 1043665 (21) 4304522/24-24 (22) 09.09,87 (46) 15.05.89. Бюл. № 18 (72) В.А.Поляков и Р.Г.Толпарев (53) 681 ° 3(088.8) (56) Авторское свидетельство СССР № 1043665, кл. С 06 F 15/36, 1982. (54) УСТРОЙСТВО ДЛЯ ОПРЕДЕЛЕНИЯ ХАРАКТЕРИСТИК СЛУЧАЙНОГО ПРОЦЕССА (57) Изобретение относится к специализированным средствам вычислительной техники и может быть использовано для определения выборочных квантилей. Цель изобретения — расшире„„SU„„1479939 A 2 ние функциональных возможностей эа счет обеспечения возможности вычисления оценок дисперсий выборочных квантилей и оценок исходной плотности вероятности в точках квантилей.

Устройство содержит регистр 1 буферной памяти, генератор 2 импульсов, переключатели 3, регистр 4, два коммутатора 5 и 9, счетчик 6 импульсов, счетчик 7 числа инверсий, схему 8 сравнения, группу элементов И 10, две группы блоков вычитания 11 и 15, группу масштабных усилителей 12, сумматоры 13, счетчик 14 циклов, группу блоков вычисления обратной величины 16, группу квадраторов 17. 1 ил.

1479939

Устройство относится к специализированным средствам вычислительной техники и может найти применение при определении статистических характеристик случайных величин. Устройство предназначено для определения выборочных квантилей, оценки их дисперсий и оценки значений плотности вероятности в точках квантилей при неизвестном исходном распределении и является усовершенствованием устройства по авт.св. Р 1043665.

Цель изобретения — расширение функциональных возможностей за счет обеспечения возможности вычисления оценок дисперсий выборочных квантилей и оценок исходной плотности вероятности в точках квантилей.

На чертеже изображена структурная схема устройства.

Устройство содержит регистр 1 буферной памяти, генератор 2 импульсов, переключатели 3, регистр

4, первый коммутатор 5, счетчик 6 импульсов, счетчик 7 числа инверсий, схему 8 сравнения, второй коммутатор 9, группу элементов И tO первую группу блоков 11 вычитания, rpvnny масштабных усилителей 12, сумматоры 13, счетчик 14 циклов, вторую группу блоков 15 вычитания, группу блоков 16 вычисления обратной величины, группу квадраторов 17.

Устройство реализует следующий . алгоритм работы.

На вход устройства поступает последовательность независимых отсчетов, закон распределения Р(х) и плотность распределения f(x) которых не известны.

Устройство периодически делает выборку объемом N членов и ранжирует ее в соответствии с х. < х.. а ...cx <» ..." х.„,(1) где — номер члена в ранжированном ряду.

Каждый член ранжированного ряда х представляет выборочный кван(>) тиль х pI с уровнем вероятности Р удовлетворяющим соотношению

В случае, если по заданной Р из (2) получается (3) (i) = P N

J л дробное число, то значение х ., обе спечиванл,е данный уровень P, нахоЪ дится с помощью интерполяции из выражения

1О

+ (4) (:(l>jg — (i) j (1) 1 — доля интерполяции; в данной формуле — операция взятия целой части; где п(я+1 х -; . — соответственно ближайший (.(11j 3 больший и ближайший мень20 ший к промежуточному между ними х<,1; члены ранжированной выборки с указанными целочисленными номерами.

Выборочные квантили, удовлетворяющие соотношению (2), а также полученные на основе (4), выводятся из устройства для использования.

Для N-2W-2 выборочных квантилей

ЗО х с номерами членов ранжированной выборки от (i) = 2 + W до (i)

= N-W-1, т,е. полученных без интерполяции, бчдут определены оценки л дисперсий 1)(хр; ) и оценки исходнОЙ плОтности верОятнОсти в тОч

Ь ках хр> квантилей f.(х ). С этой целью заранее, с помощью выражения (5) Для каждого из квантилей оценка дисперсии производится с помощью

45 л л й

Р ) (Q(11 У 1 (1 щ1

Pj (1 — P>) N где К—

Учитывая, что в данном случае (j) = i+W, из (6) следует, что операции, необходимые для определения дисперсии, надо проделать для каждой пары членов выборки с номером (i+2W) и (i), где i = 2 + (N-2Ы-i) 40 где j J — - операция взятия целой части, определяется значение W.

1479939

Для каждого из квактилей оценка исходной плотности в точках квантилей производится с помощью п

1()) Ф W j (()I 11

Генератор 2 импульсов, закончив 40 выдачу пачки иэ N тактовых импульсов, начинает. выдавать непрерывную последовательность синхроимпульсов, поступающих на синхровходы первого коммутатора 5, второго коммутатора 9 и на вход счетчика 6 импульсов.

Под воздействием последовательности синхроимпульсов и в такт с ними начинается процедура упорядочения (ранжирования) выборки. Частота сле- 50 дования синхроимпульсов выбирается такой, чтобы максимально сократить время ранжировки выборки.

Процедура ранжирования основания на алгоритме попарной перестановки.

При первом просмотре каждое значение выборки х, сравнивается с соседним ( значением х;+, . Пара значений выбор55

Алгоритм реализован в устройстве следующим образом.

Регистр 1 буферной памяти после заполнения выдает сигнал в генератор 2 импульсов. По этому сигналу генератор 2 импульсов выдает на тактовый вход регистра 4 пачку из N импульсов. Под воздействием тактовых импульсов значения последовательности из регистра 1 буферной памяти последовательно поступают через информационный вход в регистр 4 объемом N разрядов. Таким образом, в основном регистре А образуется неупорядоченная выборка из значений х, Х у ° ° ° Х у ° ° ° ухц ° 25

Частота следования тактовых им пульсов выбирается таким образом, чтобы согласовать скорость поступления исходной последовательности на вход регистра буферной памяти со скоростью записи выборки в регистр 4 и скоростью ее обработки. Значения выборки одновременно с записью в регистре 4 через информационные выходы каждого разряда регистра 4 устанавливаются на соответствующих им входах первого коммутатора 5 переключателей 3 и соответствующих им входах каждого из элементов И 10 группы.

1 ки для которых справедливо неравенство х, > х;,, меняется местами, Лопарное сравнение и перестановка членов выборки продолжается до тех пор, пока не закончится процесс ранжировки выборки, т.е. х, с х,16

° ° ° х(1 е ° ° х () °

«С ах

Факт выполнения неравенства х; > х;+, при попарном просмотре выборки является инверсией, а. функция суммы инверсий определяется выражением

N-E т и (xx,, х;, ), (8)

1э х;> x„1з где u(x<, х;+, )

О, х; а х;+<

Значение с вычисляется при каждом цикле попарного сравнения и перестановок членов выборки. Признаком окончания ранжировки выборки является значение 6 1, сосчитанное за цикл просмотра.

Алгоритм ранжировки выборки реализуется следующим образом. В каждом цикле просмотра первый коммутатор 5 в такт синхроимпульсам производит последовательное подсоединение информационных выходов каждой соседней пары разрядов регистра 4 к двум входам схемы 8 сравнения. В результате сравнения каждой пары соседних членов выборки схема 8 сравнения выдает значение инверсии

u(x; х, ) в виде "0" или "1".Эти значения поступают в счетчик 7 числа инверсий, который в соответствии с выражением (8), определяет сумму инверсий а за каждый цикл просмотра. Одновременно каждое значение инверсии поступает на вход второго коммутатора 9.

Работая синхронно с первым коммутатором 5, второй коммутатор 9 последовательно соединяет выход схемы 8 сравнения с входом того переключателя 3, который своими другими двумя входами, а также двумя выходами соединен соответственно с информацион-,— ными выходами и информационными входами каждого иэ пары разрядов регистра 4, содержимое которых анализировалось схемой 8 сравнения. Одновременно второй коммутатор 9, чтобы не усложнялась конструкция переключате5 14 лей 3, преобразует поступившее яа е-;а вход значения ц Ф 0 в значение на своем выходе 1, оставляя неизменным значение u =

Переключатель 3 соединяет подключеннные к нему информационные выходы сравниваемой пары соседних разрядов к своим выходам таким образом, чтобы при и Ф 0 на входе второго коммутатора 9 содержимое сравниваемых разряцов регистра оставалось неизменным, а при u = 1 поменялось; местами путем переписывания иэ одного разряда в другой. Затем переключатель 3 возвращается в нейтраль-. ное положение.

Конец каждого цикла просмотра определяется с помощью счетчика 6 импульсов и счетчика 14 циклов. В первом цикле просмотра счетчик 6 импульсов отсчитывает М-1 импульс из последовательности синхроимпульсов, поступающих на его вход из генератора 2 импульсов (N — объем выборки). Этим определяется граница первого цикла. В результате будет осуществлено однократное попарное сравнение и перестановка в соответствии с принятым алгоритмом всех членов выборки, содержащихся в регистре 4.

Если при этом за цикл счетчик

7 числа инверсий получит 1, то осуществляется переход к второму циклу просмотра. При о 1 потенциал, подаваемый на управляющие входы генератора 2 импульсов и группы элементов И 10 с выхода счетчика

7 числа инверсий, не меняется. Поэтому генератор 2 импульсов будет продолжать выдавать последовательность синхроимпульсов для управления коммутаторами 5 и 9, а группа элементов И 10 продолжает оставаться запертой, не разрешая выдачу информации, установленной на их входах, для дальнейшего использования. Счетчик 6 импульсов, закончив счет до

N-1, своим выходным сигналом сбрасывает в исходное состояние счетчик

7 числа инверсий и возвращает в исходное состояние первый коммутатор

5 и второй коммутатор 9, так что каждый очередной цикл просмотра будет начинаться с одних и тех же разрядов регистра 4. Одновременно счетчик 6 импульсов своим выходным сигналом записывает "1 в счетчик

14 циклов.

/9939

Ç0

6

Счетчик 14 циклов передает записанную в нем сулему (в данном случае "1") на второй вход счетчика 6 импульсов, где эта сумма вычитается из записанного там установочного числа, до которого должен считать счетчик 6 импульсов. В результате количество попарных сравнений членов вборки во втором цикле просмотра уменьшится на "1", до N-2, и самый большой член выборки в ран-. жированном ряду, уже занявший в первом цикле просмотра крайнее место, в попарном сравнении во втором и последующих циклах просмотра принимать участие не будет, т.к. первый коммутатор 5 и второй коммутатор 9 произведут до своего возвращения на исходные позиции меньшее количество переключений.

Циклы просмотра с уменьшением числа попарных сравнений на 1 в каждом очередном цикле продолжаются до тех пор, пока выполняется усло— вие Р 1. При С < 1 счетчик 7 числа инверсий своим выходным сигналом прекращает подачу синхроимпульсов генератором 2 импульсов, сбрасывает в исходное состояние счетчик 14 циклов и отпирает элементы

И 10. Этот момент соответствует формированию на входах элементов И 10 ранжированной выборки, Искомые значения квантилей, удовлетворяющих выражению (2), поступают на выходы устройства непосредственно с выходов тех элементов И 10, номера которых соответствуют (2).

Уточняющая интерполяция производится в соответствии с (4) следующим образом. Информация с каждой пары элементов И 10 с известными номерами (() ) и ((д) 3 + 1 поступает на два входа блока 11 вычитания, соответствующего этой паре.

Кроме того, выход элемента И 10 с номером ((х) j соединен с одним из входов сумматора 13, входящего в данную цепочку интерполяции. В блоке

11 определяется разность х <, — х (,1,g которая поступает в масштабный усилитель 12, входящий в эту же цепочку. В усилителе 12 полученная разность умножается на известную долю интерполяции „. Произведение поступает на второй вход сумматора

13, где складывается с выходной информацией ((i); -го элемента И 10, 1479939 для использования. Затем процесс повторяется.

Формула из обретения!

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

30

С выходов устройства значения л л л х, 1)(хр ) и f (х р ), определенные

Э У по данной вь.борке, поступают на входы других технических устройств

Составитель Е.Хуртин

Редактор М.Циткина Техред А.Кравчук Корректор С.Черни

Заказ 2544/47 Тираж 669 Подписное

ВНИИХИ Государственного комитета по.изобретениям и открытиям при ГКНТ СССР

113035, Москва, Ж-35, Раушская наб., д. 4/5

Производственно-издательский комбинат "Патент", r.Óæroðoä, ул. Гагарина,!01 равной х, „, в результате полуг 11 чается искомый квантиль х (,1, поступающий ка соответствующий выход устройства.

Оценка дисперсии квантилей производится в соответствйи с (6). Для этого информация с выходов каждой пары элементов И 10 с номерами (i) и (i+2W) поступает на два входа вто- 10 рого блока 15 вычитания,соответствующего данной паре элементов. В блоке

15 определяется разность х <11+® х(1 „,1, поступающая в соответствующие этой паре элементов И 10 блок

16 вычисления обратной величины и квадратор 17. Выходы блоков 17 являются выходами оценки дисперсии квантилей устройства.

Устройство для определения характеристик случайного процесса по авт.св. Р 1043665, о т л и ч а ю— щ е е с я тем, что, с целью расширения функциональных возможностей за счет обеспечения возможности вы-, числения оценок дисперсий выборочных квантилей и оценок исходной плотности вероятности в точках квантилей, в Hего введены вторая группа блоков вычитания, группа блоков вычисления обратной величины и группа квадраторов, причем вход уменьшаемого блока вычитания второй группы подключен к выходу i-ro элемента И группы (i = 2,...,(N — 2W—

l), W — целая часть от величины (Я, И1 — объем выборки)» а вход вычитаемого — к выходу (i+2W)-ro элемента И группы, выход i-го блока вычитания подключен к входу одноименного блока вычисления обратной величины группы, выход которого является выходом оценки плотности вероятности устройства, и подключен к .входу одноименного квадратора группы, выход которого является выходом оценки дисперсии квантилей устройства.