捷联惯导总结--初始对准,位置标定,INS姿态更新,GPS/INS组合

article/2025/9/29 3:18:18

惯导及组合导航回顾  2018.09.16

今天和17系的同学一起把惯导的流程捋了一遍,为了加深自己的记忆,这里在前面把心得大致列出来。

我们这里只考虑捷联式惯导及松组合

首先拿到惯性传感器(加速度计和陀螺仪)需要对其进行标定,是为了求出零偏,比例因子,交轴耦合这些误差。这些误差有的虽然在出厂时会给出,但是考虑到实际环境不同,还是要做标定。一般采用的是六位置旋转法,针对加速度计和陀螺仪同时进行标定。

其次需要求姿态转移矩阵C_{b}^{n},这里就是初始对准了。

这里针对激光陀螺,这种高精度陀螺仪可以敏感地球自转角速度以及重力加速度,所以可以用w_{b}^{ib}f_{sf}^{b}用双矢量定姿来求姿态阵。而MEMS这种低精度IMU,不能准确敏感w_{b}^{ib}f_{sf}^{b},则需要借助磁力计,卫导双天线定姿等来对准。

初始对准中,双矢量定姿方法得到粗略的姿态阵,这时通过静基座下输出的速度和失准角误差(估计的姿态阵与真实值之间存在的小量的数学平台失准角\o),利用误差做间接卡尔曼滤波,用姿态更新算法得到精对准姿态阵。

有了初始姿态阵,同时通过卫导得到初始位置和初始速度(多普勒),然后就是惯导的过程了。

由于惯导是主要使用积分,因此需要对每个惯导历元输出的结果进行积分。这里就是加速度积分得到速度,速度积分得到位置的过程了。这里位置是东北天坐标系下,后面和卫导松组合时要转换成经纬高(相对应位置在北和东方向的增量分别除子午圈和卯酉圈)。

在惯导和卫导松组合时,以惯导为主系统,惯导和卫导的姿态,位置,速度的差值作为状态量(更多维数时包括加速度,陀螺仪零偏,杆臂误差,时间不同步误差等),每个卫导历元输出一次组合结果。组合是通过卡尔曼滤波求出误差值,将其补偿到INS的结果中。

多传感器融合时,可以采用联邦滤波或者集中滤波,INS和GPS,高度计,磁力计分别组成误差方程进行卡尔曼滤波。这是采用同样的状态量。

0.惯导概述

导航定位原理分类

航位推算(Dead-Reckoning, DR)

 磁罗盘+里程计
 惯性导航

空间交汇测量(距离、 角度)
 经纬仪
 GPS
数据库匹配导航
 地标定位
 地形匹配

惯性传感器得到的所有观测量都是相对于惯性坐标系。
惯性传感器的基本原理是牛顿第二定律
(加速度计 & 陀螺仪, Accelerometer & Gyroscope)

惯性传感器——加速度计

 f = a – g (惯导比力方程)
 f =加速度计输出(比力, Specific Force)
 a =相对于惯性空间的运动加速度
 g =地球万有引力

惯性传感器 —— 陀螺仪

定义:一种用于测量相对于惯性参考系的角速率的传感器
机械陀螺 vs.光学陀螺
 机械陀螺 (转子陀螺,振动陀螺, …)
 光学陀螺: Sagnac 效应

  •  激光陀螺 (RLGs)
  •  光纤陀螺 (FOGs)

振动陀螺的原理——哥氏效应(Coriolis Effect)
光学陀螺的原理——萨格纳克效应 (Sagnac Effect): 反向传播光束之间的相移正比于旋转速度

ISA, IMU, INS

Inertial Sensor
Assembly (ISA)
3轴陀螺 + 3轴加速度计;
输出原始传感器数据.
Inertial
Measurement
Unit (IMU)
ISA经误差标定补偿(零偏,比例因子,etc)
和数据转换;
输出补偿后的数据.
Inertial
Navigation
System (INS)
IMU +惯性导航算法
(惯导机械编排);
输出位置、 速度、 姿态角.

惯性测量单元(IMU)= 3轴加速度计 + 3轴陀螺

平台式与捷联式对比

 Features of INS

 平台式捷联式
