Метод построения динамической частотной характеристики лазерного гирометра со знакопеременной частотной подставкой типа меандр Текст научной статьи по специальности «Физика»

УДК 621.378.33 DOI: 10.18698/0236-3933-2016-4-129-141

МЕТОД ПОСТРОЕНИЯ ДИНАМИЧЕСКОЙ ЧАСТОТНОЙ ХАРАКТЕРИСТИКИ ЛАЗЕРНОГО ГИРОМЕТРА

СО ЗНАКОПЕРЕМЕННОЙ ЧАСТОТНОЙ ПОДСТАВКОЙ ТИПА МЕАНДР

В.Ф. Судаков [email protected]

МГТУ им. Н.Э. Баумана, Москва, Российская Федерация

Аннотация

Сложные (комбинированные) знакопеременные частотные подставки в лазерных гирометрах применяются для приближения формы их частотных характеристик к идеальной прямой. Известные аналитические методы расчета частотных характеристик в этих случаях практически непригодны, а непосредственное использование персональных компьютеров приводит к неприемлемо большим затратам времени. В целях преодоления этих трудностей разработан специальный метод расчета частотных характеристик при наличии знакопеременной частотной подставки. Предложено приращение фазы сигнала биений за период подставки рассматривать как угол поворота на фазовой плоскости за такое же время вектора-решения вспомогательного линейного векторного дифференциального уравнения. Тем самым вместо нелинейного фазового уравнения решена линейная га-мильтонова система дифференциальных уравнений второго порядка. Наиболее трудоемкая часть расчета сформулирована в аналитическом виде. Вычисления (рекуррентного типа) выполнены в среде Mathcad 15 за предельно малое время. Предлагаемый метод применен для расчета частотных характеристик лазерных гирометров при подставке меандр

Ключевые слова

Лазерный гирометр, фазовое уравнение, частотная подставка, частотная характеристика, подставка типа меандр

Поступила в редакцию 30.11.2015 © МГТУ им. Н.Э. Баумана, 2016

Введение. Частотная характеристика (ЧХ) — это зависимость частоты биений (приращение фазы сигнала биений за известный интервал времени, отнесенное к величине этого интервала) от разности собственных частот кольцевого резонатора (КР). Именно так в настоящей работе и будем рассматривать входную величину и выходное значение ЧХ. Переход к угловой скорости от разности частот можно считать отдельной задачей (об определении коэффициента линейной связи между ними).

Многочисленные работы по лазерным гирометрам (ЛГ), начиная с конца 1960-х г. (например, [1-8]), в той или иной степени способствовали созданию ком-

пьютерной (имитационной) модели ЛГ. Ценность упомянутых работ была и остается до настоящего времени в проведении количественного описания ЧХ, в частности, с использованием ПК. Исторически первой количественно была описана динамическая ЧХ при наличии простой периодической подставки. На этом этапе в основном применялись аналитические методы исследования. В конечном счете, ЧХ каждого образца ЛГ определяли экспериментально. Интерес к получению ЧХ на ПК [9-12] объясняется тем, что в машинном эксперименте можно получить ЧХ виртуального ЛГ, не имея при этом реального прибора.

Использование комбинированной подставки, когда на периоде основной (медленной) подставки действует несколько десятков периодов быстрого меандра, радикально изменило подход к построению динамической ЧХ. Стало принципиально невозможным описать ЧХ аналитически. В данном случае необходимо ориентироваться на создание метода, по которому ЧХ будет строить ПК. В то же время наиболее принципиальную часть расчета целесообразно выполнить аналитически, и только рутинную часть расчетов возложить на ПК. Также следует исходить из того, что допустимое время построения ЧХ должно быть достаточно мало, поскольку предполагается оперативность получения и использования результата. Именно так и ставится задача: необходимо разработать такой метод расчета ЧХ с комбинированной подставкой, который в среде Mattaad 15 может быть реализован в пределах 8...10 с на один вариант ЧХ. Работоспособность метода целесообразно сначала проверить путем построения динамической ЧХ при наличии простой знакопеременной подставки, например, в форме меандра. Такая подставка многократно использовалась ранее (например, [10-12]), что и позволяет считать ее эталоном для проверки предлагаемого алгоритма.

Вспомогательная система линейных дифференциальных уравнений. Рассмотрим ЛГ в режиме, когда интенсивности встречных волн практически равны и постоянны. В этом случае динамика фазы ^(t) сигнала биений двух встречно распространяющихся волн может быть описана, как известно [4], одним нелинейным дифференциальным уравнением

