卡尔曼滤波实验报告 捷联惯导静基座初始对准实验 一、实验目的 ①掌握捷联惯导的构成和基本工作原理； ②掌握捷联惯导静基座对准的基本工作原理； ③了解捷联惯导静基座对准时的每个系统状态的可观测性； ④了解双位置对准时系统状态的可观测性的变化。 二、实验原理 选取状态变量为：[]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）是在平台式惯导系统基础上发展而来的，它是一种无框架系统，由三个速率陀螺、三个线加速度计和微型计算机组成。平台式惯导系统和捷联式惯导系统的主要区别是：前者有实体的物理平台，陀螺和加速度计置于陀螺稳定的平台上，该平台跟踪导航坐标系，以实现速度和位置解算，姿态数据直接取自于平台的环架；后者的陀螺和加速度计直接固连在载体上作
捷联惯导程序，依据加表和陀螺仪的输出数据来求解飞行器的姿态 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表示，经度的初始值划为弧度形式
捷联惯导系统快速罗经初始对准方法研究1 严恭敏1，严卫生1,2，徐德民1,2 1西北工业大学航海学院，西安（710072） 2水下信息处理与控制国家级重点实验室，西安（710072） E-mail：yangongmin@https://www.doczj.com/doc/7e5143467.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)=[ %====本程序为捷联惯导的解算程序(由惯性器件的输出解算出飞行器的位置、速度、姿态信息)====== 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); %初始位置 捷联式惯导系统初始对准 惯性技术是惯导(惯性导航与惯性制导)技术、惯性仪表技术、惯性测量技术以及有关设备和装置技术的统称。惯性导航与惯性制导是当今非常重要的综合技术之一，它广泛用于航空、航海、航天及陆地各领域。惯性导航系统是和用陀螺与加速度计通过最初的方向基准和位置信息来确定运载体在一特定坐标系内的姿态、位置、速度和加速度的自主式导航系统。惯性制导系统是利用运载体内部的陀螺、加速度计测量其运动参数，经过计算机发出控制指令，从而把运载体按照预定的路线准确地引导到目的地的制导系统。自主性是惯性系统最重要的特点。确定运动对象导航参数的方法和仪器有许多，例如磁、天文、无线电、水声、全球卫星定位系统等等，然而它们都有一个致命的弱点，即不是自主的，不是要向外界发出信息，就是要依赖对外观测信息，而惯性系统与上述诸方法的基本区别就在于是完全自主的，即导弹、潜艇、飞船等可以在一个完全与外界条件以及电磁波隔绝的假想“封闭”空间内实现精确导航。因此，惯导系统具有隐蔽性好、抗干扰、不受任何气象条件限制的优点，且数据更新速率高，可以提供连续实时的导航参数。 惯性系统在国防科学技术中占有非常重要的地位，因而是世界各工业强国重点发展的技术领域之一。随着惯性技术的不断发展，许多国家已将其应用领域扩大到现代化交通运输，海洋开发，大地测量与勘探，石油钻井，矿井、隧道的掘进与贯通，机器人控制，现代化医疗器械，摄影技术以及森林防护，农业播种、施肥等民用领域。惯性技术的发展表明：从传统的机械转子型陀螺向固态陀螺仪(激光、光纤陀螺仪)转移，并进一步向以半导体硅为基本材料的微机械振动陀螺发展；从框架式平台系统向捷联系统转移，从纯惯性捷联系统向以惯性系统为基础的多体制组合导航系统发展，成为今后惯性技术发展的总趋势。 捷联式惯性导航系统，导航用的加速度计是直接捆绑在运载体上，它测量的是运载体坐标系轴向比力，只要把这个比力转换到惯性坐标系上，则其他计算就和空间稳定的平台式惯性导航系统一样，而比力转换的关键就是要实时地进行姿态基准计算来提供数学平台，即实时更新姿态矩阵b C，有些资料上称姿态矩阵 g 为捷联矩阵或方向余弦矩阵b C。一般选择地理坐标系为导航坐标系，那么捷联 g 矩阵也可表示为i C，其导航原理图如下所示： g 第二章 捷联惯导系统的初试对准 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为捷联式惯导系统的原理图，陀螺仪和加速度计输出分别送入姿态矩阵计算和由载体坐标系至平台坐标系的方向余弦矩阵的计算。有了姿态矩阵，其一可以实现把载体坐标系轴向加速度信息变换到导航坐标系轴，进而可以进行所需的导航参数计算，其二利用姿态矩阵的元素，提取方位和姿态信息。 图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) 车载捷联惯导系统基本原理 一、捷联惯导系统基本原理 捷联惯导系统基本原理如图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)可得四元数姿态更新算法为： 光纤陀螺捷联惯导多位置系统级标定方法 标定技术是捷联惯导系统领域的一项关键技术,本质上是一种误差补偿技术。按照观测量的不同可以分为分立式标定和系统级标定,由于前者在研究方法上比后者更加趋于成熟,因此常常采用分立式标定来确定惯组的数学模型。分立式标定依靠高精度三轴转台来进行一系列实验来得出陀螺仪与加速计模型中的各项 参数,实验步骤多、标定时间长,且标定精度受转台限制无法进一步提升。为使标定更有效率,越来越多的学者转向系统级标定方法的研究。 本文介绍了两种系统级标定的一般方法,分别是滤波法和解析法,代表了系 统级标定的两种不同思路。相较于解析法,运用滤波的方法来求解标定模型不仅原理简单,而且不涉及复杂的数学推导,精度高,标定路径编排较少。因此本课题以卡尔曼滤波为基础来进行系统级标定方案的设计。直接将惯组的待标定参数扩充到状态量中,得到的是一个30维的大系统。 为降低试验设计难度,提高滤波精度,本文对其进行降维处理。通过引入了一种基于行初等变换的可观测性分析方法,运用该方法可以直观地识别出状态方程中哪些状态独立可观测、哪些状态不可观测、哪些状态非独立可观测。基于可观测性的分析结果,在原有的状态方程中,保留独立可观测量,剔除不可观测量、合并非独立可观测量,从而达到简化系统方程,降低系统维数的目的。最后依据降维后的系统,设计出了6个静态位置与6个单轴旋转位置,并对此进行了仿真试验,仿真试验结果表明该方案具有理论可行性。 之后进行了实际转台试验,将其与分立式标定结果进行比较。通过对比四组五级海况实验的姿态变化曲线,可以清晰的看出系统级方案的补偿效果整体上优于分立式标定,故而这种基于降维后的系统级标定滤波方案具有一定的工程实践价值。 哈尔滨工业大学《导航原理》课程 捷联惯导程序，依据加表和陀螺仪的输出数据来求解飞行器的姿态 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值 捷联惯导系统从20世纪60年代初开始发展起来，在1969年，捷联惯导系统作为"阿波罗"-13号登月飞船的应急备份装臵，在其服务舱发生爆炸时将飞船成功地引导到返回地球的轨道上时起到了决定性作用，成为捷联式惯导系统发展中的一个里程碑。 捷联式惯性导航(strap-downinertialnavigation)，捷联（strap-down）的英语原义是“捆绑”的意思。因此捷联式惯性导航也就是将惯性测量元件（陀螺仪和加速度计）直接装在导弹需要诸如姿态、速度、航向等导航信息的主体上，用计算机把测量信号变换为导航参数的一种导航技术。 一、捷联惯导系统工作原理及特点 惯导系统基本工作原理是以牛顿力学定律为基础，通过测量载体在惯性参考系的加速度，将它对时间进行积分，之后将其变换到导航坐标系，得到在导航坐标系中的速度、偏航角和位臵信息等。 捷联惯导系统（SINS）是一种无框架系统，由三个速率陀螺、三个线加速度计和微型计算机组成。由于惯性元器件有固定漂移率，会造成导航误差，因此导弹通常采用指令、GPS或其组合等方式对惯导进行定时修正，以获取持续准确的位臵参数。如采用指令+捷联式惯导 捷联惯导系统能精确提供载体的姿态、地速、经纬度等导航参数，是利用惯性敏感器、基准方向及最初的位臵信息 来确定运载体的方位、位臵和速度的自主式航位推算导航系统。在工作时不依赖外界信息，也不向外界辐射能量，不易受到干扰破坏。它完全是依靠载体自身设备独立自主地进行导航，它与外界不发生任何光、声、磁、电的联系，从而实现了与外界条件隔绝的假想的“封闭”空间内实现精确导航。所以它具有隐蔽性好，工作不受气象条件和人为的外界干扰等一系列的优点。 除此以外捷联惯导系统的最大特点是没有实体平台，即将陀螺仪和加速度计直接安装在机动载体上，在计算机中实时的计算姿态矩阵，通过姿态矩阵把导航加速度计测量的载体沿机体坐标系轴向的加速度信息变换到导航坐标系，然后进行导航计算。同时，从姿态矩阵的元素中提取姿态和航向信息．由此可见，在捷联惯导系统中平台的作用已由计算机及其软件的作用代替了，捷联式惯导系统采用的是数学平台。力学编排就是按照合适的数学模型由观测量计算出导航定位参数。具体地讲，利用陀螺仪测得的载体相对于惯性参照系的旋转角速度，计算出载体坐标系至导航计算坐标系之问的坐标转换矩阵；将测量的比力(加速度计测量载体相对于惯性空间的线加速度)变换至导航坐标系，并经过两次积分得到所需的速度位臵信息。 二、捷联惯导系统有以下独特优点： （1）去掉了复杂的平台机械系统，系统结构极为简单，捷联惯导的解算程序
捷联式惯导系统初始对准
2捷联惯性导航系统初始对准原理
捷联惯导作业
车载捷联惯导系统基本原理
光纤陀螺捷联惯导多位置系统级标定方法
Matlab捷联惯导
捷联惯导详细讲解
相关主题
文本预览