RU2526585C2 - Способ определения угла ориентации стоячей волны в твердотельном волновом гироскопе - Google Patents

Abstract

Изобретение относится к области точного приборостроения и может быть использовано при создании твердотельных волновых гироскопов и систем ориентации и навигации на их основе. Способ заключается в том, что предварительно располагают гироскоп на платформе поворотного стола таким образом, чтобы его входная ось совпадала по направлению с осью вращения платформы и была направлена вертикально. Задают платформе последовательно два эталонных значения угловой скорости и в угловом диапазоне ориентации стоячей волны ±90° измеряют показания гироскопа, определяют разницу в ширине трубок изменений показаний гироскопа в заданном угловом диапазоне ее ориентации при разных угловых скоростях платформы. Подбирают корректирующие коэффициенты для вырабатываемых приборных значений синфазных и квадратурных составляющих синусного и косинусного каналов датчика угла гироскопа, обеспечивающих минимизацию этой разницы в ширине трубок изменения показаний гироскопа, а в рабочем режиме определяют угол ориентации стоячей волны относительно резонатора с помощью аналитического выражения, параметры которого скорректированы в результате предварительной операции. Изобретение обеспечивает повышение точности выработки угла ориентации стоячей волны твердотельного волнового гироскопа относительно его резонатора. 3 ил.

Description

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

В современном гироскопическом приборостроении известны способы определения стоячей волны относительно резонатора твердотельного волнового гироскопа, в котором измеряют сигналы ω₁=Acosωtcos2υ,

ω₂=Acosωtsin2υ с датчиков угла гироскопа,

где А - амплитуда колебаний резонатора,

υ - угол ориентации стоячей волны,

ω - круговая частота колебаний рабочей формы, t - время. Используя эти сигналы ,угол ориентации волны определяют из выражения 2υ=arctgω₂/ω₁ (см., например, книгу В.А.Матвеева, В.И.Липатникова, А.В.Алехина «Проектирование волнового твердотельного гироскопа». - М., 1998, стр.50).

В этом способе погрешности определения сигналов ω₁ и ω₂ обуславливают гармонические ошибки определения угла ориентации волны в зависимости от самого угла υ. Для компенсации этих погрешностей требуются трудоемкие калибровочные операции. При этом также необходимо применение алгоритмов определения четверти, в которой находится волна, обеспечивающих требуемую точность нахождения υ, когда ω₁ стремится к нулю.

За прототип взят способ определения ориентации стоячей волны в твердотельном волновом гироскопе, использующий приращение на к-том такте вычислений угла ориентации волны к вычисленному углу на предыдущем к-1 такте (см. отчет ИПМ РАН №1274 «Теоретические исследования, разработка методик, программно-методического обеспечения, регулировки и испытаний твердотельного волнового гироскопа». - М., 1998).

В этом способе для определения ориентации стоячей волны в твердотельном волновом гироскопе по «измерительным» каналам косинуса «С» и синуса «S» измеряют сигналы датчиков углов гироскопа, выделяют синфазные С_ф, S_ф и квадратурные С_кв, S_кв составляющие этих сигналов. Определяют тригонометрические функции синуса и косинуса угла ориентации стоячей волны, используя синфазные и квадратурные составляющие сигналов датчиков углов

$$\cos 4{\varphi }_0=\left(C_{ф}^2+{С}_{кв}^2+S_{ф}^2+S_{кв}^2\right)\cdot {(A\cdot B)}^{-1}$$

;

$$\sin 4{\varphi }_0=2\left(C_{ф}^{}{S}_{ф}^{}+C_{кв}{S}_{кв}^{}\right)\cdot {(A\cdot B)}^{-1},\text{ (1)}$$

где

$$A=\sqrt{C_{ф}^2+C_{кв}^2+S_{ф}^2+S_{кв}^2+2\cdot \left(C_{ф}{S}_{кв}-{С}_{кв}{S}_{ф}\right)}$$

;

$$B=\sqrt{C_{ф}^2+C_{кв}^2+S_{ф}^2+S_{кв}^2+2\cdot \left(C_{ф}{S}_{кв}-{С}_{кв}{S}_{ф}\right)}$$

.

Для расчета угла 4φ₀ запоминают его значение и значения его тригонометрических функций на предыдущем шаге съема показаний к-1 4φ_0,к-1, соз4φ_0,к-1, sin4φ_0,к-1.

Вычисляют новые значения тригонометрических функций на к-м шаге cos4φ₀, sin4φ_0,к.

Новое значение угла 4φ_0,к рассчитывают по формуле

$$4{\varphi }_{\mathrm{0,}к}=4{\varphi }_{\mathrm{0,}к-1}+\Delta {(4{\varphi }_0)}_{к},\text{ (2)}$$

где Δ(4φ₀)_к=sinΔ(4φ₀)_к=sin4φ_0,кcos4φ_0,к-1-cos4φ_0,кsin4φ_0,к-1.

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

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

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

$$\begin{array}{cc}{С}_{фп,к}={С}_{ф,к}+\delta {С}_{ф,к};& {С}_{фп,к-1}={С}_{ф,к-1}+\delta {С}_{ф,к-1};\\ {С}_{квп,к}={С}_{кв,к}+\delta {С}_{кв,к};& {С}_{квп,к-1}={С}_{кв,к-1}+\delta {С}_{кв,к-1};\\ S_{фп,к}=S_{ф,к}+\delta S_{ф,к};& S_{фп,к-1}=S_{ф,к-1}+\delta S_{ф,к-1};\\ S_{квп,к}=S_{кв,к}+\delta S_{кв,к};& S_{квп,к-1}=S_{кв,к-1}+\delta S_{кв,к-1}\end{array}\text{ (3)}$$

где

С_ф,к, S_ф,к, С_ф,к-1, S_ф,к-1 - истинные значения синфазных составляющих по каналам косинуса «С» и синуса «S» на к и к-1 тактах вычислений;

С_кв,к, S_кв,к, С_кв,к-1, S_кв,к-1 - истинные значения квадратурных составляющих по каналам «С» и «S» на к и к-1 тактах;

