способ демодуляции сигнала волоконного интерферометра

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

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



где

a=u⁽⁰⁾(C₁ -C₂)+u⁽¹⁾(C₂-C₀)+u ⁽²⁾(C₀-C₁), b=u⁽⁰⁾(S ₁-S₂)+u⁽¹⁾(S₂-S₀ )+u⁽²⁾(S₀-S₁),

где u ⁽⁰⁾, u⁽¹⁾, u⁽²⁾ - значения отсчетов в тройке, С₀=cos[ ⁽⁰⁾], C₁=cos[ ⁽¹⁾], C₂=cos[ ⁽²⁾] и S₀=sin[ ⁽⁰⁾], S₁=sin[ ⁽¹⁾], S₂=sin[ ⁽²⁾],

где ⁽⁰⁾= _msin[ ₀], ⁽¹⁾= _msin[ ₁], ⁽²⁾= _msin[ ₂],

где _m - амплитуда модулирующего сигнала, ₀ - фазовая задержка первого отсчета относительно начала периода модуляции, ₁= ₀+2 /3, ₂= ₀+4 /3.

2. Способ по п.1 отличающийся тем, что формируют тройки отсчетов из всех трех отсчетов одного периода, с пропуском М-1 периодов, где М - целое число >1.

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



где

a=u⁽⁰⁾(C₁ -C₂)+u⁽¹⁾(C₂-C₀)+u ⁽²⁾(C₀-C₁), b=u⁽⁰⁾(S ₁-S₂)+u⁽¹⁾(S₂-S₀ )+u⁽²⁾(S₀-S₁),

где u ⁽⁰⁾, u⁽¹⁾, u⁽²⁾ - значения отсчетов в тройке, С₀=cos[ ⁽⁰⁾], C₁=cos[ ⁽¹⁾], C₂=cos[ ⁽²⁾] и S₀=sin[ ⁽⁰⁾], S₁=sin[ ⁽¹⁾], S₂=sin[ ⁽²⁾],

где для первой тройки - ⁽⁰⁾= _msin[ ₀], ⁽¹⁾= _msin[ ₁], ⁽²⁾= _msin[ ₂], для второй тройки ⁽⁰⁾= _msin[ ₁], ⁽¹⁾= _msin[ ₂], ⁽²⁾= _msin[ ₀], для третьей тройки ⁽⁰⁾= _msin[ ₂], ⁽¹⁾= _msin[ ₀], ⁽²⁾= _msin[ ₁],

где _m - амплитуда модулирующего сигнала, ₀ - фазовая задержка первого отсчета относительно начала периода модуляции, ₁= ₀+2 /3, ₂= ₀+4 /3, после чего процесс вычисления фаз повторяют.

Описание изобретения к патенту

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

Известно «Опросное устройство для волоконно-оптических линий с использованием двух склонов» [Патент US № 6778720]. Способ демодуляции сигнала, используемый в данном устройстве, включает вспомогательную модуляцию разности фаз световых волн в интерферометрах гармоническим сигналом, преобразование выходного сигнала фотоприемника, установленного на выходе линии, в последовательность цифровых отсчетов, выбор по пять отсчетов на каждом из двух склонов периодической зависимости сигнала от времени, вычисление по части выбранных отсчетов искомой разности фаз световых волн в интерферометре, вычисление амплитуды вспомогательной модуляции по остальным отсчетам и ее последующую регулировку. Недостатками данного способа являются сложность его реализации, использование в каждом цикле вычислений не менее десяти отсчетов сигнала, необходимость точного поддержания амплитуды вспомогательной модуляции и одинаковости разности плеч всех интерферометров, входящих в линию.

Известен «Четырехступенчатый дискретный способ демодуляции фазовых сдвигов для линий из волоконно-оптических датчиков» [Патент US № 6122057]. Способ включает вспомогательную фазовую модуляцию интерферометра гармоническим сигналом, формирование четырех сигналов, являющихся интегралами выходного сигнала фотоприемника в течение разных промежутков времени, преобразование интегральных сигналов в последовательности цифровых отсчетов, вычисление искомой разности фаз световых волн в интерферометре. Недостатками данного способа являются сложность его реализации, необходимость формирования четырех интегральных сигналов, дополнительного вычисления поправочного коэффициента для корректировки вычисленных значений разности фаз.