dy = кпр (t)+aL cos (у+х), (1)

dt

где K = 1--— коэффициент, учитывающий уровень накачки (®0d) и покоил/ п

луширину доплеровского контура усиления активной среды (k0wVn); Qp = &cw ~&ccw — расщепление собственных частот КР; полуширина статической зоны захвата (Ql ) и постоянная фаза х определяются амплитудами и фазами коэффициентов обратного рассеяния встречных волн.

Комбинированная подставка — это переменное частотное расщепление

вида

Qd (t) = Qr sign [sin(vr t)] + Qs sign [ sin(v5t)], (2)

т. е. представляет собой суперпозицию медленного меандра с амплитудой О) периода Т и быстрого меандра с амплитудой QR » периода ^ ^ Т). Кроме того, расщепление имеет также постоянную составляющую О. Таким образом,

K QP ^) = QD (t) + О.

Общих приемов аналитического решения фазового уравнения (1) не существует, однако, можно предложить специальный метод определения фазы в конце периода модуляции ) при начальной фазе ¥(0), что позволяет по-

ПГз) -¥(0)

строить частотную характеристику Qbeat = ■

Ts

(Qbeat — частота

биений по определению).

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

^ = At) У dt

(3)

с начальным условием у(0). Пусть матричный коэффициент имеет специальный вид

A(t) =

( 0 -y(t) ^ a(t) 0

(4)

Поскольку БрЛ = 0, то система дифференциальных уравнений (3) каноническая (гамильтонова). Следовательно, ее можно переписать в стандартном виде

с гамильтонианом H(t) = ( 0 1 ^

4 = н (t) У dt

(a(t) 0 ^ 0 y(t)

(5)

и симплектической матрицей

J =

-1 0

Введем в рассмотрение угол 9 (£), образованный вектором у^) = в соответствующей фазовой плоскости с осью уь Очевидно, что

Уг(1)

0(t) = arctg

yAt)

( y1(t)Л y2(t) ,

(6)

Угловая скорость вращения вектор-решения у (£) вспомогательного дифференциального уравнения (3) есть

dö d f y2 — = —\ arctg — |.

dt dt V yi

Вычислим эту производную как

d0 1 (dy2

dyi

\ y1 — y2

dt y12 + y22 v dt dt

(7)

Очевидно (см. (5)), что

dy f 0 1 ^

Hy__

f dyi ^ f dy2 >

v"1 0 У

dt dy2 dt

dt dy1 dt

Воспользуемся этим выражением, чтобы представить (7) в следующем виде:

d9

1

(Hy, y)

,4*. У1 - ^ У2 |=-2 — 2 1 dt П dt у J (y, y)

dt У12 + У22

Здесь использовано определение скалярного произведения векторов

(у,у) = У12 + У22, (Ну,у) = ^У1 -^У2.

Из матрицы Н () и уравнения (5) следует, что скалярное произведение (Ну, у) = ау12 + уу22. Тогда угловую скорость вектора-решения запишем так:

de dt

Поскольку ео8 0 =

y1

2 + y22

( Hy, y ) _ay12 +yy22 (y, y)

г, sinö _

y12 + y22

y2

2 + y22

(8)

то последней формуле

можно придать следующий вид:

de

— _ a cos2 ö + Y sin2 9, dt

или

d(2ö) dt

_ (a + y) + (a - y)cos (29).

(9)

Это уравнение совпадает с фазовым уравнением (1) теории ЛГ

dt

_Qd(t) + Q-Ql cos(Ч' + х + л),

/

если положить, что

а(г)- у (г) = -О, а(г) + у(г) = ОD(г) + О.

Отметим, что при этом фазе сигнала биений ¥ соответствует угол поворота 0 вектора-решения вспомогательного векторного дифференциального уравнения:

29=¥+х + л. (10)

Таким образом, вместо определения фазы сигнала биений из фазового уравнения с заданными параметрами О^, ОD(г) можно найти решение 20 уравнения (9) с заданными переменными коэффициентами

ов (г)+о-оь ов (г)+о+оь а(г) =-2-, у(г) =-2-. (11)

Фактически решать уравнение (9) не нужно, так как можно искать угловую скорость по формуле (8). Правая часть этого выражения зависит от решения линейной задачи (3), которая проще, чем эквивалентные нелинейные уравнения (9) или (1).

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

Общий вид вектора-решения вспомогательного дифференциального уравнения (3) в случае меандрподставки. Частное решение (3) можно записать в явном виде

у(г)=х(г) у(0).

Здесь использована матричная функция (матрициант) Х(г), удовлетворяющая матричному уравнению

^ = Л(г )х (12)

йг

(1 0 >

с начальным условием Х(0) = ^ ^ .

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

QD (t) = qs (t) =

Qs , 0 < t < —;

T 2 (13)

-Qs , — < t < Ts . 2

Постоянная О в (11) может изменяться только как независимая перемен-2л

ная. Постоянные Оз, Т = —, Оь считаются заданными параметрами. Учиты-

V)

вая (11) и (13), выявили, что матрица Л()) кусочно-постоянная и периодиче-^ 2л

ская с периодом Т = —. При таких матричных элементах возможно аналити-V)

ческое решение задачи, в частности — определение матрицианта. Чтобы упростить формулы, введем следующие обозначения при q = 1, 2:

ач = а ((-1)«-1 О) ), уя =у((-1)«-1 О) ), = ц( (-1)«-1 О) ) = 70^.

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

