基于捷联惯导的运动轨迹跟踪技术研究

目录

摘要 ............................................................................................................................... I ABSTRACT ..................................................................................................................... I I 第1章绪论 .. (1)

1.1课题研究的背景与意义 (1)

1.2国内外研究现状及分析 (2)

1.2.1传统捷联式惯性导航的研究现状与分析 (2)

1.2.2运动轨迹的研究现状与分析 (4)

1.3本文主要内容和结构安排 (7)

第2章捷联式惯性导航系统 (8)

2.1引言 (8)

2.2常用关系式 (8)

2.3姿态更新的四元数算法 (12)

2.3.1四元数 (12)

2.3.2四元数微分方程 (13)

2.3.3毕卡更新算法 (14)

2.4姿态更新算法的补偿 (15)

2.4.1互补滤波算法 (15)

2.4.2卡尔曼滤波器 (17)

2.4.3梯度下降法 (19)

2.5本章小结 (20)

第3章运动轨迹跟踪算法 (21)

3.1引言 (21)

3.2预处理 (21)

3.2.1差分法截取数据 (21)

3.3.2多采样率处理 (23)

3.3速度更新 (25)

3.3.1速度更新算法 (25)

3.3.2速度误差补偿 (26)

3.4安装在鞋上的惯性导航 (28)

3.4.1零速度检测 (29)

3.4.2步态周期判别 (30)

3.4.3运动轨迹校正 (32)

3.5本章小结 (34)

第4章运动轨迹跟踪系统实现 (35)

4.1引言 (35)

4.2生成模拟数据 (35)

4.2.1模拟运动姿态 (35)

4.2.2模拟运动轨迹 (36)

4.3姿态解算算法仿真 (36)

4.3.1毕卡法更新 (37)

4.3.2有误差补偿的姿态解算 (38)

4.4数据的采集与接收 (42)

4.4.1惯导硬件采集模块 (42)

4.4.2蓝牙接收 (44)

4.5运动轨迹跟踪实验 (45)

4.6本章小结 (48)

结论 (49)

参考文献 (50)

攻读硕士学位期间发表的专利及其它成果 (54)

哈尔滨工业大学学位论文原创性声明和使用权限 (55)

致谢 (56)

北航卡尔曼滤波课程-捷联惯导静基座初始对准实验

卡尔曼滤波实验报告 捷联惯导静基座初始对准实验 一、实验目的 ①掌握捷联惯导的构成和基本工作原理； ②掌握捷联惯导静基座对准的基本工作原理； ③了解捷联惯导静基座对准时的每个系统状态的可观测性； ④了解双位置对准时系统状态的可观测性的变化。 二、实验原理 选取状态变量为：[]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)

捷联惯导系统粗对准方法比较

捷联惯导系统粗对准方法比较 魏春岭 张洪钺 北京航空航天大学自动化科学与电气工程学院 北京 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

第六章 捷联惯导

第六章捷联惯导

6-1捷联惯导的原理?捷联惯导系统概述 ?捷联惯性技术的发展过程 ?捷联惯导系统与平台惯导系统的对比 ?捷联惯导系统的基本力学编排方程?捷联惯导系统的算法概述 ?捷联惯导系统原理框图的说明 ?姿态方程的解算 （1）姿态和航向角的计算 （2）姿态矩阵的微分方程 （3）四元数的运动学微分方程 （4）等效旋转矢量法及其微分方程 （5）位移角速率方程 （6）速度方程

?导航位置方程 （1）游动方位系与地球系之间的方向余弦矩阵 （2）载体位置计算 （3）方向余弦矩阵计算 ?垂直通道阻尼 ?捷联惯性器件的余度技术?单自由度陀螺仪的配置方案 （1）四陀螺仪配置方案 （2）六陀螺仪系统 ?二自由度陀螺仪的配置方案

?捷联惯导的数值计算方法?数值积分法 （1）欧拉法 （2）四阶龙格－库塔法 ?角速率信息的提取

“ 捷联（Strapdown）”这一术语的英文原义就是“捆绑”的意思。因此，所谓捷联惯性系统也就是将惯性敏感元件（陀螺与加速度计）直接“捆绑”在载体上，从而完成制导和导航任务的系统。 V-2导弹 “阿波罗-13”宇宙飞船 “海盗”火星降落器

