Научный журнал
Современные наукоемкие технологии
ISSN 1812-7320
"Перечень" ВАК
ИФ РИНЦ = 0,899

АЛГОРИТМЫ ВЫЧИСЛЕНИЙ В СИСТЕМЕ АНАЛИТИЧЕСКОГО ГИРОКОМПАСИРОВАНИЯ С ВРАЩАЮЩИМСЯ ДАТЧИКОМ УГЛОВОЙ СКОРОСТИ

Белянин Л.Н. 1 Дмитриев В.С. 1 Ву Д.К. 1
1 ФГАОУ ВО «Национальный исследовательский Томский политехнический университет»
Предлагаются алгоритмы вычислений параметров ориентации – углов Эйлера – Крылова – в системе аналитического гирокомпасирования, содержащей гироскопический датчик угловой скорости, установленный на одноосной вращающейся платформе, два акселерометра и вычислитель. По оси вращения платформы установлен точный преобразователь угла в цифровой код (энкодер). Платформа приводится во вращение с помощью привода, обеспечивающего постоянную, строго заданную угловую скорость. Прибор может быть наклонён по отношению к вертикали места в любом направлении на угол до 15–20 градусов. Основное внимание уделено поиску приёмов и алгоритмов вычисления азимутального угла, которые позволили бы исключить или уменьшить влияние погрешностей датчика угловой скорости на точность определения этого угла. К таким погрешностям относятся: смещение нуля, изменение крутизны характеристики (масштаба, чувствительности), динамические погрешности (амплитудные и фазовые) при условии их постоянства в пуске, а также шум в выходном сигнале датчика угловой скорости. Предлагается сначала вычислять угол, характеризующий ориентацию проекции вектора угловой скорости вращения Земли на экваториальную плоскость прибора, а затем вычислять азимутальный угол на основе информации об указанном угле, а также углах, характеризующих ориентацию прибора по отношению к вертикали места (плоскости горизонта). Угол, характеризующий ориентацию проекции вектора угловой скорости вращения Земли на экваториальную плоскость прибора, определяется по результатам отсчётов углов с энкодера с учётом знака производной в те моменты времени, когда сигнал с датчика угловой скорости при вращении платформы становится равным нулю.
определение ориентации неподвижного прибора
система аналитического гирокомпасирования
вращающийся датчик угловой скорости
алгоритмы вычисления
повышение точности
1. Опыт создания гироскопического инклинометра / Л.Н. Белянин, А.Н. Голиков, В.М. Мартемьянов, С.Н. Самойлов // Автоматизация и информационное обеспечение технологических процессов в нефтяной промышленности, том 2: Сб. трудов ОАО НПФ «Геофит» ВНК. – Томск: Изд-во Том. ун-та, 2002. – С. 34–49.
2. Попов Г.В. Быстрое определение азимута // Авиакосмическое приборостроение. – 2014. – № 1. – С. 3–9.
3. Редькин С.П. Азимутальное ориентирование с помощью датчиков угловой скорости // Авиакосмическое приборостроение. – 2014. – № 4. – С. 10–22.
4. Биндер Я.И. Аналитическое компасирование в инклинометрии скважин малого диаметра // Гироскопия и навигация. – 2003. – № 2. – С. 38–46.
5. Patent US 8099876B1, МПК G01C 19/38. Azimuth determination system and method therefor / Angelo Truncale, Robert J. Atmur, Joseph E. Justin, Brian Barsamian; assignee The boeing company, Chicago, IL, US. – Appl. No. 12/276,056; filled Nov. 21, 2008; pub. date Jan. 24, 2012.
6. Цыбряева И.В., Гуськов А.А., Кривошеев С.В., Стрелков А.Ю. Методы повышения точности инклинометрии скважин гироскопическим инклинометром ИГН73-100/80 // Научно-технический вестник «Каротажник». – 2013. – № 4. – С. 81–89.
7. Patent US20100088063 A1, МПК G06F 15/00, G01B 1/00. Method and apparatus for precision azimuth measurement / Darren R. Laughlin; assignee A-tech corporation, Albuquerque, NM, US. – Appl. No. 12/572,181; filled Oct. 1, 2009; pub. date Apr. 08, 2010.
8. Patent US7412775 B1, МПК G01C 19/38, G01C 19/02. Gyroscope north seeker system and method / Drew A. Karnick, Timothy J. Hanson; assignee Honeywell International Inc., Morristown, NJ, US. Appl. No. 11/833,398; filled Aug. 3, 2007; pub. date Aug. 19, 2008.
9. Ву Д.К. Алгоритмы вычисления параметров в системе ориентации и навигации горнопроходческого комбайна. НАУКА. ТЕХНОЛОГИИ. ИННОВАЦИИ // Сборник научных трудов в 10 ч. / под ред. д.т.н. Гуськова А.В. – Новосибирск: Изд-во НГТУ, 2017. – Часть 10. – С. 102–106.
10. Белянин Л.Н. Калибровка модуля акселерометров в условиях неопределенности направления их осей чувствительности / Л.Н. Белянин, Е.В. Якимова // Контроль. Диагностика. – 2011. – Специальный выпуск. – С. 64–70.
11. Лазарева Т.Я. Основы теории автоматического управления: Учебное пособие. 2-е изд., перераб. и доп. / Т.Я. Лазарева, Ю.Ф. Мартемьянов. – Тамбов: Изд-во Тамб. гос. техн. ун-та, 2004. – 352 с.
12. Лукьянов Д.П. Прикладная теория гироскопов / Д.П. Лукьянов, В.Я. Распопов, Ю.В. Филатов. – СПб.: ГНЦ РФ ОАО «Концерн «ЦНИИ «Электроприбор», 2015. – 316 с.
13. Подчезерцев В.П. Компоненты модели погрешностей динамически настраиваемого гироскопа / В.П. Подчезерцев, Тан Синюань, Цинь Цзыхао // Авиакосмическое приборостроение. – 2015. – № 1. – С. 8–18.
14. Белянин Л.Н. Алгоритмы вычислений в непрерывном гироскопическом инклинометре. Автоматизация и информационное обеспечение технологических процессов в нефтяной промышленности: Сб. статей / Под ред. А.К. Хорькова. – Томск: Изд-во Том. ун-та, 2002. – Т. 2. – С. 50–63.

