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

На правах рукописи

Линник Михаил Александрович

ИНФОРМАЦИОННО-ИЗМЕРИТЕЛЬНЫЙ КОМПЛЕКС ДЛЯ РЕГИСТРАЦИИ ГИДРОАКУСТИЧЕСКИХ СИГНАЛОВ

05.11.16 - Информационно-измерительные и управляющие системы (техника и технология)

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук

2 3 ОКГ 2014

Хабаровск - 2014

005553739

005553739

Работа выполнена в федеральном государственном бюджетном образовательном учреждении высшего профессионального образования «Тихоокеанский государственный университет»

Научный руководитель: кандидат технических наук, доцент

Бурдянский Игорь Николаевич

Официальные оппоненты: Моргунов Юрий Николаевич, доктор технических

наук, старший научный сотрудник, федеральное государственное бюджетное учреждение науки «Тихоокеанский океанологический институт им. В.И. Ильичева» Дальневосточного отделения Российской академии наук, заведующий отделом технических средств исследования океана

Калинов Геннадий Алексеевич, кандидат технических наук, Дальневосточный филиал федерального государственного унитарного предприятия «Всероссийский научно-исследовательский институт физико-технических и радиотехнических измерений», начальник отдела вторичного эталона единиц времени и частоты

Ведущая организация: Федеральное государственное бюджетное учрежде-

ние науки «Специальное конструкторское бюро средств автоматизации морских исследований» Дальневосточного отделения Российской академии

Защита состоится «26» ноября 2014 г. в 16-00 на заседании диссертационного совета Д 212.294.05 при ФГБОУ ВПО «Тихоокеанский государственный университет» по адресу: 680035, Хабаровск, ул. Тихоокеанская, 136, ауд. 315л.

С диссертацией можно ознакомиться в библиотеке и на сайте ФГБОУ ВПО «Тихоокеанский государственный университет» http://pnu.edu.ru.

Автореферат разослан «10» октября 2014 г.

Ученый секретарь

диссертационного сове^?^ Бурдинский Игорь Николаевич

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность работы. В течение последних десятилетий наблюдается постоянный рост интереса к созданию подводных гидроакустических (ГА) сенсорных сетей, способных обеспечить решение задач сбора океанографических данных, мониторинга загрязнений окружающей среды, геологической разведки, сейсмологического мониторинга и предупреждения стихийных бедствий. Элементами ГА сенсорных сетей могут быть подводные стационарные системы сбора данных, автономные подводные роботы (ПР), подводные маяки и суда обеспечения, совместно выполняющие задачу сбора данных и реагирования на изменение внешней обстановки в рамках целевой миссии.

Основными техническими средствами подводных сенсорных сетей, обеспечивающими их надежную эксплуатацию, являются ГА информационно-измерительные комплексы (ИИК). Корректное функционирование ИИК напрямую зависит от эффективной регистрации ГА сигналов в реальном масштабе времени. Решение задачи регистрации сигналов возложено на специальные программно-аппаратные средства передачи, приема и цифровой обработки ГА сигналов.

Элементы и части ГА ИИК должны проектироваться с учетом ограниченного ресурса батарей, массогабаритных показателей и характеристик ГА канала (ограниченной полосы пропускания, многолучевого распространения, сильного затухания и рассеивания ГА сигнала).

Исследования по проблемам разработки ГА ИИК отражены в работах таких российских и зарубежных ученых, как М.Д. Агеев, Ю.В. Матвиенко, JI.B. Киселев, Б.А. Касаткин, С.Г. Яковлев, А.Ж. Ермоленко, В.А. Акуличев, К.Г. Кебкал, М. Сто-янович, И. Дж. Мил, Р. Банаш, Дж. А. Катиповик, Р.Д. Крист и др.

Наибольший опыт в разработке элементов и частей ГА ИИК имеют научно-исследовательские и промышленные организации Link Quest (США), Evologics (Германия), IXSEA (Бельгия), Kongsberg Maritime (Норвегия), ИПМТ ДВО РАН, ОАО «НИИ гидросвязи «Штиль», ФГУП ОКБ ОТ РАН (Россия). Основными недостатками существующих комплексов являются низкие показатели точности обнаружения момента прихода сигнала и достоверности приема данных при соотношении сигнал/шум (SNR) на входе приемника менее 6 дБ. Повышение помехоустойчивости ГА ИИК позволяет уменьшить энергопотребление передающего оборудования и увеличить дальность действия комплекса.

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

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

Для достижения указанной цели в диссертационной работе были поставлены и решены следующие задачи:

1. Анализ современных алгоритмов и средств цифровой обработки ГА сигналов для решения задачи передачи информационно-измерительных данных (ИИД) в водной среде.

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

3. Программная реализация разработанных алгоритмов обработки ГА сигналов, разработка компьютерной модели ИИК и проведение вычислительных экспериментов.

4. Анализ современных информационных технологий для имплементации основных блоков проектируемого ГА ПИК.

5. Разработка оптимальной структуры ИИК с точки зрения энергоэффективности, массогабаритных показателей и разработка аппаратно-программного обеспечения ГА ИИК.

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

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

Научная новизна:

1. Предложен алгоритм порогового детектирования ГА шумоподобных фазоманипулированных сигналов (ФМШ), позволяющий вести асинхронный прием с высокой точностью обнаружения момента прихода сигнала за счет оценки относительной энергии шума на выходе коррелятора.

2. Разработан и исследован метод передачи ИИД по ГА каналу для ПР, обеспечивающий повышенную помехоустойчивость и снижение энергопотребления тракта излучения за счет использования ФМШ и алгоритма сжатия информации.

3. Разработай алгоритм автоматической адаптации к характеристикам ГА канала для высокоскоростной передачи данных методом многопозиционной фазоча-стотной модуляции (МФЧМ), обеспечивающий снижение относительного числа ошибок (ВЕК) за счет добавления в пакет данных пилот-сигналов.

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