从捷联技术的发展过程中我们已经看到捷联系统的优越性已越来越突出的显示出来，并在许多方面已日渐代替平台系统。为什么会出现这种情况呢？为了回答这一问题，这里从生产与使用的角度将捷联系统与平台系统做一对比。 （1）硬件和软件的复杂程度 由于捷联系统没有平台框架及相连的伺服装置，因而简化了硬件；代价是增加了计算机的负担，需要一个比较复杂的实时程序。 （2）可靠性 捷联系统的可靠性要比平台系统高，其原因是它的机械构件少，加之容易采用多敏感元件配置，实现余度技术。 （3）成本与可维护性 由于平台系统在机械结构上要复杂得多，而对于捷联系统只是算法复杂些，因而从制造成本上看捷联系统的成本要比平台系统低。从市场供应的情况来看，数字计算机的价格一直在下降，而平台系统的价格一直在上升。 此外，捷联系统比平台系统具有较长的平均故障间隔时间，加之模块设计简化了维修，从而捷联系统的可维护性比平台系统大为提高了。

捷联惯导系统初始对准技术的研究

捷联惯导系统初始对准技术的研究 摘要:初始对准是捷联惯导系统关键技术之一。初始对准精度直接影响捷联惯导系统的工作精度，初始对准时间也是反映武器系统快速反应能力的重要战术指标。 捷联惯性导航系统是将惯性器件陀螺仪、加速度计构成的惯性测量单元直接与载体固联，测量得到的载体角速度与线运动参数是沿载体固联的坐标轴上的分量。导航计算机通过计算“姿态矩阵”可以将加速度信息转换到惯性坐标系或当地地理坐标系，从而实现“数学平台”，然后再进行速度及位置计算。图1即为捷联式惯性导航系统原理框图。 捷联惯导系统的关键技术包括初始对准问题、有害加速度的消除及引力修正、惯性元件误差模型的建立和实时补偿、捷联矩阵的更新等。捷联惯性导航系统初始对准的目的是建立捷联矩阵的初始值。 1、捷联惯导系统初始对准基本概念 按对准阶段来分，初始对准一般分为两个阶段:第一阶段为粗对准，第二阶段为精对准。捷联系统粗对准的任务是得到粗略的捷联矩阵，为后续的精对准提供基础，此阶段精度可以低一些，但要求速度快。精对准是在粗对准的基础上进行的，通过处理惯性敏感元件的输出信息，精确校正真实导航坐标系与计算的导航坐标系之间的失准角，使之趋于零，从而得到精确的捷联矩阵。按照捷联惯性导航系统初始对准时载体的运行状态来分，可分为静基座对准和动基座对准。按

照初始对准时是否取得外部信息，可分为自对准和非自对准.惯性导航系统的自对准是利用重力矢量和地球自转角速率矢量通过解析的方法实现的初始对准，这种对准方法的优点是自主性强，缺点是所需的对准时间长。非自主式对准可以通过机电或光学方法将外部参考坐标系引入系统，实现惯性系统的初始对准.在捷联惯性导航系统的粗对准阶段，可以通过引入主惯导系统的航向姿态信息，通过传递对准，迅速将数学平台对准导航坐标系，减小初始失准角.在精对准阶段，可以通过组合导航的方法，利用其它导航设备(如GPS，计程仪)等提供的信息(如速度和位置)作为观测信息，通过卡尔曼滤波实现精确对准。目前有关初始对准问题的研究主要集中在误差模型的建立、模型求解方法和误差模型的可观性分析三个方面，本文正是针对这三个方面，对初始对准的国内外研究状况进行综述。 2、初始对准的误差模型 捷联惯导系统初始对准的误差模型及常用算法研究的基础模型有Ψ角误差模型和Φ角误差模型。前者基于导航坐标系，是指导航坐标系到计算坐标系(通常取计算的地理坐标系)的小角度误差；后者基于真实坐标系，是指导航坐标系到真实坐标系(通常取地理坐标系)的小角度误差。 3、求解误差模型的方法 3. 1古典方法 从频率域角度设计对准回路，并将对准分为水平对准和方位对准，对准过程首先是水平粗对准，然后是方位粗对准;在粗对准之后再精对准，首先是水平精对准，然后是方位精对准.由于对准回路频带低，响应慢，因此整个对准时间长。针对在静基座捷联惯导系统初始对准中，东向陀螺漂移没有估计效果，使得估计的方位失准角存在常值误差的问题，上海交大提出了，建立东向陀螺漂移估计位的修正方程，并对其进行修正，从而大大提高了方位失准角的估计精度。 3.2 Kalman滤波在初始对准中的应用 提高惯性导航系统初始对准精度的最佳途径之一是利用Kalman滤波这一重要数学工具，对于采用自主对准方式的惯导系统，一般采用卡尔曼滤波技术估计出系统的失调角和惯导系统误差源，然后采用一定的控制技术设计出控制角速率使失调角达到规定的要求。卡尔曼滤波的对象是用状态方程来描述随机线性系统，它按照估计误差方差最小的原则，从被污染的观测值中，实时估计出系统的各个