Принцип аналитического гирокомпасирования широко применяется для определения ориентации неподвижных объектов по отношению к связанной с Землей горизонтальной, географически ориентированной системе координат, в частности, в инклинометрии скважин [1–3]. В общем случае в состав системы аналитического гирокомпасирования входят три датчика линейного ускорения (акселерометра), оси чувствительности которых направлены вдоль осей связанной с объектом системы координат, три гироскопических датчика угловой скорости (ДУС) с осями чувствительности, направленными параллельно тем же осям, и вычислитель. Если углы отклонения объекта по отношению к плоскости горизонта (вертикали места) ограничены, то возможен усечённый вариант системы, содержащий два акселерометра, два ДУС и вычислитель.

При использовании в системе акселерометров навигационного класса обеспечивается достаточно высокая точность определения параметров, характеризующих ориентацию объекта по отношению к плоскости горизонта. Проблемным является обеспечение необходимой точности определения ориентации объекта в азимуте. Для того, чтобы погрешность определения азимутального угла не превышала 1 ° на средних широтах, погрешности измерения ДУС не должны превышать десятых долей град/час. Для обеспечения погрешности определения азимутального угла на уровне 10 угловых минут погрешности ДУС не должны превышать сотых долей град/час. Для сравнения: тело, вращающееся вокруг оси с постоянной скоростью и делающее один оборот за год, вращается со скоростью 0,0416 град/час. Поэтому актуальной является разработка способов, приёмов и алгоритмов, позволяющих исключить либо уменьшить влияние погрешностей ДУС на точность определения азимутального угла.

Эффективным в этом отношении является метод, основанный на обработке результатов двух замеров каждым ДУС, у которых пространственное положение осей чувствительности изменяется на противоположное [4, 5]. В этом случае исключается влияние на результат гирокомпасирования смещения нулей ДУС и изменения их крутизны характеристики при условии, что в течение всего цикла измерения эти параметры остаются неизменными. Иногда этот метод называют методом двойного гирокомпасирования. Однако устройство для точной переориентации осей чувствительности ДУС делает прибор более сложным и дорогим.