体积相对较大
重量
成本
性能可达最高精度最高到导航级
自标定能力
环境适应性对冲击和振动敏感抗冲击和振动

惯性导航系统特性

 误差随时间和运动距离累积
 需要初始信息
 成本高、 笨重
 完全自主性和高可靠性(军用和航空航天)
 导航信息完备、 连续

INS精度等级

 战略级
Strategic-Grade
导航级
Navigation-Grade
战术级
Tactical Grade
微机械级
MEMS
?
定位误差< 30 m/hr0.5 – 2 nmi/hr
(70-100k USD)
10-20 nmi/hr
(10-20K USD)
 
陀螺零偏0.0001 deg/hr大约地球自转的
1/1000, 0.015 deg/hr
(1 nmi/hr)
1- 10 deg/hr 
加速度计零偏1 ug50 – 100 ug100 – 1000 ug 
应用领域 洲际弹道导弹
 潜艇
通用航空
高精度测绘
 短时间应用(战
术导弹)
与GPS组合使用
 

陀螺零偏是关键

惯性器件的相对测量能力 

以导航级惯导的陀螺为例:陀螺零偏 0.01 deg/hr    陀螺量程 1000 deg/s = 3600,000 deg/hr
相对测量能力= 量程 / 零偏 = 3.6e8 !
 等效于测量 武汉—北京 距离, 精度3mm!

MEMS惯导的演化

 

传感器误差分类

零偏(Bias, 静态误差、 加性误差)

  • 包含噪声(Noise)

比例因子(Scale factor, 动态误差、 乘性误差)

  • 包含非线性(non-linearity)
  • 轴偏移/轴交叉(Axes misalignment / cross-axis)

基本的误差类型
 常值误差(Constant error)
 重复性(Repeatability) --multiple runs
 稳定性(Stability) -- within one run
 噪声和带宽(Noise & Bandwidth)
 热敏感度(Thermal sensitivity)


1.静态解析粗对准

 

2.加速度计位置标定


标定方法总结

  • 标定的精度依赖于各轴相对于参考坐标对准的准确性.
  • 为了获得准确的标定结果, 需要一些专业设备(如转台或规则的立方体)来获得IMU的精确姿态和旋转角.
  • 由于对专业设备的依赖性, 这些标定方法主要设计用于在实验室测试、 厂家校准和对相对较高精度的IMU的标定.

3.INS姿态解算

INS测量

在使用惯导系统做测量应用之前
1. 标定系统, 用来估计加速度计和陀螺仪的零偏和比例因子误差的确定性部分
2. 采集长时间的设备数据, 用来估计加速度计和陀螺的噪声特性(Gauss Marko参数)
3. 利用陀螺的随机游走参数估计初始对准的时长, 用来得到一定的航向角精度(be realistic)
 动态模式下进行测试来预测系统性能
1. 检查导航算法实现的正确性, 改善标定参数, 估计测量更新所需要的频率来获
得一定的测量精度.
2. 对于这个测试通常需要一个良好的参考轨迹 (通常是 GPS).
3. 测试通常是沿着一个L形的轨迹.
 

4.GPS/INS组合导航

松组合采用的是两套独立的系统,惯导系统与卫星接收机系统。其中利用惯导的误差方程作为系统方程,利用惯导输出的位置及速度  与  卫星接收机输出的位置及速度做差,作为系统的量测量。然后进行相应的卡尔曼滤波,利用卡尔曼滤波输出的姿态误差,位置误差以及速度误差对惯导输出的载体位置,速度,姿态进行更新。

惯导误差状态可用以下一阶状态方程表示

为估计和补偿惯导误差, 需要更高精度的外部测量值

假设可获得持续的GPS位置观测值作为量测更新: 观测方程可通过辅助信息PGPS(GPS位置观测值) 与INS位置输出PINS求差得到;用于卡尔曼滤波更新的观测向量Z表示为:

观测向量Z可写成状态量x的函数

如果速度和位置观测值均可用, 则

——————————————————————————————————————————————————————

21维

摘自武汉大学牛小骥老师讲义

