СПОСОБ ОПРЕДЕЛЕНИЯ ПОГРЕШНОСТЕЙ ОСНОВНЫХ ХАРАКТЕРИСТИК БЛОКА ИНЕРЦИАЛЬНЫХ ИЗМЕРИТЕЛЕЙ

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

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

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

Известен способ калибровки инерциального измерительного модуля по каналу акселерометров из патента РФ №2477864 с датой приоритета 30.11.2011 г., включающий оценку параметров математической модели погрешностей при задании различных ориентаций модуля относительно вектора ускорения свободного падения на грубом поворотном столе. При этом определяют нулевые сигналы и матрицу, описывающую относительное расположение измерительных осей, перекрестные связи, масштабные коэффициенты акселерометров, затем проводят привязку матрицы, описывающей относительное расположение измерительных осей, перекрестные связи, масштабные коэффициенты акселерометров к осям инерциального измерительного модуля, для чего точно определяют ориентацию модуля в двух различных положениях относительно вектора ускорения свободного падения, за исключением положений, получающихся путем разворота инерциального измерительного модуля на 180° вокруг любой оси.

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

В качестве прототипа принят способ калибровки бесплатформенных инерциальных навигационных систем (патент РФ №2406973 с датой приоритета 05.02.2009 г.), посредством которого возможно производить оценку основных характеристик измерительных каналов бесплатформенных инерциальных навигационных систем (далее - БИНС).

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

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

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

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

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

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

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

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

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

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

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

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

Решение данной задачи осуществляется в три этапа.

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

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

где ΔW_Ai - приращение интеграла проекции кажущегося ускорения на ось чувствительности i-го акселерометра за такт опроса (i=1, 2, 3, 4), м/с;

M_i - значение масштабного коэффициента i-го акселерометра, м/с;

ΔN_i- - приращение выходного кода i-го акселерометра, ед. кода;

τ_i - смещение нулевого сигнала i-го акселерометрического канала (i=1, 2, 3, 4), м/с²;

t_ИНТ - длительность интервала опроса, с.

На втором этапе формируются рабочие тройки некомпланарных векторов по показаниям наборов троек акселерометров (исходя из наличия четырех акселерометров).

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

С помощью математической модели погрешностей блока инерциальных измерителей для каждой рабочей тройки рассчитываются проекции в прямоугольную систему координат. По этим проекциям рассчитывается общее значение действующего ускорения на блок инерциальных измерителей БИНС по данным с рассматриваемой рабочей тройки. Переход от проекций на оси чувствительности акселерометров к проекциям на оси прямоугольной системы координат происходит с помощью матрицы перехода М_А (2), которая содержит углы, характеризующие погрешность ориентации оси чувствительности i-го акселерометра по отношению к ее номинальному положению (углы α_i=1..4, β_i=1..4):

где М_А - матрица, связывающая оси чувствительности акселерометров с осями прямоугольной системы координат;

ΔW_in - проекция кажущегося ускорения на оси прямоугольной системы координат (i=X, Y, Z).

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

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

где - модуль вектора кажущегося ускорения, действующего на прибор.

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

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

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

П_i - оцениваемые параметры (i=1, 2, …n).

Для решения системы (4) составляется функция рассогласования (5):

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

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

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

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