Находит применение также метод, при котором один ДУС вращается с постоянной, достаточно малой скоростью вокруг вертикальной оси, в результате чего ось чувствительности ДУС непрерывно меняет свою ориентацию по отношению к вектору угловой скорости вращения Земли. Искомый угол, характеризующий ориентацию прибора в азимуте, определяется путём минимизации функции невязки между сигналом ДУС и вычисленным значением гармонической функции угла в определённые моменты времени. Такой метод определения азимута реализован в гироскопическом инклинометре ИГН73-100/80 [6]. Недостатком метода является необходимость выполнения большого объёма вычислений. Предлагаются также другие варианты определения направления на Север с помощью вращающегося вокруг вертикальной оси ДУС [7, 8].

Настоящая статья посвящена разработке приёмов и алгоритмов вычислений в системе аналитического гирокомпасирования, в которой используется один ДУС, установленный на платформе, вращающейся с постоянной скоростью. Причём жёстких требований к вертикализации прибора не предъявляется. Прибор может быть наклонён по отношению к вертикали места в любом направлении на угол до 15–20 градусов. Прибор входит в состав разрабатываемой системы ориентации и навигации горнопроходческого комбайна [9]. Составные части прибора и их положение относительно опорной и связанной с корпусом прибора системы координат показаны на рис. 1.

bel1.tif

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

На рис. 1 обозначено:

-О0ξ0η0ζ0 – горизонтальная, географически ориентированная, правая опорная система координат (СК), ось О0ζ0 которой направлена по вертикали вниз; оси О0ξ0 и О0η0 расположены в плоскости горизонта, причем ось О0ξ0 направлена на географический Север;

- О0X0Y0Z0 – правая СК, связанная с корпусом прибора; ось О0Z0 направлена вдоль оси вращения платформы в сторону нижней части прибора; оси О0X0 и О0Y0 лежат в экваториальной плоскости прибора;

- О0XПYПZП – правая СК, связанная с платформой; ось О0ZП совпадает с осью О0Z0; оси О0XП и О0YП в исходном положении совпадают с осями О0X0 и О0Y0 соответственно;

- Ψ, Θ, Φ – углы Эйлера – Крылова, характеризующие ориентацию СК О0X0Y0Z0 по отношению к СК bel01.wmf; угол Ψ будем называть азимутальном углом прибора;

- ε – угол поворота платформы вокруг оси О0ZП;

- bel02.wmf, bel03.wmf – акселерометры, установленные по осям О0X0 и О0Y0 соответственно.

Принцип работы прибора

При вращении платформы с постоянной скоростью сигнал с ДУС изменяется по гармоническому закону в зависимости от угла ε, принимая нулевое значение дважды за один оборот. Определяются два первых последовательных угла (ε1, ε2), характеризующие положения платформы относительно корпуса в моменты времени, когда сигнал ДУС равен нулю. Вычисляется среднее арифметическое значение этих углов и таким образом находится угловое положение платформы, в котором ось чувствительности ДУС направлена либо на условный Север (εN), если в положении платформы, характеризующемся углом ε1, производная сигнала ДУС положительна, либо на условный Юг (εS), если в положении платформы, характеризующемся углом ε1, производная сигнала ДУС отрицательна (рис. 2).

Под условным Севером здесь понимается такое значение угла ε, при котором проекция вектора угловой скорости вращения Земли на ось чувствительности ДУС положительна и максимальна по модулю. Под условным Югом – значение угла ε, при котором проекция вектора угловой скорости вращения Земли на ось чувствительности ДУС отрицательна и максимальна по модулю. При строго вертикальном положении оси вращения платформы понятия «Север условный» и «Юг условный» совпадают соответственно с понятиями «Север» и «Юг» истинные. Зная угол εN (или εS) и ориентацию платформы относительно вертикали, можно определить азимутальный угол корпуса прибора. Ниже предлагается алгоритм вычисления угла Ψ по известным значениям углов εN, Θ и Φ. Алгоритмы вычисления углов Θ и Φ по информации с акселерометров здесь не рассматриваются. С ними можно ознакомиться, например, в [10].

bel2.tif

Рис. 2. График, поясняющий метод определения ориентации проекции вектора угловой скорости вращения Земли на экваториальную плоскость прибора Nусл. – условный Север; Sусл. – условный Юг

Вычисление угла Ψ по известным Θ, Φ, ε и φ0

Обозначим проекции вектора абсолютной угловой скорости вращения трёхгранника осей СК О0ξ0η0ζ0: bel04.wmf, bel05.wmf, bel06.wmf соответственно. Поскольку вращения трёхгранника происходят по причине вращения Земли с угловой скоростью bel07.wmf, то с учётом ориентации СК О0ξ0η0ζ0 относительно Земли имеем