δС_ф,к, δS_ф,к, δС_ф,к-1, δS_ф,к-1 - погрешности выработки синфазных составляющих по каналам «С» и «S» на к и к-1 тактах;

δС_кв,к, δS_кв,к, δ_Скв,к-1, δS_кв,к-1 - погрешности выработки квадратурных составляющих по каналам «С» и «S» на к и к-1 тактах.

Вырабатываемые приборные значения углов на к-м и к-1 такте и приращение угла на к-м такте:

$$\begin{array}{l}4{\varphi }_{0п,к}=4{\varphi }_{\mathrm{0,}к}+\delta 4{\varphi }_{\mathrm{0,}к};\\ 4{\varphi }_{0п,к}=4{\varphi }_{\mathrm{0,}к}+\delta 4{\varphi }_{\mathrm{0,}к};\text{ (4)}\\ \Delta {(}_4=\Delta {(}_4+\delta \Delta {(}_4,\end{array}$$

где

4φ_0,к,4φ_0,к-1 - истинные значения углов на к и к-1 тактах;

Δ(4φ₀)_к - истинное значение приращения угла на к такте;

δ4φ_0,к; Δ4φ_0,к-1 - погрешности выработки углов на к и к-1 тактах;

δΔ(4φ₀)_к - погрешность выработки угла на к-м такте.

С учетом (4) погрешность определения угла ориентации волны на к-м такте можно представить

$$\delta 4{\varphi }_{\mathrm{0,}к}=\delta 4{\varphi }_{\mathrm{0,}к-1}+\delta \Delta {(4{\varphi }_0)}_{к}.\text{ (5)}$$

Рассмотрим погрешность определения угла, вызванную погрешностью приращения угла δΔ(4φ₀)_к.

При точном задании начального угла эта погрешность обуславливает погрешность определения угла ориентации волны.

В выражениях (3) истинные значения синфазных и квадратурных составляющих по измерительным каналам косинуса и синуса можно представить:

$$\begin{array}{cc}{С}_{ф,к}=a\cos 2{\varphi }_{0.к},& {С}_{ф,к-1}=a\cos 2{\varphi }_{0.к-1},\\ {С}_{кв,к}=b\sin 2{\varphi }_{\mathrm{0,}к},& {С}_{кв,к-1}=b\sin 2{\varphi }_{\mathrm{0,}к-1},\\ S_{ф,к}=a\cos 2{\varphi }_{0.к},& S_{ф,к-1}=a\cos 2{\varphi }_{0.к-1},\\ S_{кв,к}=b\sin 2{\varphi }_{\mathrm{0,}к},& S_{кв,к-1}=b\sin 2{\varphi }_{\mathrm{0,}к-1},\end{array}\text{ (6)}$$

где

a - большая полуось эллипсовидной траектории точки в плоскости CS, изображающей состояние резонатора, b - малая полуось эллипсовидной траектории точки в плоскости CS, изображающей состояние резонатора.

Используя выражения (1), (2), (3), (6), погрешность приращения угла можно представить в виде:

$$\delta \Delta {(4{\varphi }_0)}_{к}=n_{к}{m}_{к-1}-m_{к}{n}_{к-1}-\Delta {(4{\varphi }_0)}_{к},\text{ (7)}$$

где

$$\begin{array}{l}{m}_{к}={[{\left({а}_{к}\cos 2{\varphi }_{\mathrm{0,}к}+\delta {С}_{ф,к}\right)}^2+{\left(-b_{к}\sin 2{\varphi }_{\mathrm{0,}к}+\delta {С}_{кв,к}\right)}^2-\left({а}_{к}\sin 2{\varphi }_{\mathrm{0,}к}+\delta S_{ф.к}\right)}^2-\\ -\left(b_{к}\cos 2{\varphi }_{\mathrm{0,}к}+\delta S_{кв,к}\right){}^2]\cdot \{[{\left({а}_{к}\cos 2{\varphi }_{\mathrm{0,}к}+\delta {С}_{ф,к}\right)}^2+{\left(-b_{к}\sin 2{\varphi }_{\mathrm{0,}к}+\delta {С}_{кв,к}\right)}^2+\\ +{\left({а}_{к}\sin 2{\varphi }_{\mathrm{0,}к}+\delta S_{ф.к}\right)}^2+{\left(b_{к}\cos 2{\varphi }_{\mathrm{0,}к}+\delta S_{кв,к}\right)}^2+2({\left({а}_{к}\cos 2{\varphi }_{\mathrm{0,}к}+\delta {С}_{ф,к}\right)}^2\cdot \\ \cdot {\left(b_{к}\cos 2{\varphi }_{\mathrm{0,}к}+\delta S_{кв,к}\right)}^2-{\left(-b_{к}\sin 2{\varphi }_{\mathrm{0,}к}+\delta {С}_{кв,к}\right)}^2\left({а}_{к}\sin 2{\varphi }_{\mathrm{0,}к}+\delta S_{ф.к}\right))]\cdot \\ \cdot [{\left({а}_{к}\cos 2{\varphi }_{\mathrm{0,}к}+\delta {С}_{ф,к}\right)}^2+{\left(-b_{к}\sin 2{\varphi }_{\mathrm{0,}к}+\delta {С}_{кв,к}\right)}^2+{\left({а}_{к}\sin 2{\varphi }_{\mathrm{0,}к}+\delta S_{ф.к}\right)}^2+\\ +{\left(b_{к}\cos 2{\varphi }_{\mathrm{0,}к}+\delta S_{кв,к}\right)}^2-2(\left({а}_{к}\cos 2{\varphi }_{\mathrm{0,}к}+\delta {С}_{ф,к}\right)\cdot \left(b_{к}\cos 2{\varphi }_{\mathrm{0,}к}+\delta S_{кв,к}\right)-\\ \left(-b_{к}\sin 2{\varphi }_{\mathrm{0,}к}+\delta {С}_{кв,к}\right){\left({а}_{к}\cos 2{\varphi }_{\mathrm{0,}к}+\delta {С}_{ф,к}\right))]\}}^{-\frac{1}{2}};\end{array}$$

$$\begin{array}{l}{n}_{к}=2[\left({а}_{к}\cos 2{\varphi }_{\mathrm{0,}к}+\delta {С}_{ф,к}\right)\left({а}_{к}\sin 2{\varphi }_{\mathrm{0,}к}+\delta S_{ф.к}\right)+\left(-b_{к}\sin 2{\varphi }_{\mathrm{0,}к}+\delta {С}_{кв,к}\right)\cdot \\ \cdot \left(b_{к}\cos 2{\varphi }_{\mathrm{0,}к}+\delta S_{кв,к}\right)]\cdot \{[{\left({а}_{к}\cos 2{\varphi }_{\mathrm{0,}к}+\delta {С}_{ф,к}\right)}^2+{\left(-b_{к}\sin 2{\varphi }_{\mathrm{0,}к}+\delta {С}_{кв,к}\right)}^2+\\ +{\left({а}_{к}\sin 2{\varphi }_{\mathrm{0,}к}+\delta S_{ф.к}\right)}^2+{\left(b_{к}\cos 2{\varphi }_{\mathrm{0,}к}+\delta S_{кв,к}\right)}^2+2({\left({а}_{к}\cos 2{\varphi }_{\mathrm{0,}к}+\delta {С}_{ф,к}\right)}^2\cdot \\ \cdot {\left(b_{к}\cos 2{\varphi }_{\mathrm{0,}к}+\delta S_{кв,к}\right)}^2-{\left(-b_{к}\sin 2{\varphi }_{\mathrm{0,}к}+\delta {С}_{кв,к}\right)}^2\left({а}_{к}\sin 2{\varphi }_{\mathrm{0,}к}+\delta S_{ф.к}\right))]\cdot \\ \cdot [{\left({а}_{к}\cos 2{\varphi }_{\mathrm{0,}к}+\delta {С}_{ф,к}\right)}^2+{\left(-b_{к}\sin 2{\varphi }_{\mathrm{0,}к}+\delta {С}_{кв,к}\right)}^2+{\left({а}_{к}\sin 2{\varphi }_{\mathrm{0,}к}+\delta S_{ф.к}\right)}^2+\\ +{\left(b_{к}\cos 2{\varphi }_{\mathrm{0,}к}+\delta S_{кв,к}\right)}^2-2(\left({а}_{к}\cos 2{\varphi }_{\mathrm{0,}к}+\delta {С}_{ф,к}\right)\cdot \left(b_{к}\cos 2{\varphi }_{\mathrm{0,}к}+\delta S_{кв,к}\right)-\\ \left(-b_{к}\sin 2{\varphi }_{\mathrm{0,}к}+\delta {С}_{кв,к}\right){\left({а}_{к}\cos 2{\varphi }_{\mathrm{0,}к}+\delta {С}_{ф,к}\right))]\}}^{-\frac{1}{2}};\end{array}$$

$$\begin{array}{l}{m}_{к-1}={[{\left({а}_{к-1}\cos 2{\varphi }_{\mathrm{0,}к-1}+\delta {С}_{ф,к-1}\right)}^2+{\left(-b_{к-1}\sin 2{\varphi }_{\mathrm{0,}к-1}+\delta {С}_{кв,к-1}\right)}^2-\left({а}_{к-1}\sin 2{\varphi }_{\mathrm{0,}к-1}+\delta S_{ф.к-1}\right)}^2-\\ -\left(b_{к-1}\cos 2{\varphi }_{\mathrm{0,}к-1}+\delta S_{кв,к-1}\right){}^2]\cdot \{[{\left({а}_{к-1}\cos 2{\varphi }_{\mathrm{0,}к-1}+\delta {С}_{ф,к-1}\right)}^2+{\left(-b_{к-1}\sin 2{\varphi }_{\mathrm{0,}к-1}+\delta {С}_{кв,к-1}\right)}^2+\\ +{\left({а}_{к-1}\sin 2{\varphi }_{\mathrm{0,}к-1}+\delta S_{ф.к-1}\right)}^2+{\left(b_{к-1}\cos 2{\varphi }_{\mathrm{0,}к-1}+\delta S_{кв,к-1}\right)}^2+2({\left({а}_{к-1}\cos 2{\varphi }_{\mathrm{0,}к-1}+\delta {С}_{ф,к-1}\right)}^2\cdot \\ \cdot {\left(b_{к-1}\cos 2{\varphi }_{\mathrm{0,}к-1}+\delta S_{кв,к-1}\right)}^2-{\left(-b_{к-1}\sin 2{\varphi }_{\mathrm{0,}к-1}+\delta {С}_{кв,к-1}\right)}^2\left({а}_{к-1}\sin 2{\varphi }_{\mathrm{0,}к-1}+\delta S_{ф.к-1}\right))]\cdot \\ \cdot [{\left({а}_{к-1}\cos 2{\varphi }_{\mathrm{0,}к-1}+\delta {С}_{ф,к-1}\right)}^2+{\left(-b_{к-1}\sin 2{\varphi }_{\mathrm{0,}к-1}+\delta {С}_{кв,к-1}\right)}^2+{\left({а}_{к-1}\sin 2{\varphi }_{\mathrm{0,}к-1}+\delta S_{ф.к-1}\right)}^2+\\ +{\left(b_{к-1}\cos 2{\varphi }_{\mathrm{0,}к-1}+\delta S_{кв,к-1}\right)}^2-2(\left({а}_{к-1}\cos 2{\varphi }_{\mathrm{0,}к-1}+\delta {С}_{ф,к-1}\right)\cdot \left(b_{к-1}\cos 2{\varphi }_{\mathrm{0,}к-1}+\delta S_{кв,к-1}\right)-\\ \left(-b_{к-1}\sin 2{\varphi }_{\mathrm{0,}к-1}+\delta {С}_{кв,к-1}\right){\left({а}_{к-1}\cos 2{\varphi }_{\mathrm{0,}к-1}+\delta {С}_{ф,к-1}\right))]\}}^{-\frac{1}{2}};\end{array}$$

$$\begin{array}{l}{n}_{к-1}=2[\left({а}_{к-1}\cos 2{\varphi }_{\mathrm{0,}к-1}+\delta {С}_{ф,к-1}\right)\left({а}_{к-1}\sin 2{\varphi }_{\mathrm{0,}к-1}+\delta S_{ф.к-1}\right)+\left(-b_{к-1}\sin 2{\varphi }_{\mathrm{0,}к-1}+\delta {С}_{кв,к-1}\right)\cdot \\ \cdot \left(b_{к-1}\cos 2{\varphi }_{\mathrm{0,}к-1}+\delta S_{кв,к-1}\right)]\cdot \{[{\left({а}_{к-1}\cos 2{\varphi }_{\mathrm{0,}к-1}+\delta {С}_{ф,к-1}\right)}^2+{\left(-b_{к-1}\sin 2{\varphi }_{\mathrm{0,}к-1}+\delta {С}_{кв,к-1}\right)}^2+\\ +{\left({а}_{к-1}\sin 2{\varphi }_{\mathrm{0,}к-1}+\delta S_{ф.к-1}\right)}^2+{\left(b_{к-1}\cos 2{\varphi }_{\mathrm{0,}к-1}+\delta S_{кв,к-1}\right)}^2+2({\left({а}_{к-1}\cos 2{\varphi }_{\mathrm{0,}к-1}+\delta {С}_{ф,к-1}\right)}^2\cdot \\ \cdot {\left(b_{к-1}\cos 2{\varphi }_{\mathrm{0,}к-1}+\delta S_{кв,к-1}\right)}^2-{\left(-b_{к-1}\sin 2{\varphi }_{\mathrm{0,}к-1}+\delta {С}_{кв,к-1}\right)}^2\left({а}_{к-1}\sin 2{\varphi }_{\mathrm{0,}к-1}+\delta S_{ф.к-1}\right))]\cdot \\ \cdot [{\left({а}_{к-1}\cos 2{\varphi }_{\mathrm{0,}к-1}+\delta {С}_{ф,к-1}\right)}^2+{\left(-b_{к-1}\sin 2{\varphi }_{\mathrm{0,}к-1}+\delta {С}_{кв,к-1}\right)}^2+{\left({а}_{к-1}\sin 2{\varphi }_{\mathrm{0,}к-1}+\delta S_{ф.к-1}\right)}^2+\\ +{\left(b_{к-1}\cos 2{\varphi }_{\mathrm{0,}к-1}+\delta S_{кв,к-1}\right)}^2-2(\left({а}_{к-1}\cos 2{\varphi }_{\mathrm{0,}к-1}+\delta {С}_{ф,к-1}\right)\cdot \left(b_{к-1}\cos 2{\varphi }_{\mathrm{0,}к-1}+\delta S_{кв,к-1}\right)-\\ \left(-b_{к}\sin 2{\varphi }_{\mathrm{0,}к}+\delta {С}_{кв,к}\right){\left({а}_{к}\cos 2{\varphi }_{\mathrm{0,}к}+\delta {С}_{ф,к}\right))]\}}^{-\frac{1}{2}};\end{array}$$

Δ(4φ₀)_к=sin4φ_0,кcos4φ_0,к-1-cos4φ_0,кsin4φ_0,к-1.

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

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

1 - скорость прецессии волны 12,4°/ч;

2 - скорость прецессии волны 24,8°/ч;

3 - скорость прецессии волны 37,2°/ч.

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

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

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

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

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

$${\varphi }_{\mathrm{0,}к}=\mathrm{0,25}(n_{к}{m}_{к-1}-m_{к}{n}_{к-1})+{\varphi }_{\mathrm{0,}к-1},\text{ (8)}$$

где

$$\begin{array}{l}{m}_{к}=[{\left({С}_{фп,к}-\Delta {С}_{ф}\right)}^2+{\left({С}_{квп,к}-\Delta {С}_{кв}\right)}^2-{\left(S_{фп,к}-\Delta S_{ф}\right)}^{{}_2}-\left(S_{квп,к}-\Delta S_{кв}\right){}^2]\cdot \\ \cdot \{[{\left({С}_{фп,к}-\Delta {С}_{ф}\right)}^2+{\left({С}_{квп,к}-\Delta {С}_{кв}\right)}^2+{\left(S_{фп,к}-\Delta S_{ф}\right)}^{{}_2}+{\left(S_{квп,к}-\Delta S_{кв}\right)}^2+\\ +2(\left({С}_{фп,к}-\Delta {С}_{ф}\right)\left(S_{квп,к}-\Delta S_{кв}\right)-\left({С}_{кв,к}-\Delta {С}_{кв}\right)\left(S_{фп,к}-\Delta S_{ф}\right))]\cdot \\ \cdot [{\left({С}_{фп,к}-\Delta {С}_{ф}\right)}^2+{\left({С}_{квп,к}-\Delta {С}_{кв}\right)}^2+{\left(S_{фп,к}-\Delta S_{ф}\right)}^{{}_2}+{\left(S_{квп,к}-\Delta S_{кв}\right)}^2-\\ -2{(\left({С}_{фп,к}-\Delta {С}_{ф}\right)\left(S_{квп,к}-\Delta S_{кв}\right)-\left({С}_{квп,к}-\Delta {С}_{кв}\right)\left(S_{фп,к}-\Delta S_{ф}\right))]\}}^{-\frac{1}{2}};\end{array}$$

$$\begin{array}{l}{n}_{к}=2[\left({С}_{фп,к}-\Delta {С}_{ф}\right)\left(S_{фп,к}-\Delta S_{ф}\right)+\left({С}_{квп,к}-\Delta {С}_{кв}\right)\left(S_{квп,к}-\Delta S_{кв}\right)]\cdot \\ \cdot \{[{\left({С}_{фп,к}-\Delta {С}_{ф}\right)}^2+{\left({С}_{квп,к}-\Delta {С}_{кв}\right)}^2+{\left(S_{фп,к}-\Delta S_{ф}\right)}^{{}_2}+{\left(S_{квп,к}-\Delta S_{кв}\right)}^2+\\ +2(\left({С}_{фп,к}-\Delta {С}_{ф}\right)\left(S_{квп,к}-\Delta S_{кв}\right)-\left({С}_{квп,к}-\Delta {С}_{кв}\right)\left(S_{фп,к}-\Delta S_{ф}\right))]\cdot \\ \cdot [{\left({С}_{фп,к}-\Delta {С}_{ф}\right)}^2+{\left({С}_{квп,к}-\Delta {С}_{кв}\right)}^2+{\left(S_{фп,к}-\Delta S_{ф}\right)}^{{}_2}+{\left(S_{квп,к}-\Delta S_{кв}\right)}^2-\\ -2(\left({С}_{фп,к}-\Delta {С}_{ф}\right)\left(S_{квп,к}-\Delta S_{кв}\right)-\left({С}_{квп,к}-\Delta {С}_{кв}\right){\left(S_{фп,к}-\Delta S_{ф}\right))]\}}^{-\frac{1}{2}};\end{array}$$

$$\begin{array}{l}{m}_{к-1}=[{\left({С}_{фп,к-1}-\Delta {С}_{ф}\right)}^2+{\left({С}_{квп,к-1}-\Delta {С}_{кв}\right)}^2-{\left(S_{фп,к-1}-\Delta S_{ф}\right)}^{{}_2}-\left(S_{квп,к-1}-\Delta S_{кв}\right){}^2]\cdot \\ \cdot \{[{\left({С}_{фп,к-1}-\Delta {С}_{ф}\right)}^2+{\left({С}_{квп,к-1}-\Delta {С}_{кв}\right)}^2+{\left(S_{фп,к-1}-\Delta S_{ф}\right)}^{{}_2}+{\left(S_{квп,к-1}-\Delta S_{кв}\right)}^2+\\ +2(\left({С}_{фп,к-1}-\Delta {С}_{ф}\right)\left(S_{квп,к-1}-\Delta S_{кв}\right)-\left({С}_{квп,к-1}-\Delta {С}_{кв}\right)\left(S_{фп,к-1}-\Delta S_{ф}\right))]\cdot \\ \cdot [{\left({С}_{фп,к-1}-\Delta {С}_{ф}\right)}^2+{\left({С}_{квп,к-1}-\Delta {С}_{кв}\right)}^2+{\left(S_{фп,к-1}-\Delta S_{ф}\right)}^{{}_2}+{\left(S_{квп,к-1}-\Delta S_{кв}\right)}^2-\\ -2(\left({С}_{фп,к-1}-\Delta {С}_{ф}\right)\left(S_{квп,к-1}-\Delta S_{кв}\right)-\left({С}_{квп,к-1}-\Delta {С}_{кв}\right){\left(S_{фп,к-1}-\Delta S_{ф}\right))]\}}^{-\frac{1}{2}};\end{array}$$

$$\begin{array}{l}{n}_{к-1}=2[\left({С}_{фп,к-1}-\Delta {С}_{ф}\right)\left(S_{фп,к-1}-\Delta S_{ф}\right)+\left({С}_{квп,к-1}-\Delta {С}_{кв}\right)\left(S_{квп,к-1}-\Delta S_{кв}\right)]\cdot \\ \cdot \{[{\left({С}_{фп,к-1}-\Delta {С}_{ф}\right)}^2+{\left({С}_{квп,к-1}-\Delta {С}_{кв}\right)}^2+{\left(S_{фп,к-1}-\Delta S_{ф}\right)}^{{}_2}+{\left(S_{квп,к-1}-\Delta S_{кв}\right)}^2+\\ +2(\left({С}_{фп,к-1}-\Delta {С}_{ф}\right)\left(S_{квп,к-1}-\Delta S_{кв}\right)-\left({С}_{квп,к-1}-\Delta {С}_{кв}\right)\left(S_{фп,к-1}-\Delta S_{ф}\right))]\cdot \\ \cdot [{\left({С}_{фп,к-1}-\Delta {С}_{ф}\right)}^2+{\left({С}_{квп,к-1}-\Delta {С}_{кв}\right)}^2+{\left(S_{фп,к-1}-\Delta S_{ф}\right)}^{{}_2}+{\left(S_{квп,к-1}-\Delta S_{кв}\right)}^2-\\ -2(\left({С}_{фп,к-1}-\Delta {С}_{ф}\right)\left(S_{квп,к-1}-\Delta S_{кв}\right)-\left({С}_{квп,к-1}-\Delta {С}_{кв}\right){\left(S_{фп,к-1}-\Delta S_{ф}\right))]\}}^{-\frac{1}{2}};\end{array}$$

где

С_фп,к - приборная синфазная составляющая сигнала датчика угла гироскопа по косинусному каналу на к-м такте вычислений;

С_квп,к - приборная квадратурная составляющая сигнала датчика угла по косинусному каналу на к-м такте;

S_фп,к - приборная синфазная составляющая сигнала датчика угла по синусному каналу на к-м такте;

S_квп,к - приборная квадратурная составляющая сигнала датчика угла по синусному каналу на к-м такте;

С_фп,к-1 - приборная синфазная составляющая сигнала датчика угла по косинусному каналу на к-1-ом такте вычислений;

С_квп,к-1 - приборная квадратурная составляющая сигнала датчика угла по косинусному каналу на к-1-ом такте вычислений;

S_фп,к-1 - приборная синфазная составляющая сигнала датчика угла по синусному каналу на к-1-ом такте вычислений;

S_квп,к-1 - приборная квадратурная составляющая сигнала датчика угла по синусному каналу на к-1-ом такте вычислений;

ΔС_ф - коэффициент, корректирующий постоянную ошибку определения

синфазной составляющей сигнала датчика угла по косинусному каналу;

ΔC_кв - коэффициент, корректирующий постоянную ошибку определения квадратурной составляющей сигнала датчика угла по косинусному каналу;

ΔS_ф- коэффициент, корректирующий постоянную ошибку определения

синфазной составляющей сигнала датчика угла по синусному каналу;

ΔS_кв - коэффициент, корректирующий постоянную ошибку определения квадратурной составляющей сигнала датчика угла по синусному каналу;

φ_0,к-1 - угол ориентации волны относительно резонатора, определенный на к-1-ом такте.

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

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

$$\begin{array}{l}{С}_{фп,к}-\Delta {С}_{ф,к}+\delta {С}_{ф,к}-\Delta {С}_{ф,к};\\ C_{квп,к}-\Delta {С}_{кв,к}={С}_{кв,к}+\delta {С}_{кв,к}-\Delta {С}_{кв,к};\\ S_{фп,к}-\Delta S_{ф,к}+\delta S_{ф,к}-\Delta S_{ф,к};\\ S_{квп,к}-\Delta S_{кв,к}=S_{кв,к}+\delta S_{кв,к}-\Delta S_{кв,к};\text{ (9)}\\ {С}_{фп,к-1}-\Delta {С}_{ф,к-1}+\delta {С}_{ф,к-1}-\Delta {С}_{ф,к-1};\\ C_{квп,к-1}-\Delta {С}_{кв,к-1}={С}_{кв,к-1}+\delta {С}_{кв,к-1}-\Delta {С}_{кв,к-1};\\ S_{фп,к-1}-\Delta S_{ф,к-1}+\delta S_{ф,к-1}-\Delta S_{ф,к-1};\\ S_{квп,к-1}-\Delta S_{кв,к-1}=S_{кв,к-1}+\delta S_{кв,к-1}-\Delta S_{кв,к-1}.\end{array}$$

Корректирующие коэффициенты подбирают экспериментально. Для этого предварительно располагают гироскоп на платформе поворотного стенда, задают платформе первое значение эталонной угловой скорости. На заданном угловом диапазоне ориентации стоячей волны относительно резонатора измеряют показания гироскопа, которые являются функцией угла ориентации φ₀.

$${\omega }_1({\varphi }_0+{\phi }_3,{\omega }_{п})={\omega }_{до}({\varphi }_0+{\varphi }_1)+{\omega }_{зв}+{\omega }_{п1}+{\omega }_{д1}({\varphi }_0+{\varphi }_2,{\omega }_{п}),\text{ (10)}$$

где

ω_до(φ₀+φ₁) - систематическая составляющая скорости дрейфа волны, не зависящая от угловой скорости платформы и зависящая от угла ориентации волны φ₀ с начальной фазой φ₁;

ω_зв - вертикальная составляющая угловой скорости вращения Земли;

ω_п1 - первое значение эталонной угловой скорости, задаваемой платформе;

ω_д1(φ₁+φ₂,φ_п) - систематическая составляющая эквивалентной скорости дрейфа, зависящая от входной угловой скорости ω_п угла ориентации волны φ₀ с начальной фазой φ₂;

φ₂ - начальная фаза зависимости показаний гироскопа.

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

$${\omega }_2({\varphi }_0+{\phi }_3,{\omega }_{п})={\omega }_{до}({\varphi }_0+{\varphi }_1)+{\omega }_{зв}+{\omega }_{п2}+{\omega }_{д2}({\varphi }_0+{\varphi }_2,{\omega }_{п}),\text{ (11)}$$

где

ω_п2 - второе значение эталонной угловой скорости, задаваемой платформе;

φ_д2(φ₁+φ₂) - систематическая составляющая эквивалентной скорости дрейфа, зависящая от входной угловой скорости ω_п и угла ориентации волны φ₀ с начальной фазой φ₂.

Экспериментальные исследования показали, что систематическая составляющая ω_до(φ₀+φ₁) имеет синусоидальную зависимость от угла φ₀ с периодичностью 180°. Составляющие ω_3В,ω_п1,ω_п2 в выражениях (10) (11) являются постоянными составляющими. Составляющие ω_д1(φ₀+φ₂), ω_д2(φ₀+φ₂), как показано на фиг.1, имеют синусоидальную зависимость от угла φ₀. Фазы φ₁ и φ₂ могут быть разными и быть обусловлены конструкционными параметрами и погрешностями выработки синфазных и квадратурных составляющих сигнала датчика угла гироскопа.

Таким образом, зависимости ω₁(φ₀+φ₃,ω_п) и ω₂(φ₀+φ₃,ω_п) являются синусоидальными и периодическими. Их изменения на заданном угловом диапазоне ориентации волны при измерениях можно характеризовать шириной трубок соответственно Δ₁ и Δ₂.

Учитывая, что ω_д1(φ₀+φ₂,ω_n) и ω_д2(φ₀+φ₂,ω_n) также зависят от погрешностей выработки синфазных и квадратурных составляющих сигналов по информационным каналам косинуса и синуса, то подбирают корректирующие коэффициенты к этим синфазным и квадратурным составляющим таким образом, чтобы обеспечить минимальную разницу между Δ₁ и Δ₂, что происходит, когда ω_д1(φ₀+φ₂,ω_n) и ω_д2(φ₀+φ₂,ω_n) стремятся к нулевой величине.

Используя соотношения (7) с учетом (6) и (9), получаем рабочую формулу (8) предлагаемого способа определения угла ориентации стоячей волны.

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

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

а) новые действия по установке гироскопа на платформу поворотного стола и заданию с помощью платформы последовательно двух значений угловых скоростей по входной оси гироскопа;

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

в) новое действие по определению разности в ширине трубок изменения показаний гироскопа при двух угловых скоростях по входной оси гироскопа;

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

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

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

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

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

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

Расчет проводился с применением соотношений (7) и использованием следующих параметров:

- частота обновления информации f=50 Гц;

- длительность такта обновления информации ΔТ=0,02 с;

- нормированное значение большой полуоси эллипсовидной траектории точки в плоскости CS, изображающей состояние резонатора, на к-ом и к-м тактах а_к-1=а_к=1;

- нормированное значение малой полуоси эллипсовидной траектории точки в плоскости CS, изображающей состояние резонатора, b_к-1=b_к=0,1а_к;

- погрешности выработки синфазных составляющих по каналам косинуса и синуса δС_ф,к, δС_ф,к-1, δS_ф,к, δS_ф,к-1 и квадратурных составляющих δС_кв,к, δС_кв,к-1, δS_кв,к, δS_кв,к-1 составляет 1% соответственно от а_к и b_к;

- собственная скорость дрейфа, не зависящая от угловой входной скорости $${\omega }_{{Д}_0}=-1\cdot \sin (2{\varphi }_0+{\varphi }_1)°/ч$$

.

На фиг.1 оси имеют обозначения: х1=φ₀, Y1=δΔ(φ₀).

Выполнено компьютерное моделирование предлагаемого способа.

При моделировании использовали представленные выше параметры и ошибки гироскопа. По входной оси гироскопа задавали первое значение угловой скорости ω_п=10°/ч.

При этом задавали также в соответствии с (10)

ω_Д0=-1·sin(2φ₀+φ₁) °/ч, где φ₁=0; ω_зв=12,4 °/ч.

Результаты расчетов представлены на фиг.2 зависимостью 1.

Результаты расчетов при втором значении входной угловой скорости

ω_п2=20 °/ч представлены на фиг.3 зависимостью 1.

На фиг.2 и фиг.3 оси имеют обозначение Х2=ХЗ=φ₀, Y2=ω₁, Y3=ω₂.

Из сравнения зависимостей 1 на фиг.2 и фиг.3 определяли разность в ширине трубок изменения показаний, которая составила 0,25 °/ч.

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

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

Погрешности выработки приращения угла ориентации волны на такте рассчитывались при разных угловых скоростях, направленных по входной оси гироскопа, которые представлены на фиг.1 в виде зависимостей 1, 2, 3 соответственно для угловых скоростей 12,4 °/ч, 24,8 °/ч, 37,2 °/ч.

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

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

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

Claims (1)

1. Способ определения угла ориентации стоячей волны в твердотельном волновом гироскопе, включающий измерение сигналов с синусного и косинусного каналов датчика угла гироскопа, выделение синфазных и квадратурных составляющих этих сигналов, определение тригонометрических функций синуса и косинуса угла ориентации стоячей волны, используя синфазные и квадратурные составляющие сигналов, вычисление угла ориентации волны относительно резонатора гироскопа на каждом такте вычислений по формуле с применением тригонометрических функций синуса и косинуса угла волны, полученных для к-1-го и к-го тактов, а также значение угла ориентации волны на к-1-м такте, отличающийся тем, что предварительно располагают гироскоп на платформе поворотного стола таким образом, чтобы его входная ось совпадала по направлению с осью вращения платформы и была направлена вертикально, задают платформе первое эталонное значение угловой скорости и в угловом диапазоне ориентации стоячей волны относительно резонатора ±90° измеряют показания гироскопа, задают второе эталонное значение угловой скорости платформы и на этом же угловом диапазоне ориентации волны измеряют показания гироскопа, определяют разницу в ширине трубок изменений показаний гироскопа в заданном угловом диапазоне ее ориентации при разных угловых скоростях платформы, подбирают корректирующие коэффициенты для вырабатываемых приборных значений синфазных и квадратурных составляющих синусного и косинусного каналов датчика угла гироскопа, минимизирующих эту разницу в ширине трубок изменения показаний гироскопа, после чего в рабочем режиме, используя подобранные корректирующие коэффициенты, определяют угол ориентации стоячей волны относительно резонатора по формуле
   φ_0,к=0,25·(n_кm_к-1-m_кn_к-1)+φ_0,к-1,
   где
   $$\begin{array}{l}{m}_{к}=[{\left({С}_{фп,к}-\Delta {С}_{ф}\right)}^2+{\left({С}_{квп,к}-\Delta {С}_{кв}\right)}^2-{\left(S_{фп,к}-\Delta S_{ф}\right)}^{{}_2}-\left(S_{квп,к}-\Delta S_{кв}\right){}^2]\cdot \\ \cdot \{[{\left({С}_{фп,к}-\Delta {С}_{ф}\right)}^2+{\left({С}_{квп,к}-\Delta {С}_{кв}\right)}^2+{\left(S_{фп,к}-\Delta S_{ф}\right)}^{{}_2}+{\left(S_{квп,к}-\Delta S_{кв}\right)}^2+\\ +2(\left({С}_{фп,к}-\Delta {С}_{ф}\right)\left(S_{квп,к}-\Delta S_{кв}\right)-\left({С}_{кв,к}-\Delta {С}_{кв}\right)\left(S_{фп,к}-\Delta S_{ф}\right))]\cdot \\ \cdot [{\left({С}_{фп,к}-\Delta {С}_{ф}\right)}^2+{\left({С}_{квп,к}-\Delta {С}_{кв}\right)}^2+{\left(S_{фп,к}-\Delta S_{ф}\right)}^{{}_2}+{\left(S_{квп,к}-\Delta S_{кв}\right)}^2-\\ -2{(\left({С}_{фп,к}-\Delta {С}_{ф}\right)\left(S_{квп,к}-\Delta S_{кв}\right)-\left({С}_{квп,к}-\Delta {С}_{кв}\right)\left(S_{фп,к}-\Delta S_{ф}\right))]\}}^{-\frac{1}{2}};\end{array}$$

   

   
   $$\begin{array}{l}{n}_{к}=2[\left({С}_{фп,к}-\Delta {С}_{ф}\right)\left(S_{фп,к}-\Delta S_{ф}\right)+\left({С}_{квп,к}-\Delta {С}_{кв}\right)\left(S_{квп,к}-\Delta S_{кв}\right)]\cdot \\ \cdot \{[{\left({С}_{фп,к}-\Delta {С}_{ф}\right)}^2+{\left({С}_{квп,к}-\Delta {С}_{кв}\right)}^2+{\left(S_{фп,к}-\Delta S_{ф}\right)}^{{}_2}+{\left(S_{квп,к}-\Delta S_{кв}\right)}^2+\\ +2(\left({С}_{фп,к}-\Delta {С}_{ф}\right)\left(S_{квп,к}-\Delta S_{кв}\right)-\left({С}_{квп,к}-\Delta {С}_{кв}\right)\left(S_{фп,к}-\Delta S_{ф}\right))]\cdot \\ \cdot [{\left({С}_{фп,к}-\Delta {С}_{ф}\right)}^2+{\left({С}_{квп,к}-\Delta {С}_{кв}\right)}^2+{\left(S_{фп,к}-\Delta S_{ф}\right)}^{{}_2}+{\left(S_{квп,к}-\Delta S_{кв}\right)}^2-\\ -2(\left({С}_{фп,к}-\Delta {С}_{ф}\right)\left(S_{квп,к}-\Delta S_{кв}\right)-\left({С}_{квп,к}-\Delta {С}_{кв}\right){\left(S_{фп,к}-\Delta S_{ф}\right))]\}}^{-\frac{1}{2}};\end{array}$$

   

   
   $$\begin{array}{l}{m}_{к-1}=[{\left({С}_{фп,к-1}-\Delta {С}_{ф}\right)}^2+{\left({С}_{квп,к-1}-\Delta {С}_{кв}\right)}^2-{\left(S_{фп,к-1}-\Delta S_{ф}\right)}^{{}_2}-\left(S_{квп,к-1}-\Delta S_{кв}\right){}^2]\cdot \\ \cdot \{[{\left({С}_{фп,к-1}-\Delta {С}_{ф}\right)}^2+{\left({С}_{квп,к-1}-\Delta {С}_{кв}\right)}^2+{\left(S_{фп,к-1}-\Delta S_{ф}\right)}^{{}_2}+{\left(S_{квп,к-1}-\Delta S_{кв}\right)}^2+\\ +2(\left({С}_{фп,к-1}-\Delta {С}_{ф}\right)\left(S_{квп,к-1}-\Delta S_{кв}\right)-\left({С}_{квп,к-1}-\Delta {С}_{кв}\right)\left(S_{фп,к-1}-\Delta S_{ф}\right))]\cdot \\ \cdot [{\left({С}_{фп,к-1}-\Delta {С}_{ф}\right)}^2+{\left({С}_{квп,к-1}-\Delta {С}_{кв}\right)}^2+{\left(S_{фп,к-1}-\Delta S_{ф}\right)}^{{}_2}+{\left(S_{квп,к-1}-\Delta S_{кв}\right)}^2-\\ -2(\left({С}_{фп,к-1}-\Delta {С}_{ф}\right)\left(S_{квп,к-1}-\Delta S_{кв}\right)-\left({С}_{квп,к-1}-\Delta {С}_{кв}\right){\left(S_{фп,к-1}-\Delta S_{ф}\right))]\}}^{-\frac{1}{2}};\end{array}$$

   

   
   $$\begin{array}{l}{n}_{к-1}=2[\left({С}_{фп,к-1}-\Delta {С}_{ф}\right)\left(S_{фп,к-1}-\Delta S_{ф}\right)+\left({С}_{квп,к-1}-\Delta {С}_{кв}\right)\left(S_{квп,к-1}-\Delta S_{кв}\right)]\cdot \\ \cdot \{[{\left({С}_{фп,к-1}-\Delta {С}_{ф}\right)}^2+{\left({С}_{квп,к-1}-\Delta {С}_{кв}\right)}^2+{\left(S_{фп,к-1}-\Delta S_{ф}\right)}^{{}_2}+{\left(S_{квп,к-1}-\Delta S_{кв}\right)}^2+\\ +2(\left({С}_{фп,к-1}-\Delta {С}_{ф}\right)\left(S_{квп,к-1}-\Delta S_{кв}\right)-\left({С}_{квп,к-1}-\Delta {С}_{кв}\right)\left(S_{фп,к-1}-\Delta S_{ф}\right))]\cdot \\ \cdot [{\left({С}_{фп,к-1}-\Delta {С}_{ф}\right)}^2+{\left({С}_{квп,к-1}-\Delta {С}_{кв}\right)}^2+{\left(S_{фп,к-1}-\Delta S_{ф}\right)}^{{}_2}+{\left(S_{квп,к-1}-\Delta S_{кв}\right)}^2-\\ -2(\left({С}_{фп,к-1}-\Delta {С}_{ф}\right)\left(S_{квп,к-1}-\Delta S_{кв}\right)-\left({С}_{квп,к-1}-\Delta {С}_{кв}\right){\left(S_{фп,к-1}-\Delta S_{ф}\right))]\}}^{-\frac{1}{2}};\end{array}$$

   

   
   С_фп,к - приборная синфазная составляющая сигнала датчика угла гироскопа по косинусному каналу на к-м такте вычислений;
   С_квп,к - приборная квадратурная составляющая сигнала датчика угла по косинусному каналу на к-м такте;
   S_фп,к - приборная синфазная составляющая сигнала датчика угла по синусному каналу на к-м такте;
   S_квп,к - приборная квадратурная составляющая сигнала датчика угла по синусному каналу на к-м такте;
   С_фп,к-1 - приборная синфазная составляющая сигнала датчика угла по косинусному каналу на к-1-ом такте;
   С_квп,к-1 - приборная квадратурная составляющая сигнала датчика угла по косинусному каналу на к-1-ом такте;
   S_фп,к-1 - приборная синфазная составляющая сигнала датчика угла по синусному каналу на к-1-ом такте;
   S_квп,к-1 - приборная квадратурная составляющая сигнала датчика угла по синусному каналу на к-1-ом такте;
   ΔС_ф - коэффициент, корректирующий постоянную ошибку определения синфазной составляющей сигнала датчика угла по косинусному каналу;
   ΔС_кв - коэффициент, корректирующий постоянную ошибку определения квадратурной составляющей сигнала датчика угла по косинусному каналу;
   ΔS_ф - коэффициент, корректирующий постоянную ошибку определения синфазной составляющей сигнала датчика угла по синусному каналу;
   ΔS_кв - коэффициент, корректирующий постоянную ошибку определения квадратурной составляющей сигнала датчика угла по синусному каналу;
   φ_0,к-1 - угол ориентации волны относительно резонатора, определенный на к-1-ом такте.

Images