5. Предложена структура унифицированного ИИК, позволяющая в реальном масштабе времени проводить многоканальную регистрацию ГА сигналов и сократить энергетические и массогабаритные показатели системы за счет использования технологии программируемой системы на кристалле (ПСнК).

Практическая ценность результатов работы. Разработан универсальный унифицированный ГА ИИК, позволяющий повысить показатели точности обнаружения момента прихода сигналов, помехоустойчивости и достоверности передачи ИИД в водной среде, уменьшить энергетические и массогабаритные показатели. Разработанный ИИК позволил сократить время создания, отладки и испытаний новых аппаратно-программных средств цифровой обработки ГА сигналов.

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

Реализация результатов работы. Основные результаты работы используются в Институте проблем морских технологий (ИПМТ) ДВО РАН при выполнении научно-исследовательских и опытно-конструкторских работ по разработке систем навигации, телеметрии и управления автономных необитаемых подводных аппаратов серии «ММТ». Результаты работы используются в учебном процессе ТОГУ в рамках дисциплин «Системы передачи данных», «Системы цифровой обработки сигналов», «Многоканальные навигационные системы», а также в курсовом и дипломном проектировании по специальностям «Вычислительные машины, комплексы, системы и сети» и «Многоканальные телекоммуникационные системы».

Связь темы исследования с научными программами. Исследования проводились при поддержке грантов в рамках следующих научно-исследовательских работ: НИР по аналитической ведомственной целевой программе «Развитие научного потенциала высшей школы» (2009-2013 годы) на тему «Высокоскоростная система ГА связи для осуществления телекоммуникационного взаимодействия между ПР» (№ гос. per. 2.1.2/2620); НИР по федеральной целевой программе «Научные и научно-педагогические кадры инновационной России (2009-2013 годы)» на темы «Исследование принципов построения систем наведения и позиционирования ПР с применением современных микроэлекгронных средств на примере матричной логики и цифровых сигнальных процессоров» (ГК № П497 от 13.06.2010), «Математическая модель ГА канала и средство комплексного анализа акустических сигналов на основе экспериментальных и теоретических исследований» (ГК № 14.740.11.1138 от 09.06.2011); НИР по программе У.М.Н.И.К. на тему «Разработка высокоточного модуля ГА навигационной системы для определения наклонной дальности» (ГК № У-2008-6/4 от 28.10.2008, ГК № У-2010-1 /11 от 01.03.2010); НИР тематического плана ФГБОУ ВПО ТОГУ на тему «Разработка учебно-лабораторного стенда для исследования систем телеуправляемых ПР» (№ 4.07, ТОГУ); НИР тематического плана ИПМТ ДВО РАН «Фундаментальные исследования, утвержденные Президиумом РАН» на темы «Исследование и разработка принципов создания многофункциональных и специализированных автономных необитаемых подводных аппаратов на основе прогрессивных технологий и роботизированных систем» (№ гос. per. 01200955364, 2009-2011 гг.); «Исследование приоритетных направлений создания «интеллектуальных» подводных робототехнических комплексов для автоматизированного обследования водной среды и рельефа дна, морских инфраструктур, физических полей океана» (№ гос. per. 0120125673, 2012 г.).

Основные положения, выносимые на защиту:

1. Алгоритм порогового детектирования ГА сигнала на основе вычисления величины отношения двух корреляционных функций по маске излучаемого и неизлу-чаемого ФМШ, обеспечивающий асинхронный прием, защиту от влияния импульсных помех и погрешность обнаружения момента прихода сигнала, не превышающую длительность одного символа ФМШ при SNR меньше 0 дБ.

2. Метод передачи ИИД по ГА каналу с использованием ФМШ и алгоритма сжатая данных, обеспечивающий при его использовании достоверный прием информации с BER не хуже 10"6 при SNR меньше 0 дБ и снижение энергопотребления тракта излучения ПР до 5 раз.

»жлт™, Алгоритм адаптации приемника к характеристикам ГА канала для МФЧМ, Обеспечивающий уменьшение BER до 10 раз при SNR, равном 6 и более децибелам.

4. Алгоритм синхронизации ИИК, обеспечивающий при его использовании снижение погрешности формирования ИИ до 100 раз для кварцевых генераторов с относительной погрешностью формирования частоты 10'5.

5. Алгоритм управления тактовой частотой работы параллельно-конвейерной структуры корреляторов в унифицированном ИИК на основе технологии ПСнК, обеспечивающий при его использовании сокращение энергопотребления комплекса до 10 раз.

Апробация результатов работы. Основные результаты диссертационной работы докладывались и обсуждались на всероссийских и международных конференциях-XXVI конференции памяти выдающегося конструктора гироскопических приборов Н.Н Острякова (г. Санкт-Петербург, 2008 г.); X Всероссийской научно-технической конференции «Теоретические и прикладные вопросы современных информационных технологий» (г. Улан-Удэ, 2009 г.); Ill, IV и V Всероссийских научно-технических конференциях «Технические проблемы освоения мирового океана - ТПОМО» (г Владивосток, 2009 2011, 2013 гг.); 12-ой, 13-ой и 14-ой международных конференциях и выставках «Цифровая обработка сигналов и ее применение DSPA» (г. Москва 2010 2011

2012 гг.); The First Russia and Pacific Conference on Computer Technology ¡nd Applications «RPC 2010» (Владивосток, 2010 г.); Международной научно-практической конференции «Информационные технологии и высокопроизводительные вычисления» (г Хабаровск, 2011 г.); XII и XIII краевых конкурсах молодых ученых и аспирантов (г Хабаровск, 2010,2011 гг.); V международной научной конференции «Проектирование инженерных и научных приложений в среде MATLAB» (Украина, г. Харьков, 2011 г.)-XVIII Санкт-Петербургской Международной конференции по интегрированным навигационным системам (г. Санкт-Петербург, 2011 г.); IEEE Santander Conference and Exhibition «OCEANS 2011» (Испания, г. Сантандер, 2011 г.); Хабаровском краевом молодежном инновационном конвенте (г. Хабаровск, 2010 г.); Международной on-line конференции «Инновационное будущее: мир ждет твою идею!» (Испания г. Мадрид г Москва, г. Хабаровск, 2011 г.); XVII и XVIII Санкт-Петербургских Международных выставках «Высокие технологии. Инновации. Инвестиции» (Hi-Tech) (г. Санкт-Петербург, 2011,2012 гг.); XI Всероссийской выставке научно-технического творчества молодежи «НТТМ-2011» (г. Москва, 2011 г.); VIII и IX Всероссийских научно-практических конференциях «Перспективные системы и задачи управления» (п. Домбай,

2013 г., г. Сочи, 2014 г.). Работа в целом обсуждалась на научных семинарах кафедры