Известен «Способ и устройство демодуляции выходных сигналов интерферометра повышенной точности», выбранный за прототип [Патент US № 6556509]. Способ демодуляции включает модуляцию разности фаз световых волн в интерферометре по гармоническому закону с амплитудой, преобразование периодического сигнала, получаемого с фотоприемника, установленного на выходе интерферометра, (период сигнала равен периоду вспомогательной модуляции) в поток цифровых отсчетов с получением двенадцати отсчетов в течение каждого периода. Отсчеты регистрируются через равные промежутки времени, частота дискретизации в двенадцать раз выше частоты модуляции. Выбор в течение каждого периода сигнала шести отсчетов из двенадцати и вычисление по этим шести отсчетам значения искомой разности фаз один раз за период с использованием формул, полученных для амплитуды модуляции. После чего осуществляют вычисление по оставшимся шести отсчетам амплитуды модуляции, проверка равенства амплитуды значению радиан и корректировка амплитуды.

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

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

Для решения поставленной задачи предложено два варианта способа демодуляции сигнала волоконного интерферометра.

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

где

a=u⁽⁰⁾ (C₁-C₂)+u⁽¹⁾(C₂-C ₀)+u⁽²⁾(C₀-C₁), b=u ⁽⁰⁾(S₁-S₂)+u⁽¹⁾(S ₂-S₀)+u⁽²⁾(S₀-S₁ ),

где u⁽⁰⁾, u⁽¹⁾, u ⁽²⁾ - значения отсчетов в тройке, С₀=cos[ ⁽⁰⁾], C₁=cos[ ⁽¹⁾], C₂=cos[ ⁽²⁾] и S₀=sin[ ⁽⁰⁾], S₁=sin[ ⁽¹⁾], S₂=sin[ ⁽²⁾],

где ⁽⁰⁾= _msin[ ₀], ⁽¹⁾= _msin[ ₁], ⁽²⁾= _msin[ ₂],

где _m - амплитуда модулирующего сигнала, ₀ - фазовая задержка первого отсчета относительно начала периода модуляции, ₁= ₀+2 /3, ₂= ₀+4 /3. Тройки отсчетов могут быть сформированы из всех трех отсчетов одного периода с пропуском М-1 периодов, где М - целое число >1.

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

где

a=u⁽⁰⁾ (C₁-C₂)+u⁽¹⁾(C₂-C ₀)+u⁽²⁾(C₀-C₁), b=u ⁽⁰⁾(S₁-S₂)+u⁽¹⁾(S ₂-S₀)+u⁽²⁾(S₀-S₁ ),

где u⁽⁰⁾, u⁽¹⁾, u ⁽²⁾ - значения отсчетов в тройке, С₀=cos[ ⁽⁰⁾], C₁=cos[ ⁽¹⁾], C₂=cos[ ⁽²⁾] и S₀=sin[ ⁽⁰⁾], S₁=sin[ ⁽¹⁾], S₂=sin[ ⁽²⁾],

где для первой тройки - ⁽⁰⁾= _msin[ ₀], ⁽¹⁾= _msin[ ₁], ⁽²⁾= _msin[ ₂], для второй тройки ⁽⁰⁾= _msin[ ₁], ⁽¹⁾= _msin[ ₂], ⁽²⁾= _msin[ ₀], для третьей тройки ⁽⁰⁾= _msin[ ₂], ⁽¹⁾= _msin[ ₀], ⁽²⁾= _msin[ ₁],

где _m - амплитуда модулирующего сигнала, ₀ - фазовая задержка первого отсчета относительно начала периода модуляции, ₁= ₀+2 /3, ₂= ₀+4 /3, после чего процесс вычисления фаз повторяют.

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

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

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

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

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

Для модуляции разности фаз в интерферометре используем гармонический сигнал

