CC BY
60УДК 681.787
М.А. Земляникин, Ю.В. Раевская
ОБРАБОТКА СИГНАЛА, ПРИНИМАЕМОГО ВОЛОКОННО-ОПТИЧЕСКИМ
ГИДРОФОНОМ
Нижегородский государственный технический университет им. Р.Е. Алексеева
Рассматриваются две схемы демодуляции сигнала, принимаемого волоконно-оптическим гидрофоном -схема гомодинного метода демодуляции на основе перекрестного перемножения и схема демодуляции на основе вычисления значений функции арктангенса. Приводятся результаты расчетов амплитудных характеристик схем демодуляции, производится их сравнение с экспериментальными характеристиками, а также делается сравнение двух схем демодуляции между собой, даются результаты работы схем демодуляции.
Ключевые слова: волоконно-оптический гидрофон, амплитудная характеристика, гомодинный метод демодуляции.
Введение
Особенностью водной среды является малая дальность распространения электромагнитных волн в ней. Это затрудняет создание приборов, применяющих беспроводную связь, подобных тем, которые используются на суше. Поэтому большинство телеизмерений в водной среде проводятся с помощью гидроакустических систем. Так, например, задача организации судоходства и осуществления связи военных и гражданских судов успешно решается с помощью гидролокации, разведки морских нефтегазоносных бассейнов, включая входящие в арктический шельф, а также осуществляется компактными и протяженными гидроакустическими системами, применяемыми в системах наблюдения акваторий портов и охраны морских границ.
На протяжении всего прошлого столетия гидроакустические измерения проводились с помощью приборов, использующих пьезокерамические чувствительные элементы [1]. В пьезоэлектрическом гидрофоне используется прямой пьезоэффект некоторых веществ (сегнето-ва соль, ниобат калия), при котором механические напряжения, приложенные к кристаллу, вызывают появление поверхностных электрических зарядов на его гранях. При переменном механическом напряжении возникает переменная электродвижущая сила на электродах-обкладках. Вещество, обладающее пьезоэлектрическими свойствами, не обязательно является монокристаллом, поэтому нашли широкое применение поликристаллические материалы, называемые пьезокерамикой (например, титанат бария, титанат-цирконат свинца). Чувствительные элементы пьезоэлектрических гидрофонов изготавливают в виде пакетов прямоугольной или цилиндрической формы. Однако такие чувствительные элементы обладают некоторыми существенными недостатками - они имеют большой вес и объем, и их достаточно сложно мультиплексировать.
Поэтому за несколько последних десятилетий количество гидроакустических систем, построенных на основе волоконно-оптических интерферометрических датчиков, увеличилось в разы. Причиной этой перемены является ряд преимуществ волоконно-оптических датчиков над пьезоэлектрическими [2, 3]:
• высокая чувствительность,
• устойчивость к электромагнитным помехам,
• электрическая пассивность,
• малые габариты,
• легкость мультиплексирования,
• низкая стоимость изготовления чувствительных элементов,
• приемники могут быть удалены на расстояние 5-10 км от блока сбора данных,
• по одному оптическому волокну можно передавать сигналы многих датчиков.
© Земляникин М.А., Раевская Ю.В., 2018.
Несмотря на значительные успехи в области построения современных волоконно-оптических измерительных гидроакустических систем, все еще существует ряд проблем, которые не являются до конца решенными. Поэтому детального рассмотрения требуют вопросы проектирования волоконно-оптического гидроакустического датчика, а также его чувствительного элемента, организации сети волоконно-оптических гидрофонов, определения оптимальных способов демодуляции интерферометрических сигналов и алгоритмов последующей их обработки, которая обеспечивает достижение целей, поставленных перед разработчиком. В настоящей статье рассматриваются схемы демодуляции сигналов, принимаемых волоконно-оптическим гидрофоном и алгоритмы их последующей обработки.
Схема гомодинного метода демодуляции на основе перекрестного перемножения
Схема гомодинного метода демодуляции на основе перекрестного перемножения [4] представлена на рис. 1. Левая часть схемы представляет собой волоконно-оптический гидрофон, а правая - схему демодуляции сигналов.
Рис. 1. Схема гомодинного метода демодуляции на основе перекрестного перемножения:
Л - лазерный источник излучения; У - волоконный разветвитель; ОВ - оптическое волокно; ФМ - фазовый модулятор; ФП - фотоприемник; ОГ - опорный генератор; ФНЧ - фильтр низких частот; ФВЧ - фильтр высоких частот, используемый для исключения влияния механических воздействий, а также температуры; / - интегратор
Выходной сигнал интерферометра, входящего в состав волоконно-оптического гидрофона, представляется синусоидальной функцией, аргументом которой является разность фаз между плечами интерферометра 9(t). В разности фаз содержится необходимая информация об измеряемом давлении. Ток фотоприемника на выходе интерферометра представляется формулой
i(t) = AcosQ(t) + В,
где А и В - параметры, относящиеся к интенсивности сигнала и чувствительности детектора.
Во многих схемах демодуляции используется также высокочастотное колебание частоты ш0, поэтому ток фотоприемника имеет вид
i(t) = Acos[C cosu0t + q(t)]+В, (2)
где С cos №0t - несущий сигнал; С - индекс модуляции; ф(0 включает себя не только интересующий нас сигнал, но и смещение рабочей точки интерферометра.
Ток i имеет значительное количество гармонических составляющих, возникающих из-за фазовой модуляции. Амплитуда n-й гармоники связана с функцией Бесселя n-го порядка с аргументом фт, используя это разложение, получим формулу:
i(t) = А
/с(С) + 2 ^(-1)fc/2fc(C) cos 2fc«ot fc=i
cos ф(0 —
2 ^(—1)fcW(C) cos(2fc + 1)^ot
fc=0
(3)
sin9(t)[ + 5.
Например, амплитуда компоненты с частотой основного тона представляется формулой:
-2Л/1(С) соз ш^зтф(0, (4)
амплитуда второй гармоники:
—2Л/2(С) cos 2w0t cos ф( t).
(5)
Из этого выражения видно, что если ф(0 = 0, то только четные гармоники присутствуют в выходном сигнале (рис. 2), в то время как при ф(0 = п/2 останутся только нечетные гармоники (рис. 3).
Рис. 2. Спектр выходного сигнала интерферометра при ф(£) = 0, /о = 10 кГц
Рис. 3. Спектр выходного сигнала интерферометра при ф(£) = тс/2, /о = 10 кГц
Похожим образом ф(0 может быть разложен на сигнал с частотой w (циклическая частота измеряемого сигнала) и сдвиг из-за влияния окружающей среды
ф(0 = D cos wt + ^(t).
(6)
При разложении в ряд Фурье получаем:
cos ф(0 =
m
/о(Я) + 2 ^(—1)*/2*(Я) cos 2fcwt
sin^(t),
fc=1
cos^(t) —
2^(—1)fcW(£)cos(2fc + 1)wt fc=0
(7)
sin ф(0 =
2 ^(—1)fc/2fc+i(^) cos(2fc + 1)wt
sin^(t).
cos^(t) —
& = 0
со
/0ф)+ 2^(-1)*/2Л(Я)соз2£шГ зтф(0. (8)
Эти уравнения, аналогично (3), показывают, что в случае ф(0 = 0 в спектре присутствуют только четные гармоники частоты ш вокруг четных гармоник несущего сигнала с частотой ш0, а также нечетные гармоники частоты ш вокруг нечетных гармоник несущего сигнала (рис. 4).
со
m
m
При = п/2 в спектре присутствуют только нечетные гармоники частоты ш вокруг четных гармоник несущего сигнала с частотой ш0, а также нечетные гармоники частоты ш вокруг нечетных гармоник несущего сигнала (рис. 5).
Чтобы выделить ф(0, сначала выделим амплитуды первой и второй гармоник частоты ш0 из ¿(С) с помощью смесителя, перемножив ¿(С) с Gcosш0t и Hcos2ш0t, соответственно. Данные сигналы подают на фильтр нижних частот, на выходе которого получаем амплитуды основной и второй гармоник, соответственно:
-2Ав]г(С) smф(t),
2АН]2(С) ^ф(г)
(9) (10)
Рис. 4. Спектр сигнала при ф(£) = О,
f0 = 10 кГц, /=1кГц
Рис. 5. Спектр сигнала при ф(£) = ъ/2,
f0 = 10 кГц, /=1кГц
Затем дифференцируем обе компоненты по времени:
-2Ав]1(С) ^ф(г)ф(г),
(11)
2АН J2(С)smф(t)ф(t). (12)
Перемножив (10) и (11), а также (9) и (12), и вычтя второе произведение из первого,
получим:
(-2А)Н]2(С) ^ ф( г)(-2А)в]1(С) ^ ф( -
-(-2А)]1(С) sш фф 2АЩ2(С^т ф(г) ф&) = = 4А2 йН/1(С)Ь(С)(^ ф( г)2 + sш ф(02)ф(0 = = А2 СНМС)МС)ф(0.
(13)
Таким образом, полученный сигнал прямо пропорционален производной фазы ф(^. После интегратора сигнал будет иметь вид:
4А2 СНМС)МС)ф(0.
(14)
Из формулы (14) становится понятно, что выходной сигнал прямо пропорционален произведению /1(С}]2(С) и для обеспечения максимальной чувствительности волоконно-оптического интерферометра необходимо максимизировать данное произведение. Результаты моделирования, представленные на рис. 6, показывают, что максимальное значение, равное 0,2246, достигается при С = 2,37 рад, что согласуется с результатами, приведенными в литературе.
Рис. 6. Сравнение результатов расчетной и экспериментальной зависимостей амплитуды выходного сигнала от индекса фазовой модуляции С
Таким образом, получим искомую фазу ф(0, из которой можно найти значение давления.
Достоинством данного метода является независимость выходного сигнала от положения рабочей точки интерферометра, что объясняет его широкое использование.
Недостаток - зависимость выходного сигнала от мощности лазера, т.е. неизбежные флуктуации мощности лазерного источника излучения приводят к отклонению демодулиро-ванного сигнала от измеряемого.
Рис. 7. Результат работы схемы демодуляции на основе перекрестного перемножения
Моделирование данной схемы демодуляции было проведено с использованием средств пакета программ математического моделирования MATLAB. В качестве фильтра нижних частот использовался фильтр Баттерворта из-за отсутствия пульсаций и гладкости его амплитудно-частотной характеристики в полосе пропускания. Фильтрация производилась с помощью функции йШШ, которая исключает фазовый сдвиг, но в два раза увеличивает порядок фильтра. Для применения численного дифференцирования была реализована функция diff, вычисляющая производную в точке по четырем узлам. Численное интегрирование производилось с помощью встроенной функции cumtrapz, использующей метод трапеций.
Результат работы реализованного алгоритма приведен на рис. 7. Измеряемый сигнал равен сумме двух гармонических сигналов с частотами 1 кГц и 2 кГц, и амплитудами, равными 0,5 рад. Также используется нормировка на/2(С) и на В2. Из рисунка видно, что входной и демодулированный сигналы практически совпадают.
Схема демодуляции на основе вычисления значений функции арктангенса
Недостатком схемы гомодинного метода демодуляции на основе перекрестного перемножения, как было изложено ранее, является зависимость выходного сигнала от мощности лазера. Был разработан ряд схем демодуляции, в которых влияние данного фактора исключено. Одной из них является схема гомодинного метода демодуляции на основе вычисления значений функции арктангенса [5], которая приведена на рис. 8.
Рис. 8. Схема гомодинного метода демодуляции на основе вычисления функции арктангенса
Элементы данной схемы идентичны соответствующим элементам схемы на основе перекрестного перемножения, Х/У - блок деления, аг^ - блок вычисления значений функции арктангенса, фазовращ - фазовращатель, ФВЧ - фильтр высоких частот.
Как было показано, в схеме на основе перекрестного перемножения амплитуды сигналов после фильтров нижних частот имеют вид
-2AGJ1(Qsmф(t), (15)
-2АН]2(С)^ф(г). (16)
\ \
\ \ \ \ -\ \
\,
01234567
С. индекс фазовой модуляции, рад
/ (Г)
Рис. 9. Значения функции * в зависимости от индекса фазовой модуляции
Тогда при делении (15) на (16) выходной сигнал будет равен -1—tg9(t) при равных
/ 2(С)
амплитудах G и H. Чтобы найти искомую фазу, необходимо разделить на ^l(C) и вычислить
функцию арктангенса. Но в случае, если ^^ = 1, процедуру деления производить не требу-
ЖС.)
/i (О
ется. Построив зависимость 1 от индекса модуляции C (рис. 9), на графике можно заме-
тО
тить, что ^ =1 при С = 2,63 и С = 6,09.
При данных индексах модуляции результат после применения функции арктангенса будет
I (С)
без искажений. При индексе модуляции С = 4,48 = -1, так как arctg(-х) = -arctg(x),
и для получения правильного результата необходимо инвертировать знак полученного значения.
Данный алгоритм был реализован с помощью пакета прикладных программ MATLAB, реализация левой части схемы аналогична реализации предыдущей схемы. Значения функции арктангенса вычислялись с помощью встроенной функции atan. Результаты работы модели представлены на рис. 10.
Рис. 10. Результат работы схемы демодуляции на основе вычисления значений функции арктангенса при амплитуде сигнала меньше —
Измеряемый сигнал равен сумме двух гармонических сигналов с частотами 1 кГц и 2 кГц, и амплитудами равными 0,5 рад.
Как известно, функция агС^ х принимает значения в промежутке (— П; П), в то время
как значения измеряемой разности фаз могут превышать это значение и достигать десятков радиан. Поэтому без применения фазовращателя выходной сигнал при амплитуде гармонического входного сигнала А = 6 рад и частоте / = 1 кГц выглядел бы так, как показано на рис. 11.
Алгоритм работы данного фазовращателя заключается в следующем: вычисляются значения функции арктангенса в двух соседних точках. Если разность > П, то из всех после-
ПП
дующих значений вычитается - иначе, если разность < — —, то ко всем последующим значениям прибавляется П.
Рис. 11. Результат работы схемы демодуля- Рис. 12. Результат работы схемы демодуляции ции на основе вычисления значений функции на основе вычисления значений функции арктангенса при амплитуде сигнала больше - арктангенса при амплитуде сигнала больше -
без использования фазовращателя
при использования фазовращателя
Сравнение схем демодуляции на основе перекрестного перемножения и вычислении значений функции арктангенса
Амплитудные характеристики данных схем гомодинного метода демодуляции на основе перекрестного перемножения и вычисления значений функции арктангенса представлены на рис. 13.
X РОС эсм РЗС агс±д - идеальная хар-ка
1
ЯГ
Амплитуда входного сигналу рад
Рис. 13. Амплитудные характеристики схем демодуляции
х Рве осм теоретическая хар-ка
- идеальная хар-ка
/С****** X \ X
Xх \ х
X
V х
Амплитуда входного сигнала, рад
Рис. 14. Сравнение теоретической и экспериментальной амплитудной характеристик метода демодуляции на основе перекрестного перемножения
Из графика видно, что амплитудная характеристика первого метода значительно отличается от линейной. В литературе [6] показано, что амплитудная характеристика имеет вид:
N-1
к=0
(17)
где N равно количеству гармоник измеряемого сигнала, укладывающихся в полосе пропускания фильтра нижних частот.
Сравнение теоретической и экспериментальной [7] амплитудных характеристик для метода перекрестного перемножения представлено на рис. 14.
Качество работы реализованных демодуляторов можно оценить по табл. 1, в которой указана величина среднеквадратического отклонения в зависимости от вида модулирующего сигнала.
Таблица 1
Величина среднеквадратического отклонения (СКО) в зависимости от вида модулирующего сигнала
Модулирующее сообщение Величина СКО при использовании PGC DCM Величина СКО при использовании PGC arctg
Двухтональное (500 Гц, 1 кГц) 0,0004 0,0014
Двухтональное (500 Гц, 2 кГц) 0,0025 0,0016
Трехтональное (500 Гц, 1 кГц, 2 кГц) 0,0452 0,0023
Результаты в табл. 1 свидетельствуют о том, что СКО в первом методе резко увеличивается при усложнении сигнала.
Заключение
Была рассмотрена реализация схем гомодинного метода демодуляции на основе перекрестного перемножения и на основе вычисления значений функции арктангенса. Реализация производилась с помощью пакета прикладных программ MATLAB. Проверялась работоспособность данных схем на типовых сигналах. Полученные при моделировании амплитудные характеристики сравнивались с теоретическими характеристиками. Также оценивалось качество алгоритмов по величине среднеквадратического отклонения.
Из приведенного анализа становится ясно, что метод на основе вычисления значений функции арктангенса обладает большей точностью и его амплитудная характеристика линейна в большем диапазоне, поэтому при реализации схемы демодуляции в цифровом виде метод гомодинной демодуляции на основе вычисления значений функции арктангенса более предпочтителен.
Библиографический список
1. Рыбянец, А.Н. Современное состояние и перспективы развития пьезоэлектрической керамики за рубежом / А.Н. Рыбянец, В.П. Сахненко // Микросистемная техника. - 2002. - № 3.
- С. 16-22.
2. Fiber Optic Hydrophones. Stockbridge: 2011. - 5 p.
3. Eisenmenger, W. Fiber optic probe hydrophone / W. Eisenmenger, J. Staudenraus, C. Wurster // Proc. of the IEEE Ultrasonics Symposium. - 1994. - V. 2. - P. 941-944.
4. Dandridge, A. Homodyne demodulation scheme for fiber optic sensors using hase generated carrier / A. Dandridge , T.G. Gialloronzi, A.B. Tveten // IEEE Journal of Quantum Electronics. - 1982.
- V. 18. - № 10. - P. 1647-1653. doi: 10.1109/JQE.1982.1071416
5. Liao, Y. The arctangent approach of digital PGC demodulation for optic interferometric sensors / Y. Liao, X. Mao, L. Wang, M. Zhang // Proceedings of SPIE - The International Society for Optical Engineering. 2006. V. 6292. Art. 62921E. doi: 10.1117/12.678455
6. Плотников, М. Ю. Волоконно-оптический гидрофон : дисс. ... канд. техн. наук : 05.11.01 / Плотников, М. Ю. - СПб.: НИУ ИТМО, 2014. - 155 c.
7. Беликин, М.Н. Экспериментальное сравнение алгоритмов гомодинной демодуляции сигналов для фазового волоконно-оптического датчика / М.Н. Беликин [и др.] // Научно-технический вестник информационных технологий, механики и оптики. - 2015. - Т. 15. - № 6. - С. 1008-1014.
Дата поступления в редакцию 01.02.2018
M.A. Zemlyanikin, Yu.V. Raevskaya PROCESSING OF SIGNALS ACCEPTED BY FIBER-OPTICAL HYDROPHONE
Nizhny Novgorod state technical university. R.E. Alekseev
Purpose: Despite significant progress in building modern fiber-optic measuring sonar systems, there are still a number of problems that are not completely resolved. Therefore, detailed consideration requires the design of a fiber-optic hy-droacoustic sensor, as well as its sensitive element, the organization of a network of fiber-optic hydrophones, the determination of optimal methods for demodulating interferometric signals and algorithms for their subsequent processing, which ensures the achievement of the goals set for the developer. The purpose of this work is to investigate the operation of possible demodulation schemes for signals coming from an interferometer included in a fiber-optic hydrophone, as well as algorithms for the subsequent processing of these signals.
Design/methodology/approach: The operation of demodulation circuits was tested using the MATLAB software package.
Findings: The analysis showed that the homodyne modulation scheme based on calculating the values of the arctangent function is more accurate and its amplitude characteristic is linear over a larger range, therefore, when the demodulation scheme is implemented digitally, the homodyne demodulation method based on calculating the values of the arctangent function is more preferable.
Research limitations/implications: Amplitude characteristics were obtained for circuits for two signal demodulation schemes (based on cross multiplying and calculating the values of the arc tangent function), a comparison was made between the theoretical and experimental amplitude characteristics.
Originality/value: The amplitude characteristics of the demodulation circuits and algorithms for processing interferometer signals obtained as a result of the investigation can be used in the future to design a fiber-optic hydrophone and organize a network of fiber-optic hydrophones.
Key words: fiber-optic hydrophone, amplitude characteristic, homodyne demodulation method.