техника» и Лаборатории интеллектуальных технологий и систем 10I У (2008-2014 гг.).

Публикации и личный вклад автора. Основные результаты исследований отражены в 49 печатных работах, среди которых 4 публикации в изданиях, входящих в перечень российских рецензируемых научных журналов, 2 патента на изобретение 20 публикаций в сборниках международных и всероссийских конференций, 19 свидетельств о регистрации программы по тематике исследований.

Личный вклад автора при проведении исследований заключается в разработке методов передачи и регистрации ГА сигналов, структуры и аппаратно-программного обеспечения ГА ИИК, проведении вычислительных экспериментов и натурных испытании.

Структура и объём работы. Диссертация состоит из введения, четырёх глав заключения, списка литературы, приложений. Работа содержит 120 страниц основного текста, 76 рисунков, 22 таблицы, библиографический список из 148 наименований.

СОДЕРЖАНИЕ РАБОТЫ

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

В первой главе определены основные требования и задачи, которые должны решаться ИИК для регистрации ГА сигналов. Под регистрацией подразумевается цифровая обработка ГА сигналов для решения задачи обнаружения и оценивания параметров сигнала при осуществлении информационного обмена по ГА каналу. Информационный обмен может включать в себя задачи обсервации, телеметрии и управления подводными объектами. В работе представлен сравнительный обзор некоторых ГА ИИК различных производителей. Описаны основные недостатки существующих комплексов: каждую задачу информационного обмена выполняет отдельное устройство; ИИК построены на низкопроизводительных при многоканальной обработке сигналов цифровых сигнальных процессорах и специализированных микроконтроллерах; снижение энергетических показателей достигается за счет перехода в спящий режим (отключение комплекса на фиксированное время); заявленные показатели точности обнаружения момента прихода ГА сигнала и BER существующих ИИК обеспечены при SNR на входе приемника, равном б и более децибелам.

В работе проанализирован классический протокол информационного обмена по ГА каналу с совмещением передачи данных обсервации, телеметрии и управления (последовательная полудуплексная передача). Объектами взаимодействия являются ПР, береговая станция, маяки-ответчики и суда обеспечения. Согласно протоколу информационный обмен осуществляется через равные временные интервалы. Каждый временной интервал, в свою очередь, логически разделен на три ИИ по характеру взаимодействия подводных объектов (обсервация, телеметрия и управление).

Проанализированы существующие технологии модуляции, применяемые при передаче ГА ИИД. Показано, что для увеличения помехоустойчивости и обеспечения кодового разделения каналов обработки целесообразно применять модуляцию с расширением спектра методом прямой последовательности.

В работе рассмотрены возможные аппаратные средства для реализации основных блоков цифровой обработки сигналов ИИК. При проектировании ИИК целесообразно использовать технологию ПСнК типа FPGA (Field-Programmable Gate Агтау - программируемая пользователем вентильная матрица), позволяющую в отличие от цифровых сигнальных процессоров реализовать распараллеленную работу нескольких равновесных вычислительных процессов внутри одного устройства. Также использование FPGA позволяет проектировать реконфигурируемую систему, которая может быть быстро наделена дополнительными функциями в зависимости от условий эксплуатации и требований потребителя. Проведен анализ рынка современных FPGA, характеристики которых удовлетворяют требованиям надежности и энергопотребления. Этим требованиям отвечают микросхемы семейста Virtex 4-7 (фирма Xilinx) и Cyclone III—V (фирма Altera).

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

обработка в реальном масштабе времени; сокращение энергетических и малогабаритных показателей; создание единого унифицированного ИИК; использование таких современных микропроцессорных технологий, как ПСнК; компенсация ошибок формирования ИИ для уменьшения погрешности обнаружения момента прихода ГА сигналов.

Вторая глава посвящена разработке компьютерной модели ГА ИИК, состоящей из следующих модулей: 1) модель устройства генерации и регистрации ГА сигналов-2) шумовая модель канала передачи; 3) модель протокола передачи ИИД; 4) модель синхронизации объектов информационного обмена; 5) блок статистической обработки данных. Особенность компьютерной модели ИИК состоит в том, что каждый модуль запрограммирован с учетом последующей имплементации на БРвА.

В работе используются ФМШ, закодированные двоичными псевдослучайными последовательностями. Лучшими по своим корреляционным свойствам (отношение максимума автокорреляционной функции к боковым лепесткам) являются псевдослучайные коды типа М-последовательности и коды Голда.