bel08.wmf; bel09.wmf; bel10.wmf,

где φ0 – широта местоположения прибора.

Найдём проекции вектора bel11.wmf на ось О0XП. С этой целью сначала запишем выражения для проекций вектора bel12.wmf на оси СК О0X0Y0Z0:

bel13.wmf;

bel14.wmf

bel15.wmf

Искомая проекция вектора угловой скорости вращения Земли на ось чувствительности ДУС О0XП определится

bel16.wmf. (1)

При ε = εN функция bel17.wmf имеет экстремум и при этом положительна. Это позволяет записать следующие уравнение и неравенство:

bel18.wmf; (2)

bel19.wmf. (3)

Поставив в уравнение (2) и неравенство (3) выражения для угловых скоростей bel20.wmf и bel21.wmf, получим

bel22.wmf (4)

bel23.wmf (5)

где

bel24.wmf;

bel25.wmf;

bel26.wmf;

bel27.wmf;

bel28.wmf;

bel29.wmf.

Уравнение (4) даже с учётом неравенства (5) не даёт однозначного решения для угла Ψ. Поступим следующим образом. Будем считать неизвестными в полученных уравнении и неравенстве значения sin Ψ и cos Ψ, добавив к этим уравнению и неравенству уравнение, связывающее между собой sin Ψ и cos Ψ

bel30.wmf. (6)

Выразим из (4) cos Ψ:

bel31.wmf. (7)

Подставив (7) в (6), получим квадратное уравнение относительно sin Ψ

bel32.wmf,

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

bel33.wmf. (8)

Дискриминант этого уравнения имеет вид

bel34.wmf.

Если bel35.wmf, то уравнение (8) имеет два решения. Первое решение будет

bel36.wmf. (9)

Подставив первое решение для sin Ψ в (7), получим

bel37.wmf. (10)

Тогда второе решение уравнения (8) будет иметь вид

bel38.wmf. (11)

Соответственно, второе решение для cos Ψ будет иметь вид

bel39.wmf. (12)

Из двух пар значений bel40.wmf и bel41.wmf выберем ту, которая удовлетворяет неравенству (5). Эта полученная пара значений синуса и косинуса искомого угла позволяет определить его единственное значение. Главное значение угла вычисляется по формуле

bel42.wmf. (13)

Здесь под главным значением понимается значение острого угла, который может быть как положительным, так и отрицательным. Действительное значение угла Ψ будет определяться с учётом квадранта, в котором лежит угол по известным формулам приведения. Квадрант определяется по знакам sin Ψ и sin Ψ. Алгоритм определения действительного значения угла Ψ представлен в виде табл. 1.

Таблица 1

Определение действительного значения угла Ψ

 

Квадрант, в котором лежит угол Ψ

Действительное значение угла Ψ

sin Ψ > 0, cos Ψ > 0

первый

Ψ = Ψ*

sin Ψ > 0, cos Ψ < 0

второй

Ψ = Ψ* + 108 °

sin Ψ < 0, cos Ψ < 0

третий

Ψ = Ψ* + 108 °

sin Ψ < 0, cos Ψ > 0

четвёртый

Ψ = Ψ* + 360 °

Достоверность полученных алгоритмов подтверждена численными экспериментами.

Рассмотрим влияние основных погрешностей ДУС на точность определения угла εN .

Смещение нуля и изменение крутизны характеристики ДУС (чувствительности, масштаба) от пуска к пуску

Если указанные параметры остаются неизменными в пуске, то они не оказывают влияния на точность определения угла εN, а, следовательно, и угла Ψ. На рис. 3 показано изменение выходного сигнала ДУС в функции от угла Ψ при наличии смещения нуля и изменения чувствительности.

bel3.tif

Рис. 3. Выходной сигнал ДУС при наличии смещения нуля и изменения крутизны его характеристики

Видно, что углы ε1 и ε2 при этом изменяются, но значение εN, как среднее арифметическое между ними, остаётся неизменным.

Фазовая погрешность ДУС

