Адаптивный статистический анализатор

 

Изобретение относится к специализированным средствам автоматики и вычислительной техники, предназначенным для вычисления плотности распределения вероятностей и ее параметров в условиях недостатка информации о свойствах сигнала. С целью повьппения точности вычислений в условиях недостатка априорной информации о частотных и амплитудных свойствах сигнала в анализатор введены второй регистр микрокоманд, последовательно соединенные блок выборки и запоминания , блок вычитания, первый коммутатор , дискр;иинатор знака, элемент 2И-Ш1И, цифроаналоговый преобразователь и второй коммутатор. При этом выходы коммутаторов подключены соответственно к информационному и опор- g ному входам аналого-цифрового преобразователя . 1 з.п. ф-лы, 3 ил. (Л

СООЗ СОВЕТСНИХ

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

РЕСПУБЛИН (51) 4 G 06 F 15/36

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

H АВТОРСКОМУ СВИДЕТЕЛЬСТВУ

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

ПО ДЕЛАМ ИЗОБРЕТЕНИЙ И ОТНРЫТИЙ (21 ) 4229291/24-24 (22) 13.04.87 (46) 30.10 ° 88. Бюл. ¹ 40 (71) Ленинградский электротехнический институт им. В.И. Ульянова (Ленина) (72) В.И. Якименко, И.Б. Столбов, А,Ф. Бульбанюк и Ц.Б. Эпштейн (53) 681.3(088.8) (56) Авторское свидетельство СССР № 1135340, кл. G 06 F 15/36, 1983.

Авторское свидетельство СССР № 1001110, кл. G 06 F 15/36, 1981. (54) АДАПТИВНЫИ СТАТИСТИЧЕСКИЙ АНАЛИЗАТОР (57) Изобретение относится к специализированным средствам автоматики и вычислительной техники, предназнаSU» 1434453 А 1 ченным для вычисления плотности распределения вероятностей и ее параметров в условиях недостатка информации о свойствах сигнала. С целью повьппения точности вычислений в условиях недостатка априорной информации о частотных и амплитудных свойствах сигнала в анализатор введены второй регистр микрокоманд, последовательно соединенные блок выборки и запоминания, блок вычитания, первый коммутатор, дискриминатор знака, элемент

2И-ИЛИ, цифроаналоговый преобразователь и второй коммутатор. При этом выходы коммутаторов подключены соответственно к информационному и опор- <о нову входам аналого-циФрового лраовразователя. 1 з.п. ф-лы, 3 ип.

1434453

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

Целью изобретения является повы- щ шеяие точности оценки плотности ðàñпределения вероятности. На фиг. 1 показана структурно-функциональная схема адаптивного анализатора, на фиг. 2 — схема первого бло- 15 ка определения размаха, на фиг. 3— плотность распределения вероятностей — нормального Ч и гамма-расн пределеиия М,, вычисленная без адаптивных режимов обработки и после ада- 20 птивной обработки входного сигнала.

Адаптивный статистический анализатор содержит последовательно соединенные блок 1 выборки и запоминания, ! усреднитель 2, блок 3 вычитания, пер- 25 вый коммутатор 4, аналого-,цифровой преобразователь 5, блок 6 памяти и суммирующий счетчик 7, последовательно соединенные генератор 8 импульсов управляющий делитель 9 частоты, первый элемент ИЛИ 10, элемент 11 задержки, второй элемент ИЛИ 12 и первый регистр 13 микрокоманд, вход записи которого подключен к выходу элемента И 14, элементы ИЛИ 15 и 16, триггеры 17 и 18, второй элемент И 19, выход которого подключен к управляющему входу дискриминатора 20 .знака, первый и второй выходы которого соединены с соответствующими входами элемента 2И-ИЛИ 21, к выходу которого подключены последовательно соединенные цифроаналоговый преобразователь 22 и второй коммутатор 23.

Информационные входы первого и вто-45 рого блоков 24 и 25 определения размаха объединены, а выходы подключены к первому и второму входам делителя

26, выход которого подключен к первому входу блока 27 сравнения, к второму входу которого подключен код порога сравнения.

Выход элемента ИЛИ 28 соединен с тактовым входом второго регистра 29 микрокоманд, выходы которого подключены к управляющим входам соответствующих элементов. Выход переполнения счетчика 20 выборок подключен к тактовому входу блока 27 сравнения и пер-, вому входу элемента И 31, второй вход которого объединен с инверсным входом элемента И 14, вторым входом элемента И 19 и подключен к инверсному выходу триггера 18. Выход элемента И 31 является управляющим выходом анали- затора (для задания начала вывода ре" зультата обработки).

Блок 24 определения размаха содержит, например, последовательно соединенные первый элемент И 32, первый регистр 33, первый элемент 34 сравнения. Второй элемент И 35, второй регистр 36 и второй элемент 37 сравяения последовательно соединены, разрядные входы элементов И 32 и 35 объединены и подключены к вторым входам элементов 34 и 37 сравнения и являются кодовым входом блока 24, а знаковый разряд кода подключен к соответствующему входу элемента И 32 и инверсному входу элемента И 35. Тактовый вход блока соединен с тактовыми входами элементов 34, 37 сравнения, выходы которых соединены с соответствующими входами элементов И 32 и

35, тактовыми входами регистров 33 и

36 и подключены к соответствующим входам элемента ИЛИ 38, выход которого соединен с тактовым входом сумматора 39, кодовые входы которого подключены соответственно к выходу элементов И 32 и 35, а выход сумматора

39 является информационным выходом блока 24 (распределение кодов разных знаков может быть также выполнено при помощи компаратора, открывающего один из элементов И 32 или

35), Входы элемента ИЛИ 40 являются входами обнуления б,пока определения размаха.

Адаптивный анализатор осуществляет вычисление плотности распределения вероятностей sa два этапа: этап обучения, при котором производится автоматический поиск рациональной частоты дискретизации входного сигнала н определение максимального, минимального уровней и размаха сигнала,и этап оценивания, во время которого вычисляется оценка функции Ъ „ при заданных на этапе обучения начальных условиях.

В исходном состоянии блоки анализатора обнулены (установлены в нулевое состоние ), в генераторе 8 импульсов установлена тактовая частота, определяемая максимальным

453

30

55 т м Н Вх

1434 быстродействием логических элементов анализатора. В управляемом делителе 9 частоты установлен коэффициент деления, определяющий интервал dt временной дискретизации входного сигнала, а в счетчике 30 выборок — коэффициент пересчета N 1, соответствующий объему анализируемых отсчетов. Время задержки с„ элемента 11 задержки несколько больше времени выполнения операцнй в усреднителе 2 и блоке 3 вычитания. К первому входу блока 27 сравнения подключен код F допустимой относительной погрешности оценки среднего значения m „

С этапа обучения начинается анализ сигнала, подключенного к входу блока 1 выборки и запомийания: одновременно включается генератор 8 импульсов, а на вход "Пуск" подается одиночный импульс запуска "1", который устанавливает в исходное состояние триггеры 17 и 18, а через элемент ИЛИ 10 поступает на тактовый вход блока 1 выборки и запоминания, осуществляя отсчет х(1 dt), через элементы 11 задержки и элемент ИЛИ 12 — на тактовый вход регистра !3 микрокоманд и одновременно через открытый элемент И 14 — на вход записи регистра 13 микрокрманд. При этом в счетчике 30 выборок устанавливается код n = 1, а с первого разряда регистра 13 микрокоманд через элемент ИЛИ 16 на тактовый вход аналого-цифрового преобразователя 5 (АЦП) подается микрокоманда

Сl разрешения квантования центрироо ванного отсчета х, = (х(1 Ь t) — ш,Д, поступившего с выхода блока 3 вычитания через открытый первый вход первого коммутатора 4.

На этапе обучения квантование отсчетов в АЦП 5 осуществляется при подключенниом к его опорному входу эталонном напряжении E (через первый вход коммутатора 23). Сформио о рованный код х, = x(l dt) центрированного сигнала подается на информационные входы блоков 24 и 25 определения размаха. Тактовый импульс с выхода генератора 8 импульсов осуществляет сдвиг единичного импульса в регистре 13 микрокоманд во второй разряд, формируя микрокоманду С2, по которой разрешаются о запись кода х, в блок 24 определения pasMaxa и срабатывание триггера

17, открывающего второй вход коммутатора 4, подключенный к вьжоду

5 усреднения 2, в котором определяется текущее среднее значение m входного сигнала.

Следующий тактовый импульс с вью хода генератора 8 импульсов формирует в регистре 13 микрокоманд следующую микрокоманду-СЗ, которая разрешает квантованиеш„, подключенного к информационному входу АЦЙ 5 с выхода усреднителя 2, а следующая микрокоманда С4 разрешает запись кода текущего среднего значения m в.

В блок 25 определения размаха и возвращает триггер 17 в исходное состояние, при котором открыт первый вход коммутатора 4. После этого по микрокоманде С5, поступающей íà управляющий вход делителя 26 кода, осуществляется деление кодов, сформированных в блоках 24 и 25 определения размаха:

8, = х

К„, После этого с выхода управляемого делителя 9 частоты через элемент HJIH

10 на управляющий вход устройства

1 выборки и запоминания подается импульс, осуществляющий сброс в нуль предыдущего отсчета сигнала и выбор-! ку и запоминание следующего отсчета х . Этот же импульс поступает однок временно на входы С и Я регистра 13 микрокоманд, вследствие чего он записывается в первый разряд, т.е. форми40, рует новую микрокоманду Cl, а в счетчике 30 выборок устанавливается код

n = 2. После этого последовательность тактовых импульсов с выхода генератора 8 импульсов, поступая на тактовый вход С, формирует в регистре

13 микрокоманд последовательно микрокоманды С2, СЗ, С4 и С5, которые разрешают преобразование центрированного о отсчета х и значения m„ в соответствующих блоках, аналогично предыдущей последовательности микрокоманд о L для отсчета xi и значения mq„

В результате обработки N выборок на вход делителя 26 кода подается код текущего отношения размахов, сформированных на интервале N анализа

5 !43445

При этом в счетчике 30 выборок образуется импульс переполнения, который поступает на управляющий вход блока 27 сравнения, разрешая провер С 5 ку выполнения условия 8 g При невыполнении этого условия (т.е. выполняется соотношение о 7 811) на первом выходе блока 27 сравнения формируется импульс, сбрасывающий в ну- 10 левое состояние содержимое блоков

24 и 25 определения размаха, а в управляемом делителе 9 частоты устанавливается новое значение коэффициента деления частоты тактовых импульсов, соответствующее увеличению шага дискретизации дС входного сигнала в блоке 1 выборки и запоминания.

После этого вновь выполняется цикл обучения для новых 111 выборок сигнала. 20

На каждом цикле обучения счетчик

30 выборок через каждБ е И отсчетов сигнала формирует импульс переполнения, являющийся микрокомандой "Сравнение" для блока 27 сравнения. 25

После нескольких циклов корректировки шага, дискретизации достигается выполнение условия некоррелированности отсчетов, при котором вычисляемая оценка k(x ) плотности распреде- gp ления является статистически достоверной с необходимой точностью (допустимой погрешностью 3„ ). При этом осуществляется адаптивный выбор не только частотного диапазона анализа

ЙР = —, но и амплитудного диадС пазона, так как при выполнении условия 3 .й 8д не осуществляется сброс информации об амплитудных свойствах 4р сигнала, накопленной в блоке 24 определения размаха В „ = 11х „1+ х„,„).

D(m„) л

Механизм осуществления такой адаптивной процедуры основан на стремлении статистического анализатора обес- 45 печить некоррелированность отсчетов сигнала. Это достигается периодической проверкой выполнения определенных математических зависимостей между параметрами входного сигнала и 50 параметрами вычисляемой плотности распределения W(x ).

Иэ соотношения между дисперсией

D1, входного сигнала и дисперсией оценки математического ожидания

0(ш „) этого сигнала

5 следуют соотношения для среднеквадратических отклонений 8,„ G„ пропорциональных им размахов В и R (фиг.3) соответственно для плотности распределения оценки математиче-! ского ожидания Жд и плотности распределения амплитуд сигнала У„

1 (2)

Я ) т

Однако соотношения (1 ) и (2) для параметров справедливы только для случая некорреллированных отсчетов сигнала. Поэтому, используя эту закономерность, можно задавать объем выборок (в счетчике 30) и уровень допустимой (пороговой) относительной погрешности (код 8„ в устройстве 27 сравнения) (3) х

Следовательно, процедура контроля текущего отношения размахов путем сравнения текущего значения 8 с пороговым 8 позволяет за несколько циклов обучения адаптивно регулировать степень коррелированности отсчетов до обеспечения их некоррелированности: адаптивное регулирование шага дискретизации д1 по микроко1 мандам "подстройка", которые подаются с первого выхода блока 27 сравнения на управляющий вход управляемого делителя 9 частоты, устанавливается такой шаг дt. = gt, при ко. тором для определяемых в блоках 24, 25 и 26 параметров (2) и (3) начинает выполняться условие т- п (4) Причем одновременно с операциями подстройки оптимального шага и определения текущего значения разл маха R „осуществляется выделение минимального и максимального уровней етсчетов, подключенных к выходам 2 и 3 блока 24 определения размаха.

Блок 24 (25) определения размаха работает, например, следующим обра. зом (фиг. 2)

Отсчеты входного сигнала поступают одновременно на элементы И 32 и 35 и на соответствующие входы элементов

34 и 37 сравнения, настроенные на о о срабатывание (при х„ х, р в элемен34453 н о м

7 14 о о те 34, а при х„а х„в элементе 37 сравнения). Результат текущего сравнения поступает на один из входов элемента И 32 (35) и на тактовый вход регистра 33 (36), в котором запоминается новое экстремальное значение отсчета, а через элемент ИЛИ

38 управляющий импульс разрешает операцию сложения в сумматоре 39, в котором„формируется новое значение размаха R „. Перед началом работы и в течение циклов обучения элементы устанавливаются в нулевое состояние по импульсам пуска или обнуления, поступающим через элемент ИЛИ 40. При выполнении условия 3 = — Я в

В т R — и х регистрах 33 и 36 остаются записанные экстремальные значения (минимальные и максимальное) кодов от.счетов, необходимые для использования на цикле оценивания.

Выполнение условия (4) в блоке

27 сравнения характеризует окончание циклов обучения: на втором выходе блока 27 сравнения формируется микрокоманда "Начало оценивания", которая подается на второй вход триггера 18.

На этапе оценивания опорный вход

АЦП 5 подключен к выходу цифроаналогового преобразователя 22, закрыт элемент И 14 и открыты элемент: И 31 и элемент И 19, через который импульс дискретизации подается для записи в первый разряд регистра 29 микрокоманд.

Одновременно первый импульс дискретизации с выхода управляемого делителя 9 частоты. поступает на вход счетчика 30 выборки, устанавливая его содержимое равным n = 1, и на уп— равляющий вход устройства 1 выборки и запоминания, на выходе которого формируется отсчет сигнала х, = х х(1 ht ). Блок 3 вычитания осуществляет центрирование с использованием оценки математического ожидания

m„:х(1Ь1

Аналогично формируется ряд отсчетов на интервале анализа 11,N>.

1.

С выхода блока 3 вычитания каждый о центрированный отсчет (2 ° x„) через, коммутатор 4 подается на информационный вход АЦП 5 и на вход дискрими10

35 ао

55 атора 20 знака, который срабатывает, ткрывая первый или второй вход элеента 2И-ИЛИ 21, в зависимости от зна. ка центрированного отсчета. При этом к опорному входу АЦП 5 через элемент

2И-ИЛИ 21, цифроаналоговый преобразователь 22 и коммутатор 23 подключается напряжение х „или х;„с выхода 2 или 3 блока 24 (фиг. 2), задающее амплитудный диапазон квантоо вания отсчета х„ в АЦП 5.

По микрокоманде Cl с первого выхода регистра 29 микрокоманд в АЦП 5 осуществляется преобразование отсчета х„ в код, который подключен к адресному входу блока 6 памяти. При сдвиге импульса в регистре 29 микрокоманд во второй разряд формируется микрокоманда С2, которая поступает на первый управляющий вход блока Ь памяти, разрешая считывание предыдущего содержимого ячейки, определяемой кодом адреса А ®„, и его запись в суммирующий счетчик 7. Следующий тактовый импульс, поступающий на вход с регистра 29 микрокоманд, формирует микрокоманду СЗ, разрешающую операцию суммирования в суммирующем счетчике ?: его содержание увеличивается на единицу и становится равным

У + 1. По микрокоманде С4 раэреша 1 ется запись этого кода с выхода суммирующего счетчика 7 в ту же и-ю ячейку блока 6 памяти, а по микрокоманде С5 осуществляется обнуление суммирующего счетчика ?.

Таким образом, дискриминатор 20 знака осуществляет подключение к опорному входу АЦП 5 определенного уровня напряжения (х,„) или !х,„ „ с выхода АЦП 22 (в зависимости от знака ь анализируемого текущего отсчета х сигнала), т.е. управляет шагом квантования по амплитуде и соответствующим распределением этих данных в определенные ячейки блока памяти 6 (фиг. Зб). Благодаря этим операциям в блоке 6 памяти накапливается оценка л функции Ж„(х ) по некоррелированным отсчетам, т.е. с высокой статистической достоверностью и в широком диапазоне амплитуд анализа вследствие накапливания отдельно левой и правой ветвей гистограммы (фиг.Зб).

После обработки N-го отсчета сигнала в счетчике 30 выборок формируется импульс переполнения, который поступает через открытый элемент И 31

9 143445 на выход адаптивного анализатора, являясь микрокомандой "Вывод - для счиывания результатов анализа функции

g 4c) на регистраторы или последующие блоки измерительно-вычислитель5 ного комплекса.

Следовательно, выполнение этапов обучения и оценивания позволяет не „ только с высокой оперативностью иск»- )p лючить погрешности возможной. коррелированности входных отсчетов (как в прототипе), но и одновременно обеспечить адаптацию к диапазону амплитуд с распределением кодов в цифровой блок 6 памяти при помощи цифровых элементов ° распределения кодов и формирования микрокоманд "Квантование"

