捷联惯导系统粗对准方法比较
捷联惯导系统粗对准方法比较
魏春岭 张洪钺
北京航空航天大学自动化科学与电气工程学院 北京 100083
摘 要 通过误差分析对三种捷联惯导系统解析粗对准方法进行了比较。指出在
相同的传感器精度条件下,利用正交向量计算捷联矩阵比传统方法有更高的对准
精度,直接计算法不仅精度高,而且计算简单,更适合工程应用。
主题词 捷联惯导系统 解析粗对准
Comparison of Analytic Coarse Alignment Methods
Wei Chunling Zhang Hongyue
Beijing University of Aeronautics and Astronautics,Beijing100083
Abstract Three analytic coarse alignment methods to strapdo wn inertial navigation system
are com pared via error analysis.The later two are superior to the traditional one because their
east level dri f t misalignment angles are not corrupted b y gyro uncertainty.Due to its high ac-
curacy and com putation e ff iciency,the direct method is more suitable for practical applica-
tions.
Subject terms Strapdown inertial navigation systems Analytic coarse alignment
作为一种航迹推算系统,惯性导航系统对初始解算条件有较高要求,初始对准误差会直接影响导航的精度。对于捷联式惯性导航系统,初始对准的目的就是要确定捷联矩阵C n b。解析粗对准就是利用加速度计和陀螺仪对重力加速度和地球自转角速度的测量值估算C n b,为精对准提供初始条件,因此选择算法简单、精度更高的粗对准方法有其实际意义。本文通过误差分析与计算机仿真比较了三种解析粗对准方法,指出直接计算法更适合工程应用。
1 解析粗对准方法
假定当地纬度 已知,地理系采用东北天坐标系,则重力加速度g和地球自转角速度
收稿日期 1999年12月
16
ie在地理系中的分量都是确定已知的,可以表示为
g n=[0 0 -g]T(1)
n ie=[0 cos sin ]T(2)其中 为地球自转角速率,g为重力加速度的幅值。
解析粗对准就是要利用已知g n, n ie和对g b, b i e的测量值来估计C n b。由于粗对准时间较短,通常采用对传感器测量值取平均的办法来提高对准精度。为表达方便,用a b和 b ie 分别表示加速度计和陀螺仪的测量平均值。
1 1 利用g、 ie和构造的g ie,计算C n b
为了直接解出正交阵C n b中的所有元素,需要构造新的向量来增加方程的数目,传统方法是构造辅助向量g ie。根据坐标转换关系,g和 ie在机体系中的投影可表示成
[g b b ie g b b ie]=C b n[g n n ie g n n ie](3) 考虑传感器的测量偏差和外界干扰,加速度计和陀螺仪的测量值可分别写成
a b=g b+ a b(4)
b= b ie+ b(5)其中, a b表示加速度计的测量误差, b表示陀螺仪的测量误差。则C n b的估计值为[1]
C^n b=(g n)T
( n ie)T
(g n n ie)T
-1
(a b)T
( b)T
(a b b)T
(6)
由于传感器测量误差的存在,使得C^n b不满足正交性要求,可按下式进行正交化[1]
(C^n b)o=C^n b[(C^n b)T C^n b]-1/2(7)其中,(C^n b)0表示正交化后的估计值,上式使得[(C^n b)o-C^n b]T[(C^n b)o-C^n b]的迹最小。
1 2 利用三个相互正交的向量g,g ie和(g ie) g,计算C n b
为了求解C n b,利用三个相互正交的向量a b,a b b和(a b b) a b计算C^n b,公式如下[2]
C^n b=(g n)T
(g n n ie)T
[(g n n ie) g n]T
-1
(a b)T
(a b b)T
[(a b b) a b]T
(8)
采用(7)式对其进行正交化处理。
1 3 由测量值直接计算C n b
为了表达方便,将导航坐标系到机体坐标系的转移矩阵C n b写成
C n b=[x y z](9)其中x,y,z为构成C b n的三个列向量,且彼此间具有正交约束。
考虑测量误差,k时刻加速度计和陀螺仪的测量值可以表示成
a b(k)=-gz+ a b(k)(10)
b(k)= cos y+ sin z+ b(k)(11) 采用正交约束下的最小二乘方法估计C b n的三个分向量[3]
17
^z=-1
N
k=1
a b(k)(12)
y^=1
N
k=1
b(k)+ ^z(13)
x^=y^ ^z=1
N
k=1
b(k) ^z(14)
其中,N为总的测量次数,
= N k=1a b(k)T N k=1a b(k)1/2(15)
=1 N
k=1
a b(k)T N k=1 b(k)(16)
= N
k=1 b(k)T N
k=1
b(k)- 21/2(17)
则捷联矩阵的估计值可以写成
C^n b=[x^ y^ ^z]T(18) 按照式(12)~(17)计算得到的C^n b已经是正交阵,不需要再做正交化处理。
2 误差分析
分析各种对准方法的误差特性对于选择更好的对准方案有着实用意义。Britting给出了方法1的误差模型[1],由于传感器误差采用的仍然是机体系中的表示,不利于物理解释;Jiang比较了前两种方法正交化前的误差特性[2],指出方法2具有更高的对准精度。这里采用Britting的方法分析传感器测量误差对三种粗对准方法精度的影响,为了便于和平台系统进行比较,将误差模型写成地理系中的表达形式。
由于传感器测量误差的影响,正交化前的估计值C^n b主要含有刻度系数误差、歪斜误差和漂移误差三项[4,5]。通过正交化可以消除前两项误差,因此,可以用漂移误差组成的反对称阵来描述正交化后的(C^n b)0与理想捷联阵C n b间的关系,即
(C^n b)0=[I- ]C n b(19)其中, 为漂移误差角 =[ E, N, U]T构成的反对称阵。
将C^n b的计算公式(6)和(8)写成如下形式
C^n b=MQ^=M(Q+ Q)(20)其中,M对应于和g n, n ie有关的矩阵,Q^对应于和测量a b, b有关的矩阵, Q R3 3为加速度计和陀螺仪测量误差构成的矩阵。将(7)式中的平方根项按级数展开,有近似的正交化公式
(C^n b)0=[I+1
2
(M QC b n-C n b Q T M T)]C n b(21)
与(19)式进行比较,可以得到由漂移误差角构成的反对称阵为
=1
2C n b Q T M T-M QC b c=1
2
C n b Q T M T-(C n b Q T M T)T(22)
18
2 1 方法1的误差分析
对于方法1,计算捷联阵所采用的三个向量并不满足正交性
M=
001/g cos
tan /g1/ cos 0
-1/g00
(23)
Q T=[ a b b (a b b)](24)定义传感器误差在地理系中的投影为
a n=C n
b a b=[ a E a N a U]T(25)
n=C n b b=[ E N U]T(26)忽略二阶小量,可得
C n b Q T=[ a n n ( a n n ie+g n n)](27)将(23)和(27)代入(22),可以得到漂移误差角的表达式为
E=1
2 a N
g+
a U
g tan +
U
cos (28) N=-
a E
g
(29)
U=- a E
g
tg -
E
cos (30)
可见:方位误差角与平台罗经对准相似,主要取决于东向陀螺漂移,所不同的是东向漂移误差角不仅与北向加速度计误差有关,还与天向加速度计误差和天向陀螺漂移有关,而且主要取决于天向陀螺漂移的大小。
2 2 方法2的误差分析
对于方法2,采用三个相互正交的向量计算C n b,M具有更简单的表达式
M=
01/g cos 0
001/g2 cos
-1/g00
(31)
Q T= a b (a b b) [(a b b) a b](32)忽略二阶小量,有
C n b Q T=[ a n ( a n n ie+g n n)( a n n ie) g n+(g n n)
g n+(g n n ie) a n)](33)将(31)和(33)代入(22),得到漂移误差角的表达式为
E= a N
g(34)
N=- a E
g
(35)
U=- a E
g
tg -
E
(36)
比较(34)和(28),此时东向漂移误差角只和北向加速度计误差有关,而与天向陀螺和加速度计的测量误差无关,这与自修正精对准方法得到的漂移误差角的表达式相
19
同[1,6]。可见采用方法2进行解析粗对准,比方法1具有更高的水平对准精度。
2 3 直接计算法与方法2的等价性
我们通过证明直接计算法与方法2表达式的等价性来分析它的误差特性。将(31)代入(8)有
C^n b=
01/g cos 0
001/g2 cos
-1/g00
(a b)T
(a b b)T
[(a b b) a b]T
(37)
仍然将C b n写成(9)的形式,由上式可以得到对三个列向量的估计分别为
x^ =a b b
g cos
(38)
y^ =(a b b) a b
g2 cos
(39) ^z =-
a b
g(40)
显然x^ ,y^ ,^z 两两正交,按(7)式对C^b n进行正交化,等于对三个列向量进行归一化,
即x^0=a b b
a b b
(41)
y^0=
(a b b) a b
(a b b) a b
(42) ^z0=-a
b
a b
(43)
其中,| |表示向量的模,y^0=^z0 x^0,捷联矩阵由(C^b n)0的转置获得。为了便于比较,可将(12)和(14)用测量的时间均值表示,有
^z=-a b
a b
(44)
x^=1
(
b ^z)=(a b b)
a b b
(45)
可见:x^=x^0,^z=^z0且y^=^z x^=y^0,因此,直接计算法和方法2的表达式是等价的,即误差特性相同。
3 仿真验证
仿真条件:纬度 =40 ;方位角 =-45 ,俯仰角 =30 ,倾斜角 =10 ;陀螺漂移 =0 1( )/h,加速度计零偏 =100 g;测量噪声强度 =0 01( )/h, =10 g。对准时间为1min,采样时间间隔为100ms。
按三种方法进行粗对准得到的漂移误差角分别为(单位:( ))
1=[18 1862 -0 3965 -34 4145]T
2=[-0 1683 -0 3965 -34 4145]T
3=[-0 1683 -0 3965 -34 4145]T
由仿真结果可以看出在相同的传感器精度条件下,后两种方法比传统方法具有更高的
20
水平对准精度。仿真结果与前面通过误差分析得到的结论是一致的。
4 结 论
通过上面的误差分析和计算机仿真验证,可见采用g, ie和g ie计算C n b时,东向漂移误差角不仅和北向加速度计误差有关,而且与天向加速度计误差和天向陀螺漂移有关, E主要取决于天向陀螺漂移。采用三个相互正交的向量g,g ie,(g ie) g 计算C n b和直接计算法是等价的,且具有更高的对准精度,此时 E只和北向加速度计误差有关。直接计算法因避免了三角函数运算以及正交化过程,公式简单,更适合于工程应用。
参 考 文 献
1 Bri tting,K R Inertial Navigation Systems Analysis Wiley Interscience,New York,1971
2 Jiang,Yeon-Fuh Error Analysis of Analytic Coarse Ali gnment Methods IEEE Trans Aerospace and Electronic
Systems,1998,34(1),334~337
3 Schi melevich,L and R Naor,New Approach to Coarse Alignmen t,PLANS 94,324~327
4 Mortensen,R E ,Strapdown Guidance Error Analysis,IEEE Trans Aerospace and Electronic Systems,1974,
10(4),451~457
5 陈哲,捷联惯导系统原理,宇航出版社,1986
6 袁信、郑谔,捷联式惯性导航原理,航空专业教材编审组,1985
上接第15页
5 结论
输入指令整形器以前馈的形式作用于转台伺服系统时,既不影响系统的稳定性,又可以很好地抑制频带拓宽时的机械谐振,而且对系统参数变化具有较高的鲁棒性。本文给出的前馈整形器与高开环增益的反馈控制相结合的控制算法,其算法简单,实现和调试也很容易,在电机功率足够的情况下,可以作为拓宽飞行仿真转台频带的一种有效的新途径。
参 考 文 献
1 N. C.Singer,W.P.seering.Preshaping Command Inputs to Reduce System Vibrati on.Tansactions of the AS ME,
1992,112:76~81
2 李智铭 三轴飞行仿真转台频带拓宽技术研究 [学位论文]中国运载火箭技术研究院,1999,4
3 Brett R.Murphy Ichiro Watanabe.Digital Shaping Filters for Reducing Machine Vibration.IEEE Transactions on
Robotics and Automation,1992,8:285~289
21
北航卡尔曼滤波课程-捷联惯导静基座初始对准实验
卡尔曼滤波实验报告 捷联惯导静基座初始对准实验 一、实验目的 ①掌握捷联惯导的构成和基本工作原理； ②掌握捷联惯导静基座对准的基本工作原理； ③了解捷联惯导静基座对准时的每个系统状态的可观测性； ④了解双位置对准时系统状态的可观测性的变化。 二、实验原理 选取状态变量为：[]T E N E N U x y x y z X V V δδεεε=ψψψ??，其
中导航坐标系选为东北天坐标系，E V δ为东向速度误差，N V δ为北向速度误差，E ψ为东向姿态误差角，N ψ为北向姿态误差角，U ψ为天向姿态误差角，x ?为东向加速度偏置，y ?为北向加速度偏置，x ε为东向陀螺漂移，y ε为北向陀螺漂移，z ε为天向陀螺漂移。则系统的状态模型为： X AX W =+ (1) 其中， 1112212211 12 1321222331323302sin 000002sin 000000000sin cos 0000sin 000000cos 0000000000000000000000000000000000000000000000000000 0L g C C L g C C L L C C C L C C C L C C C A Ω-? ? ??-Ω????Ω-Ω? ?-Ω????Ω=? ?????? ?????????? ? [00000]E N E N U T V V W W W W W W δδψψψ=，E D V W W δψ 为零均值高斯 白噪声，分别为加速度计误差和陀螺漂移的噪声成分，Ω为地球自转角速度，ij C 为姿态矩 阵n b C 中的元素，L 为当地纬度。 量测量选取两个水平速度误差：[ ]T E N Z V V δδ=，则量测方程为： 10000000000100000000E E N N V X V δηδη???? ??=+???????????? (2) 即Z HX η=+ 其中，H 为量测矩阵，[]T E N ηηη=为量测方程的随机噪声状态矢量，为零均值高 斯白噪声。 要利用基本卡尔曼滤波方程进行状态估计，需要将状态方程和量测方程进行离散化。 系统转移矩阵为： 2323/1111102!3!! n n k k k k k k n T T T I TA A A A n ∞ -----=Φ=++++=∑ (3)
第六章 捷联惯导
第六章捷联惯导
6-1捷联惯导的原理?捷联惯导系统概述 ?捷联惯性技术的发展过程 ?捷联惯导系统与平台惯导系统的对比 ?捷联惯导系统的基本力学编排方程?捷联惯导系统的算法概述 ?捷联惯导系统原理框图的说明 ?姿态方程的解算 （1）姿态和航向角的计算 （2）姿态矩阵的微分方程 （3）四元数的运动学微分方程 （4）等效旋转矢量法及其微分方程 （5）位移角速率方程 （6）速度方程
?导航位置方程 （1）游动方位系与地球系之间的方向余弦矩阵 （2）载体位置计算 （3）方向余弦矩阵计算 ?垂直通道阻尼 ?捷联惯性器件的余度技术?单自由度陀螺仪的配置方案 （1）四陀螺仪配置方案 （2）六陀螺仪系统 ?二自由度陀螺仪的配置方案
?捷联惯导的数值计算方法?数值积分法 （1）欧拉法 （2）四阶龙格－库塔法 ?角速率信息的提取
“ 捷联（Strapdown）”这一术语的英文原义就是“捆绑”的意思。因此，所谓捷联惯性系统也就是将惯性敏感元件（陀螺与加速度计）直接“捆绑”在载体上，从而完成制导和导航任务的系统。 V-2导弹 “阿波罗-13”宇宙飞船 “海盗”火星降落器
从捷联技术的发展过程中我们已经看到捷联系统的优越性已越来越突出的显示出来，并在许多方面已日渐代替平台系统。为什么会出现这种情况呢？为了回答这一问题，这里从生产与使用的角度将捷联系统与平台系统做一对比。 （1）硬件和软件的复杂程度 由于捷联系统没有平台框架及相连的伺服装置，因而简化了硬件；代价是增加了计算机的负担，需要一个比较复杂的实时程序。 （2）可靠性 捷联系统的可靠性要比平台系统高，其原因是它的机械构件少，加之容易采用多敏感元件配置，实现余度技术。 （3）成本与可维护性 由于平台系统在机械结构上要复杂得多，而对于捷联系统只是算法复杂些，因而从制造成本上看捷联系统的成本要比平台系统低。从市场供应的情况来看，数字计算机的价格一直在下降，而平台系统的价格一直在上升。 此外，捷联系统比平台系统具有较长的平均故障间隔时间，加之模块设计简化了维修，从而捷联系统的可维护性比平台系统大为提高了。
捷联惯导姿态算法中的圆锥误差与量化误差
第27卷 第1期航 空 学 报 Vol 127No 11 2006年 1月ACTA A ERONAU TICA ET ASTRONAU TICA SIN ICA J an. 2006 收稿日期:2004209230;修订日期:2005204227基金项目:国家自然科学基金(60234030)、国家杰出青年科学基金 (60225015)和教育部高校青年教师奖资助项目 文章编号:100026893(2005)0120098206 捷联惯导姿态算法中的圆锥误差与量化误差 练军想,胡德文,胡小平,吴文启 (国防科技大学机电工程与自动化学院自动控制系,湖南长沙 410073) R esearch on Coning E rror and Q uantization E rror of SINS Attitude Algorithm L IAN J un 2xiang ,HU De 2wen ,HU Xiao 2ping ,WU Wen 2qi (Department of Automatic Control ,College of Mechatronics and Automation ,National University of Denfense Technology ,Changsha 410073,China ) 摘 要:对捷联惯导系统的误差源进行了研究,利用几何方法分析了不可交换性误差和量化误差的形成机理,以及它们的相互影响。针对工程应用中激光陀螺输出脉冲采样量化条件,就多子样算法进行了讨论,并设计了基于MA TL AB/Simulink 的仿真。研究结果表明,当考虑量化误差的影响时,选取适当的量化因子,三子样等效旋转矢量算法比其它算法具有更好的综合性能。 关键词:捷联惯导;姿态算法;不可交换性误差;量化误差;多子样中图分类号:V249.3 文献标识码:A Abstract :In this paper ,the error sources of strapdown inertial navigation system (SINS )are researched.The mechanism of noncommutativity error and quantization error is illustrated f rom the geometrical point of view ,and the mutual influence between them is analyzed.Considering the sample condition of the output pulses of the ring laser gyroscope (RL G ),the multi 2sample attitude algorithm of SINS is discussed.The simu 2lation is carried out using MA TL AB/Simulink.It is concluded that the three 2sample rotation vector algorithm with the proper quantization factor outperforms others when the impact of quantization errors is taken into ac 2count. K ey w ords :SINS ;attitude algorithm ;noncommutativity error ;quantization error ;multi 2sample 对于捷联惯导系统姿态更新算法的研究,目前国内外同行们有很多成果。普遍的观点认 为[1,2]:四元数姿态更新算法比方向余弦和欧拉角等其它方法的计算量小、适应性好;双子样、三子样等效旋转矢量算法在抑制不可交换性误差方面比单子样有显著改善;将姿态更新分成快慢两个回路,即快速回路更新旋转矢量,慢速回路更新姿态参数,既能有效抑制不可交换性误差,还可以减少算法的计算量。但上述观点很少有结合捷联惯导系统工程实际来考虑的。本文针对工程中激光陀螺输出脉冲采样量化条件,就多子样算法展开研究,并得出了相应的结论。 1 误差源分析 评价算法优劣的最重要的标准之一是算法的 精度。要提高算法的精度就必须研究各种误差及其补偿、抑制方法。 在对惯导系统作原理性分析时,尚可以理想化地认为平台坐标系(平台式或捷联式)准确、无误差地模拟地理坐标系,但实际的惯导系统不可避免地受到各种误差的影响[3]。通常,可将误差分为源信息误差和解算方法误差。对捷联惯导系统,解算方法误差是指数值求解姿态矩阵和比力方程时产生的原理性误差,比如不可交换性误差、舍入误差、截断误差等。源信息误差可分为解析源信息误差和测量源信息误差。测量源信息误差就是惯性传感器的误差,此类误差因惯性器件的物理原理和仪表的精度不同而各不相同。比如陀螺的常值漂移和随机漂移、刻度因子误差、非线性误差、输出值的量化误差等。 下面重点对捷联姿态解算时不可交换性误差和陀螺采样的量化误差进行分析。111 不可交换性误差 在利用方向余弦矩阵微分方程或四元数微分方程进行姿态更新时,都会遇到角速度矢量的积分问题。由于刚体有限转动的不可交换性,当转
捷联惯导系统初始对准技术的研究
捷联惯导系统初始对准技术的研究 摘要:初始对准是捷联惯导系统关键技术之一。初始对准精度直接影响捷联惯导系统的工作精度，初始对准时间也是反映武器系统快速反应能力的重要战术指标。 捷联惯性导航系统是将惯性器件陀螺仪、加速度计构成的惯性测量单元直接与载体固联，测量得到的载体角速度与线运动参数是沿载体固联的坐标轴上的分量。导航计算机通过计算“姿态矩阵”可以将加速度信息转换到惯性坐标系或当地地理坐标系，从而实现“数学平台”，然后再进行速度及位置计算。图1即为捷联式惯性导航系统原理框图。 捷联惯导系统的关键技术包括初始对准问题、有害加速度的消除及引力修正、惯性元件误差模型的建立和实时补偿、捷联矩阵的更新等。捷联惯性导航系统初始对准的目的是建立捷联矩阵的初始值。 1、捷联惯导系统初始对准基本概念 按对准阶段来分，初始对准一般分为两个阶段:第一阶段为粗对准，第二阶段为精对准。捷联系统粗对准的任务是得到粗略的捷联矩阵，为后续的精对准提供基础，此阶段精度可以低一些，但要求速度快。精对准是在粗对准的基础上进行的，通过处理惯性敏感元件的输出信息，精确校正真实导航坐标系与计算的导航坐标系之间的失准角，使之趋于零，从而得到精确的捷联矩阵。按照捷联惯性导航系统初始对准时载体的运行状态来分，可分为静基座对准和动基座对准。按
照初始对准时是否取得外部信息，可分为自对准和非自对准.惯性导航系统的自对准是利用重力矢量和地球自转角速率矢量通过解析的方法实现的初始对准，这种对准方法的优点是自主性强，缺点是所需的对准时间长。非自主式对准可以通过机电或光学方法将外部参考坐标系引入系统，实现惯性系统的初始对准.在捷联惯性导航系统的粗对准阶段，可以通过引入主惯导系统的航向姿态信息，通过传递对准，迅速将数学平台对准导航坐标系，减小初始失准角.在精对准阶段，可以通过组合导航的方法，利用其它导航设备(如GPS，计程仪)等提供的信息(如速度和位置)作为观测信息，通过卡尔曼滤波实现精确对准。目前有关初始对准问题的研究主要集中在误差模型的建立、模型求解方法和误差模型的可观性分析三个方面，本文正是针对这三个方面，对初始对准的国内外研究状况进行综述。 2、初始对准的误差模型 捷联惯导系统初始对准的误差模型及常用算法研究的基础模型有Ψ角误差模型和Φ角误差模型。前者基于导航坐标系，是指导航坐标系到计算坐标系(通常取计算的地理坐标系)的小角度误差；后者基于真实坐标系，是指导航坐标系到真实坐标系(通常取地理坐标系)的小角度误差。 3、求解误差模型的方法 3. 1古典方法 从频率域角度设计对准回路，并将对准分为水平对准和方位对准，对准过程首先是水平粗对准，然后是方位粗对准;在粗对准之后再精对准，首先是水平精对准，然后是方位精对准.由于对准回路频带低，响应慢，因此整个对准时间长。针对在静基座捷联惯导系统初始对准中，东向陀螺漂移没有估计效果，使得估计的方位失准角存在常值误差的问题，上海交大提出了，建立东向陀螺漂移估计位的修正方程，并对其进行修正，从而大大提高了方位失准角的估计精度。 3.2 Kalman滤波在初始对准中的应用 提高惯性导航系统初始对准精度的最佳途径之一是利用Kalman滤波这一重要数学工具，对于采用自主对准方式的惯导系统，一般采用卡尔曼滤波技术估计出系统的失调角和惯导系统误差源，然后采用一定的控制技术设计出控制角速率使失调角达到规定的要求。卡尔曼滤波的对象是用状态方程来描述随机线性系统，它按照估计误差方差最小的原则，从被污染的观测值中，实时估计出系统的各个
捷联惯导算法心得
捷联惯导算法心得 1、四个概念：“地理”坐标系、“机体”坐标系、他们之间换算公式、换算公式用的系数。 地理坐标系：东、北、天，以下简称地理。在这个坐标系里有重力永远是（0,0,1g），地磁永远是（0,1,x）（地磁的垂直不关心）两个三维向量。 机体坐标系：以下简称机体，上面有陀螺、加计、电子罗盘传感器，三个三维向量。 换算公式：以下简称公式，公式就是描述机体姿态的表达方法，一般都是用以地理为基准，从地理换算到机体的公式，有四元数、欧拉角、方向余弦矩阵。 换算公式的系数：以下简称系数，四元数的q0123、欧拉角的ROLL/PITCH/YAW、余弦矩阵的9个数。系数就是描述机体姿态的表达方法的具体数值。 姿态，其实就是公式+系数的组合，一般经常用人容易理解的公式“欧拉角”表示，系数就是横滚xx度俯仰xx度航向xx度。 2、五个数据源：重力、地磁、陀螺、加计、电子罗盘，前两个来自地理，后三个来自机体。 3、陀螺向量：基于机体，也在机体上积分，因为地理上无参考数据源，所以很独立，直接在公式的老系数上积分，得到新系数。 狭义上的捷联惯导算法，就是指这个陀螺积分公式，也分为欧拉角、方向余弦矩阵、四元数，他们的积分算法有增量法、数值积分法（X阶龙格-库塔）等等 4、加计向量、重力向量：加计基于机体，重力基于地理，重力向量（0,0,1g）用公式换算到机体，与机体的加计向量算出误差。理论上应该没有误差，这误差逆向思维一下，其实就是换算公式的系数误差。所以这误差可用于纠正公式的系数（横滚、俯仰），也就是姿态。 5、电子罗盘向量、地磁向量：同上，只不过要砍掉地理上的垂直向量，因为无用。只留下地理水平面上的向量。误差可以用来纠正公式的系数（航向）。 6、就这样，系数不停地被陀螺积分更新，也不停地被误差修正，它和公式所代表的姿态也在不断更新。 如果积分和修正用四元数算法（因为运算量较少、无奇点误差），最后用欧拉角输出控制PID（因为角度比较直观），那就需要有个四元数系数到欧拉角系数的转换。常用的三种公式，它们之间都有转换算法。 再搞个直白一点的例子： 机体好似一条船，地理就是那地图，姿态就是航向（船头在地图上的方位），重力和地磁是地图上的灯塔，陀螺/积分公式是舵手，加计和电子罗盘是瞭望手。舵手负责估计和把稳航向，他相信自己，本来船向北开的，就一定会一直往北开，
捷联惯导算法心得重点讲义资料
1、四个概念：“地理”坐标系、“机体”坐标系、他们之间换算公式、换算公式用的系数。 地理坐标系：东、北、天，以下简称地理。在这个坐标系里有重力永远是（0,0,1g），地磁永远是（0,1,x）（地磁的垂直不关心）两个三维向量。 机体坐标系：以下简称机体，上面有陀螺、加计、电子罗盘传感器，三个三维向量。 换算公式：以下简称公式，公式就是描述机体姿态的表达方法，一般都是用以地理为基准，从地理换算到机体的公式，有四元数、欧拉角、方向余弦矩阵。 换算公式的系数：以下简称系数，四元数的q0123、欧拉角的ROLL/PITCH/YAW、余弦矩阵的9个数。系数就是描述机体姿态的表达方法的具体数值。 姿态，其实就是公式+系数的组合，一般经常用人容易理解的公式“欧拉角”表示，系数就是横滚xx度俯仰xx度航向xx度。 2、五个数据源：重力、地磁、陀螺、加计、电子罗盘，前两个来自地理，后三个来自机体。 3、陀螺向量：基于机体，也在机体上积分，因为地理上无参考数据源，所以很独立，直接在公式的老系数上积分，得到新系数。 狭义上的捷联惯导算法，就是指这个陀螺积分公式，也分为欧拉角、方向余弦矩阵、四元数，他们的积分算法有增量法、数值积分法（X阶龙格-库塔）等等 4、加计向量、重力向量：加计基于机体，重力基于地理，重力向量（0,0,1g）用公式换算到机体，与机体的加计向量算出误差。理论上应该没有误差，这误差逆向思维一下，其实就是换算公式的系数误差。所以这误差可用于纠正公式的系数（横滚、俯仰），也就是姿态。 5、电子罗盘向量、地磁向量：同上，只不过要砍掉地理上的垂直向量，因为无用。只留下地理水平面上的向量。误差可以用来纠正公式的系数（航向）。 6、就这样，系数不停地被陀螺积分更新，也不停地被误差修正，它和公式所代表的姿态也在不断更新。 如果积分和修正用四元数算法（因为运算量较少、无奇点误差），最后用欧拉角输出控制PID（因为角度比较直观），那就需要有个四元数系数到欧拉角系数的转换。常用的三种公式，它们之间都有转换算法。 再搞个直白一点的例子： 机体好似一条船，地理就是那地图，姿态就是航向（船头在地图上的方位），重力和地磁是地图上的灯塔，陀螺/积分公式是舵手，加计和电子罗盘是瞭望手。 舵手负责估计和把稳航向，他相信自己，本来船向北开的，就一定会一直往北开，觉得转了90度弯，那就会往东开。 当然如果舵手很牛逼，也许能估计很准确，维持很长时间。不过只信任舵手，肯定会迷路，所以一般都有地图和瞭望手来观察误差。 瞭望手根据地图灯塔方位和船的当前航向，算出灯塔理论上应该在船的X方位。然而看到实际灯塔在船的Y方位，那肯定船的当前航向有偏差了，偏差就是ERR=X-Y。 舵手收到瞭望手给的ERR报告，觉得可靠，那就听个90%*ERR，觉得天气不好、地图误差大，那就听个10%*ERR，根据这个来纠正估算航向。。
捷联惯导系统快速罗经初始对准方法研究
捷联惯导系统快速罗经初始对准方法研究1 严恭敏1，严卫生1,2，徐德民1,2 1西北工业大学航海学院，西安（710072） 2水下信息处理与控制国家级重点实验室，西安（710072） E-mail：yangongmin@https://www.360docs.net/doc/01201580.html, 摘要：在分析平台罗经初始对准原理基础上，提出了捷联罗经初始对准的原理并推导了适合于软件编程的算法。将捷联罗经对准的具体实现划分为四个阶段：方位角未知情况下的水平对准、粗略方位自对准、重新水平对准和罗经方位对准，通过对大方位误差角捷联惯导非线性误差方程的简化，推导了粗略方位自对准的算法公式。如果导航计算机存储容量足够大并且计算能力足够强，根据捷联惯导系统数学平台多样性和可进行逆向姿态控制的特点，设计了一种用于缩短捷联罗经初始对准时间的具体步骤。最后，试验表明快速捷联罗经对准方案是有效的。 关键词：捷联惯导系统，罗经效应，初始对准，逆向控制 中图分类号：V249.3 1. 引言 平台惯导系统罗经初始对准过程通常可分为两步，先是水平调平，然后是方位对准。方位对准在水平调平的基础上进行，一般采样罗经方位对准方法。方位罗经对准利用的是罗经效应，也就是，在正确的平台跟踪当地地理坐标系的角速率控制指令下，如果平台存在方位轴向的偏差角，平台将产生绕东向轴的倾斜，该倾斜能由北向加速度计感测到，利用北向加速度计的输出并设计适当的控制规律，控制平台方位轴朝减小方位偏差方向转动，实现平台自动寻北。捷联惯导系统初始对准通常可分为粗对准和精对准两个阶段：在粗对准阶段，利用地球自转角速度和重力加速度作为参考量，通过惯性器件的测量输出建立粗略的导航计算坐标系；在精对准阶段，通过现代控制理论最优估计方法估计出失准角，获得准确的姿态矩阵[1,2]。 捷联惯导系统经典解析式粗对准方法难以适应晃动干扰环境，有不少文献研究了晃动基座下的初始对准问题并且也出现一些应用实例，激光陀螺和光纤陀螺的发展和不断成熟为捷联罗经的研究注入了新的活力[3-6]。从本质上说，捷联惯导系统与平台惯导系统是相同的，前者以数学平台（利用姿态矩阵、四元数或欧拉角等数学工具）模拟后者的实体平台，描述捷联惯导系统相对于参考坐标系的空间方位。平台惯导系统中实体平台具有隔离外界干扰的作用，因而平台罗经能够实现晃动基座下的初始对准，同理，在捷联惯导系统初始对准中也可以根据平台罗经初始对准的特点，建立相应的数学平台隔离晃动影响。经典控制理论与现代最优估计方法相比，前者的优点之一是勿需精确的数学模型与噪声模型，应用经典控制理论进行罗经对准的设计方法已经非常成熟，为捷联罗经对准方案设计提供了大量的参考，然而初始对准时间长是平台罗经的一大缺点。快速初始对准是国内在捷联罗经对准方法研究中亟待解决的一个主要问题，该问题在某些西方国家已得到较好解决，例如法国iXSea公司的OctansIII型光纤陀螺罗经在动态环境下，能在3min内完成初始对准，达到0.2o×sec(L)的精度[5]，成为捷联罗经研究与应用中的佼佼者，它为我们的研究和工程开发目标提供了参考。 本文从分析平台罗经初始对准的原理出发，提出了捷联罗经初始对准的原理并推导了便于软件编程的算法，通过对大方位误差角捷联惯导非线性误差方程的简化，推导了粗略方位 1本课题得到水下信息处理与控制国家级重点实验室基金（9140C230206070C2306）的资助。
捷联式惯导系统初始对准
捷联式惯导系统初始对准 惯性技术是惯导(惯性导航与惯性制导)技术、惯性仪表技术、惯性测量技术以及有关设备和装置技术的统称。惯性导航与惯性制导是当今非常重要的综合技术之一，它广泛用于航空、航海、航天及陆地各领域。惯性导航系统是和用陀螺与加速度计通过最初的方向基准和位置信息来确定运载体在一特定坐标系内的姿态、位置、速度和加速度的自主式导航系统。惯性制导系统是利用运载体内部的陀螺、加速度计测量其运动参数，经过计算机发出控制指令，从而把运载体按照预定的路线准确地引导到目的地的制导系统。自主性是惯性系统最重要的特点。确定运动对象导航参数的方法和仪器有许多，例如磁、天文、无线电、水声、全球卫星定位系统等等，然而它们都有一个致命的弱点，即不是自主的，不是要向外界发出信息，就是要依赖对外观测信息，而惯性系统与上述诸方法的基本区别就在于是完全自主的，即导弹、潜艇、飞船等可以在一个完全与外界条件以及电磁波隔绝的假想“封闭”空间内实现精确导航。因此，惯导系统具有隐蔽性好、抗干扰、不受任何气象条件限制的优点，且数据更新速率高，可以提供连续实时的导航参数。 惯性系统在国防科学技术中占有非常重要的地位，因而是世界各工业强国重点发展的技术领域之一。随着惯性技术的不断发展，许多国家已将其应用领域扩大到现代化交通运输，海洋开发，大地测量与勘探，石油钻井，矿井、隧道的掘进与贯通，机器人控制，现代化医疗器械，摄影技术以及森林防护，农业播种、施肥等民用领域。惯性技术的发展表明：从传统的机械转子型陀螺向固态陀螺仪(激光、光纤陀螺仪)转移，并进一步向以半导体硅为基本材料的微机械振动陀螺发展；从框架式平台系统向捷联系统转移，从纯惯性捷联系统向以惯性系统为基础的多体制组合导航系统发展，成为今后惯性技术发展的总趋势。 捷联式惯性导航系统，导航用的加速度计是直接捆绑在运载体上，它测量的是运载体坐标系轴向比力，只要把这个比力转换到惯性坐标系上，则其他计算就和空间稳定的平台式惯性导航系统一样，而比力转换的关键就是要实时地进行姿态基准计算来提供数学平台，即实时更新姿态矩阵b C，有些资料上称姿态矩阵 g 为捷联矩阵或方向余弦矩阵b C。一般选择地理坐标系为导航坐标系，那么捷联 g 矩阵也可表示为i C，其导航原理图如下所示： g
捷联惯导姿态算法中的圆锥误差与量化误差
第27卷 第1期航 空 学 报 Vo l 27No 1 2006年 1月ACT A A ERON A U T ICA ET A ST RO N AU T ICA SIN ICA Jan. 2006 收稿日期:2004 09 30;修订日期:2005 04 27基金项目:国家自然科学基金(60234030)、国家杰出青年科学基金 (60225015)和教育部高校青年教师奖资助项目 文章编号:1000 6893(2005)01 0098 06捷联惯导姿态算法中的圆锥误差与量化误差 练军想,胡德文,胡小平,吴文启 (国防科技大学机电工程与自动化学院自动控制系,湖南长沙 410073) Research on C oning Error and Quantization Error of SINS Attitude Algorithm LIAN Jun x iang ,H U De w en,H U Xiao ping ,WU W en qi (Department of A utomatic Contr ol,Colleg e of M echat ronics and A utomation,N atio nal U niver sity of Denfense T echno log y,Chang sha 410073,China) 摘 要:对捷联惯导系统的误差源进行了研究,利用几何方法分析了不可交换性误差和量化误差的形成机理,以及它们的相互影响。针对工程应用中激光陀螺输出脉冲采样量化条件,就多子样算法进行了讨论,并设计了基于M A T L A B/Simulink 的仿真。研究结果表明,当考虑量化误差的影响时,选取适当的量化因子,三子样等效旋转矢量算法比其它算法具有更好的综合性能。 关键词:捷联惯导;姿态算法;不可交换性误差;量化误差;多子样中图分类号:V 249.3 文献标识码:A Abstract:In this paper,the er ror so urces o f str apdow n iner tial nav igat ion system (SIN S)are r esear ched.T he mechanism of nonco mmut ativity er ro r and quantizatio n er ro r is illustr ated fr om the g eometrical point of v iew,and the mutual influence betw een them is analy zed.Considering the sample quantizing co ndition of the o utput pulses of the r ing laser g yro sco pe (R LG ),the multi sample att itude alg or ithm o f SIN S is discussed.T he simu latio n is car ried out using M A T L A B/Simulink.It is co ncluded that the three sam ple rotat ion vector alg or ithm wit h t he pro per quant izat ion factor outper for ms o thers when the impact of quant izat ion err ors is taken into ac count. Key words:SI NS;attitude algo rithm;no nco mmutat ivit y er ror ;quantizatio n er ro r;multi sample 对于捷联惯导系统姿态更新算法的研究,目 前国内外同行们有很多成果。普遍的观点认为[1,2]:四元数姿态更新算法比方向余弦和欧拉角等其它方法的计算量小、适应性好;双子样、三子样等效旋转矢量算法在抑制不可交换性误差方面比单子样有显著改善;将姿态更新分成快慢两个回路,即快速回路更新旋转矢量,慢速回路更新姿态参数,既能有效抑制不可交换性误差,还可以减少算法的计算量。但上述观点很少有结合捷联惯导系统工程实际来考虑的。本文针对工程中激光陀螺输出脉冲采样量化条件,就多子样算法展开研究,并得出了相应的结论。 1 误差源分析 评价算法优劣的最重要的标准之一是算法的 精度。要提高算法的精度就必须研究各种误差及其补偿、抑制方法。 在对惯导系统作原理性分析时,尚可以理想化地认为平台坐标系(平台式或捷联式)准确、无误差地模拟地理坐标系,但实际的惯导系统不可避免地受到各种误差的影响[3]。通常,可将误差分为源信息误差和解算方法误差。对捷联惯导系统,解算方法误差是指数值求解姿态矩阵和比力方程时产生的原理性误差,比如不可交换性误差、舍入误差、截断误差等。源信息误差可分为解析源信息误差和测量源信息误差。测量源信息误差就是惯性传感器的误差,此类误差因惯性器件的物理原理和仪表的精度不同而各不相同。比如陀螺的常值漂移和随机漂移、刻度因子误差、非线性误差、输出值的量化误差等。 下面重点对捷联姿态解算时不可交换性误差和陀螺采样的量化误差进行分析。1 1 不可交换性误差 在利用方向余弦矩阵微分方程或四元数微分 方程进行姿态更新时,都会遇到角速度矢量的积分问题。由于刚体有限转动的不可交换性,当转
2捷联惯性导航系统初始对准原理
第二章 捷联惯导系统的初试对准 2.1引言 惯导系统是一种自主式导航系统。它不需要任何人为的外部信息，只要给定导航的初始条件（例如初始速度、位置等），便可根据系统中的惯性敏感元件测量的比力和角速率通过计算机实时地计算出各种导航参数。由于“平台”是测量比力的基准，因此“平台”的初始对准就非常重要。对于平台惯导系统，初试对准的任务就是要将平台调整在给定的导航坐标系的方向上。若采用游动方位系统，则需要将平台调水平---称为水平对准，并将平台的方位角调至某个方位角处---称为方位对准。对于捷联惯导系统，由于捷联矩阵T 起到了平台的作用，因此导航工作一开始就需要获得捷联矩阵T 的初始值，以便完成导航的任务。显然捷联惯导系统的初始对准就是确定捷联矩阵的初始值。在静基座条件下，捷联惯导系统的加速度计的输入量为---b g ，陀螺的输入量为地球自转角速率b ie ω。因此b g 与 b ie ω就成为初始对准的基准。将陀螺与加速度计的输入引出计算机，通过计算机 就可以计算出捷联矩阵T 的初始值。 由以上的分析可以看出，陀螺与加速度计的误差会导致对准误差；对准飞行器的干扰运动也是产生对准误差的重要因素。因此滤波技术对捷联系统尤其重要。由于初始对准的误差将会对捷联惯导系统的工作造成难以消除的影响，因此研究初始对准的误差传播方程也是非常必要的。 2.2 捷联惯导系统的基本工作原理 捷联式惯性导航系统，陀螺仪和加速度计直接与载体固联，加速度计测量是载体坐标系轴向比力，只要把这个比力转换到导航坐标系上，则其它计算就与平台式惯性导航系统一样，而比力转换的关键就是要实时地进行姿态基准计算来提供数学平台，即实时更新姿态矩阵n b C ，姿态矩阵也称为捷联矩阵。一般选择地理坐标系为导航坐标系，那么捷联矩阵n b C 也可表示为t b C ， 其导航原理图如图2.1所示。
捷联式惯性导航系统
1 绪论 随着计算机和微电子技术的迅猛发展，利用计算机的强大解算和控制功能代替机电稳定系统成为可能。于是，一种新型惯导系统--捷联惯导系统从20世纪60年代初开始发展起来，尤其在1969年，捷联惯导系统作为"阿波罗"-13号登月飞船的应急备份装置，在其服务舱发生爆炸时将飞船成功地引导到返回地球的轨道上时起到了决定性作用，成为捷联式惯导系统发展中的一个里程碑。 捷联式惯性导航(strap-down inertial navigation) ，捷联（strap-down）的英语原义是“捆绑”的意思。因此捷联式惯性导航也就是将惯性测量元件（陀螺仪和加速度计）直接装在飞行器、舰艇、导弹等需要诸如姿态、速度、航向等导航信息的主体上，用计算机把测量信号变换为导航参数的一种导航技术。现代电子计算机技术的迅速发展为捷联式惯性导航系统创造了条件。惯性导航系统是利用惯性敏感器、基准方向及最初的位置信息来确定运载体的方位、位置和速度的自主式航位推算导航系统。在工作时不依赖外界信息，也不向外界辐射能量，不易受到干扰破坏。它完全是依靠载体自身设备独立自主地进行导航，它与外界不发生任何光、声、磁、电的联系，从而实现了与外界条件隔绝的假想的“封闭”空间内实现精确导航。所以它具有隐蔽性好，工作不受气象条件和人为的外界干扰等一系列的优点，这些优点使得惯性导航在航天、航空、航海和测量上都得到了广泛的运用[1] 1.1 捷联惯导系统工作原理及特点 惯导系统主要分为平台式惯导系统和捷联式惯导系统两大类。惯导系统（INS）是一种不依赖于任何外部信息、也不向外部辐射能量的自主式导航系
统，具有隐蔽性好，可在空中、地面、水下等各种复杂环境下工作的特点。 捷联惯导系统（SINS）是在平台式惯导系统基础上发展而来的，它是一种无框架系统，由三个速率陀螺、三个线加速度计和微型计算机组成。平台式惯导系统和捷联式惯导系统的主要区别是：前者有实体的物理平台，陀螺和加速度计置于陀螺稳定的平台上，该平台跟踪导航坐标系，以实现速度和位置解算，姿态数据直接取自于平台的环架；后者的陀螺和加速度计直接固连在载体上作为测量基准，它不再采用机电平台，惯性平台的功能由计算机完成，即在计算机内建立一个数学平台取代机电平台的功能，其飞行器姿态数据通过计算机计算得到，故有时也称其为"数学平台"，这是捷联惯导系统区别于平台式惯导系统的根本点。由于惯性元器件有固定漂移率，会造成导航误差，因此，远程导弹、飞机等武器平台通常采用指令、GPS或其组合等方式对惯导进行定时修正，以获取持续准确的位置参数。如采用指令+捷联式惯导、GPS+惯导（GPS/INS）。美国的战斧巡航导弹采用了GPS+INS +地形匹配组合导航。 惯导系统基本工作原理是以牛顿力学定律为基础，通过测量载体在惯性参考系的加速度，将它对时间进行积分，之后将其变换到导航坐标系，得到在导航坐标系中的速度、偏航角和位置信息等。对捷联惯导系统而言，平台的作用和概念体现在计算机中，它是写在计算机中的方向余弦阵。直接安装在载体上的惯性元件测得相对惯性空间的加速度和角加速度是沿载体轴的分量，将这些分量经过一个坐标转换方向余弦阵，可以转换到要求的计算机坐标系内的分量。如果这个矩阵可以描述载体和地理坐标系之间的关系，那么载体坐标系测得的相对惯性空间的加速度和角速度，经过转换后便可得到沿地理坐标系的加速度和角速度分量，有了已知方位的加速度和角速度分量之后，导航计算机便
车载捷联惯导系统基本原理
车载捷联惯导系统基本原理 一、捷联惯导系统基本原理 捷联惯导系统基本原理如图2-1所示： 图中陀螺和加速度计直接与载体系b固联，用来测量载体的角运动信息和线运动信息。导航解算的本质是根据初值进行积分的过程，通过求解姿态微分方程完成对姿态和航向角的积分，通过求解比力微分方程完成对速度的积分，通过求解位置微分方程实现对位置的积分。捷联惯导的姿态矩阵C n 相当于“数学平台”，取代了平台惯导中的实体平台，而ω?相当于对数学平台“施矩”的指令角速率。
二、捷联惯导微分方程 （一）姿态微分方程 在捷联惯导系统中，导航坐标系n 和载体坐标系b 之间的角位置关系通常用姿态矩阵、四元数和欧拉角表示，相应也存在姿态矩阵微分方程、四元数微分方程和欧拉角微分方程三种形式。 姿态矩阵微分方程的表达式为：
在欧拉角微分方程式(2.2-7)中，当俯仰角θ趋于90o时，cosθ趋于0，tanθ趋于无穷，方程存在奇异性，所以这种方法不能在全姿态范围内正常工作；姿态矩阵微分方程式(2.2-1)可全姿态工作，但姿态矩阵更新相当于求解包含9个未知量的线性微分方程组，计算量大；四元数微分方程式(2.2-6)同样可以全姿态工作，且更新算法只需求解4个未知量的线性微分方程组，计算量小，算法简单，是较实用的工程算法。 （二）速度微分方程 速度微分方程即比力方程，是惯性导航解算的基本关系式： 三、捷联惯性导航算法 捷联惯导解算的目的是根据惯性器件输出求解载体姿
态、速度和位置等导航信息，实际上就是求解三个微分方程的过程，相应存在姿态更新算法、速度更新算法和位置更新算法。 （一）姿态更新算法 求解微分方程式(2.2-6)可得四元数姿态更新算法为：
导航原理_捷联惯导系统
导航原理作业（惯性导航部分）
一枚导弹采用捷联惯性导航系统，三个速率陀螺仪Gx, Gy, Gz 和三个加速度计Ax, Ay, Az 的敏感轴分别沿着着 弹体坐标系的Xb, Yb, Zb轴。初始时刻该导弹处在北纬 45.75度，东经126.63度。 第一种情形：正对导弹进行地面静态测试（导弹质心相对地面静止）。 初始时刻弹体坐标系和地理坐标系重合，如图所示，弹体的Xb轴指东，Yb轴指北，Zb轴指天。此后弹体坐标系Xb-Yb-Zb 相对地理坐标系的转动如下： 首先，弹体绕Zb（方位轴）转过-10 度； 接着，弹体绕Xb（俯仰轴）转过15 度； 然后，弹体绕Yb（滚动轴）转过20 度； 最后弹体相对地面停止旋转。 请分别用方向余弦矩阵和四元数两种方法计算：弹体经过三次旋转并停止之后，弹体上三个加速度计Ax, Ay, Az的输出。取重力加速度的大小g = 9.8m/s2。 第二种情形：导弹正在飞行中。 初始时刻弹体坐标系仍和地理坐标系重合；且导弹初始高度200m，初始北向速度1800 m/s，初始东向速度和垂直速度都为零。 陀螺仪和加速度计的输出都为脉冲数形式，陀螺输出的每个脉冲代表0.00001弧度的角增量。加速度计输出的每个脉冲代表1μg，1g = 9.8m/s2。陀螺仪和加速度计输出的采样频率都为10Hz，在200秒内三个陀螺仪和三个加速度计的输出存在了数据文件gaout.mat中，内含一矩阵变量ga，有2000行，6列。每一行中的数据代表每个采样时刻三个陀螺Gx, Gy, Gz和三个加速度计Ax, Ay, Az 的输出的脉冲数。格式如下表（前10行）
将地球视为理想的球体，半径6371.00公里，且不考虑仪表误差，也不考虑弹体高度对重力加速度的影响。选取弹体的姿态计算周期为0.1秒，速度和位置的计算周期为1秒。 （1）请计算200秒后弹体到达的经纬度和高度，东向和北向速度； （2）请计算200秒后弹体相对当地地理坐标系的姿态四元数； （3）请绘制出200秒内导弹的经、纬度变化曲线（以经度为横轴，纬度为纵轴）； （4）请绘制出200秒内导弹的高度变化曲线（以时间为横轴，高度为纵轴）。 二、程序设计说明及代码 1.第一种情形 （1）方向余弦矩阵法 1)程序代码 clear;clc; thetax=15*pi/180;thetay=20*pi/180;thetaz=(-10)*pi/180; A0=[0;0;-9.8]; Theta=[0,-thetaz,thetay;thetaz,0,-thetax;-thetay,thetax,0]; theta0=sqrt(thetax^2+thetay^2+thetaz^2); S=(sin(theta0))/theta0;C=(1-cos(theta0))/theta0^2; CT=eye(3)+S*Theta+C*(Theta^2); CTN=inv(CT); A1=CTN*A0 2)输出结果 （2）四元数法 1)程序代码
捷联惯导系统性能分析
捷联惯导系统性能分析 摘要 本文简要介绍捷联式惯性导航系统的各种分析技术,对捷联惯导系统算法验证的过程进行了讨论。封闭形式分析仿真驱动程序，可以用来锻炼/验证捷联算法方程。分析的精度捷联络筒，划船和分析方法位置融合算法（包括位置算法的折叠效果）函数的算法重复率和系统振动输入。包括的是一个简化的分析模型，该模型描述可用于翻译系统到惯性传感器作为传感器组件的函数的输入的振动安装失衡。捷联系统静态漂移和的旋转测试程序/方程描述捷联式传感器的校准系数确定。该文件概述了卡尔曼滤波器的设计和协方差分析技术，并介绍了验证辅助捷联惯导系统的一般步骤卡尔曼滤波器的配置。最后，论述了系统集成测试的一般过程验证所有的硬件，软件，系统功能操作进行正确和准确和界面元素。 坐标框架 在本文中，使用一个坐标系是一个分析性的抽象定义的三个相互垂直单位矢量。一个坐标系可以看作一组三个相互垂直的线（轴）通过一个共同的点（原点）与来自沿着坐标轴的原点上的单位向量。在本文中，每个坐标系的原点的物理位置是任意的。主坐标帧利用有以下几种： B帧=“身体”捷联式惯性传感器轴平行的坐标系。 N帧=“导航”的坐标系在当地具有Z轴平行地垂直向上定位。A“漂移方位”N帧有水平的X，Y轴旋转相对于非旋转的惯性空间在本地地球的速度的垂直分量绕Z轴。“自由方位”N帧的转动惯量为零率X，Y轴的Z轴周围。“地理”N帧的X，Y轴绕Z轴旋转，以维持当地的真北Y轴平行。 E型=“地球”的引用与固定角几何相对于地球坐标系。 I帧=“惯性”非旋转坐标系。 符号 V =向量没有特定的坐标系指定。向量是一个有长度的参数和方向。纸中使用的载体，被分类为“免费的矢量”，因此，没有坐标框架中，他们分析描述的理想地点。 V A =列矩阵的元素等于V的坐标系A轴的投影。“ 投影V对每个Frame A轴等于与坐标框架的V的点积A轴的单位向量。 V A×=斜对称的（或跨产品）的形式表示的正方形矩阵的V A 0 - VZA VYA VZA 0 - VXA - VYA VXA 0 VXA，VYA，VZA是V A的组成部分。“ 另一种框架向量V A×矩阵乘积等于跨产品的V A用该载体在A帧。 CA2 A1 =方向余弦矩阵，将一个向量的坐标系A2投影 形成坐标系A1投影形式。 ωA1A2的坐标系相对坐标系A1 A2 =角速率。当A1是非 旋转，ωA1A2是将角速度传感器测量的角速率安装在A2帧。 = e DT =相对于时间的导数。 t =时间。 一、简介
捷联惯导详细讲解
捷联惯导系统从20世纪60年代初开始发展起来，在1969年，捷联惯导系统作为"阿波罗"-13号登月飞船的应急备份装臵，在其服务舱发生爆炸时将飞船成功地引导到返回地球的轨道上时起到了决定性作用，成为捷联式惯导系统发展中的一个里程碑。 捷联式惯性导航(strap-downinertialnavigation)，捷联（strap-down）的英语原义是“捆绑”的意思。因此捷联式惯性导航也就是将惯性测量元件（陀螺仪和加速度计）直接装在导弹需要诸如姿态、速度、航向等导航信息的主体上，用计算机把测量信号变换为导航参数的一种导航技术。 一、捷联惯导系统工作原理及特点 惯导系统基本工作原理是以牛顿力学定律为基础，通过测量载体在惯性参考系的加速度，将它对时间进行积分，之后将其变换到导航坐标系，得到在导航坐标系中的速度、偏航角和位臵信息等。 捷联惯导系统（SINS）是一种无框架系统，由三个速率陀螺、三个线加速度计和微型计算机组成。由于惯性元器件有固定漂移率，会造成导航误差，因此导弹通常采用指令、GPS或其组合等方式对惯导进行定时修正，以获取持续准确的位臵参数。如采用指令+捷联式惯导 捷联惯导系统能精确提供载体的姿态、地速、经纬度等导航参数，是利用惯性敏感器、基准方向及最初的位臵信息
来确定运载体的方位、位臵和速度的自主式航位推算导航系统。在工作时不依赖外界信息，也不向外界辐射能量，不易受到干扰破坏。它完全是依靠载体自身设备独立自主地进行导航，它与外界不发生任何光、声、磁、电的联系，从而实现了与外界条件隔绝的假想的“封闭”空间内实现精确导航。所以它具有隐蔽性好，工作不受气象条件和人为的外界干扰等一系列的优点。 除此以外捷联惯导系统的最大特点是没有实体平台，即将陀螺仪和加速度计直接安装在机动载体上，在计算机中实时的计算姿态矩阵，通过姿态矩阵把导航加速度计测量的载体沿机体坐标系轴向的加速度信息变换到导航坐标系，然后进行导航计算。同时，从姿态矩阵的元素中提取姿态和航向信息．由此可见，在捷联惯导系统中平台的作用已由计算机及其软件的作用代替了，捷联式惯导系统采用的是数学平台。力学编排就是按照合适的数学模型由观测量计算出导航定位参数。具体地讲，利用陀螺仪测得的载体相对于惯性参照系的旋转角速度，计算出载体坐标系至导航计算坐标系之问的坐标转换矩阵；将测量的比力(加速度计测量载体相对于惯性空间的线加速度)变换至导航坐标系，并经过两次积分得到所需的速度位臵信息。 二、捷联惯导系统有以下独特优点： （1）去掉了复杂的平台机械系统，系统结构极为简单，