Основным методом для обнаружения ФМШ является вычисление классической корреляционной функции. Использование правила максимума корреляционной функции позволяет определить момент прихода сигнала. Часто возникает потребность обнаружения момента прихода сигнала в реальном масштабе времени при неизвестном периоде излучения сигнала. Для этого используются пороговые способы детектирования. Обнаружение сигнала по фиксированному значению корреляционной функции является самым простым в реализации пороговым способом, но не обладает адаптивностью, т.е. не подходит для среды с изменяющимися характеристиками (ГА канал). Использование посимвольной корреляционной функции позволяет адаптивно обнаруживать ФМШ, но данный способ имеет погрешность определения момента прихода сигнала, составляющую как минимум длительность одного символа ФМШ.

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

(1)

где & - значение корреляционной функции для искомого сигнала; Ак - значение корреляционной функции по маске неизлучаемого сигнала (ортогонального ФМШ той же длительности, что и искомый сигнал, который не используемого при передаче сигналов); рис ~ настраиваемый параметр.

Для сокращения ложных срабатываний перед нахождением параметра цис необходимо провести медианную фильтрацию над Як и Ак в интервале ЛЬ, равном по длительности двум символам ФМШ. После срабатывания адаптивного порога необходимо найти максимум Л* на интервале ЛЪ. Использование двух корреляционных функций позволяет оценить ЯШ на выходе коррелятора по выражению

Структурная схема предложенного способа детектирования ФМШ по двум корреляционным функциям представлена на рис. 1, где ЫК - число каналов корреляционной обработки, щ - последовательность входных отсчетов обрабатываемого сигнала; тц ~ маска искомого сигнала; тг! - маска неизлучаемого сигнала.

Для исследования способов порогового детектирования сигналов была использована апробированная шумовая модель ГА канала, при которой к сигналу добавлялся шум реального ГА канала с параметрами, соответствующими характеристикам ГА канала в шельфовой зоне Японского моря Приморского края. С помощью данной модели была проведена серия вычислительных экспериментов.

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

р(ЯУЙг) = |, (3)

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

Полное количество вычислительных экспериментов Ь рассчитывалось для каждого ЯУЛ на основе теоремы Муавра-Лапласа, которая позволяет оценить близость частоты р и вероятности Р в схеме Бернулли. В результате моделирования была выбрана относительная частота приема сигнала /о //й = 4, где^/Ь - частота дискретизации; /я - несущая частота сигнала. Данная величина позволяет минимизировать аппаратные затраты при реализации приемника и не снижает точности обнаружения.

На рис. 2 представлены зависимости р от БЫЯ при обнаружении момента прихода сигнала, не превышающего один символ ФМШ с параметрами С? = 127, Г5 = 8 для разных методов порогового детектирования ГА сигнала, где 7у - число периодов ^ на символ последовательности; Сх -число символов в последовательности. Согласно рис. 2, прием сигнала по фиксированному значению корреляционной функции возможен при ЯУЛ* = 7 дБ, для посимвольной корреляции - при Ж/?* = -9 дБ, для способа по двум корреляционным функциям - при ¡¡Ш* = -11 дБ, где - минимальное ¿ТУЯ при р = 1. Вычислительные эксперименты показали, что предлагаемый способ порогового детектирования по двум корреляционным функциям дает выигрыш в £№?* на 2 и более децибел по сравнению с посимвольной корреляцией.

В работе была проанализирована максимальная погрешность обнаружения момента прихода ФМШ в зависимости от способом порогового детектирования на основе двух корреляционных функций. Для ФМШ с Су = 127 и Тэ = 8 погрешность обнаружения момента прихода при £ЛЖ* = -11 дБ не превышает 40 мкс.

Рис. 1

-20.0 -15.0 -10.0 -5.0 0,0 5.0 10.0 эт, дБ

—•—Две корреляционные функции -» -Посимвольная корреляция —Фиксированное значение

Рис. 2

При помощи вычислительных экспериментов для разных ФМШ были определены параметр /л!ес и (табл. 1). Пороговое значение детектора для ФМШ, закодированных М-последовательностями или кодами Голда, определяется числом, кратным степени 2, что уменьшает вычислительную сложность данного способа.

Для классического протокола информационного обмена был разработан метод передачи МИД с использованием ГА ФМШ. В предлагаемом методе информационное сообщение состоит из п подпакетов, каждый подпакет состоит из к информационных слов, между которыми формируется защитный интервал для уменьшения влияния межсимвольной интерференции. В начале каждого подпакета добавлены синхронизирующие слова, используемые для синхронизации и передачи служебной информации (тип пакета, тип абонента и т.д.). В конце информационного сообщения добавлена контрольная сумма, вычисляемая сложением по модулю два всех информационных слов в сообщении. На рис. 3 представлен формат информационного сообщения, где ЗИ - защитный интервал; КС - контрольная сумма; - синхронизирующие слова, ¡е [1, п]\ Ц -слова данных, )е[1,к].

Синхронизирующие слова представляют собой М-последовательности, а данные кодируются кодами Голда. Объем данных, который можно передать с помощью одного информационного слова, определяется как Ши = /о^М^, где № - количество различных последовательностей одной длины, используемых в сообщении. Скорость передачи данных в предложенном способе определяется выражением

г- И« - (4)

Ттл

где ¡д= Су • длительность информационного слова; ¡зи- длительность защитного

интервала; Тизл - длительность излучения; Ыд - количество передаваемых данных в информационном сообщении. При определении момента прихода синхронизирующих слов имеется возможность при известной скорости звука оценивать наклонную дальность и тем самым осуществлять обсервацию.

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

Таблица I

С? Гя=8 7^ = 4

8Ж*, дБ МшС, ОТН. ед. ЯЛЖ*, дБ ОТН. ед.

31 -6 3,98 -4 1,95

63 -8 3,99 -6 1,98

127 -11 4,0 -8 1,99

255 -13 4,0 -10 2,0

511 -15 4,0 -13 2,0

1023 -18 4,0 -15 2,0

Информационное сообщение

Симе ол Подпакет 1 Подпакет