"Считывание", "Суммирование", "Запись". При этом максимально эффектив- gp но используется, информация, полученная на этапе обучения. формула изобретения

l. Адаптивный статистический ана- 25 лизатор, содержащий генератор импульсов, выход которого подключен к информационному входу управляемого делителя частоты, счетчик выборок, выход которого подключен к первому Зр ! входу первого элемента И, выход которого является выходом окончания оценки распределения анализатора, 1 аналого-цифровой преобразователь, выход которого подключен к информационному входу первого блока определения размаха и к адресному входу блока памяти, информационный вход которого соединен с выходом суммирующего счетчика, информационный вход 4р которого подключен к выходу блока памяти и является информационным выходом анализатора, выход размаха первого блока определения размаха соединен с входом делителя блока деле- 4б ния, вход делимого которого подключен к выходу размаха кода второго блока определения размаха, выход блока деления соединен с первым входом блока сравнения, второй вход ко- 5р торого является входом задания первого анализатора, о т л и ч а ю— шийся тем, что, с целью повьппения точности оценки плотности распределения вероятностей, в него введены два регистра микрокоманд, усредкитель, первый и второй триггеры, первый, второй и третий элементы И, с первого по пятый элементы ИЛИ, элеl0 мент задержки, блок выборки и запоминания, блок вычитания, первый и второй коммутаторы, дискриминатор знака, элемент 2И-ИЛИ, цифроаналоговый преобразователь, при этом выход первого элемента ИЛИ подключен к счетному входу счетчика выборок и к тактовому входу блока выборки и запоминания, информационный вход которого является информационным входом анализатора, выход блока выборки и запоминания соединен с входом уменьшаемого блока вычитания и с входом усреднителя; выход которого подключен к входу вычитаемого блока вычитания и к первому информационному входу первого коммутатора, второй инфор1 мационный вход которого соединен с выходом блока вычитания, выход первого коммутатора подключен к информационному входу дискриминатора знака и к информационному входу аналогоцифрового преобразователя, опорнйй вход которого соединен с выходом второго коммутатора, первый информационный вход которого является входом эталонного источника анализатора, первый и второй адресные входы первого коммутатора подключены соответственно к прямому и инверсному выходам первого триггера, вход установки в "О" которого соединен с выходом второго элемента ИЛИ, первый вход которого объединен с входом установки в "О" второго триггера и является пусковым входом анализатора, выход элемента задержки подключен к первому входу третьего элемента ИЛИ и к прямому входу второго элемента И, выход которого соединен с входом записи первого регистра микрокоманд, тактовый вход которого подключен к вьжоду генератора импульсов, инверсный вход второго элемента И объединен с первым адресным входом второго коммутатора, с первым входом третьего элемента И и соединен с инверсным выходом второго триггера, прямой выход которого подключен к второму адресному входу второго коммутатора, второй информационный вход которого соединен с выходом цифроаналогового преобразователя, тактовый вход аналого-цифрового преобразователя подключен к вьжоду четвертого элемента ИЛИ,первый вход которого соединен с первым вьжодом первого регистра микрокоманд, ка.