捷联式惯性导航系统

1 绪论 随着计算机和微电子技术的迅猛发展，利用计算机的强大解算和控制功能代替机电稳定系统成为可能。于是，一种新型惯导系统--捷联惯导系统从20世纪60年代初开始发展起来，尤其在1969年，捷联惯导系统作为"阿波罗"-13号登月飞船的应急备份装置，在其服务舱发生爆炸时将飞船成功地引导到返回地球的轨道上时起到了决定性作用，成为捷联式惯导系统发展中的一个里程碑。 捷联式惯性导航(strap-down inertial navigation)，捷联（strap-down）的英语原义是“捆绑”的意思。因此捷联式惯性导航也就是将惯性测量元件（陀螺仪和加速度计）直接装在飞行器、舰艇、导弹等需要诸如姿态、速度、航向等导航信息的主体上，用计算机把测量信号变换为导航参数的一种导航技术。现代电子计算机技术的迅速发展为捷联式惯性导航系统创造了条件。惯性导航系统是利用惯性敏感器、基准方向及最初的位置信息来确定运载体的方位、位置和速度的自主式航位推算导航系统。在工作时不依赖外界信息，也不向外界辐射能量，不易受到干扰破坏。它完全是依靠载体自身设备独立自主地进行导航，它与外界不发生任何光、声、磁、电的联系，从而实现了与外界条件隔绝的假想的“封闭”空间内实现精确导航。所以它具有隐蔽性好，工作不受气象条件和人为的外界干扰等一系列的优点，这些优点使得惯性导航在航天、航空、航海和测量上都得到了广泛的运用[1] 1.1 捷联惯导系统工作原理及特点 惯导系统主要分为平台式惯导系统和捷联式惯导系统两大类。惯导系统（INS）是一种不依赖于任何外部信息、也不向外部辐射能量的自主式导航系统，具有隐蔽性好，可在空中、地面、水下等各种复杂环境下工作的特点。 捷联惯导系统（SINS）是在平台式惯导系统基础上发展而来的，它是一种无框架系统，由三个速率陀螺、三个线加速度计和微型计算机组成。平台式惯导系统和捷联式惯导系统的主要区别是：前者有实体的物理平台，陀螺和加速度计置于陀螺稳定的平台上，该平台跟踪导航坐标系，以实现速度和位置解算，姿态数据直接取自于平台的环架；后者的陀螺和加速度计直接固连在载体上作

捷联惯导Matlab程序求解飞行器的姿态

捷联惯导程序，依据加表和陀螺仪的输出数据来求解飞行器的姿态 clc; clear; format long; %设置数据精度为15位小数 Data=importdata(''); % 导入实验所采集的数据，以矩阵形式赋给Data变量，必须与该M文件在同一个文件夹中 Px=Data(:,3); % Px,Py,Pz为陀螺仪的输出值 Py=Data(:,4); Pz=Data(:,5); Nx=Data(:,6); % Nx,Ny，Nz为加速度计的输出值 Ny=Data(:,7); Nz=Data(:,8); % 陀螺仪模型参数标定如下： Sx = ; Sy = ; Sz = ; Mxy = ; Mxz = ; Myx = ; Myz = ; Mzx = ; Mzy = ; Dx = ; Dy = ; Dz = ; GyroCali_A = [ 1 -Mxy -Mxz ; -Myx 1 -Myz ; -Mzx -Mzy 1 ]; % 加速度计模型参数标定如下： Kx = ; Ky = ; Kz = ; Ixy = ; Ixz = ; Iyx = ; Iyz = ; Izx = ; Izy = ; Bx = ; By = ; Bz = ; AccCali_A = [1 -Ixy -Ixz ; -Iyx 1 -Iyz ; -Izx -Izy 1 ]; Delta_t = ; %采样时间为秒 Delta_Theta_x = 0; Delta_Theta_y = 0; Delta_Theta_z = 0; %定义陀螺仪输出的角度增量 Delta_Vx = 0; Delta_Vy = 0; Delta_Vz = 0; %定义加速度计输出的速度增量 L = zeros(1,12001); L(1)= *pi/180 ; %纬度用L表示，纬度的初始值划为弧度形式，因为后面计算位置矩阵更新 L(2)= *pi/180 ; %时需要用到前两次的L值来计算当前L值，所以在此定义2个初始L值Lamda = *pi/180 ; %经度用Lamda表示，经度的初始值划为弧度形式 h = 136 ; %高度用h表示