старт ЗИ ¡1 ЗИ ЗИ, ЗИ ( КС стоп

Тизл

Рис.3

послан запрос о повторной передаче всего информационного сообщения. В работе с помощью разработанной модели был проанализирован выигрыш по энергопотреблению для подвижных объектов, у которых ИИД меняются часто. Например, при исследовании подводной акватории ПР движется по траектории типа «галс» с постоянной скоростью и передачей следующих параметров: координат, курса, скорости, высоты, глубины, дальности, пеленга. В работе получена зависимость энергопотребления тракта излучения Рим от разной скорости движения ПР и длительности информационного интервала Тии■ Так, при скорости движения, равной 1 м/с, и при Тии= 10 с энергопотребление тракта излучения сокращается до S раз. Предложенный метод передачи ИИД по ГА каналу защищен патентом на изобретение.

В работе был оценен параметр BER при приеме всех ФМШ, составляющих информационное сообщение с параметрами Ts = 8, fs = 12 кГц, Cs е [63, 127, 255, 511] (рис. 4).

Если SNR на входе приемника равняется 6 и более децибелам, предлагается к каждому синхронизирующему слову информационного сообщения добавлять пакет данных с сигналами МФЧМ, используемый для высокоскоростной передачи данных по ГА каналу. Применение предлагаемого порогового способа детектирования по двум корреляционным функциям, позволяющего оценить SNR на входе корреляционного приемника (2), дает возможность автоматического перехода на МФЧМ. Пакет данных разбивается на сегменты, в начало каждого сегмента добавляется пилот-сигнал, используемый для адаптации параметров приемного тракта к характеристикам ГА канала в полосе рабочих частот. Скорость передачи данных определяется длительностью символа в сегменте и его информационной емкостью. Использование алгоритма адаптации позволяет увеличить показатель BER с 10~3 до 10г4 для SNR > 6 дБ.

Согласно применяемому в работе протоколу последовательной полудуплексной передачи, информационный обмен по ГА каналу осуществляется через равные временные интервалы. Временные интервалы задаются высокостабильными генераторами, а ИИ - стандартными широко применяемыми кварцевыми генераторами. Для уменьшения ошибки обнаружения момента прихода сигналов необходимо как можно точнее формировать каждый ИИ. Относительная погрешность формирования частоты широко применяемых кварцевых генераторов (MEC, Epson) составляет 10"5. В работе проводилось исследование по оценке стабильности работы кварцевых генераторов МЕС и Epson. Распределение ошибки формирования системной тактовой частоты /г исследуемых кварцевых генераторов является близким к нормальному. Рассчитанное значение математического ожидания M(At) составляет порядка 100 мкс, среднеквадратического отклонения a(\t) - около В мкс. Характер погрешности мультипликативный, т.е. возрастает с увеличением длительности ИИ.

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

SNR, дБ

-20 -18 -16 -14 -12 -10 -8 -6 -4 -2 0

---Cs = 511 Cs = 255 -~-Cs=127 .....Cs = 63

~РшГ4

в корректировании счетчика текущего временного интервала cnt в зависимости от рассчитанной длительности предыдущего интервала по выражению

cnt + 1, при cnt ф Д,

cnt + 2, при cnt = Д, sign = О cnt=\ , (5)

cnt, при cnt = Д, sign = 1

О, при cnt — Tn_x

где Д =

-r.J+i

- корректирующий коэффициент; sign -

Ш0„„-тл< о

знаковый коэффициент; Тcarat - константная длина временного интервала, отн. ед.; Tn-i - длина предыдущего временного интервала, отн. ед.

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

Третья глава посвящена синтезу элементов и частей ГА ИИК на базе технологии ПСнК типа FPGA. Синтез схемотехнического дизайна, верификация элементов и частей комплекса производились в САПР Xilinx ISE Design и Altera Quartus на языках программирования VHDL и Verilog.

Рассматриваемые способы обработки ГА сигналов реализованы в блоке регистрации ГА сигналов, который представляет собой специализированный процессор. Основным методов обработки сигналов является способ адаптивного порогового детектирования. При вычислении корреляционной функции с приходом каждого значения обрабатываемого сигнала выполняется скалярное произведение массива (вектора) отсчетов цифровой копии сигнала (эталонный сигнал) с вектором, содержащим новое значение

иN-1 предыдущих значений обрабатываемого сигнала, где N = CS- -Ts -длина об-

fs

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

Для определения момента прихода сигнала с точностью до одного периода /о необходимо вычислить значение корреляционной функции Rk за время 7>, которое должно удовлетворять условию Тр < 1 //о. Следует отметить, что время обработки 7> непосредственно зависит от TV и /г и определяется выражением Тр = N/fi. Исходя из этого, /г > N ■ fD. Например, для обработки сигнала с N = 4064 (Cs = 127, Ts = 8) и /d = 48 кГц необходимо fi > 195 МГц.

Для уменьшения/г предлагается реализация con-уровневой вертикальной каскадной параллельной структуры построения канала корреляционной обработки, представленная на рис. 5. При использовании данной схемы можно за счет увеличения уровней каскада добиться получения значений корреляционной функции каждый такт /Ь вне зависимости от значения/г. Реализация устройства многоканальной обработки достигается путем повторения данных схем с изменением значений эталонного сигнала в памяти RAM.

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

Для выполнения функций корреляционной обработки на микросхемах FPGA Xilinx и Altera использовались специализированные блоки конвейерного умножения и сложения MAC и блоки быстрой памяти для хранения отсчетов входного сигнала и масок. Выражение для определения количества МАС-блоков следующее:

N

Nx+l

Рис. 5

2ышг-2ки*'

(6)

где кьн - разрядность обрабатываемых данных; Ышс - выходная разрядность МАС-блока.

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

(7)

N

где Nmem- размер блока быстрой памяти в битах.

Использование МАС-блоков и блоков быстрой памяти в FPGA при построении блока корреляционной обработки не накладывает жестких ограничений на разрядность входных, выходных и промежуточных данных и тем самым позволяет создавать системы, в которых минимизируются ошибки, вызванные округлением.

В работе показано, что в современных микросхемах FPGA энергопотребление в большей степени зависит от тактовой частоты и в меньшей степени от объема занятой логики. Аппроксимирующая полиноминальная кривая зависимости потребляемой энергии от fi, полученная экспериментальным способом (величина достоверной аппроксимации 0,9989), определяется выражением

y(fi) = -0,0009-fi3 + 0,1415 ■fir2 + 11,739 ■ fi + 125. (8)

Исходя из данных оценки энергопотребления FPGA в зависимости от fi с использованием реализации построения многоканального коррелятора на базе многоуровневой вертикальной каскадной параллельной структуры для более низкой fir, в работе предлагается алгоритм уменьшения энергопотребления ИИК, сущность которого заключается в следующем:

1. В исходном состоянии работают только одноканальный приемник (регистратор стартовой последовательности), сторожевой таймер и блок управления/Ь Одноканальный приемник реализован по параллельно-конвейерной схеме для работы на более низкой /г, позволяющей обрабатывать сигнал в реальном времени.

2. При регистрации запускающей стартовой последовательности система переходит на повышенную/г, запускаются все блоки.

3. По стоповой последовательности или по сигналу таймаута от сторожевого таймера система возвращается в исходное состояние.

Коэффициент энергоэффекгивности (отношение энергопотребления комплекса при работе всех блоков к энергопотреблению ИИК, работающего по предложенному алгоритму) приемного устройства определяется выражением

уУР-Т

ЭЭ -

- + 0,5

(9)

/Гс)-тп+у(/Г)-(т-гп) где Г - интервал оценки энергопотребления; Тп- время работы устройства при приеме информационного сообщения;//*™ - минимальная fi работы одноканального приемника во время ожидания прихода информационного сообщения; frmKC - максимальная /г работы устройства во время приема информационного сообщения. Например если Г = 60 с, Тп = 6 с,fr* = 100 МГц/т«- = 1 МГц, то кээ = б. Предложенный алгоритм является реализацией защищенного патентом на изобретение способа работы маяка гидроакустической навигационной системы.