11

143445 второй выход которого подключен к входу установки в "1" первого триггера и к тактовому входу первого блока определения размаха, выход максимального кода которого соединен с первым входом элемента 2И-ИЛИ, выход которого подключен к входу цифроаналогового преобразователя, вто- . рой вход элемента 2И-ИЛИ соединен с 0 положительным выходом дискриминатора знака, отрицательный выход которого соединен с третьим входом элемента

2И-ИЛИ, четвертый вход которого подключен к выходу минимального кода пер-1Б вого блока определения размаха, входы установки в "0" первого и второго блоков определения размаха объединены и подключены к выходу "Больше" блока сравнения, тактовый вход которого соединен с выходом переполнения счетчика выборок, второй вход четвертого элемента ИЛИ подключен к третьему выходу первого регистра микрокоманд, четвертый выход которого соеди- 25 нен с вторым входом второго элемента

ИЛИ и тактовым входом второго блока определения размаха, информационный вход которого подключен к выходу аналого-цифрового преобразователя, выход первого элемента ИЛИ соединен с входом элемента задержки и с вторым входом третьего элемента И, выход которого подключен к тактовому входу дискриминатора знака, к перво35 му входу пятого элемента ИЛИ и входу записи второго регистра микрокоманд, тактовый вход которого соединен с выходом пятого элемента ИЛИ, второй вход которого подключен к выходу генератора импульсов,третий вход четвертого элемента ИЛИ соединен с первым выходом второго регистра микрокоманд, второй выход которого подключен к входу считывания блока па- 4 мяти и тактовому входу суммирующего счетчика, счетный вход которого соединен с третьим выходом второго регистра микрокоманд, четвертый выход которого подключен к входу считывания блока памяти, пятый выход второго регистра микрокоманд соединен с входом установки в "0" суммируюшего счетчика, выход "Меньше или