捷联惯导系统快速罗经初始对准方法研究

捷联惯导系统快速罗经初始对准方法研究1 严恭敏1，严卫生1,2，徐德民1,2 1西北工业大学航海学院，西安（710072） 2水下信息处理与控制国家级重点实验室，西安（710072） E-mail：yangongmin@https://www.sodocs.net/doc/3217001236.html, 摘要：在分析平台罗经初始对准原理基础上，提出了捷联罗经初始对准的原理并推导了适合于软件编程的算法。将捷联罗经对准的具体实现划分为四个阶段：方位角未知情况下的水平对准、粗略方位自对准、重新水平对准和罗经方位对准，通过对大方位误差角捷联惯导非线性误差方程的简化，推导了粗略方位自对准的算法公式。如果导航计算机存储容量足够大并且计算能力足够强，根据捷联惯导系统数学平台多样性和可进行逆向姿态控制的特点，设计了一种用于缩短捷联罗经初始对准时间的具体步骤。最后，试验表明快速捷联罗经对准方案是有效的。 关键词：捷联惯导系统，罗经效应，初始对准，逆向控制 中图分类号：V249.3 1. 引言 平台惯导系统罗经初始对准过程通常可分为两步，先是水平调平，然后是方位对准。方位对准在水平调平的基础上进行，一般采样罗经方位对准方法。方位罗经对准利用的是罗经效应，也就是，在正确的平台跟踪当地地理坐标系的角速率控制指令下，如果平台存在方位轴向的偏差角，平台将产生绕东向轴的倾斜，该倾斜能由北向加速度计感测到，利用北向加速度计的输出并设计适当的控制规律，控制平台方位轴朝减小方位偏差方向转动，实现平台自动寻北。捷联惯导系统初始对准通常可分为粗对准和精对准两个阶段：在粗对准阶段，利用地球自转角速度和重力加速度作为参考量，通过惯性器件的测量输出建立粗略的导航计算坐标系；在精对准阶段，通过现代控制理论最优估计方法估计出失准角，获得准确的姿态矩阵[1,2]。 捷联惯导系统经典解析式粗对准方法难以适应晃动干扰环境，有不少文献研究了晃动基座下的初始对准问题并且也出现一些应用实例，激光陀螺和光纤陀螺的发展和不断成熟为捷联罗经的研究注入了新的活力[3-6]。从本质上说，捷联惯导系统与平台惯导系统是相同的，前者以数学平台（利用姿态矩阵、四元数或欧拉角等数学工具）模拟后者的实体平台，描述捷联惯导系统相对于参考坐标系的空间方位。平台惯导系统中实体平台具有隔离外界干扰的作用，因而平台罗经能够实现晃动基座下的初始对准，同理，在捷联惯导系统初始对准中也可以根据平台罗经初始对准的特点，建立相应的数学平台隔离晃动影响。经典控制理论与现代最优估计方法相比，前者的优点之一是勿需精确的数学模型与噪声模型，应用经典控制理论进行罗经对准的设计方法已经非常成熟，为捷联罗经对准方案设计提供了大量的参考，然而初始对准时间长是平台罗经的一大缺点。快速初始对准是国内在捷联罗经对准方法研究中亟待解决的一个主要问题，该问题在某些西方国家已得到较好解决，例如法国iXSea公司的OctansIII型光纤陀螺罗经在动态环境下，能在3min内完成初始对准，达到0.2o×sec(L)的精度[5]，成为捷联罗经研究与应用中的佼佼者，它为我们的研究和工程开发目标提供了参考。 本文从分析平台罗经初始对准的原理出发，提出了捷联罗经初始对准的原理并推导了便于软件编程的算法，通过对大方位误差角捷联惯导非线性误差方程的简化，推导了粗略方位 1本课题得到水下信息处理与控制国家级重点实验室基金（9140C230206070C2306）的资助。

捷联式惯导系统误差解析解研究