На рис. 6 представлена структура ГА ИИК, состоящая из ГА преобразователя, аналогового приемо-дередатчика, портативного компьютера или планшета, системы единого времени (GPS приемник и/или высокостабильный кварцевый генератор), системы тактового питания и FPGA с программным или аппаратным процессором. Особенность разработанного ИИК заключается в том, что структура является модульной и остается открытой для наращивания аппаратно-программного обеспечения и оптимизации методов обработки данных.

Коми>,шкашо1шый блок

ЕЕ

Программ отправки и сбора шфориашн

Программа управления

Программа обработки информации

Портативный компьютер

HDD SATA SD

Интерфейс USB

Интерфейс Elhcmct

Интерфейс RS-232

Интерфейс SATA

Интерфейс АЦП/ ЦАП

Классическая корреляция

Блок управлеит

Посимвольная коррелляшм

ИДДдДД Ñ,

Пороговый детектор

Многоканальный приемник/демодулятор

Передатчик/ модулятор

Синхронизацм

Устройство управления

НЕ

Одноканальный приемник

Управление тактовой частотой

Блок уменьшения энергопотребления

Систс ма такто вого питания

Сторожевой

Гндровкустмссю1Й преобразователь

А1алоговый приемо-передатчиь

Преду силитсяь с игната

Усилитесь мощности

Полосовой фильтр

АЦП

Рис. 6

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

В четвертой главе представлены результаты стендовых, лабораторных и натурных испытаний разработанного ГА ИИК.

Была проведена оценка погрешности формирования ИИ ИИК на базе кварцевого генератора с относительной погрешностью формирования частоты Ю-5 с помощью сертифицированного государственного вторичного эталона единиц времени и частоты ВЭТ 1-7, суммарная погрешность которого составляет не более 10"14. Оценивался разбег фаз Дг между сигналами ИИ от ВЭТ 1-7 и от внутреннего кварцевого генератора ИИК (с алгоритмом коррекции и без него). Число испытаний - 6 ООО, период ИИ - 30 с. Для нескорректированного сигнала ИИ максимальное значение разбега фаз ЛЬлшг = 115 мкс, = 99 мкс, с(А= 9 мкс. Для скорректированного сигнала ИИ

Мтах= 1,05 мкс, М(М) = 0,002 МКС, о(Д^ = 0,64 мкс.

С помощью ВЕТ 1-7 проводилась оценка погрешности формирования сигналов синхронизации платы точного времени КАС подводных роботов ИПМТ ДВО РАН. В основе платы КАС лежит высокоточный кварцевый генератор с относительной погрешностью формирования частоты 10'7. Полученные параметры: ДГмя = 2,5 мкс, М(Ы) = 2,2 мкс, а('Ы) = 0,4 мкс (число испытаний - 6 000).

Для оценки временной синхронизации двух ИИК была проанализирована погрешность формирования сигналов единого времени разнесенных вРв приемников, как разбег фаз Дг между ними. На рис. 7 представлен график распределения значения Дг для разного количества фиксируемых спутников 5аг на вРЭ приемниках (число испытаний - 19 000). При 7 и более фиксируемых спутниках

Мтах = 100 НС, М(М) = 1 НС, а(Дг) = 5 не. Проведенные экспериментальные исследования показали, что в реальных условиях эксплуатации погрешность временной синхронизации не превышала 100 не, что на порядок меньше погрешности формирования ИИ.

