СПОСОБ КОРРЕКЦИИ БЕСПЛАТФОРМЕННОЙ ИНЕРЦИАЛЬНОЙ НАВИГАЦИОННОЙ СИСТЕМЫ

Изобретение относится к измерительной технике и может быть использовано для морских, воздушных и наземных объектов. Задачей изобретения является повышение точности бесплатформенной инерциальной навигационной системы (БИНС) путем создания способа непрерывной коррекции курсовертикали.

Классическим алгоритмом вычисления углов ориентации является пересчет показаний ДУС в угловые скорости , , с последующим их интегрированием. Недостатком такой системы является накапливаемая во времени погрешность и, как следствие, ограниченное время работы. Для устранения указанного недостатка в систему необходимо вводить дополнительную информацию, характеризующую реальную угловую ориентацию объекта, в частности беспилотного летательного аппарата (БПЛА). Источником такой информации служат датчики линейного ускорения (ДЛУ). Основные погрешности системы маятниковой коррекции возникают в результате действия постоянных или медленно меняющихся ускорений. В настоящий момент данная проблема решается путем отключения маятниковой коррекции на высокоманевренных участках полета или путем комплексирования блока гироскопов (БГ) с другими системами ориентации (магнитометрическая, видеосистема и др.).

Проблема коррекции курсовертикали БИНС заключается в том, что при маневрировании ЛА моменты времени, когда статические оценки крена и тангажа обладают достаточной точностью, могут возникать недопустимо редко. В связи с этим предлагается способ адаптивной коррекции углов крена и тангажа, в котором коррекция выполняется непрерывно.

Известен способ коррекции инерциальной навигационной системы (ИНС), описанный в патенте RU 2345326 С1, МПК G01C 21/06, опубл. 27.01.2009, автора Прохорцова А.В., принятый нами за прототип.

Сущность способа заключается в следующем. По показаниям акселерометров, входящих в состав ИНС, определяют абсолютное ускорение, действующее на объект, на котором установлена ИНС, по формуле:

,

где g_x - показания акселерометра, измеряющего ускорение по продольной оси объекта, на котором установлена ИНС; g_y - показания акселерометра, измеряющего ускорение по вертикальной оси;

g_z - показания акселерометра, измеряющего ускорение по поперечной оси объекта, на котором установлена ИНС.

В момент времени, когда абсолютное ускорение равно ускорению силы тяжести для местности, где находится ИНС, производится коррекция ИНС по углам тангажа и крена. Этот момент времени соответствует равномерному движению объекта. Для коррекции ИНС по формуле υ=-arcsin(g_x/g) находится истинное значение угла тангажа, а по формуле γ=-arctg(g_z/g_y) находится истинное значение угла крена. Далее показания ИНС по углам тангажа и крена заменяются на вычисленные.

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

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

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

Благодаря предлагаемому способу коррекции БИНС происходит подавление влияния кажущегося ускорения, достаточное для обеспечения необходимой точности оценивания крена и тангажа. В результате зависимость маятниковой коррекции от вида движения ЛА ослабляется до уровня, позволяющего использовать датчики ДУС и ДЛУ средней и низкой точности, в том числе микромеханического типа.

Суть работы способа излагается ниже.

По измерениям датчиков ДУС и ДЛУ оценивается вектор состояния:

Здесь ϑ - угол тангажа; γ - угол крена; V - модуль вектора скорости летательного аппарата относительно земли; i - номер дискретного момента времени измерений датчиков; - априорное нормальное распределение вектора состояния; - априорное математическое ожидание; - априорная ковариационная матрица.

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

.

Дискретная динамическая модель объекта наблюдения имеет вид.

Адаптивный фильтр Калмана строится относительно вектора состояния (1). Вектор наблюдений, обозначаемый далее Z, содержит измерения акселерометров.

Здесь v_i - вектор ошибок измерений с заданной постоянной ковариационной матрицей R.

Функции f_x, f_y, f_z определяют связь измерений перегрузок с параметрами полета. Точные соотношения для этих функций имеют вид.

Здесь V_x, V_y, V_z - проекции вектора земной скорости на связанные оси ЛА.

Полный учет соотношений (5) при ограничении состава датчиков только акселерометрами и гироскопами не представляется возможным, поэтому принимается упрощающее допущение о малости углов атаки и скольжения, а также допущение о постоянстве модуля путевой скорости на интервале дискретизации Δt.

При этом имеют место соотношения: , и уравнения (5) упрощаются.

С учетом (6) матрица Якоби вектора наблюдений (5) имеет вид.

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

где

Текущий вектор состояния (1) рассчитывается по уравнениям Эйлера (2), после чего полагается, что он известен с точностью до погрешностей датчиков. При этом уравнения объекта принимаются в виде:

Здесь x_iq - вектор, в котором компоненты крена и тангажа рассчитаны по соотношениям (3), а компонента скорости принимается равной ее априорному значению на момент текущих измерений;

w_i - вектор возмущений с ковариационной матрицей Q_i:

Вектор возмущений q_ϑ, q_γ задается с учетом точности ДУС, q_V - учитывает изменение скорости.

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

Адаптивный алгоритм фильтра Калмана имеет следующий вид:

Способ коррекции БИНС проверялся при помощи обработки полетных данных вертолета Robinson (фиг. 1 и 2).

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

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

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

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

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

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

Таким образом, с помощью приведенных соотношений решается задача оценивания вектора (1) по наблюдениям (4) с учетом одношагового алгоритма ориентации (2). Получаемые при этом оценки вектора (1) на каждом шаге подставлялись в уравнения Эйлера (2).

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

На основе результатов обработки полетных данных вертолета:

1. Вычислялась ориентация ЛА по алгоритму комплексирования измерений приемника СНС с датчиками ДУС и ДДУ на скользящем интервале наблюдений.

2. Вычислялись углы крена и тангажа по способу адаптивной коррекции БИНС. Полученные данные подтвердили, что заявляемое изобретение обеспечивает повышение точности и непрерывности коррекции углов тангажа и крена курсовертикали БИНС в условиях маневрирования в полете.

Предложенный способ коррекции БИНС позволяет использовать датчики ДУС и ДЛУ средней и низкой точности, в том числе микромеханического типа, так как из-за непрерывной коррекции ошибки не накапливаются.