3 1 2 равно" блока сравнения подключен к входу установки в "1" второго счетчика, инверсный выход которого соединен с вторым входом первого элемента И, выход управляемого делителя частоты соединен с первым входом первого элемента ИЛИ, второй вход которого является входом пуска анализатора, пуска блока определения размаха, выход "Больше" блока сравнения подключен к тактовому входу управляемого делителя частоты.

2. Анализатор по п. 1, о т л и— ч а ю шийся тем, что блок определения размаха содержит два элемента И, два регистра, два элемента сравнения, два элемента ИЛИ и сумматор, выход которого является выходом размаха кода блока, тактовыми входами которого являются входы элементов сравнения, выход Меньше™ первого элемента сравнения соединен с входом синхронизации первого регистра и с первыми входами первого элемента ИЛИ и первого элемента И, выход которого соединен с информационным входом первого регистра, выход которого соединен с первыми входами сумматора и первого элемента сравнения и является выходом максимального кода блока, входом установки в 0" которого является первый вход второго элемента ИЛИ, выход которого подлючен к входам сброса регистров и сумматоров, выход второго регистра соединен с первым входом второго элемента сравнения и с вторым входом сумматора и является выходом минимального кода блока, информационный вход которого подключен к вторым информационным входам первого и второго элементов сравнения и к второму входу первого и первому входу второго элементов И, выход которого соединен с информационным входом второго регистра, выход "Больше или равно" втЬрого элемента сравнения подключен к второму входу первого элемента ИЛИ, к входу синхронизации второго регистра и к второму входу второго элемента И, второй вход второго элемента ИЛИ является входом пуска бло1434453 1434453

Составитель Л.Григорьян- Чтенц

Техред М.яндык КорректорВ. Бутяга

Редактор Л. Веселовская

Тираж 704

Подписное

Заказ 5557/52

В Государственного комитета СССР по делам изобретений и открытий

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

Производственно-полиграфическое предприятие, г, Ужгород, ул. Проектная, 4