где _m и f_м - амплитуда и частота модулирующего сигнала.

Известно, что напряжение, получаемое с фотоприемника, установленного на выходе двухплечного волоконного интерферометра, при наличии модуляции (t) имеет вид

где _С - искомая разность фаз, которую необходимо определить в результате демодуляции, U₀ и U_m - постоянная составляющая и амплитуда интерференционного сигнала.

Предполагается медленное изменение _С, U₀ и U_m по сравнению с колебаниями модулирующего сигнала, и в течение одного периода модуляции эти величины можно считать постоянными. Данное условие может быть формально выражено следующими неравенствами:

Преобразуем полученный сигнал (2) в последовательность отсчетов u_n с частотой дискретизации f_d , в три раза превышающей частоту модуляции f_d=3·f _M,

где n - номер отсчета от начала первого периода модуляции (n=0, 1, 2, ), ₀=2 f_М t - фазовая задержка первого отсчета относительно начала периода модуляции, t_n - момент времени n-го отсчета интерференционного сигнала, _n - значения модулирующего сигнала (t) в моменты времени t_n.

Для описания формируемых сигналов удобно ввести в рассмотрение значения аргумента модулирующего фазу гармонического колебания (1) в моменты регистрации отсчетов t_n

Для вычисления значения искомой разности фаз _C выбираем три отсчета интерференционного сигнала. Это условие определяется тем, что функция (2), описывающая модулированный интерференционный сигнал, имеет три неизвестных параметра, которые исходя из условия (3) полагаются постоянными в течение периода модуляции: искомая разность фаз _С, U₀ и U_m.

Обозначив тройку отсчетов сигнала в течение одного периода модуляции {u⁽⁰⁾, u⁽¹⁾, u⁽²⁾} и записав их значения по формуле (4), получим систему уравнений, решение которой приводит к выражениям для искомой разности фаз _C в диапазоне [- , ]:

где

a=u⁽⁰⁾ (C₁-C₂)+u⁽¹⁾(С₂-C ₀)+u⁽²⁾(C₀-C₁), b=u ⁽⁰⁾(S₁-S₂)+u⁽¹⁾(S ₂-S₀)+u⁽²⁾(S₀-S₁ ),

где C₀=cos[ ⁽⁰⁾], C₁=cos[ ⁽¹⁾], C₂=cos[ ⁽²⁾] и S₀=sin[ ⁽⁰⁾], S₁=sin[ ⁽¹⁾], S₂=sin[ ⁽²⁾],

где ⁽⁰⁾= _msin[ ₀], ⁽¹⁾= _msin[ ₁], ⁽²⁾= _msin[ ₂].

Тройки отсчетов могут быть сформированы из всех трех отсчетов одного периода, с пропуском М-1 периодов, где М - целое число >1.

Второй вариант способа демодуляции аналогичен первому, но тройки {u ⁽⁰⁾, u⁽¹⁾, u⁽²⁾} формируют из трех соседних отсчетов со сдвигом в каждой последующей тройке на один отсчет.

Формула (6), по которой проводят вычисления значений искомой разности фаз _C, не меняются, но при вычислениях используют три варианта коэффициентов, соответствующих трем разным наборам { ⁽⁰⁾, ⁽¹⁾, ⁽²⁾}:

для первой тройки { ⁽⁰⁾, ⁽¹⁾, ⁽²⁾}={ _msin[ ₀], _msin[ ₁], _msin[ ₂]},

для второй тройки { ⁽⁰⁾, ⁽¹⁾, ⁽²⁾}={ _msin[ ₁], _msin[ ₂], _msin[ ₀]},

для третьей тройки { ⁽⁰⁾, ⁽¹⁾, ⁽²⁾}={ _msin[ ₂], _msin[ ₀], _msin[ ₁]}.

Каждому из этих трех наборов { ⁽⁰⁾, ⁽¹⁾, ⁽²⁾} соответствует свой набор коэффициентов {С₀, С₁, С₂, S₀, S₁, S₂}, вычисляемый по тем же формулам: C₀=cos[ ⁽⁰⁾], C₁=cos[ ⁽¹⁾], С₂=cos[ ⁽²⁾] и S₀=sin[ ⁽⁰⁾], S₁=sin[ ⁽¹⁾], S₂=sin[ ⁽²⁾].

