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

 

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

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

РЕСПУБЛИК

ЗС59 Н 03 К 13 02

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

ПО ДЕЛАМ ИЗОБРЕТЕНИЙ И ОТНРЫТИЙ (21 ) 3427512/21-06 (22) 14 ° 04. 82 (46 ) 30.05. 84. Бюл. Р 20 (72) В.A.Çàéöåâ, С.К.Турченкова, Н.В.Попенко и A.Á.Ùàäðèí (71) Ленинградский ордена Ленина электротехнический институт им. В. И.Ульянова (Ленина).и СевероЗападный заочный политехнический институт (53) 681.325(088.8). (56) 1. Авторское свидетельство СССР

М 660239, кл. Н 03 К 13/17, 1979.

Авторское свидетельство СССР а 767962, кл. H 03 К 13/02, 1980. (5 4 ): (.57 ) 1. АДАПТ ИВНЫИ ВРЕМЕННОИ

ДИСКРЕТИЗАТОР, состоящий из L каналов- сжатия, каждый из которых содержит первый преобразователь напряжение-код, вторые . входы-выходы которого соединены с .первыми выходом-входом шестого управляюцего регистра, первый блок шин, пятые вход-выход кото-. рого соединены через пятый управляющий регистр с первыми выходом-входом, первой микроЭВМ, вторые входы-выходыо каналов сжатия соединены через тринадцатые управляющие регистры соответственно с первыми L --тыми выходомвходом второго блока шин,(Ь +11 -е вход-выход которого соединены через двенадцатый управляющий регистр с выходом-входом пятой микро ЭВМ (блок программного управления), о т л ич а ю шийся тем, что, с целью повышения точности и быстродействия многопараметрического адаптивного преобразования формы случайных сигналов, в каждый из каналов сжатия дополнительно введенй управляемый делитель напряжения, управляемый усилитель, первый и второй блоки дискриминаторов уровней, второй преобразователь напряжение-код, преобразователь код-напряжение, первый и второй вентили, первый и второй ре„„SU„„1095386 А вер ивные счетчики, первый-четвертый, седьмой-одиннадцатый управляющие регистры, вторая, третья, четвертая микроЭВМ (процессорыj, причем первый вход управляемого делителя напряжения является первым входом канала сжатия, второй вход управляемого делителя напряжения соединен с первым выходом первого реверсивного счетчика, выход управляемого дели" теля напряжения параллельно соединен с первым входом управляемого усилителя и входом первого блока дискриминаторов уровней, выход которого соединен параллельно с вторыми входами первого вентиля и первого реверсивного счетчика, первый вход которого Я соединен с. выходом первого вентиля, второй вход управляемого усилителя . фф соединен с выходом преобразователя M4 код-напряжение, третий вход соединен с первым выходом второго реверсивного счетчика, выход соединен параллельно с входами первого и второгр преобразователей напряжениекод и второго блока дискриминаторов уровней, выход которого паралельно соединен с вторыми входами второго вентиля и второго реверсивного счет чика, первый вход которого соединен с выходом второго вентиля, выход второго преобразователя напряжениекод соединен параллельно с первыми входами первого и второго вентилей ф"» и третьего управляющего регистра, вторые выходы первого и второго реверсивных счетчиков и вход преобразсвателя код-напряжение соединены соответственно с первыми входами 3 первого, второго и выходом четвертого управляющих регистров, вторые входывыходы первого-четвертого управляющих регистров соединены соответственнно с первыми-четвертыми выходамивходами первого блока шин, первые и вторые входы-выходы второй микроЭВФ соединены через седьмой и восьмой

1095386 управляющие регистры соответственно с вторыми выходами-входами первой н третьей микроЭ1 М,, первые вход-выход которой соедийены с вторыми выходомI входом шестого управляющего регистра, первые вторые и третьи входы-выходы четвертой микроЭВМ соединены через

Изобретение относится к средстваь. электроиэмерительной техники, в частности к измерительно-вычислительным устройствам комплексной статистической обработки, предназначено для многопараметрического адаптивного преобразования формы кваэистационарного широкополосного случайного сигнала с произвольной формой графиков текущих оценок спектральной плотности и характеристической функции на участках стационарности, и может быть использовано в распределенных системах сбора и обработки телеметрической информации, системах испытания и диагностики сложных объектов в машиностроении, самолетостроении, системах управления сложными технологическими процессами в химической и мета . ургической промышленности для решения задачи повыаения эффек- 20 тивности (точности и быстродействия1 адаптивного преобразования формы сложных сигналов.

Известен адаптивный временной дискретиэатор, содержащий канал сжа- 25 тия, состоящий из преобразователя напряжение-код, вход которого соединен с входной клеммой канала сжатия, арифметический блок и блок программного управления (1) . 30

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

Наиболее близким по технической сущности к предлагаемому является адаптивный временной дискретизатор. содержащий Ь каналов сжатия. вторые входы-выхожу которых соединены через управляющие регистры с первыми

1 -тыми выходами-входами второго блока шин.(Ь +1) -е вход-выход которого соединены с выходом-входом второго блока программного управления через его управляющий регистр. причем каждый канал сжатия содержит преобразо- 50 ватель напряжение-код, первый вход девятый, десятыи и одиннадцатый управляющие регистры соответственно с третьими выходами-входами первой, второй и третьей микроЭВМ, а чет. вертые вход-выход являются вы ходом и вторым входом канала сжатия. которого является первым входом канала сжатия, первый выход является первым выходом канала сжатия, третий вход соединен с первым входом преобразователя интервал-код, первый выход которого является четвертым выходом канала сжатия, третьим выходом которого является третий выход преобразователя код-интервал, вторые входвыход преобразователя напряжение-код, преобразователя код-интервал, преобразователя интервал-код, входы-выходы арифметического блока и первого блока программного управления соединены через их управляющие регистры соответственно с первыми-пятыми входами-выходами первого блока шин, шестые вход-выход которого являются вторыми входом-выходом канала сжатия (2).

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

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

Поставленная цель достигается тем, что в адаптивный временной дискриминатор, состоящий из Ь каналов сжатия, каждый иэ которых содержит первый преобразователь напряжениекод, вторые входы-выходы которого соединены с первыми выходом-входом шестого управляющего регистра, первый блок шин, пятые вход-выход которого соединены через пятый управляющий регистр с первыми выходом-входом первой микроЭВМ,. вторые входы-выходы каналов сжатия соединены через тринадцатые управляющие регистры соответственно с первыми Ь -тыми выходомвходом второго блдка шин,(L, +1) -e вход-выход которого соединены через двенадцатый управляющий регистр с выходом-входом. пятой микроЭВМ, . (блок программного управления), до1095386 полнительно введены в каждый канал сжатия управляемый делитель напряжения, управляемый усилитель, первый и второй блоки дискриминаторов уровней, второй преобразователь напряже)ние-код, преобразователь код-напря- 5 жение, первый и орой вентили, пер.вый и второй реверсивные счетчики, первый-четвертый, седьмой-одиннадцатый управляющие регистры, вторая, третья, четвертая микроЭВМ (про- 10 цессоры), причем первый вход управляемого делителя напряжения является первым входом канала сжатия, второй вход управляемого делителя напряжения соединен с первым выходом перво- 15

ro реверсивного счетчика, выход управляемого делителя напряжения параллельно соединен с первым входом управляемого усилителя и входом первого блока дискриминаторов уровней, выход которого соединен параллельно с вторыми входами первого вентиля и первого реверсивного счетчика, первый вход которого соединен с выходом первого вентиля, второй вход управяемого усилителя соединен с выходом преобразователя код-напряжение, третий вход соединен с первым выходом второго реверсивного счетчика, выход соединен параллельно с входами перво го и второго преобразователей напряжение-код и второго блока дискриминаторов уровней, выход которого параллельно соединен с вторыми входами второго вентиля и второго реверсивного счетчика, первый вход которого соединен с выходом второго вентиля, выход второго преобразователя напряжение-код соединен параллельно с первыми входами первого и второго вентилей и третьего управляющего 40 регистра, вторые выходы первого и второго реверсивных счетчиков и вход преобразователя код-напряжение соединены соответственно с первыми входами первого, второго и выходом чет- 45 вертого управляющих регистров, вторые входы-выходы первого-четвертого управляющих регистров соединены соответственно с первыми-четвертыми выходами-входами первого блока шин, йервые и вторые входы-выходы второй микроЭВМ - соединены через седьмой. и восьмой управляющие регистры соответственно с вторыми выходами-входами первой и тРетьей микРоЭВМ, первые вход-выход которой соединены с вторыми выходом-входом шестого управляющего регистра, первые, вторые и третьи входы-выходы четвертой микроЭВМ соединены через девятый, десятый и одиннадцатый управляющие регистры соответственно с третьими выходамивходами первой, второй и третьей микроЭВМ, а четвертые вход-выход являются выходом и вторым входом канала сжатия. 65

На чертеже представлена электрическая структурная схема адаптивного временного дискретизатора.

Адаптивный временный дискретизатор состоит из L каналов сжатия 1 —

1, каждый из которых содержит управляемый делитель напряжения 2, первый вход которого является первым входом канала сжатия, первый блок 3 дискриминаторов уровней, управляемый усилитель 4, второй блок 5 дискриминаторов уровней, первый реверсивный счетчик 6, второй реверсивный счетик 7, второй преобразователь 8 напряжение-код, первый 9 и второй 10 вентили, преобразователь 11 коднапряжение, первый преобразователь

12 напряжение-код, первый 13, второй

14, третий 15 и четвертый 16, управляющие регистры, первый блок 17 шин, пятый 18 и шестой 19 управляющие регистры, первую микроЭВМ 20, седьмой управляющий регистр 21, вторую микроЭВМ 22, восьмой управляющий регистр 23, третью микроЭВМ 24, девятый 25, десятый 26 и одиннадцатый

27 управляющие регистры, четвертую .микроЭВм 28, четвертые входы-выходы которой являются соответственно вторым входом и выходом канала сжатия каждого из Ь тринадцатых управляющих регистров 29„ -29„, пятой микроЭВМ

30, выполняющей функции блока программного управления, двенадцатого управляющего регистра 31 и второго блока шин 32. Причем второй вход и выход каждого из L каналов сжатия

1 -1 соединены через соответствующие тринадцатые управляющие регист- " ры 29 -29 соответственно с первыми

L-тыми выходами-входами второго блока 32 шин,(L +1) -е вход-выход которого соединены через двенадцатый управляющий регистр 31 с выходом-входом пятой микроЭВМ, (блок программного управления5 30. Второй вход управ-, ляемого делителя напряжения 2 соединен с первым выходом первого реверсивного счетчика 6, выход соединен параллельно с первым входом управляемого усилителя 4 и входом первого блока 3 дискриминаторов уровней, выход которого соединен с вторыми входами первого вентиля 9 и первого реверсивного счетчика 6, первый вход которого соединен с выходом первого вентиля 9, второй вход управляемого усилителя 4 соединен с выходом преобразователя 11 код-напряжение, третий вход соединен с первым выходом второго реверсивного счетчика 7, выход параллельно соединен с входами первого 12 и второго 8 преобразователей напряжение-код и второго блока 5 дискриминаторов уровней, выход которого соединен с вторыми входами вто рого вентиля 10 и второго реверсивного счетчика 7, первый вход которо1095 386

ro соединен с выходом второго вентиля 10, вторые вход-выход первого преобразователя 12 напряжение-код .оединены через шестой управляющий регистр 19 с первыми выходом-входом третьей микроЭВМ 24, выход второго преобразователя 8 напряжение-код соединен параллельно с первыми входами первого 9 и второго 10 вентилей, вторые выходы первого б и второго 7 реверсивных счетчиков, выход второго преобразователя 8 и вход преобразователя 11 код-напряжение соединены соответственно с первыми входами первого 13, второго 14, третьего 15 и первым выходом четвертого 16 управ-IS ляющих регистров, вторые вход-выход которых соединены соответственно с первыми-четвертыми выходами-входами, первого блока шин 17, пятые вход.выход которого подключены через пя- 2О тый управляющий регистр 18 к первым

Ъ. выходу-входу первой микроЭВМ 20, первые, вторые и третьи входывыходы второй микроЭВМ 22. соединены соответственно через седьмой 21, 25 восьмой 23 и десятый 26 управляющие регистры с вторыми выходами-входами первой 20, третьей 24 и четвертой

28 микроЭВМ., первые и третьи входывыходы которой соединены через де- 30 вятый 25 и одиннадцатый 27 управляющие регистры соответственно с третьими в:;ходами-входами первой 20 и третьей 24 микроЭВМ °

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

Сигнал (Опоступает на вход управляемого делителя напряжения 2, на выходе которого формируется сигнал

i)(t) =h„> х(1) где Ь„, =1 1/2 1/3 40

1/15 — масштабный коэффициент.

Сигнал (1) с выхода управляемого делителя напряжения 2 поступает на входы управляемого усилителя 4 и перво- го блока 3 дискриминаторов уровней, которые срабатывают при превышении сигналом у(4}границ диапазона(=ф, При этом первый блок 3 дискриминаторов уровней формирует импульс

"+1" на втором суммирующем входе первого реверсивного счетчика 6. а также блокирует Io BTQDQMv BxoIIv первый вентиль 9. На выходе управляемого усилителя 4 формируется сигнал

z() =h«.q(O =h„"h2 У(Ца где

II,) — — 1. 2, ..., 15 — масштабный ко- 55

2". эффициент. Сигнал Я (О с выхода управляемого усилителя 4 поступает на входы второго преобразователя 8 напряжение-код, пебвого преобразователя 12 напряжение-код. второго 60 блока 5 дискриминаторов уровней, которые срабатывают при превышении сигналом и (t) границы диапазона (3 г, а,.), и формируют импульс "1" на втором вычитающем входе второго De- 65 версивного счетчика 7 и ыа это время блокируют по второму входу второй вентиль 10.

Второй преобразователь 8 напряжение-код преобразует с грубым шагом квантования по амплитуде Ы сигнал

2 (1} =h h2q < (tj в с„щественные отсчеты i2; f с шагом, минимальное значение которого ь с ® 10 мкс. Величина

Ь t определяется суммой переходных процессов в управляемом делителе напряжений 2 (измерительных ключей, управляющих "сеткой" образцовых сопротивлений ), в управляемом усилителе 4, в первом 3 и втором 5 блоках дискриминаторов уровней (определяется временем "гашения" измерительных усилителей, работающих в момент компарирования в режиме возбуждения), а также во втором преобразователе 8, который также содержит "сетку" параллельно работающих измерительных усилителей, работающих в момент KOMпарирования в режиме возбуждения.

При передаче отсчетов 1a; ) с выхода второго преобразователя 8 напряжение-код на вход третьего управляюще(го регистра 15 формируются сопровож дающие импульсы "-1", поступающие через первый вентиль 9 на первый (вычитающий } вход первого реверсивного счетчика б, и "+1", поступающие с выхода второго вентиля 10 на первый (суммирующий) вход второго реверсивного счетчика 7, масштабные коэффи циенты и h д — в первый 13 и второй 14 уйравляющие регистры. Такое построение обеспечивает ослабление или усиление текущего размаха сигнала

8 И) в следящем режиме за интервал времени At, причем адаптация к выходу за граничные значения(уш „= 1) выполняется за интервал времени порядка 0,3 64 . Масштабные коэффициенты

h„, h „1 и отсчеты (2„ ) с первоготретьего управляющих регистров 1315 через первый блок шин 17 и пятый4 управляющий регистр 18 поступают с интервалом а1 в память данных первой микроЭВМ 20. В каждом цикле длительностью /д1 по программам, поступающим с памяти программ первой микроЭВМ 20, в ней вычисляются текущие оценки средних значений существенных отсчетов (г; ) и их интервалов I B. по рекурентным алготС) ритмам î 1-1+2 (-1с 1-1) ье,=ье,. „+2 (ae; -ьО.,„) где m;, tn;< и ь 8; и а 8ÄÄ вЂ” о енки текущих средних значений в i-м цикле

1095386 по множествам (z ), (д8

2 (» — объем выборки отсчетов при усреднении, с =5-14; лЬ =а8

=4((-1 .

1 о °

На выходе второго преобразователя 8 напряжение-код формируются

5 только существенные отсчеты (z„c) в моменты 1;, когда сигнал z(t) пересекает любой из уровней, равномерно распределенных с шагом b z по всему диапазону r f-2, + а„) . Интервал ь9;g между существенными отсчетами z „1,, (2„.с формируется с помощью первой микроЭВМ: 20. Разрядность кодов коэффициентов Ь„, Sz и отсчетов z — по

ic

4 бит каждый, а кодов интервалов 15 д8; и средйего интервала йО; — но

8 бит ка><дый. Текущая оценка среднего значения rrr; через пятый управляющий регистр 18, первый блок шин 17, четвертый управляющий регистр 16 посту- 70 пает с.интервалом сглаживания порядка

1 50-100) at . на вход преобразователя

11 код-напряжение, на выходе которого формируется сигнал ;, поступающий на второй вход управляемого усилителя 4. Это обеспечивает эффективное использование рабочего диапазона управляемого усилителя 4, на выходе которого получается сигнал

Ъ2М(%1 Х(Ц " ) h2gz(t), 30 где z(t) — центрированный, предварительно масштабированный сигнал. Разрядность оценок rrr; 4 бит. Оценка .д7;=, =at.д6; среднего интервала появления существенных отсчетов ; } позволяет косвенно контролировать текущее значение среднеквадратической частоты спектра квазистационарного сигнала х(1) и обеспечивает адаптивный выбор интервала рекуррентного усреднения * „= 2".д1 a e i (50-100 } a t i, В памяти данных первой микроЭВМ 20 по программам экспресс-обработки размещают сообщение о "гру- 45 бых" отсчетах и через седьмой управляющий регистр 21 передается в память данных второй микроЭВМ 22, где оно используется для вйчисления норми-.ванной функции спектральной 5р плот "сти по алгоритму Фурье-преобразования

П1

rm

S (ju} = — z(iat}cos iatdt-) 3(iat} (+ - 2 2 т (55 о

s n iat at где Г =a — число отсчетов амплитуд

70 функции спектральной плотности 5"(jrrr) 60

Поскольку точность вычисления оценок отсчетов функции плотности распределения W(i) зависит от точности равномерной расстановки по рабочемУ диапазо"У (zr-,z„) yPoaHeA KaaH- 65 тования по амплитуде с шагом az «az, т то первый преобразователь 12 напряжение-код должен быть времяимпульсного типа с фиксатором мгновенного значения сигнала 2(t;I на входе и иметь число уровней квантования 2z„/az не менее 28

Преобразователь 12 напряжение-код преобразует сигнал z(t) в точные отО счеты z;т с постоянным шагомais»5at, который в 5 раз более превышает время преобразования во втором более грубом преобразователе 8 напряжение -код. Отсчеты z; T через шестой управляющий регистр 19 поступают в память данных третьей микроЭВМ 24, в которой вычисляются оценкй %(2; } функции 4+ (z} плотности распределе ния по аМплитуде ф„(z. ) =%„(z;„ )<2 j Ä(;, )- "„(l- т) () где Wg(zir} «((K(zi „т} (и ((-1) — текущие оценки отсчетов функции плотности распределения;

4к(2; ) =i при iC= i - вероятность попадания в интервал 2„ -ат

Для остальных отсчетов функцйи

% „(2; .) в i -том цикле выполйяется процедура "утечки", т.е. в алгоритме (3) значение Ч(Й;,,) =0 при К=(Вычисленная функция %(k;,т) с помощью восьмого управляющего регистра 25 передается в память данных второй микроЭВМ 22, в которой вычисляют характеристическую функцию в Crrr

8 (j(r(;} g % (az)cosy(az dz-) ) % т. ( ф ю т т )

rrn

f re . (Caz,) si 2az, g8z (+) где = m-a g — число отсчетов функции (jy)(re = 2e.

На основе Фурье — преобразования (ункции плотности распределения

% (z )=%+(1дх. ) алгоритмы (2) и (4) вычисления S () ) и 8 () g) аналогичны, причем число отсчетов для обоих равно

2 . Вычисленные функции 5 (j,ù), 9" (), e) через десятый управляющий регистр 26 поступают в память данных четвертой микроЭВМ(28, которой производится вычисление граничных параметров ((йр и у> функций 5 (1 } и e+(jy), при которых "хвосты" этих функций по уровню не превышают az2 и д 2, т. е. г

Inorg )5" (jrd„))(432 рмх) 9 (+ь ) что обеспечивает повышение точности оптимального адаптивного выбора шагов точного квантования сигнала K(t) по амплитуде Ьг =(2 „}"и по времени

s 8 =2(и„;М1." . С учетом того, что

1095386 для алгоритмов вычисления моментных статистических характеристик важно сокращение времени адаптивного выбора интервала сглаживания T„=2 ЫЬО;, а также, что указанные алгоритмы слабо критичны к точности контроля граничных параметров функций 5 (j(Q) и

6+ (j q) исследуемого сигнала z (t) в предложенном устройстве формируется экспресс-сообщение о "грубых" отсчетах P ;) в форме коров: (h ), 10 (1, ), p2„), (й„), fh 8„) (е„, хранящиеся в памятй данных первой микроЭВМ. 20 и передающееся для дальнейшем обработки в память данных четвертой микроЭВМ 28. Точные отсчеты (z „,q) хранятся в памяти данных третьей микроЭВМ 24.

С помощью четвертой микроЭВМ 28 также вычисляются оценки бптимального числа „= (ь2 ) — отсчетов различно ных (дифференциальной, интегральной, условной } функций распределения по вероятности, а также оптимальный интервал Т =2 .Д 6 „.а+ рекуррентного

Qf, усреднения при их вычислении по алго- 5 ритмам типа (3 ), при этом to и Т запоминаются в дальнейших вычислейиях.

В адаптивном временном .дискретизаторе выбор шага rlBС, выполняется прямо по 1г, определяемой по вычисленной во второй микроЭВМ оценке спектра

5(r) входного сигнала, и определяет}

П 1l ся погрешностью вычисления хвоста спектра, т.е., в основном, числом циклов усреднения Й . При К =10„ d" =

=(e„- e}порядка 2 "о =O,1a, при этом 35 инвариантна в отличие от известного устройства в виду распределения входного сигнала.

Выбор pB+ по Хгр в заявляемом устройстве обеспечивает возможность 40 автоматического выбора шага по времени по задержке a@+ = h4<, шага по

1 частоте д1 = - + — в алгоритмах восст4 новления сигнала вычисления корреляционной и спектральной функций сигнала. При этом обеспечивается погрешность от дискретизации порядка 0,1Ъ.

Кроме того, в адаптивном дискретизаторе вычисление плотности распределения вероятности %(} в третьей микроЭВМ и на ее основе характеристической функции 9(jq) во второй микроЭВМ позволяет вычислять граничный параметр д 2д= в четвертой микро2@г

ЭВМ. Оценка z„ o3 er (аналогично оценке bB ) автоматически поддерживать составляющую погрешность от дискретизации по уровню и по вероятности для восстановления функций интегрального, дифференциального и условного распределений вероятностей.

Использование в предлагаемом устройстве автоматического выбора диапазона по амплитуде (1/15-15В} в соответствии с текущим размахом исследуемого сигнала в следящем режиме автоматическое центрирование процесса

1 о усиления обеспечивает оптимальный

ыбор рабочего диапазона управляемого усилителя и первого и второго преоб-! разователей напряжение-код, а также уменьшение общей погрешности адаптивного кодирования в 4-5 раз по сравнению с известным устройством.

Повышение быстродействия в заявляемом устройстве, по крайней мере, в 3 раза обеспечивается за счет сокращения числа участков пробных измерений до одного длительностью по времени 2"© д1 .Распараллеливание операций вычисления оценок %(x} и

5(jf) позволяет дополнительно сократить время выбора в 2,5 раза. Общий выигрыш по быстродействию можно оценить, по крайней мере, в 7,5 раза.

1095386

Составитель В. Ткачев

Редактор М. Келемеш Техред Ж. Кастелевич

Корректор Н.Лазарева

Филиал ППП "Патент", г.Ужгород, ул.Проектная, 4

Заказ 3621/41 Тираж 862 Подписное

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

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