第22卷　第11期计　算　机　仿　真2005年11月 文章编号:1006-9348(2005)11-0042-04 捷联式惯导系统误差解析解研究 张宾,刘藻珍 (北京理工大学机电工程学院,北京100081) 摘要:该文在一定的假设条件下利用捷联惯导系统的三维误差状态模型求解出了单通道误差状态方程的解析解,列表给出 了各误差源对于某一特定误差状态的动态影响。然后利用某型导弹的弹道数据通过对两种误差模型在同一条件下进行仿 真的方法验证了单通道误差状态方程解析解的正确性。单通道误差模型对分析各种误差源对系统的影响,确定在满足系统 精度要求的条件下主要误差源的选择范围,进行系统精度分配提供了十分方便直观的方法。 关键词:捷联;误差模型;误差分析 中图分类号:V249.32 文献标识码:A Research on the Error Ana lyti c Soluti on of Strapdown I nerti a l Nav i ga ti on System ZHANG B in,L I U Zao-zhen (School of Mechanical Electr onic Engineering,Beijing I nstitute of Technol ogy,Beijing100081,China) ABSTRACT:I n this paper,err or state model of strapdown inertial navigati on syste m(SI N S)is educed and analytic s oluti on t o monochannel err or state equati on is worked out under the conditi on of certain hypotheses.The lists of the effect of each err or s ource t o a given err or status are p r ovided when SI N S is in the state of moving.The correctness of analytic s oluti on t o monochannel err or state equati on is validated by the means of t w o err or models’si m ulati on excer p2 ting the same actual traject ory data of a certain type m issile.Monochannel err or model gives a convenient and intu2 iti onistic way t o analyze the effect of all kinds of err or s ources t o the system,deli m it the selective range of main err or s ource which can meet the requirement of the syste m accuracy and all ot syste m accuracy. KE YWO RD S:Strapdown;Err or model;Err or analysis 1　引言 在导航过程中,希望惯导系统能准确地提供各种导航信 息。但各种误差源的存在,使导航信息具有一定的误差。本 文在一定的假设条件下利用捷联惯导系统的三维误差状态 模型求解出了单通道误差状态方程的解析解,列表给出了各 误差源对于某一特定误差状态的动态影响。然后利用某型 导弹的弹道数据通过对两种误差模型在同一条件下进行仿 真的方法验证了单通道误差状态方程解析解的正确性。 单通道误差模型对分析各种误差源对系统的影响,确定 在满足系统精度要求的条件下主要误差源的选择范围,进行 系统精度分配提供了十分方便直观的方法。 2　捷联惯导误差模型 当地水平坐标系(L)中,捷联惯导系统力学编排方程计 算输出的状态变量包括:大地坐标(φ,λ,h),运动速度(V e , V n,V u)及姿态信息(r,p,y)等量。此时相应的误差状态向量 δX(t)=[

捷联惯导作业

一、原理分析： 捷联式惯导系统是将惯性器件（陀螺仪和加速度计）直接固连在载体上的系统。图1为捷联式惯导系统的原理图，陀螺仪和加速度计输出分别送入姿态矩阵计算和由载体坐标系至平台坐标系的方向余弦矩阵的计算。有了姿态矩阵，其一可以实现把载体坐标系轴向加速度信息变换到导航坐标系轴，进而可以进行所需的导航参数计算，其二利用姿态矩阵的元素，提取方位和姿态信息。 图1. 捷联式惯导系统的原理图 姿态速率微分方程为: 12b tb ωΛ=Λ (1) 其中； () b b b t t tb ib t ie et C ωωωω=-- (2) b ib ω为陀螺仪测量经补偿后的值；

0cos sin t iex t t ie iey ie t ie iez L L ωωωωωω?? ? ????? ??==???????????? ，为地球自转角速率； tan t ety t yt etx t t t etx et ety xt t etz t etx xt V R V R V L R ωωωω??-?? ? ???? ????? ??==?????????????????? ，为地理坐标系相对地球坐标系的转动角速率； 导航坐标系到载体坐标系的姿态矩阵为： cos cos sin sin sin sin cos cos sin sin cos sin sin cos cos cos sin cos sin sin sin cos sin sin cos sin cos cos cos t t C ψ?ψθ? ψ?ψθ? θ?ψθ ψθ θ ψ?ψθ? ψ?ψθ? θ?-+-?? ?? =-?? ??+-?? (3) 对应的四元素初值为: 0123cos cos cos sin sin sin 2 2 2 22 2 cos sin cos sin cos sin 2 2 2 2 2 2 cos cos sin sin sin cos 2 2 2 2 2 2 cos sin sin sin cos cos 2 2 2 2 2 2ψ θ ? ψ θ ? λψ θ ? ψθ ? λψ θ ? ψ θ ? λψ θ ? ψ θ ? λ? =-???=-???=+???=+? (4) 四元素姿态矩阵为： 22220123120313022 2 2 2 12030123 230122221302230101232() 2()2() 2()2() 2() b t C λλλλλλλλλλλλλλλλλλλλλλλλλλλλλλλλλλλλ?? +++++?? =--+-+????----+? ? (5) 将姿态速率微分方程展开成矩阵形式： 0112233001020b b b tbx tby tbx b b b tbx tbz tby b b b tby tbz tbx b b b tbz tby tbx λλωωωλλωωωωωωλλωωωλλ???? ??---??????-??????=??????-??????-???????????? (6)

捷联式惯导系统初始对准

捷联式惯导系统初始对准 惯性技术是惯导(惯性导航与惯性制导)技术、惯性仪表技术、惯性测量技术以及有关设备和装置技术的统称。惯性导航与惯性制导是当今非常重要的综合技术之一，它广泛用于航空、航海、航天及陆地各领域。惯性导航系统是和用陀螺与加速度计通过最初的方向基准和位置信息来确定运载体在一特定坐标系内的姿态、位置、速度和加速度的自主式导航系统。惯性制导系统是利用运载体内部的陀螺、加速度计测量其运动参数，经过计算机发出控制指令，从而把运载体按照预定的路线准确地引导到目的地的制导系统。自主性是惯性系统最重要的特点。确定运动对象导航参数的方法和仪器有许多，例如磁、天文、无线电、水声、全球卫星定位系统等等，然而它们都有一个致命的弱点，即不是自主的，不是要向外界发出信息，就是要依赖对外观测信息，而惯性系统与上述诸方法的基本区别就在于是完全自主的，即导弹、潜艇、飞船等可以在一个完全与外界条件以及电磁波隔绝的假想“封闭”空间内实现精确导航。因此，惯导系统具有隐蔽性好、抗干扰、不受任何气象条件限制的优点，且数据更新速率高，可以提供连续实时的导航参数。 惯性系统在国防科学技术中占有非常重要的地位，因而是世界各工业强国重点发展的技术领域之一。随着惯性技术的不断发展，许多国家已将其应用领域扩大到现代化交通运输，海洋开发，大地测量与勘探，石油钻井，矿井、隧道的掘进与贯通，机器人控制，现代化医疗器械，摄影技术以及森林防护，农业播种、施肥等民用领域。惯性技术的发展表明：从传统的机械转子型陀螺向固态陀螺仪(激光、光纤陀螺仪)转移，并进一步向以半导体硅为基本材料的微机械振动陀螺发展；从框架式平台系统向捷联系统转移，从纯惯性捷联系统向以惯性系统为基础的多体制组合导航系统发展，成为今后惯性技术发展的总趋势。 捷联式惯性导航系统，导航用的加速度计是直接捆绑在运载体上，它测量的是运载体坐标系轴向比力，只要把这个比力转换到惯性坐标系上，则其他计算就和空间稳定的平台式惯性导航系统一样，而比力转换的关键就是要实时地进行姿态基准计算来提供数学平台，即实时更新姿态矩阵b C，有些资料上称姿态矩阵 g 为捷联矩阵或方向余弦矩阵b C。一般选择地理坐标系为导航坐标系，那么捷联 g 矩阵也可表示为i C，其导航原理图如下所示： g

捷联惯导详细讲解

捷联惯导系统从20世纪60年代初开始发展起来，在1969年，捷联惯导系统作为"阿波罗"-13号登月飞船的应急备份装臵，在其服务舱发生爆炸时将飞船成功地引导到返回地球的轨道上时起到了决定性作用，成为捷联式惯导系统发展中的一个里程碑。 捷联式惯性导航(strap-downinertialnavigation)，捷联（strap-down）的英语原义是“捆绑”的意思。因此捷联式惯性导航也就是将惯性测量元件（陀螺仪和加速度计）直接装在导弹需要诸如姿态、速度、航向等导航信息的主体上，用计算机把测量信号变换为导航参数的一种导航技术。 一、捷联惯导系统工作原理及特点 惯导系统基本工作原理是以牛顿力学定律为基础，通过测量载体在惯性参考系的加速度，将它对时间进行积分，之后将其变换到导航坐标系，得到在导航坐标系中的速度、偏航角和位臵信息等。 捷联惯导系统（SINS）是一种无框架系统，由三个速率陀螺、三个线加速度计和微型计算机组成。由于惯性元器件有固定漂移率，会造成导航误差，因此导弹通常采用指令、GPS或其组合等方式对惯导进行定时修正，以获取持续准确的位臵参数。如采用指令+捷联式惯导 捷联惯导系统能精确提供载体的姿态、地速、经纬度等导航参数，是利用惯性敏感器、基准方向及最初的位臵信息

来确定运载体的方位、位臵和速度的自主式航位推算导航系统。在工作时不依赖外界信息，也不向外界辐射能量，不易受到干扰破坏。它完全是依靠载体自身设备独立自主地进行导航，它与外界不发生任何光、声、磁、电的联系，从而实现了与外界条件隔绝的假想的“封闭”空间内实现精确导航。所以它具有隐蔽性好，工作不受气象条件和人为的外界干扰等一系列的优点。 除此以外捷联惯导系统的最大特点是没有实体平台，即将陀螺仪和加速度计直接安装在机动载体上，在计算机中实时的计算姿态矩阵，通过姿态矩阵把导航加速度计测量的载体沿机体坐标系轴向的加速度信息变换到导航坐标系，然后进行导航计算。同时，从姿态矩阵的元素中提取姿态和航向信息．由此可见，在捷联惯导系统中平台的作用已由计算机及其软件的作用代替了，捷联式惯导系统采用的是数学平台。力学编排就是按照合适的数学模型由观测量计算出导航定位参数。具体地讲，利用陀螺仪测得的载体相对于惯性参照系的旋转角速度，计算出载体坐标系至导航计算坐标系之问的坐标转换矩阵；将测量的比力(加速度计测量载体相对于惯性空间的线加速度)变换至导航坐标系，并经过两次积分得到所需的速度位臵信息。 二、捷联惯导系统有以下独特优点： （1）去掉了复杂的平台机械系统，系统结构极为简单，

捷联惯导的解算程序

%====本程序为捷联惯导的解算程序(由惯性器件的输出解算出飞行器的位置、速度、姿态信息)====== clear all; close all; clc; deg_rad=pi/180; %由度转化成弧度 rad_deg=180/pi; %由弧度转化成度 %-------------------------------从源文件中读入数据 ---------------------------------- fid_read=fopen('IMUout.txt','r'); %path1_Den.dat 是由轨迹发生器产生的数据 [AllData NumofAllData]=fscanf(fid_read,'%g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g',[17 inf]); AllData=AllData'; NumofEachData=round(NumofAllData/17); Time=AllData(:,1); longitude=AllData(:,2); %经度单位：弧度 latitude=AllData(:,3); %纬度单位：弧度 High=AllData(:,4); %高度单位：米 Ve=-AllData(:,6); % 东向、北向、天向速度单位：米/妙 Vn=AllData(:,5); Vu=AllData(:,7); fb_x=AllData(:,9); %比力（fx，fy，fz） fb_y=AllData(:,8); %指向右机翼方向为x正方向，指向机头方向为y正向，z轴与x轴和y轴构成右手坐标系单位：米/秒2 fb_z=-AllData(:,10); %右前上 pitch=AllData(:,11); %俯仰角（向上为正）单位:弧度 head=-AllData(:,13); %偏航角（偏西为正） roll=AllData(:,12); %滚转角（向右为正） omigax=AllData(:,15); %陀螺输出（单位：弧度/秒，坐标轴的定义与比力的相同） omigay=AllData(:,14); omigaz=-AllData(:,16); %-------------------------------程序初始化 -------------------------------------- latitude0=latitude(1); longitude0=longitude(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所示。

捷联惯性技术的发展及与平台惯导系统的对比

捷联惯性技术的发展及与平台惯导系统的对比 [2009-06-20] 作者：admin 来源： 1.惯性技术与惯性导航的概述 惯性技术是惯性导航技术、惯性制导技术、惯性仪表技术、惯性测量技术以及惯性测试设备和装置技术的统称。它已有四十多年的发展历史了。由于惯性技术的自主性等特点，它不需要引人外界信息便可实现制导于导航。所以，它在国防科技中占有非常重要的地位，广泛的运用于航天、航空、航海等军事领域；随着惯性技术和计算机技术的不断发展以及成本降低，许多国家将其应用领域扩大到现代化交通运输、海洋开发、大地测量与勘探、石油钻井、矿井、隧道的掘进与贯通、机器人控制、现代化医疗器械、摄影技术以及森林防护、农业播种、施肥等民用领域。 惯性导航系统(Inertial Navigation System)，简称惯导，是利用惯性敏感元件、基准方向及最初的位置信息来确定运载体的方位、姿态和速度的自主式航位推算系统。惯性导航系统可以分为平台式惯导系统和捷联式惯导系统两大类：平台式惯导系统是将陀螺仪和加速计安装在一个稳定平台上，以平台坐标系为基准，测量运载体运动参数的惯性导航系统；捷联式惯导系统(Strapdown Inertial Navigation System , SI )是将惯性敏感元件(陀螺仪和加速计)直接安装在运载体上，是一种不再需要稳定平台或常平架系统的惯性导航系统。 导航的目的就是为了得到运载体的实时的方位、姿态和速度。在工程运用中，能够测定物体运动参数的方法很多：如测量位移可以用里程计，还可以用无线电定位技术、天文定位技术和卫星定位技术等；要测速度可以用测速计；要测转角可用角位置传感器(电位计、光电码盘等等)；要测角速度可以用转速表、测速电机等等。但是，以上各种测量手段还没有一种能够在同一时刻单独实时而又高精度地测量运载体的线运动和角运动，而惯性技术恰是测量这些运动参数的最理想的手段。

捷联惯导系统性能分析

捷联惯导系统性能分析 摘要 本文简要介绍捷联式惯性导航系统的各种分析技术,对捷联惯导系统算法验证的过程进行了讨论。封闭形式分析仿真驱动程序，可以用来锻炼/验证捷联算法方程。分析的精度捷联络筒，划船和分析方法位置融合算法（包括位置算法的折叠效果）函数的算法重复率和系统振动输入。包括的是一个简化的分析模型，该模型描述可用于翻译系统到惯性传感器作为传感器组件的函数的输入的振动安装失衡。捷联系统静态漂移和的旋转测试程序/方程描述捷联式传感器的校准系数确定。该文件概述了卡尔曼滤波器的设计和协方差分析技术，并介绍了验证辅助捷联惯导系统的一般步骤卡尔曼滤波器的配置。最后，论述了系统集成测试的一般过程验证所有的硬件，软件，系统功能操作进行正确和准确和界面元素。 坐标框架 在本文中，使用一个坐标系是一个分析性的抽象定义的三个相互垂直单位矢量。一个坐标系可以看作一组三个相互垂直的线（轴）通过一个共同的点（原点）与来自沿着坐标轴的原点上的单位向量。在本文中，每个坐标系的原点的物理位置是任意的。主坐标帧利用有以下几种： 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 =时间。 一、简介

捷联惯导Matlab程序求解飞行器的姿态

捷联惯导程序,依据加表与陀螺仪的输出数据来求解飞行器的姿态 clc; clear; format long; %设置数据精度为15位小数 Data=importdata('temp、txt'); % 导入实验所采集的数据,以矩阵形式赋给Data变量,temp、txt 必须与该M文件在同一个文件夹中 Px=Data(:,3); % Px,Py,Pz为陀螺仪的输出值 Py=Data(:,4); Pz=Data(:,5); Nx=Data(:,6); % Nx,Ny,Nz为加速度计的输出值 Ny=Data(:,7); Nz=Data(:,8); % 陀螺仪模型参数标定如下: Sx = -4、085903e-006 ; Sy = -4、085647e-006 ; Sz = -4、085170e-006 ; Mxy = 5、059527e-003 ; Mxz = -1、031103e-003 ; Myx = -3、355451e-003 ; Myz = 3、508468e-003 ; Mzx = -1、266671e-003 ; Mzy = -2、318244e-004 ; Dx = -2、009710e-006 ; Dy = 8、156346e-007 ; Dz = -5、749059e-007 ; GyroCali_A = [ 1 -Mxy -Mxz ; -Myx 1 -Myz ; -Mzx -Mzy 1 ]; % 加速度计模型参数标定如下: Kx = 9、272930e-004 ; Ky = 9、065544e-004 ; Kz = 9、443748e-004 ; Ixy = 6、533872e-003 ; Ixz = 9、565992e-004 ; Iyx = -6、319376e-003 ; Iyz = -6、902339e-004 ; Izx = -1、144549e-003 ; Izy = -3、857963e-004 ; Bx = -3、400847e-002 ; By = -8、916341e-003 ; Bz = -9、947414e-003 ; AccCali_A = [1 -Ixy -Ixz ; -Iyx 1 -Iyz ; -Izx -Izy 1 ]; Delta_t = 0、05; %采样时间为0、05秒 Delta_Theta_x = 0; Delta_Theta_y = 0; Delta_Theta_z = 0; %定义陀螺仪输出的角度增量 Delta_Vx = 0; Delta_Vy = 0; Delta_Vz = 0; %定义加速度计输出的速度增量 L = zeros(1,12001); L(1)= 45、7328*pi/180 ; %纬度用L表示,纬度的初始值划为弧度形式,因为后面计算位置矩阵更新 L(2)= 45、7328*pi/180 ; %时需要用到前两次的L值来计算当前L值,所以在此定义2个初始L值 Lamda = 126、6287*pi/180 ; %经度用Lamda表示,经度的初始值划为弧度形式 h = 136 ; %高度用h表示 V = [ 0 ; 0 ; 0 ]; %导航坐标系中的东北天初始速度都为0 Vx = 0; %方便后面的速度计算与速度更新 Vy = 0; Vz = 0; Theta = 0; Gama = 0;