本发明公开了一种动态变形下基于双重滤波器的传递对准方法，动态变形下产生的动态变形角以及动态变形和机体运动之间的耦合角会降低传递对准的精度；将传递对准滤波器分为两部分，第一部分针对弯曲变形角和耦合角进行估计，采用姿态匹配方法；第二部分针对动态杠杆臂进行估计，采用“速度+角速度”匹配方法。相比现有的传递对准滤波器的设计，本发明将耦合角作为滤波器的状态量，在滤波器一中对耦合角进行估计，并用于滤波器二补偿杠杆臂误差，两个滤波器同步进行，在最后一步进行信息融合，估计出子惯导系统的初始姿态误差，进行初始姿态校准。本发明不仅提高了传递对准系统的精度，同时缩短了传递对准过程的时间。
1.一种动态变形下基于双重滤波器的传递对准方法，应用于飞机机翼变形测量系统中，其中主惯导系统安装在机舱，子惯导系统安装在机翼，其特征在于，包括以下步骤：(1)用轨迹发生器产生主惯导系统的姿态、速度和位置信息以及陀螺仪和加速度计的输出，用二阶马尔科夫模拟主惯导系统、子惯导系统之间的弯曲变形角 和弯曲变形角速度 对弯曲变形进行几何分析，推导出由载体动态变形和载体运动所引起的耦合角 表达式；
(2)将弯曲变形角、弯曲变形角速度和耦合角作为状态量，采用姿态匹配方法，建立滤波器一模型；
(3)利用步骤(2)中估计的弯曲变形角和耦合角建立动态杠杆臂模型，推导速度误差表达式和角速度误差表达式；
(4)利用步骤(3)推导的速度误差表达式和角速度误差表达式，采用“速度+角速度”匹配方法，建立滤波器二的模型，估计子惯导系统的初始姿态误差，并将此误差用于子惯导系统的初始姿态校准，完成传递对准过程。
2.根据权利要求1所述的一种动态变形下基于双重滤波器的传递对准方法，其特征在于，所述步骤(1)中，对弯曲变形进行几何分析，推导出由载体动态变形和载体运动所引起的耦合角 表达式为：其中， M表示为：
其中， 分别表示东、北、天三个方向下子惯导系统角速度理想值。
3.根据权利要求2所述的一种动态变形下基于双重滤波器的传递对准方法，其特征在于，所述步骤(2)中，将弯曲变形角、弯曲变形角速度和耦合角作为状态量，采用姿态匹配方法建立滤波器一模型，具体如下：选取滤波器一的状态量为：
其中， 表示姿态误差， 表示子惯导系统陀螺仪测量零漂， 表示主惯导系统、子惯导系统之间的初始安装角误差；
其中，F1表示滤波器一状态转移矩阵，G1表示滤波器一系统噪声分配矩阵，w1表示滤波器一系统噪声，状态转移矩阵F1表示为：其中， 表示导航系相对于惯性系的旋转， 表示反对称矩阵， 表示子惯导系统理想坐标系与导航坐标系之间的转换矩阵， βi(i＝x,y,z)表示东、北、天三个方向上二阶马尔科夫模型的系数，F64＝MB2，F65＝MB1；
其中，y1表示姿态真实值与滤波器估计值的差值，H1表示滤波器一量测矩阵，μ1表示滤波器一量测噪声。
4.根据权利要求2所述的一种动态变形下基于双重滤波器的传递对准方法，其特征在于，所述步骤(3)中，推导出的速度误差表达式和角速度误差表达式，具体如下：角速度误差表达式为：
其中， 表示主惯导系统坐标系下主惯导系统的角速度，表示主惯导系统、子惯导系统之间的理想误差角， 表示幅值矩阵， 表示 方向上的单位矩阵， 表示子惯导系统坐标系下子惯导系统的角速度，U＝[1 1 1]T；
其中， 为地球自转引起的导航系旋转， 为子系统在地球表面移动因地球表面弯曲引起的导航系的旋转， 分别表示主惯导系统、子惯导系统在导航坐标系下的速度矢量，为主惯导系统、子惯导系统之间的弯曲变形角， 为主惯导系统、子惯导系统之间的耦合角， 表示导航系下主惯导系统的角速度， 表示子惯导系统与导航坐标系之间的转换矩阵， 表示子惯导系统加速度计测量零偏， 表示子惯导系统在导航坐标系下的比力， 表示动态杠杆臂， 表示静止状态下杠杆臂，x0 y0 z0分别表示东、北、天三个方向的静止状态下杠杆臂，R0可表示为：
5.权利要求4所述的一种动态变形下基于双重滤波器的传递对准方法，其特征在于，所述步骤(4)中，滤波器二采用“速度+角速度”匹配，利用步骤(3)推导的速度误差表达式和角速度误差表达式建立量测量方程，建立卡尔曼滤波器模型，具体如下：选取卡尔曼滤波器二的状态量为：
其中， 表示速度误差， 表示子惯导系统加速度计测量零偏；
其中，G2表示滤波器二系统噪声分配矩阵，w2表示滤波器二系统噪声，状态转移矩阵F2表示为：其中， F13＝R0B2+R0M(B1B2+B2),
其中，y2表示速度、角速度真实值与滤波器估计值的差值，μ2表示滤波器二量测噪声，
一种动态变形下基于双重滤波器的传递对准方法
技术领域
[0001] 本发明属于惯性导航技术领域，利用惯性导航系统测量飞机的机翼变形，其中涉及高精度主惯导系统对低精度子惯导系统进行校准的过程，具体涉及一种动态变形下基于双重滤波器的传递对准方法。
背景技术
[0002] 飞机的承载能力有限，特别是机翼部分，因此飞机机翼动态变形测量对测量设备的重量和尺寸有非常严格的要求，而IMU单元的测量精度与重量和尺寸成正比，每个负载处无法同时安装高精度的IMU。
[0003] 目前飞机机翼变形测量采用机身安装高精度的POS，而机翼部分则采用低精度IMU单元，通过主、子系统间传递对准获取各定位点的高精度位置、姿态信息。但是主、子之间挠曲变形产生的附加速度、角速度和角度是影响其精度的主要因素，现有的飞机机翼动态变形测量将机翼视为刚体，不考虑挠曲变形，其传递对准精度难以达到所需要的精度。
发明内容
[0004] 发明目的：针对现有技术的不足，本发明目的在于提供一种动态变形下基于双重滤波器的传递对准方法，对飞机机翼动态变形测量传递对准过程中机体运动和动态变形之间的耦合所引起的误差角度和角速度进行几何建模和数学分析，推导出耦合角度和角速度的表达式，并将传递对准滤波器分为两部分，第一部分针对弯曲变形角和耦合角进行估计，采用姿态匹配方法；第二部分针对动态杠杆臂进行估计，采用“速度+角速度”匹配方法，此设计在提高传递对准精度的前提下，同时缩短了传递对准过程的时间。
[0005] 技术方案：为实现上述发明目的，本发明采用的技术方案如下：
[0006] 一种动态变形下基于双重滤波器的传递对准方法，应用于飞机机翼变形测量系统中，其中主惯导系统安装在机舱，子惯导系统安装在机翼，该方法包括以下步骤：
[0007] (1)用轨迹发生器产生主惯导系统的姿态、速度和位置信息以及陀螺仪和加速度计的输出，用二阶马尔科夫模拟主、子惯导之间的弯曲变形角 和弯曲变形角速度 对弯曲变形进行几何分析，推导出由载体动态变形和载体运动所引起的耦合角度 表达式；
[0008] (2)将弯曲变形角、弯曲变形角速度和耦合角作为状态量，采用姿态匹配方法，建立滤波器一模型；
[0009] (3)利用步骤(2)中估计的弯曲变形角和耦合角建立动态杠杆臂模型，推导速度误差表达式和角速度误差表达式；
[0010] (4)利用步骤(3)推导的速度误差表达式和角速度误差表达式，采用“速度+角速度”匹配方法，建立滤波器二的模型，估计子惯导系统的初始姿态误差，并将此误差用于子惯导系统的初始姿态校准，完成传递对准过程。
[0011] 进一步地，所述步骤(1)中，对弯曲变形进行几何分析，推导出由载体动态变形和载体运动所引起的耦合角度 表达式为：
[0012]
[0013] 其中， M表示为：
[0014]
[0015] 其中， 分别表示东、北、天三个方向下子惯导系统角速度理想值。
[0016] 进一步地，所述步骤(2)中，将弯曲变形角、弯曲变形角速度和耦合角作为状态量，采用姿态匹配方法建立滤波器一模型，具体如下：
[0017] 选取滤波器一的状态量为：
[0018]
[0019] 其中， 表示姿态误差， 表示子系统陀螺仪测量零漂， 表示主、子系统之间的初始安装角误差；
[0020] 滤波器一的状态方程为：
[0021]
[0022] 其中，F1表示滤波器一状态转移矩阵，G1表示滤波器一系统噪声分配矩阵，w1表示滤波器一系统噪声，状态转移矩阵F1表示为：
[0023]
[0024] 其中， 表示导航系相对于惯性系的旋转， 表示反对称矩阵， 表示子系统理想坐标系与导航坐标系之间的转换矩阵，βi(i＝x,y,z)表示东、北、天三个方向上二阶马尔科夫模型的系数，F64＝MB2，F65＝MB1；
[0025] 系统量测方程为：
[0026] y1＝H1x1+μ1
[0027] 其中，y1表示姿态真实值与滤波器估计值的差值，H1表示滤波器一量测矩阵，μ1表示滤波器一量测噪声。
[0028] 进一步地，所述步骤(3)中，推导出的速度误差表达式和角速度误差表达式，具体如下：
[0029] 角速度误差表达式为：
[0030]
[0031] 其中， 表示主系统坐标系下主系统的角速度，表示主、子系统之间的理想误差角， 表示幅值矩阵， 表示 方向上的单位矩阵， 表示子系统坐标系下子系统的T角速度，U＝[1 1 1]；
[0032] 速度误差表达式为：
[0033]
[0034] 其中， 为地球自转引起的导航系旋转， 为子系统在地球表面移动因地球表面弯曲引起的导航系的旋转， 分别表示主、子系统在导航坐标系下的速度矢量，为主、子惯导之间的弯曲变形角， 为主、子惯导之间的耦合角， 表示导航系下主系统的角速度， 表示子系统与导航坐标系之间的转换矩阵， 表示子系统加速度计测量零偏， 表示子系统在导航坐标系下的比力， 表示动态杠杆臂， 表示静止状态下杠杆臂，x0 y0 z0分别表示东、北、天三个方向的静止状态下杠杆臂，R0可表示为：
[0035]
[0036] 进一步地，所述步骤(4)中，滤波器二采用“速度+角速度”匹配，利用步骤(3)推导的速度误差表达式和角速度误差表达式建立量测量方程，建立卡尔曼滤波器模型，具体如下：
[0037] 选取卡尔曼滤波器二的状态量为：
[0038]
[0039] 其中， 表示速度误差， 表示子系统加速度计测量零偏；
[0040] 滤波器的状态方程为：
[0041]
[0042] 其中，G2表示滤波器二系统噪声分配矩阵，w2表示滤波器二系统噪声，状态转移矩阵F2表示为：
[0043]
[0044] 其中， F13＝R0B2+R0M(B1B2+B2),F53＝R0MB2,F54＝R0+R0MB1，系统量测方程为：
[0045] y2＝H2x2+μ2
[0046] 其中，y2表示速度、角速度真实值与滤波器估计值的差值，μ2表示滤波器二量测噪声，
[0047] 有益效果：与现有技术相比，本发明考虑了载体运动主、子系统之间刚体运动和动态弹性形变耦合误差，对动态弹性变形下主、子系统之间的角度和角速度误差进行空间几何建模和数学分析，得出动态形变下主、子系统之间的耦合角度误差，由此推导出动态形变下主、子系统之间的角速度误差表达式，并采用双重滤波器的方法，两滤波器同步进行，在最后一步进行融合；传统的传递对准过程，不考虑动态形变和机体运动之间的耦合误差，传递对准精度无法达到高精度传递对准的要求，另一方面，采用24维的滤波器，计算量大，本发明对主、子系统之间的耦合角度进行几何分析，得出耦合角度的表达式，并将状态量分为两组，分别在两个滤波器中同步进行，此设计在提高传递对准精度的前提下，同时缩短了传递对准过程的时间。
附图说明
[0048] 图1为本发明基于双重滤波器的传递对准流程图；
[0049] 图2为角速度矢量与附加动态弯曲角速度矢量之间的空间关系示意图；
[0050] 图3为动态变形下主、子惯导之间耦合角度(投影到yoz平面)示意图；
[0051] 图4为主、子系统相对位置关系示意图。
具体实施方式
[0052] 以下结合具体的实施方案和附图对本发明作进一步详细说明：
[0053] 如图1所示，本发明实施提出的一种动态变形下基于双重滤波器的传递对准方法，用轨迹模拟器模拟飞机主系统的姿态、速度、位置和惯性器件的输出数据，同时采用二阶马尔科夫模拟输出主、子系统之间的弯曲变形角 和弯曲变形角速度 对载体运动和挠曲变形解耦合，得到耦合角并将其作为滤波器一的状态量，采用姿态匹配；滤波器二利用滤波器一的结果对杠杆臂误差进行补偿，并采用速度+角速度匹配方法。下面进行详细的分析：
[0054] 步骤1：轨迹发生器产生主惯导系统的姿态、速度和位置信息以及惯性器件(陀螺仪和加速度计)的输出，用二阶马尔科夫模拟主惯导与子惯导之间的弯曲变形角 和弯曲变形角速度 并对弯曲变形进行几何分析，推导出由主、子系统之间动态变形所引起的主、子系统之间的耦合角度 主、子系统之间弯曲变形角 可用二阶马尔科夫表示为：
[0055]
[0056] 其中，β＝2.146/τ，τ表示相关时间， 表示高斯白噪声，主、子系统之间的理想误差角矢量 表示为：
[0057]
[0058]
[0059] 其中， 表示理想状态下子系统陀螺仪的输出， 表示主系统陀螺仪的输出，表示主、子系统之间的姿态矩阵， 表示主系统初始安装误差角矢量， 表示弯曲变形角，由于动态弯曲变形的作用，产生了附加的角速度 其可表示为 则如图2所示，实际状态下子系统的角速度输出 可表示为：
[0060]
[0061] 取弯曲变形耦合角速度所引起的主、子系统之间的耦合误差角矢量为[0062]
[0063] 下标x,y,z分别表示东、北、天三个方向， 为弯曲变形耦合角速度所引起的主、子系 统之 间的 耦合 误差 角 ，即 与 的 夹 角 ，取则有：
[0064]
[0065] 如图3所示，由几何关系有：
[0066]
[0067] 用泰勒级数将反正切函数arctan展开，并略去高次项，可得：
[0068]
[0069] 其中，M可表示为：
[0070]
[0071] 步骤2：将弯曲变形角、弯曲变形角速度和耦合角作为状态量，采用姿态匹配方法建立滤波器一的模型，具体如下：
[0072] 选取滤波器一的状态量为：
[0073]
[0074] 其中， 表示姿态误差， 表示子系统陀螺仪测量零漂， 表示主、子系统之间的初始安装角误差；
[0075] 滤波器一的状态方程为：
[0076]
[0077] 其中，F1表示状态转移矩阵，G1表示系统噪声分配矩阵，w1表示系统噪声，根据步骤(1)中得到的耦合角模型，状态转移矩阵F1可表示为：
[0078]
[0079] 其中， 表示导航系相对于惯性系的旋转， 表示反对称矩阵， 表示子系统理想坐标系与导航坐标系之间的转换矩阵，βi(i＝x,y,z)表示x,y,z三个方向上二阶马尔科夫模型的系数，F64＝MB2，F65＝MB1；
[0080] 系统量测方程为：
[0081] y1＝H1x1+μ1
[0082] 其中，y1表示姿态真实值与滤波器估计值的差值，H1表示量测矩阵，H矩阵的具体表达式见文献“Multi-node Transfer Alignment based on Mechanics Modeling for Airborne DPOS”，μ1表示滤波器一量测噪声；
[0083] 步骤3：推导速度误差表达式和角速度误差表达式，具体如下：
[0084] 主、子系统之间的角速度之差 可表示为：
[0085]
[0086]
[0087] 其中，主、子系统之间的耦合误差角矢量为 则主、子系统之间的转换矩阵可以表示为 主、子系统之间的误差角速度可以表示为：
[0088]
[0089] 其中， 为 在 上的投影，由于主系统到子系统之间的旋转矢量 为小量，则有
[0090]
[0091] 故有：
[0092]
[0093] 其中， 表示反对称矩阵，取 为：
[0094]
[0095] 其中， 表示幅值矩阵， 表示 方向上的单位矩阵， 表示 与 之间的夹角矢量， 表示由 到 之间的转换矩阵， 其中，U＝[1 1 1]T，且有：
[0096]
[0097] 符号| |表示求模，将 代入 的表达式，则有：
[0098]
[0099] 将 代入 的表达式，则有：
[0100]
[0101] 主、子系统之间的位置关系如图4所示，表达式为：
[0102]
[0103] 其中， 分别表示主、子节点到地心的矢量，表示主、子节点之间的动态杠杆臂矢量，则在惯性系下可表示为：
[0104]
[0105] 其中， 表示主系统到惯性系的转换矩阵， 表示主系统坐标系下动态杠杆臂矢量，根据牛顿第二定律，有：
[0106]
[0107]
[0108] 其中， 表示子系统在惯性系下的比力， 表示主系统在惯性系下的比力，g表示重力加速度， 表示地球自转角速度，联立上式有：
[0109]
[0110] 由于，主、子系统速度矢量微分方程可表示为：
[0111]
[0112]
[0113] 其中， 分别表示主、子系统在导航坐标系下的速度矢量， 为地球自转引起的导航系旋转， 为子系统在地球表面移动因地球表面弯曲引起的导航系的旋转， 表示导航系下主系统的角速度， 表示子系统在导航坐标系下的比力， 表示子系统加速度计测量零偏，速度误差矢量方程表示为：
[0114]
[0115] 将上式两边微分有：
[0116]
[0117] 动态杠杆臂矢量可表示为：
[0118]
[0119]
[0120] 其中， 表示静止状态下杠杆臂，x0 y0 z0分别表示东、北、天三个方向的静止状态下杠杆臂, 将上式两边微分有：
[0121]
[0122]
[0123] 将动态杠杆臂表达式代入速度误差矢量表达式中，有：
[0124]
[0125] 步骤4：传递对准滤波器二采用“速度+角速度”匹配，利用步骤3推导的速度误差表达式和角速度误差表达式建立量测量方程，建立卡尔曼滤波器模型，估计子节点的初始姿态误差，并将此误差用于子节点的初始姿态校准，完成传递对准过程。
[0126] 本步骤中选取滤波器二的状态量为：
[0127]
[0128] 其中， 表示速度误差， 表示子系统加速度计测量零偏；
[0129] 滤波器的状态方程为：
[0130]
[0131] 其中，G2表示滤波器二的系统噪声分配矩阵，w2表示滤波器二的系统噪声，F2表示为：
[0132]
[0133] 其中， F13＝R0B2+R0M(B1B2+B2),F53＝R0MB2,F54＝R0+R0MB1，系统量测方程为：
[0134] y2＝H2x2+μ2
[0135] 其中，y2表示速度、角速度真实值与滤波器估计值的差值，μ2表示滤波器二量测噪声，