После третьей тройки наборы { ⁽⁰⁾, ⁽¹⁾, ⁽²⁾} и процесс вычисления разности фаз повторяют для 4-й, 5-й и 6-й троек и т.д.

Для реализации способа демодуляции необходимо регистрировать колебания _C(t), вызываемые воздействием на интерферометр в частотном диапазоне, не превышающем полосу 5 кГц и амплитудой не более 1 радиана (например, акустические или вибрационные воздействия).

По первому варианту способа модулируем разность фаз в интерферометре гармоническим сигналом с частотой f_M =50 кГц и амплитудой _m=2 радиана. При этом для максимальной частоты и амплитуды измеряемого сигнала F_C=5 кГц неравенства условия (3) принимают вид

f_М/F _С=10>>1 и (d _С/dt)=2 ·1·5=10 <<2 f_М=100 ,

т.е. условия выполняются.

Преобразование полученного периодического сигнала в последовательность цифровых отсчетов производим с частотой f_d=3f _М=150 кГц, т.е. в три раза больше, чем частота модуляции. При этом допустим, что сдвиг начального отсчета относительно начала периода соответствует фазовому сдвигу ₀=0.2 радиан (т.е. t=0.637 мкс).

Пусть в течение некоторого периода модуляции параметры интерференционного сигнала имеют значения U₀=1, U_m=1, а искомая фаза составляет _C=2 радиана. Тогда, согласно (4) для этого периода будет сформирована тройка отсчетов

{u⁽⁰⁾ , u⁽¹⁾, u⁽²⁾}={0.264, 0.063, 1.995}

Для вычисления искомой разности фаз используем фиксированные коэффициенты:

⁽⁰⁾= _msin[ ₀]=0.397, ⁽¹⁾= _msin[ ₁]=1.499, ⁽²⁾= _msin[ ₂]=-1.896 радиан;

C₀ =cos[ ⁽⁰⁾]=0.922, C₁=cos[ ⁽¹⁾]=0.072, C₂=cos[ ⁽²⁾]=-0.32;

S₀=sin[ ⁽⁰⁾]=0.387, S₁=sin[ ⁽¹⁾]=0.997, S₂=sin[ ⁽²⁾]=-0.948.

В результате расчета по формуле (6) для указанной тройки получим:

a=u ⁽⁰⁾(C₁-C₂)+u⁽¹⁾(С ₂-C₀)+u⁽²⁾(C₀-C₁ )=0.264 (0.072+0.32)+0.063 (-0.32-0.922)+

1.995 (0.922-0.072)=1.721;

b=u⁽⁰⁾(S₁ -S₂)+u⁽¹⁾(S₂-S₀)+u ⁽²⁾(S₀-S₁)=0.264 (0.997+0.948)+0.063 (-0.948-0.387)+

1.995 (0.387-0.997)=-0.788;

a/b=-2.185; arctg(a/b)=-1.42;

_С=arctg(a/b)+ =2 радиана.

Таким образом, получили правильное значение искомой разности фаз _C.

Данный пример иллюстрирует вычисление искомой разности фаз _C и в случае формирования троек из отсчетов каждого периода, и в случае формирования троек из всех отсчетов одного периода с пропуском М-1 периодов, где М - целое число >1, например при М=2 пропускаем каждый второй период.

Во втором варианте способа все действия выполняем аналогично первому, но формируем тройки отсчетов, объединяющие каждые три соседних отсчета со сдвигом на один отсчет. Например, пусть в течение некоторых двух периодов модуляции параметры интерференционного сигнала имеют значения U₀=1, U_m=1, а искомая фаза составляет _C=2 радиана. Тогда, согласно (4) для этих двух периодов будут зафиксированы следующие значения отсчетов

0.264; 0.063; 1.995; 0.264; 0.063; 1.995.

По этому набору отсчетов будут сформированы четыре тройки отсчетов:

первая {u⁽⁰⁾, u⁽¹⁾, u⁽²⁾ }={0.264, 0.063, 1.995}

вторая {u⁽⁰⁾ , u⁽¹⁾, u⁽²⁾}={0.063, 1.995, 0.264}

третья {u⁽⁰⁾, u⁽¹⁾, u⁽²⁾ }={1.995, 0.264, 0.063}

четвертая {u⁽⁰⁾ , u⁽¹⁾, u⁽²⁾}={0.264, 0.063, 1.995}

При вычислении искомой разности фаз используем три набора фиксированных коэффициентов. Для первой (4-й, 7-й и т.д.) тройки:

⁽⁰⁾= _msin[ ₀]=0.397, ⁽¹⁾= _msin[ ₁]=1.499, ⁽²⁾= _msin[ ₂]=-1.896 радиан;

C₀ =cos[ ⁽⁰⁾]=0.922, C₁=cos[ ⁽¹⁾]=0.072, C₂=cos[ ⁽²⁾]=-0.32;

S₀=sin[ ⁽⁰⁾]=0.387, S₁=sin[ ^(l)]=0.997, S₂=sin[ ⁽²⁾]=-0.948.

Для второй (5-й, 8-й и т.д.) тройки:

⁽⁰⁾= _msin[ ₁]=1.499, ⁽¹⁾= _msin[ ₂]=-1.896, ⁽²⁾= _msin[ ₀]=0.397 радиан;

C₀ =cos[ ⁽⁰⁾]=0.072, C₁=cos[ ⁽¹⁾]=-0.32, C₂=cos[ ⁽²⁾]=0.922;

S₀=sin[ ⁽⁰⁾]=0.997, S₁=sin[ ⁽¹⁾]=-0.948, S₂=sin[ ⁽²⁾]=0.387.

Для третьей (6-й, 9-й и т.д.) тройки:

⁽⁰⁾= _msin[ ₁]=1.499, ⁽¹⁾= _msin[ ₂]=0.397, ⁽²⁾= _msin[ ₀]=-1.896 радиан;

C₀ =cos[ ⁽⁰⁾]=-0.32, C₁=cos[ ⁽¹⁾]=0.922, С₂=cos[ ⁽²⁾]=0.072;

S₀=sin[ ⁽⁰⁾]=-0.948, S₁=sin[ ⁽¹⁾]=0.387, S₂=sin[ ⁽²⁾]=0.997.

Расчет коэффициентов а и b по формулам

a=u⁽⁰⁾(C₁ -C₂)+u⁽¹⁾(С₂-C₀)+u ⁽²⁾(C₀-C₁); b=u⁽⁰⁾(S ₁-S₂)+u⁽¹⁾(S₂-S₀ )+u⁽²⁾(S₀-S₁); для каждой тройки даст следующие результаты.

Для первой тройки:

а=0.264(0.072+0.32)+0.063(-0.32-0.922)+1.995(0.922-0.072)=1.721;

b=0.264(0.997+0.948)+0.063(-0.948-0.387)+1.995(0.387-0.997)=-0.788;

Для второй тройки:

a=0.063(-0.32-0.922)+1.995(0.922-0.072)+0.264(0.072+0.32)=1.721;

b=0.063(-0.948-0.387)+1.995(0.387-0.997)+0.264(0.997+0.948)=-0.788;

Для третьей тройки:

a=1.995(0.922-0.072)+0.264(0.072+0.32)+0.063(-0.32-0.922)=1.721;

b=1.995(0.387-0.997)+0.264(0.997+0.948)+0.063(-0.948-0.387)=-0.788;

Для четвертой тройки:

а=0.264(0.072+0.32)+0.063(-0.32-0.922)+1.995(0.922-0.072)=1.721;

b=0.264(0.997+0.948)+0.063(-0.948-0.387)+1.995(0.387-0.997)=-0.788;

Соответственно по всем тройкам согласно (6) получим одинаковые значения искомой разности _С=arctg(a/b)+ =2 радиана, т.е. правильное значение искомой фазы.

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