Хя = X ((-1)«-1 О) ), Ня = Н ((-1)«-1 О) ).

Рассмотрим полупериод

0, Z

2

, где матрица A(t) постоянна с элемента-

ми ах, ух. Собственные значения этой матрицы находим как корни уравнения

(

det

| -У1 ^ a1

= 0.

Отсюда получаем два корня Цх = ±^/ухах. Каждому из собственных значений соответствует собственный вектор. Два таких вектора-столбца образуют матрицу

T=

(i -i ^ a1 a1

1 1

С ее помощью можно привести матрицу Л к диагональному виду

(ь

A = T

i£1

0 ^

0 -¿|

T -

(14)

Формальное решение матричного дифференциального уравнения (12) с постоянной матрицей коэффициентов есть матричная экспонента Х()) = еЛ). По

аналогии с (14) найдем явный вид матрицианта на полупериоде

0, Ts

2

X1(t) = eAt = T ( exp (i£1t)

0 exp (-i£1t)

T-1 =

cos (|0,1t) - sin (|0,1t)

Mi

a1

sin (|M1t)

-1

cos(|Ц))

(15)

Матрициант (15) позволяет найти решение вспомогательного векторного дифференциального уравнения (3) в пределах указанного полупериода

у(") = Х^) у(0). (16)

В конечной точке рассматриваемого полупериода получим

y\f j_ -1 f fl y (»).

На следующем полупериоде

, 2 TL

2 2

решение вспомогательного век-

торного дифференциального уравнения (3) находим аналогично

Ts

Ts

y(t) _ X2 | t " f j y if J _ * f t " f j * f T

y1 (0)' y2 (0 )

На этом полупериоде коэффициенты а2, у2 и матрициант имеют другой вид

*2\ t "TL |_

С08

I , Ts

81П

М t

Ts

М2

a2

-1

- 81П

С08

i , Ts

i , Ts Ц2\ t-—

£2

a2

Вектор-решение на полупериоде

TL, 2 TL

2 2

можно представить так:

y«)_Л2И"T] yVT

В конечной точке t = Т получаем у ( Т ) = Х21 ^ у ^ |. Таким образом,

y(t) _

*1(t) y(0),

0 < t < — 2

*2\t - TL j y \ T) • T < t < 2 f-

(17)

Формула (17) является основной в последующем.

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

1 Т&

=—| ) dt. Реально эту зависимость часто заменяют более простой: Т

S 0

^ ( т ) — ¥(0)

Оьш =— —-• Мгновенная частота сигнала биений (подчеркнем, что это

Тз

не частота биений) есть производная фазы сигнала биений, т. е. определяется при нашем подходе выражением (8)

w = 2 (Hy, y) .

dt (y,y) '

(18)

По приведенному определению частота биений пропорциональна приращению фазы, т. е. в силу (18) имеем

^beat ~

V(Ts)-¥(0) = 2_T (H(t)y(t), y(t)) t = ^ Ts 0 (y(t), y(t))

Ts

= 2_Ts (H(t)X(t)y (0), X(t)y(0)) ^ Ts 0 (X(t)y(0), X(t)y(0)) '

(19)

начальным вектором /(0) =

т. е. начальным значением угла

Пусть ¥(0) = —л. Отсюда можно установить естественную связь между