Натурные испытания ГА ИИК проводились в период 2007-2013 гг. в акваториях Японского моря (бухты Витязь, Рында). На базе разработанных ИИК были смонтированы стационарные и подвижные пункты приема и излучения, расположенные на пирсах, судне обеспечения и лодке. В экспериментах использовались ГА преобразователи с резонансными частотами от 1,7 до 41,7 кГц. Обработке подвергались информационные сообщения, состоящие из ФМШ длительностью от 31 до 1023 символов. Расстояние между пунктами приема и излучения изменялось от десятков метров до 20 км. За время проведения натурных экспериментов было проанализировано около 50 000 ФМШ.

Оценивалась инструментальная погрешность разработанного комплекса при обнаружения ГА сигналов как разница Дг в моментах прихода соседних ФМШ в информационном сообщении (рис. 8). Максимальный разброс от среднего значения Дг составил ± 1 //о. Благодаря разработанному алгоритму коррекции, ошибка формирования ИИ, составляющая не более 1 мкс, незначительно влияет на погрешность обнаружения момента прихода ГА сигналов.

м, не

Бм = 7.....Баг = 6— - Ба1 = 5

Бас = 4-«- Ба! = 3 Рис. 7

ФМШ зи

ГА сигнал

симв. бит Скорость передачи, бит/с ¿Ж* дБ ВЕЛ, отн. ед. Вшах, МКС

31 37 (5) 968 -5 -4 10 160

63 69(6) 571 -7 -4 10 80

127 276 (8) 378 -10 -6 10 40

255 282 (8) 189 -12 -6 10 40

511 1 082(10) 118 -14 - 20

1 023 3 159(11) 65 -17 - 20

В табл. 2 представлены обобщенные данные по 5Ж*, ВЕЯ и максимальной погрешности етах обнаружения момента прихода ГА сигналов при/я = 41,666 кГц и/в = 192 кГц. Расхождение значения минимального БИЯ* в натурном эксперименте от данных моделирования на 1 - 2 дБ позволяет судить об адекватности используемой компьютерной модели комплекса.

Были проведены испытания высокоскорстной передачи данных при помощи МФЧМ с 4 несущими частотами и 4 фазами и алгоритмом адаптации к параметрам ГА канала. Испытания показали, что при ЯУЯ > 6 дБ возможен прием данных с ВЕЯ 10"4 при длительности символа более 1,5 мс.

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

Параметр Значение

Разрядность АЦП, ЦАП до 18 разрядов

Тип выходного сигнала ШИМ, аналоговый

Предусилитель 10 дБ

Полоса пропускания Д/ от 1 кГц до 50 кГц

Максимальная /Ь 200 кГц

Период синхронизации от 1 до 120 с

Погрешность формирования ИИ не более 1 мкс

Инструментальная погрешность обнаружения момента прихода ГА сигнала не более 1 /УЬ

Системная тактовая частота от 1 до 200 МГц

Максимальное количество каналов обработки до 256 шт.

Энергопотребление в режиме сна / ожидания / приема не более 0,002 /0,1 / 1 Вт

Массогабаритные показатели 150 х 150 х 15 мм, 120 г.

В приложении приводятся документы о внедрении результатов работы, акт и программа испытаний блока формирования ИИ, патенты и свидетельства о регистрации программ для ЭВМ, фотография опытного образца ГА ИИК.

ОСНОВНЫЕ РЕЗУЛЬТАТЫ ДИССЕРТАЦИОННОЙ РАБОТЫ

1. Разработан алгоритм порогового детектирования ГА сигнала на основе вычисления величины отношения двух корреляционных функций по маске излучаемого и неизлучаемого ФМШ, обеспечивающий асинхронный прием и погрешность обнаружения момента прихода сигнала не превышающую длительность одного символа ФМШ при SNR меньше 0 дБ.

2. Разработан и исследован метод передачи ИИД по ГА каналу, обеспечивающий достоверный прием информации с BER не хуже Ю-6 при SNR меньше 0 дБ и снижение энергопотребления тракта излучения ПР до 5 раз за счет использования ФМШ и алгоритма сжатия данных.

3. Разработан алгоритм автоматической адаптации к характеристикам ГА канала для МФЧМ, используемой при высокоскоростной передачи данных, позволяющий уменьшить BER до 10 раз при SNR, равном 6 и более децибелам за счет добавления в пакет данных пилот-сигналов.

4. Разработан алгоритм синхронизации ИИК, обеспечивающий снижение погрешности формирования ИИ до 100 раз для кварцевых генераторов с относительной погрешностью формирования частоты 10"5 за счет компенсации нестабильности системной тактовой частоты.

5. Разработан алгоритм управления тактовой частотой параллельно-конвейерной структуры корреляторов в ИИК на основе технологии ПСнК типа FPGA, обеспечивающий сокращение энергопотребления комплекса до 10 раз.

6. Разработан универсальный унифицированный ИИК на базе технологии ПСнК типа FPGA, позволяющий в реальном масштабе времени производить многоканальную регистрацию ГА сигналов, повысить показатели точности обнаружения момента прихода сигналов, помехоустойчивости и достоверности передачи ИИД в водной среде, уменьшить энергетические и массогабаритные показатели.

7. Разработанный ГА ИИК используются в ИПМТ ДВО РАН при выполнении научно-исследовательских и опытно-конструкторский работ по разработке систем навигации, телеметрии и управления автономных необитаемых подводных аппаратов серии «ММТ».

ОСНОВНЫЕ ПУБЛИК АЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ

1. Линник, М.А. Гидроакустическая система связи для осуществления телекоммуникационного взаимодействия мевду подводными роботами / М.А. Линник, И.Н. Бурдинский, A.C. Миронов // Цифровая обработка сигналов. - 2010. - № 2. - С. 53-56.

2. Линник, М.А. Мобильный измерительный комплекс для проведения гидроакустических исследований / М.А. Линник, И.В. Карабанов, И.Н. Бурдинский, A.C. Миронов, C.B. Власьевский // Информатика и системы управления. - 2012. - № 2 (32). - С. 82-91.