GNSS/INS松组合常采用误差状态卡尔曼滤波(间接卡尔曼滤波)进行数据融合, 以解决系统的非线性问题。 根据惯导误差微分方程和传感器误差模型, 将惯性传感器主要误差参数(加速度计和陀螺的零偏及比例因子误差) 增广到卡尔曼滤波的系统状态中, 可得卡尔曼滤波的状态向量及连续时间系统状态方程

________________________________________________________________________________________________

17维

摘自西工大严恭敏老师《捷联惯导算法与卡尔曼滤波原理讲义》


http://chatgpt.dhexx.cn/article/AiFRZyYM.shtml

相关文章

捷联式惯导系统初始对准

1 初始对准简介 1.1 初始对准任务 所谓对准指的是确定惯导系统各坐标轴相对于参考坐标系指向的过程。 捷联式惯导系统初始对准的任务有两项:第一&#xff0c;机体起飞前将初始速度和初始位置引人惯导系统;第二&#xff0c;机体坐标系与导航坐标系的初始变换短阵。另外&#…

c语言函数指针 的定义方法,C语言 函数指针一(函数指针的定义)

//函数指针 #include #include #include //函数指针类型跟数组类型非常相似 //函数名就是函数的地址&#xff0c;函数的指针&#xff0c;对函数名进行&取地址操作&#xff0c;还是函数名本身&#xff0c;这是C语言编译器的特殊处理 void test(int a){ printf("a%d\n&q…

C语言函数指针调用函数

C语言可以定义指针变量指向函数&#xff0c;通过指针来调用函数。 使用&#xff1a; 1、定义函数指针变量&#xff1a; 函数返回值类型 (*指针变量名)(); 2、将已有函数入口赋值给函数指针变量&#xff1a;fnPt avg; 3、使用指针变量调用函数&#xff1a;(*指针变量名)(参…

C语言函数指针与调用

C语言函数指针与调用 平时我们使用函数传递的参数一般为数据变量&#xff0c;那么是否可以传递函数呢&#xff1f; 答案是不但可以&#xff0c;而且习惯以后&#xff0c;会用上瘾的。通过传递不同的函数指针&#xff0c;我们可以实现在函数中调用不同的子函数。 下面就举个栗子…

C语言函数指针与指针函数

在大家刚开始学习c语言的时候&#xff0c;总是分不清函数指针和指针函数&#xff0c;就算是知道了它们之间的区别&#xff0c;也不了解它们的使用场景&#xff0c; 我写此博客帮大家缕一缕&#xff0c;也帮我自己缕一缕 1、函数指针与指针函数的概念以及区别 指针函数 从名…

C语言 函数指针做函数参数(即回调函数)

文章目录 函数指针做函数参数(回调函数)回调函数概念一般有三种调用方式回调函数的作用回调函数调用时刻回调函数的语法&#xff1a;1.简单的函数类型为&#xff1a;无参数、无返回值的函数。2.完全形式的回调函数注&#xff1a; 代码案例&#xff1a;模拟计算器 函数指针做函数…

C语言 函数指针

一.前言 最近学了一点函数指针的东西 感觉还是比较有意思的 在某先方面用起来也确实会方便些 给大家在这里分享一下。 二.定义 函数指针就是指向代码段中函数入口地址的指针。 从上述这句话就可以看出 函数指针的本质是一个指针&#xff0c;只不过是指向函数的指针。 三.声明形…

C语言函数指针用法

C语言函数指针用法 函数指针本质上是指针&#xff0c;它指向一个函数。 例如int (*p)(int); – 从 p 处开始, 先与指针结合, 说明 p 是一个指针, 然后与()结合, 说明指针指向的是一个函数, 然后再与()里的 int 结合, 说明函数有一个int 型的参数, 再与最外层的 int 结合, 说明…

C语言函数指针详解

1、概念 函数指针就是一个指针&#xff0c;指针指向某个内存区域&#xff0c;函数指针就是指向函数入口地址的这么一个指针变量&#xff0c;在.c文件中编写一个函数&#xff0c;将.c编译为可执行程序后&#xff0c;在.c文件中编写的函数会存放在可执行程序的代码段中&#xff…

C语言——函数指针