( /i(0) W cos 0(0)Л J2(0))~{ sin 0(0) J'

поворота, и начальной фазой у. Она следует из равенства (11)

у

20(0) = ¥(0) + у + л: 0(0) = — = 5. Следовательно, вычислять интеграл необходимо при условии

y (0 ) =

cos

sin 5

cos^ 2

v

sin ^

2 j

Таким образом, частота биений оказывается зависящей от коэффициентов обратного рассеяния не только через О^ , но и через х.

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

Ts (H(t)X(t)y(0), X(t)y(0)) t = Ts¡2

0 (X(t)y(0), X(t)y(0)) 0

(H1X1 (c) y(0), Xi (c) y(0)) (Xi (C) y(0), X1 (C) y(0))

Ts } v Ts

H2X2 (C)y[-2j, X2 (C)y^-f X2 (C) y (f), X2 (C) y (f j)

dC

или (в другом представлении) к виду

| ( HX(t)y(0), Х(") у(0) ) ^ 2 ^ ( НдХд (С) у (Тд ) , Хд (^ у ( Тд ))

I (()у(0), Х(")у(0)) ¿1 I (Хд(С)у(Т),Хд (С)у(Тд)) ^ Т

где Тд = (д-1)у.

Таким образом, исходный интеграл распался на два частных интеграла (по

числу полупериодов протяженностью Т- на периоде Т$). Каждый из этих

частных интегралов может быть вычислен аналитически.

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

(((С) y (T ) (С) y (T ))_

c0s (Ц?С)У1 (Tq sin (n?C)y2 ( Tq )

aa

4'

-1

v aq J

sin (|q C)y1 ( Tq ) + c0s (^qC)y2 (Tq ) _ A0q + ACq cos (2|J,qC)+ ASq sin (2|qC),

где

A0q = 2

ff \-2 >

vaq J

f a q ^ 1+-1

f

y1 (Tq )2

f|q ^

vaq J

2

y2 (Tq )2

V

/

ACq _■

У q

f f Л-2 ^

1J lq ^

y1 (Tq )2 +1

f у >

1+iq

v aq J

va q J

y2 ( Tq ) ;

f fl У^

1—

aq

yq

y1 (Tq )2 y1 (Tq )2 +1

1 -

lq

vaq J

f 1 -Ii '

v aq J

y2 (Tq )2 y2(Tq)2;

ASq _ y1 (Tq )y2 (Tq )

f|q V

vaq J

aq

_ y1 ( Tq )y2 ( Tq )

aq У q lq

(HqXq (С) y (Tq ), Xq (С) y (Tq ))_y qy2 (Tq )2 +aqy1 (Tq )2 _ B0q

В этих обозначениях интеграл (20) примет следующий вид:

Ts

2 2 (HqXq (С)y (Tq) Xq (С)y (Tq)) df _

Й 0 (((0y (Tq) Xq (С) y (Tq))

_ ¿Г

q_1 I у0 q

arctg

■J2 q f TS ^ У1 q

"0 q

q

arctg

f у

у1 q

? J

(21)

где

Еея =. ;

V Уq aq

B0q У^ , \2 ая , \2

я = = у2 (Т) +—у (Т) ;

Мя Мя Мя

Е2я = ^у (Тя)2 У2 (Тя)2;

У я ая

Е1я = А5я =-1^у1 (Тя)у2 (Тя).

Мя

Формула (21) определяет приращение фазы сигнала биений на периоде подставки, т. е. частоту биений. Особенность вычисления в том, что предварительно надо найти начальные условия у (Т9) для каждого полупериода. Соответствующие вычисления запишем в виде:

у(Т1 ) = у(0), у(Т2 ) = у(^) = X (Тг)У(0).

Таким образом, на ПК не возлагается задача непосредственного интегрирования, требующая больших затрат времени. Предлагаемый метод позволяет организовать вычисления в виде цикла, что упрощает запись программы вычислений. Общее время расчета одного варианта даже в среде МаШсас! 15 не превышает 10 с.

С помощью предложенного в метода была построена ЧХ (рис. 1) и выявлено ее отличие от идеальной ЧХ (рис. 2) при = 2л • 1 • 103, = 2л • 4, = 2л • 800. На рисунках введено обозначение: 0.ьеап — массив, соответствующий функции при дискретных значениях аргумента 0„. Эти графики, продолженные в левую полуплоскость, симметричны относительно центра.

0 0,8 1,6 П„/у5

-0,75 -

-1,50 -41

Рис. 2. Отличие рассчитанной частотной характеристики от идеальной

abeatJvS

Рис. 1. Рассчитанная частотная характеристика

Заключение. Предложен и реализован с использованием ПК метод расчета динамической ЧХ в ЛГ с меандрподставкой. В теории кольцевых лазеров предложенные подходы ранее не использовались и аналогов не имеют. Физическое содержание действия подставки выявляется в самой структуре предложенного метода. Метод допускает различные обобщения. Его можно использовать в задачах, описанных, например, в [7-12]. Это позволяет при необходимости учесть различные особенности реальной подставки (отклонения от меандра, запаздывания и т. д.). Неизбежные упрощения при построении ЧХ были минимальными.

Доказана принципиальная возможность моделирования динамической ЧХ даже в такой простой вычислительной среде как MathCad 15. В то же время выявлены некоторые важные недостатки (применительно к решению поставленной задачи) данного пакета, например, отсутствие встроенной обратной функции Arctg (есть только ее главное значение аг^), что в данном конкретном случае приводит к дополнительным неточностям вычисления.

ЛИТЕРАТУРА

1. Климонтович Ю.Л., Курятов В.Н., Ланда П.С. О синхронизации волн в газовом лазере с кольцевым резонатором // ЖЭТФ. 1967. Т. 51. Вып. 1. С. 3-12.

2. Азарова В.В., Голяев Ю.Д., Савельев И.И. Зеемановские лазерные гироскопы // Квант. электрон. 2015. Т. 45. № 2. С. 171-179.

3. Судаков В.Ф. К теории кольцевого генератора с изменяющейся разностью частот резонатора // ЖПС. 1975. Т. 23, В. 5. С. 811-819.

4. Бирман А.Я., Наумов П.Б., Савушкин А.Ф., Тропкин Е.Н. Анализ динамической частотной характеристики кольцевого лазера на основе теории Флоке // Квант. электрон. 1986. T. 13. № 8. С. 1638-1644.

5. Хромых А.М. Динамическая характеристика кольцевых лазеров с периодической частотной подставкой // Электронная техника. Сер. 11. Лазерная техника и оптоэлектро-ника. 1990. Вып. 2 (54). С. 44.

6. Филатов Ю.В. Оптические гироскопы. М.: Изд-во ГНЦ РФ ЦНИИ «Электроприбор», 2005. 139 с.

7. Schreiber K.U., Gebauer A., Wells J.-P.R. Closed-loop locking of an optical frequency comb to a large ring laser // Optics Letters. 2013. 38 (18). Р. 3574-3577.

8. Shahriar S.M., Yablon J., Tseng S., Salit M. An inhomogeneously broadened superluminal ring laser for rotation sensing and accelerometry // Frontiers in Optics. 2011. Paper FWZ6.

9. Long Xingwu, Yuan Jie. Method for eliminating mismatching error in monolithic triaxial ring resonators // Chinese Optics Letters. 2010. 8 (12). Р. 1135-1138.

10. Fan Zhenfang, Luo Hui, Lu Guangfeng, Hu Shaomin. Dynamic lock-in compensation for mechanically dithered ring laser gyros // Chinese Optics Letters. 2012. 10 (6). Р. 061403.

11. Reanalysis of generalized sensitivity factors for optical-axis perturbation in nonplanar ring resonators / Yuan Jie, Chen Meixiong, Li Yingying, Tan Zhongqi, Wang Zhiguo // Optics Express. 2013. 21 (2). Р. 2297-2306.

12. Wen Dandan, Li Dong, Zhao Jianlin. Generalized sensitivity factors for optical-axis perturbation in nonplanar ring resonators // Optics Express. 2011. 19 (20). Р. 19752-19757.

Судаков Владимир Федорович — д-р техн. наук, профессор кафедры «Электротехника и промышленная электроника» МГТУ им. Н.Э. Баумана (Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5).

Просьба ссылаться на эту статью следующим образом:

Судаков В.Ф. Метод построения динамической частотной характеристики лазерного гирометра со знакопеременной частотной подставкой типа меандр // Вестник МГТУ им. Н.Э. Баумана. Сер. Приборостроение. 2016. № 4. С. 129-141. БО!: 10.18698/0236-3933-2016-4-129-141

METHOD FOR CONSTRUCTING THE DYNAMIC FREQUENCY CHARACTERISTICS OF LASER GYROMETER WITH ALTERNATING FREQUENCY MEANDER STAND

V.F. Sudakov [email protected]

Bauman Moscow State Technical University, Moscow, Russian Federation

Abstract

Complex (combined) stands with alternating frequency are used in laser gyrometers (LH) to approximate the shape of their frequency characteristics (CH) to an ideal straight line. The known analytical methods for calculating of CH cannot be applied in these cases. It takes unacceptably long to use personal computers (PC). This article offers a method to calculate the dynamic frequency characteristics (VX) for alternating frequency stands. In particular, the increment of the beat signal phase during the period of the stand can be regarded as the rotation angle from the phase plane for the same time. Vector is the solution of an auxiliary linear vector differential equation. Thus, instead of solving nonlinear phase equation the article offers to solve a linear Hamilto-nian system of differential equations of the second order. As a result, the most time-consuming part of the calculation is formulated in analytical form. The calculations of the recurrent type take little time and are executed in Mathcad 15. Meander type of stand has been studied repeatedly and its characteristics are well known. Therefore, it is a good option for a test of the proposed method. This article applies the method for the calculation of the dynamic frequency characteristics (VX) for alternating frequency meander stand

Keywords

Laser gyrometer, the phase equation, stand frequency, frequency response, stand type meander

REFERENCES

[1] Klimontovich Yu.L., Kuryatov V.N., Landa P.S. Wave synchronization in a gas laser with a ring resonator. J. Exp. Theor. Phys., 1967, vol. 24, no. 1, pp. 1-7.

[2] Azarova V.V., Golyaev Yu.D., Savel'ev I.I. Zeeman laser gyroscopes. Quantum Electronics, 2015, vol. 45, no. 2, pp. 171-179.

[3] Sudakov V.F. Theory of a ring laser with varying difference of resonator frequencies. Journal of Applied Spectroscopy (JAS), 1975, vol. 23, iss. 5, pp. 1454-1460.

DOI: 10.1007/BF00615854

[4] Birman A.Ya., Naumov P.B., Savushkin A.F., Tropkin E.N. Analysis of the dynamic frequency characteristic of a ring laser based on the Floquet theory. Soviet Journal of Quantum Electronics, 1986, vol. 13, no. 8, pp. 1069-1073.

[5] Khromykh A.M. The dynamic characteristics of ring lasers with a periodic frequency support. Elektron. Tekh. Ser. 11. Lazernaya Tekh. Optoelektron. [Electron. Eng. Laser Technol. Optoelectron.], 1990, vol. 53, no. 1, p. 44 (in Russ.).

[6] Filatov Yu.V. Opticheskie giroskopy [Optical gyroscopes]. Moscow, TsNII Elektropribor Publ., 2005. 139 p.

[7] Schreiber K.U., Gebauer A., Wells J.-P.R. Closed-loop locking of an optical frequency comb to a large ring laser. Optics Letters, 2013, 38 (18), pp. 3574-3577.

[8] Shahriar S.M., Yablon J., Tseng S., Salit M. An inhomogeneously broadened superluminal ring laser for rotation sensing and accelerometry. Frontiers in Optics, 2011, Paper FWZ6.

[9] Long Xingwu, Yuan Jie. Method for eliminating mismatching error in monolithic triaxial ring resonators. Chinese Optics Letters, 2010, 8 (12), pp. 1135-1138.

[10] Fan Zhenfang, Luo Hui, Lu Guangfeng, Hu Shaomin. Dynamic lock-in compensation for mechanically dithered ring laser gyros. Chinese Optics Letters, 2012, 10 (6), p. 061403.

[11] Jie Yuan, Meixiong Chen, Yingying Li, Zhongqi Tan, Zhiguo Wang. Reanalysis of generalized sensitivity factors for optical-axis perturbation in nonplanar ring resonators. Optics Express, 2013, 21 (2), pp. 2297-2306.

[12] Wen Dandan, Li Dong, Zhao Jianlin. Generalized sensitivity factors for optical-axis perturbation in nonplanar ring resonators. Optics Express, 2011, 19 (20), pp. 19752-19757.

Sudakov V.F. — Dr. Sci. (Eng.), Professor of Electrical Engineering and Industrial Electronics Department, Bauman Moscow State Technical University (2-ya Baumanskaya ul. 5, Moscow, 105005 Russian Federation).

Please cite this article in English as:

Sudakov V.F. Method for Constructing the Dynamic Frequency Characteristics of Laser Gyrometer with Alternating Frequency Meander Stand. Vestn. Mosk. Gos. Tekh. Univ. im. N.E. Baumana, Priborostr. [Herald of the Bauman Moscow State Tech. Univ., Instrum. Eng.], 2016, no. 4, pp. 129-141. DOI: 10.18698/0236-3933-2016-4-129-141