При измерении ДУС угловой скорости, изменяющейся по гармоническому закону, возникают известные динамические погрешности, которые относятся к категории методических и делятся на амплитудные и фазовые. При постоянной скорости вращения платформы эти погрешности постоянны. Выше показано, что изменение εN. Наличие фазовой погрешности (фазового сдвига) ДУС приводит к одинаковым погрешностям определения углов ε1 и ε2, следовательно, и угла εN.

Для ДУС, построенных на основе двухстепенного гироскопа с механическим упругим элементом, и компенсационного типа при условии безынерционности всех звенев, образующих цепь компенсации (цепь обратной связи), передаточная функция представляет собой передаточную функцию устойчивого колебательного звена. Сдвиг фазы такого звена определяется формулой [11]

bel43.wmf, (14)

где Δφ – сдвиг фазы звена;

ξ – коэффициент затухания звена;

η – безразмерная частота, определяемая отношением круговой частоты изменения угловой скорости (скорости вращения платформы) к собственной частоте ДУС [12].

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

Таблица 2

Результаты расчёта сдвига фаз по формуле (14)

п/п

Коэффициент затухания звена ξ

Безразмерная частота η

Фазовая погрешность Δφ, угл. мин.

1

0,6

0,000003

–0,012600

2

0,6

0,000030

–0,123600

3

0,6

0,000300

–12,37800

4

0,6

0,003000

–12,37620

5

0,6

0,030000

–123,8202

В ДУС, построенном на основе трёхстепенного гироскопа, в том числе на основе динамически настраиваемого гироскопа, в случае безынерционных цепей обратных связей, передаточная функция представляет собой передаточную функцию апериодического звена первого порядка. Фазовый сдвиг этого звена определяется формулой

bel44.wmf, (15)

где η – безразмерная частота – отношение круговой частоты изменения угловой скорости к частоте изгиба асимптотической амплитудно-частотной характеристики звена.

Результаты расчёта сдвига фаз по формуле (15) представлены в табл. 3.

Таблица 3

Результаты расчёта сдвига фаз по формуле (15)

п/п

Безразмерная частота ζ

Фазовая погрешность Δφ, угл. мин.

1

0,000003

–0,0103200

2

0,000030

–0,1031400

3

0,000300

–1,0313400

4

0,003000

–10,313520

5

0,030000

–103,10454

Из приведенных таблиц следует, что при определённых условиях величина ошибки в определении угла εN, вызванной фазовой погрешностью, может принимать недопустимо большие значения. Например, при скорости вращения платформы 10 град/с и собственной частоте колебаний ДУС на основе двухстепенного гироскопа, равной 10 Гц, при ξ = 0,6 имеем η = 0,00278, Δφ = 11,46 угл. мин.

Для исключения влияния фазовой погрешности ДУСа на результат определения угла εN предлагается после снятия отсчётов углов bel45.wmf и bel46.wmf при вращении платформы в одну сторону изменить направление вращения на противоположное и снять отсчёты bel47.wmf и bel48.wmf. В обозначениях углов ε1 и ε2 стрелки сверху означают направление вращения платформы (направление изменения угла ε). Стрелки после индекса означают знак производной при данном значении угла. Стрелка вверх – производная положительна, стрелка вниз – производная отрицательна.

Для определения значения угла εN вычисляются значения углов bel49.wmf и bel50.wmf и находится среднее значение между ними.

Устранение влияния фазовой погрешности в этом случае поясняется рис. 4.

На рис. 4 изображено:

- кривая 1 – график, отражающий изменение сигнала ДУС при отсутствии сдвига фаз;

- кривая 2 – график, отражающий изменение сигнала ДУС при вращении платформы в положительном направлении при наличии сдвига фаз;

- кривая 3 – график, отражающий изменение сигнала ДУС при вращении платформы в отрицательном направлении при наличии сдвига фаз.

bel4.tif

Рис. 4. К пояснению способа устранения влияния фазовой погрешности ДУС

При вращении платформы в положительном направлении будем иметь

bel51.wmf. (16)

При вращении платформы в отрицательном направлении будем иметь

bel52.wmf. (17)

Вычислив среднее значение углов, определяемых выражениями (16) и (17), найдём

bel53.wmf. (18)

Значение εN, вычисляемое по формуле (18), не зависит от запаздывания, вносимого ДУС.

Запаздывание в определении угла ε по показаниям энкодера

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