函数指针 C语言中的指针类型有很多种&#xff0c;但是函数指针可能有些小伙伴没听过&#xff0c;下面我们引入一段定义。 函数指针是指向函数的指针变量。 因此“函数指针”本身首先应是指针变量&#xff0c;只不过该指针变量指向函数。这正如用指针变量可指向整型变量、字符型…

《信号与系统》小结

参考书目&#xff1a;《信号与系统(第四版) 徐亚宁 苏启常 编著》 信号与系统 信号是消息的表现形式&#xff0c;消息是信号的具体内容。也就是说任何可以承载某种消息的物理量都可以是信号。电信号只是生活中最快捷&#xff0c;也是最常用的一种。我们可以从变化的电流或电压…

信号与系统(一) 信号与系统的基本概念

1.1信号的描述 信号的描述方式采用的有3种——函数、图形和数据 1.2信号的大小 对于一个信号x(t)&#xff0c;它的能量大小为&#xff0c;但是当&#xff0c;而信号的强度不趋向于0时&#xff0c;信号能量为无限大&#xff0c;这时候就需要使用信号的功率来判断大小&#xf…

连续系统分析【信号与系统四】

连续系统分析 一、已知描述连续系统的微分方程&#xff0c;计算该系统的响应并与理论结果比较二、研究具有以下零极点的连续系统2.1 1个极点s-0.1&#xff0c;增益k12.2 1个极点s0&#xff0c;增益k12.3 2个共轭极点 &#xff0c;增益k12.4 2个共轭极点 &#xff0c;增益k12.5 …

信号与系统_第二章 连续系统的时域分析

第二章 连续系统的时域分析 ( 连续系统在时域上进行分析 ) 2.1 引言 2.2 LTI系统的微分方程表示及响应 2.3 零输入响应与零状态响应 2.4 单位冲激响应 2.5 卷积积分 2.6 卷积积分的性质 2.7 奇异函数 2.8 连续时间系统的模拟 习题 2.01_连续系统的描述: 电路图建立微分…

《信号与系统》解读 第1章 信号与系统概述-1:信号与系统的描述和分析方法

目录 1. 什么是广义的信号和系统 2. 什么是狭义的电信号和电系统 3. 两个研究对象&#xff1a;信号系统 4. 两种分析方法&#xff1a;几何图形法数学函数法 5. 两个维度与视角&#xff1a;时域频域 6. 两种坐标&#xff1a;实平面复平面 7. 目标&#xff1a;模拟信号处理…

【信号与系统学习笔记 3】—— 系统,以及系统的性质以及判断方法解析

文章目录 一、系统的分类二、系统的性质2.1 记忆性和无记忆性2.2 可逆性与可逆系统2.3 稳定性2.4 因果性2.5 时不变性2.7 线性性 一、系统的分类 其实&#xff0c;系统的分类无非就是两大类&#xff1a;离散时间系统、一个是连续时间系统 那么&#xff0c;如何判断呢&#xf…

信号与系统知识点记录(P2-P3)

P2: 连续余弦信号和离散余弦信号的不同。 周期不同。 时移和相位关系不同。 还有个有趣的性质&#xff0c;离散信号里面&#xff0c;频率相差2Π的整数倍&#xff0c;得到的离散信号是一样的。 后面又实数的指数信号和复数的指数信号的推导过程 离散脉冲函数在信号和系统里…

信号与系统_第2章 连续系统的时域分析

2.1 LTI连续系统的响应 前言 LTI连续系统的时域分析&#xff0c;归结为&#xff1a;建立并求解线性微分方程。 由于在其分析过程涉及的函数变量均为时间t&#xff0c;故称为时域分析法。这种方法比较直观&#xff0c;物理概念清楚&#xff0c;是学习各种变换域分析法的…

数字信号处理:重要知识点整理

文章目录 0 最重要&#xff01;DFT和FFT的区别1 连续时间信号频域分析2 通过离散时间信号的Z变换表达式X(z)直接写出时域离散信号&#xff08;序列&#xff09;x(n)的方法3 部分分式法的MATLAB实现&#xff08;求X(z)的部分展开式&#xff09;4 稳定系统5 求频响特性&#xff0…