3. Линник, М.А. Модель распространения сложного фазоманипулированного сигнала в гидроакустическом канале / М.А. Линник, И.В. Карабанов, И.Н. Бурдинский, A.C. Миронов // Известия ЮФУ. Технические науки. - 2013. - № 3 (140). - С. 103-112.

4. Линник, М.А. Методика определения порогового значения корреляционной функции для обнаружения сложного гидроакустического сигнала / М.А. Линник, И.В. Карабанов, И.Н. Бурдинский, A.C. Миронов, Л.Г. Вайнер // Информатика и системы управления. - 2014. - № 2 (40). - С. 45-55.

5. Маяк-ответчик гидроакустической навигационной системы : пат. 2426142 Рос. Федерация : МПК G01S 1/72 / Карабанов И.В., Бурдинский И.Н., Линник М.А., Миронов А.С., Безручко Ф.В., Чемерис Д.С. ; заявитель и патентообладатель Тихоокеанский государственный университет. -№ 2010104816/28 ; заявл. 11.02.2010 ; опубл. 10.08.2011, Бюл. № 22. -5с.: ил.

6. Способ передачи телеметрической информации для гидроакустической навигационной системы : пат. 2452976 Рос. Федерация : МПК G01S 15/00, G08C 23/02 / Карабанов И.В., Бурдинский И.Н., Линник М.А., Миронов А.С., Безручко Ф.В., Чемерис Д.С.; заявитель и патентообладатель Тихоокеанский государственный университет - К« 2010142029/28 ; заявл. 13.10.2010 ; опубл. 10.06.2012, Бюл. №6.-6 е.: ил.

7. Линник, М.А. Высокоточный модуль гидроакустической навигационной системы для определения наклонной дальности / М.А. Линник, И.В. Карабанов, И.Н. // Рефераты докладов XXVI конференции памяти выдающегося конструктора гироскопических приборов Н.Н. Острякова. - С.-Петербург: ЦНИИ «Электроприбор», 2008. - С. 57.

8. Linnik, М.А. Acoustic communication system for underwater vehicle telecommunication / M.A Linnik, I.V. Karabanov, I.N. Burdinskiy, A.S. Mironov // The First Russia and Pacific Conference on Computer Technology and Applications (Russia Pacific Computer 2010) 6 - 9 September, 2010 Russian Academy of Sciences, Far Eastern Branch. - Владивосток, 2010 - С 317-320.

9. Линник, M.A. Улучшение качественных показателей гидроакустической системы связи подводного робота / М.А. Линник, И.Н. Бурдинский // 1П Всероссийская научно-техническая конференция «Технические проблемы освоения мирового океана» : материалы конференции. - Владивосток : Дальнаука, 2009. - С. 366-371.

10. Линник, М.А. Проектирование и реализация цифровых модулей высокоскоростной гидроакустической системы связи И Материалы III конкурса-конференции научных работ молодых учёных Тихоокеанского государственного университета, Хабаровск, 1-4 дек. 2009 г. - Хабаровск : Изд-во Тихоокеан. гос. ун-та, 2010. - С. 29-34.

11. Линник, М.А. Беспроводная подводная система передачи телеметрической информации [Электронный ресурс] / М.А. Лчнник // Ученые заметки ТОГУ. - 2010. - Т. 1. - № 2, С. 14-20. - Режим доступа: http://ejournal.khstu.ru/

12. Линник, М.А. Модем гидроакустической системы связи / М.А. Линник, И.Н. Бурдинский, А.С. Миронов // 12-ая Международная конференция и выставка «Цифровая обработка сигналов и ее применение - DSPA-2010» : доклады. - М„ 2010. - Вып XII-2 -С.244-247.

13. Линник, М.А. Гидроакустическая система навигации и телеметрии подводного робота / М.А. Линник // Молодые ученые - Хабаровскому краю : в 2 т. : материалы XIII краевого конкурса молодых ученых и аспирантов, Хабаровск, 14-25 января 2011 г. - Хабаровск : Изд-во Тихоокеан. Гос. Ун-та, 2011. - Т. 2. - С. 189-193

14. Linnik, М.А. Estimation of different acoustic signals usage efficiency for navigation systems of underwater autonomous vehicles / M.A Linnik, I.V. Karabanov, I.N. Burdinskiy, A.S. Mironov // Proceedings of OCEANS 2011 IEEE Santander Conference and Exhibition - Santander, 2011.

15. Linnik, M.A. Underwater acoustic telecommunication system / M.A Linnik, I.V. Karabanov, I.N. Burdinskiy, A.S. Mironov // '.8th Saint Petersburg International Conference on Integrated Navigation System : proceedings. - Saint Petersburg, 2011. - P. 298-305

16. Линник, M.A. Информационно-измерительный комплекс для регистрации гидроакустических сигналов / М.А. Линник, И.В. Карабанов, И.Н. Бурдинский, А.С. Миронов, Ю.Г. Ларионов // Пятая Всероссийская научно-техническая конференция «Технические проблемы освоения Мирового океана» : материалы конференции. - Владивосток ■ Даль-наука, 2013. - С. 480-486.

Линник Михаил Александрович

ИНФОРМАЦИОННО-ИЗМЕРИТЕЛЬНЫЙ КОМПЛЕКС ДЛЯ РЕГИСТРАЦИИ ГИДРОАКУСТИЧЕСКИХ СИГНАЛОВ

Автореферат диссертации на соискание ученой степени кандидата технических наук

Подписано в печать 25.09.14. Формат 60x84 1/16. Бумага писчая. Гарнитура Times New Roman. Печать цифровая.

Усл. печ. л. 1,05. Тираж 120 экз. Заказ 273. Отдел оперативной полиграфин издательства Тихоокеанского государственного университета 680035, Хабаровск, ул. Тихоокеанская 136.