Однако в действительности за счёт наличия порога чувствительности компаратора и некоторой задержки выдачи кода схемой управления энкодера возникнет запаздывание. Если принять величину указанного запаздывания постоянной и не зависящей от направления вращения, то, очевидно, что при вычислении угла εN по формуле (18) это запаздывание на точность определения угла εN влияния также оказывать не будет.

Помехи (шум) в выходном сигнале ДУС

Реальный сигнал на выходе ДУС отличается от того, что изображён на рис. 3 и 4. Он содержит шум с широким частотным спектром. Характер спектра зависит от типа используемого ДУС, внешних механических воздействий на ДУС в процессе измерений, а также характера и уровня внутренних и внешних наводок на электронную часть прибора [13]. Например, в ДУС на основе роторных гироскопов присутствие сравнительно низкочастотных гармоник помех в выходном сигнале обусловлено воздействием на прибор внешних вибраций и ударов, шумом, создаваемым шарикоподшипниками подвеса ротора, дисбалансом последнего и другими дефектами механической части ДУС. Высокочастотные гармоники, как правило, обусловлены помехами электрического происхождения, возникающими в электронной части ДУС, а также внешними наводками.

Наличие помех в выходном сигнале ДУС приведёт к разбросу моментов срабатывания компаратора относительно момента равенства нулю сигнала ДУС при отсутствии помех. Возникает погрешность Δε, которая может быть как положительной, так и отрицательной. Возникновение погрешности Δε при наличии в выходном сигнале ДУС одной высокочастотной гармонической составляющей помехи показано на рис. 5.

bel5a.tif bel5b.tif

а) б)

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

Если принять, что погрешность Δε является случайной величиной с нормальным законом распределения относительно нулевого значения, то для снижения влияния этого фактора на точность определения угла εN можно рекомендовать следующий приём. Наблюдения и вычисления ведутся в течение промежутка времени, за который платформа делает n оборотов в одну и n оборотов в другую сторону. При этом угол bel54.wmf вычисляется как среднее между углами bel55.wmf

bel56.wmf, (19)

где n – число пар отсчётов (число оборотов платформы).

И угол bel57.wmf вычисляется как среднее между углами bel58.wmf

bel59.wmf. (20)

Проведённый анализ позволяет записать итоговую формулу, вычисление по которой угла εN обеспечит наибольшую точность

bel60.wmf. (21)

Влияние шума в выходном сигнале ДУС на точность определения азимутального угла будет снижаться по мере увеличения числа n, а это напрямую связано со временем определения азимутального угла.

Последовательность действий по определению углов Ψ, Θ, Φ

Приведём последовательность действий по определению углов Эйлера – Крылова Ψ, Θ, Φ, характеризующих ориентацию прибора по отношению к горизонтальной, географически ориентированной системе координат.

1. В вычислитель вводится значение широты места φ0, а также значение ускорения силы тяжести на экваторе g0.

2. Вычисляется величина ускорения силы тяжести g места [14].

3. Вычисляются углы Θ и Φ по сигналам с акселерометров с учётом найденного значения ускорения силы тяжести.

4. Платформа приводится во вращение в положительном направлении и с энкодера снимаются и записываются в память вычислителя отсчёты bel61.wmf и bel62.wmf в течение промежутка времени, необходимого для снятия n пар отсчётов.

5. Меняется направление вращения платформы на противоположное, снимаются и записываются в память вычислителя n пар отсчётов bel63.wmf и bel64.wmf.

6. Вычисляется значение угла εN по формуле (21).

7. Вычисляются значения bel65.wmf и bel66.wmf по формулам (9–12).

8. Определяется пара значений sin Ψ и cos Ψ, которая удовлетворяет неравенству (5).

9. Вычисляется главное значение угла Ψ по формуле (13).

10. Определяется действительное значение угла Ψ с использованием алгоритмов, представленных в табл. 1.

Заключение

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


Библиографическая ссылка

Белянин Л.Н., Дмитриев В.С., Ву Д.К. АЛГОРИТМЫ ВЫЧИСЛЕНИЙ В СИСТЕМЕ АНАЛИТИЧЕСКОГО ГИРОКОМПАСИРОВАНИЯ С ВРАЩАЮЩИМСЯ ДАТЧИКОМ УГЛОВОЙ СКОРОСТИ // Современные наукоемкие технологии. – 2018. – № 6. – С. 25-33;
URL: https://top-technologies.ru/ru/article/view?id=37028 (дата обращения: 18.09.2021).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.074