[发明专利]传递对准中考虑杆臂补偿的变形角测量方法无效
申请号: | 201310156752.2 | 申请日: | 2013-04-26 |
公开(公告)号: | CN103256942A | 公开(公告)日: | 2013-08-21 |
发明(设计)人: | 高伟;鲍桂清;奔粤阳;李倩;张勇刚;黄平;吴磊;王罡;陈世同;程建华 | 申请(专利权)人: | 哈尔滨工程大学 |
主分类号: | G01C25/00 | 分类号: | G01C25/00 |
代理公司: | 暂无信息 | 代理人: | 暂无信息 |
地址: | 150001 黑龙江*** | 国省代码: | 黑龙江;23 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | 本发明提供的是一种在传递对准中杆臂补偿下变形角的测量方法。将主、子惯导系统分别安装在船上;启动主惯导系统,使其完成初始对准后进入导航状态;测量主、子惯导系统之间的距离;启动子惯导系统,用计算机采集陀螺和加速度计的输出;将主惯导系统的速度、姿态信息传递给子惯导系统,进行对准;建立船体变形的模型,计算变形后的主、子惯导系统之间的距离;建立杆臂效应误差模型;建立卡尔曼滤波的状态方程和观测方程;用卡尔曼滤波估计出主、子惯导系统之间的安装误差角和船体的变形角。本发明中变形角的测量是在杆臂补偿下进行测量的,具有一定的实际意义。 | ||
搜索关键词: | 传递 对准 考虑 补偿 变形 测量方法 | ||
【主权项】:
1.一种在传递对准中杆臂补偿下变形角的测量方法,其特征包括以下步骤:(1)将主、子惯导系统安装在船上,并用线连接好;(2)启动主惯导系统,完成主惯导的初始对准,并使其处于导航状态;(3)测量主、子惯导系统之间的距离;(4)启动子惯导系统,用计算机采集陀螺和加速度计的输出;(5)将主惯导系统的速度信息和姿态信息传递给子惯导系统,完成对准;(6)建立船体变形的模型,计算变形之后主、子惯导系统之间的距离;船在航行时,船体会发生船体变形,其变形的模型为:θ · · + 2 β θ · + β 2 θ = η ]]> 其中:θ为船体的变形角,和分别为θ的一阶导数和二阶导数,β=1.073,η为白噪声。由于船体变形的影响,主、子惯导系统之间的距离就会发生变化,设变化之后的距离为r,则:r = 1 θ z - θ y - θ z 1 θ x θ y - θ x 1 · r p ]]> 其中:θx、θy、θz为船体变形角在三轴上的分量,rp为原来主、子惯导系统之间的距离。(7)建立杆臂效应的误差模型;a f = r · · + 2 ω ib × r · + ω · ib × r + ω ib × ( ω ib × r ) ]]> 其中:af为杆臂加速度,r为主、子惯导系统之间的距离,和分别是r的一阶导数和二阶导数,,ωib为陀螺的输出。(8)上述所建立的卡尔曼滤波器的状态方程和观测方程;选取系统的状态变量为:X = δ v E δ v N φ mx φ my φ mz φ ax φ ay φ az θ x θ y θ z θ · x θ · y θ · z T ]]> 其中:δvE、δvN分别为东向速度误差和北向速度误差,φmx、φmy、φmz分别为主、子惯导系统之间的失准角在三个轴上的分量,φax、φayφaz分别为主、子惯导系统之间的安装误差角在三个轴上的分量,θx、θy、θz为船体变形角在三轴上的分量,分别是θx、θy、θz的一阶导数,上标T表示取矩阵的转置。故该系统的滤波模型的状态方程为:X · = AX + W ]]> 则:A = A 1 A 2 - A 2 N 1 M 1 0 3 × 2 A 3 - A 3 I 3 × 3 0 3 × 3 0 3 × 2 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 2 0 3 × 3 0 3 × 3 0 3 × 3 I 3 × 3 0 3 × 2 0 3 × 3 0 3 × 3 A 4 A 5 , ]]>A 4 = - β x 2 0 0 0 - β y 2 0 0 0 - β z 2 , ]]>A 2 = c 13 f y - c 12 f z - c 13 f x + c 11 f z c 12 f x - c 11 f z c 23 f y - c 22 f z - c 23 f x + c 21 f z c 22 f x - c 21 f z , ]]>A 5 = - 2 β x 0 0 0 - 2 β y 0 0 0 - 2 β z , ]]>f s s = f x f y f z T , ]]>A 3 = 0 ω nsz s - ω nsy s - ω nsz s 0 ω nsx s ω nsy s - ω nsx s 0 , ]]>C s ′ n = c 11 c 12 c 13 c 21 c 22 c 23 c 31 c 32 c 33 , ]]>N 1 = β x 2 r z · β y 2 - r y · β z 2 - r z · β x 2 β y 2 r x · β z 2 , ]]>M 1 = - 2 β x + r z · ω ibz + r y · ω iby r z · β y - ω ibz - r x · ω iby - r y · β z - r x · ω ibz + ω iby - r z · β x + ω ibz - r y · ω ibx β y + r z · ω ibz + r x · ω ibx r x · β z - r y · ω ibz - ω ibx . ]]> 其中是状态向量X的导数,ωie为地球的自转角速度,为子惯导系统的姿态矩阵,为子惯导三轴的陀螺输出,Re为地球卯酉曲率半径,为当地的纬度,A4和A5中的系数均是挠曲变形角的相关值,为当地的纬度,W为系统噪声阵,设为白噪声。系统选取速度误差和姿态为观测量,则:Z=[δvE δvN φmx φmy φmz]T其中:δvE、δvN分别为东向速度误差和北向速度误差,φmx、φmy、φmz分别为主、子惯导系统之间的失准角在三个轴上的分量,上标T表示取矩阵的转置。故该系统的观测方程为:Z=HX+V其中:观测矩阵H = I 2 × 2 0 2 × 3 0 2 × 3 0 2 × 3 0 3 × 2 I 3 × 3 0 3 × 3 0 3 × 3 , ]]> I2×2、I3×3分别表示二阶和三阶单位矩阵,V为观测噪声阵,设为白噪声。(9)用卡尔曼滤波估计出主、子惯导系统之间的安装误差角φa和船体变形角θ。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于哈尔滨工程大学,未经哈尔滨工程大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/patent/201310156752.2/,转载请声明来源钻瓜专利网。
- 上一篇:一种碱式氟化铜的制备方法
- 下一篇:一种用于搪瓷的耐酸碱铸铁锑面釉