发明名称 一种基于滤波方法和曲线拟合方法相结合的协同导航方法
摘要 本发明属于应用曲线拟合方法的协同导航领域,具体为一种基于滤波方法和曲线拟合方法相结合的协同导航方法。本发明包括:建立从艇运动的系统方程和量测方程;滤波估计从艇轨迹;拟合数据保存;选取拟合向量;曲线拟合;利用拟合曲线预测下一个4秒从艇的位置。当从艇的轨迹近似为直线或机动性较弱时,由于MEMS陀螺误差较大且难以建立准确误差模型,这种方法相比滤波估计能提高定位精度。由于采用从艇当前位置之前的5个数据点进行拟合,充分利用了从艇的惯性,避免了使用之前全部数据点进行拟合多项式阶数高,难以准确估计从艇位置的缺点。
申请公布号 CN104034328B 申请公布日期 2017.03.29
申请号 CN201410216152.5 申请日期 2014.05.21
申请人 哈尔滨工程大学 发明人 徐博;金辰;刘杨;董海波;邱立民;贺浩;高伟;白金磊;单为
分类号 G01C21/00(2006.01)I 主分类号 G01C21/00(2006.01)I
代理机构 代理人
主权项 一种基于滤波方法和曲线拟合方法相结合的协同导航方法,其特征在于:(1)建立从艇运动的系统方程和量测方程<img file="FDA0000508589600000011.GIF" wi="648" he="257" />x<sub>k</sub>、y<sub>k</sub>表示从艇在k时刻的位置,v<sub>k</sub>为从艇速度,t表示推位时间间隔,<img file="FDA0000508589600000012.GIF" wi="63" he="71" />表示MEMS陀螺测得的航向,<img file="FDA0000508589600000013.GIF" wi="87" he="66" />表示MEMS陀螺误差,用估计值修正航向测量值<img file="FDA0000508589600000014.GIF" wi="56" he="69" />的偏差:X<sub>k+1</sub>=f(X<sub>k</sub>,u<sub>k</sub>,t)+w<sub>k</sub>线性化后得:X<sub>k+1</sub>=F<sub>k</sub>X<sub>k</sub>+B<sub>k</sub>u<sub>k</sub>+w<sub>k</sub>式中,<img file="FDA0000508589600000015.GIF" wi="1714" he="244" /><img file="FDA0000508589600000016.GIF" wi="1113" he="243" />系统噪声w<sub>k</sub>~N(0,Q<sub>k</sub>);根据主从艇的距离建立量测方程:<img file="FDA0000508589600000017.GIF" wi="1057" he="98" />观测量Z<sub>k</sub>表示主从艇的距离,x<sub>a</sub>、y<sub>a</sub>表示主艇的位置,x<sub>b</sub>、y<sub>b</sub>表示从艇位置,<img file="FDA0000508589600000018.GIF" wi="865" he="130" />V<sub>k</sub>~N(0,R<sub>k</sub>)为测量噪声;(2)滤波估计从艇轨迹时间更新:<img file="FDA0000508589600000019.GIF" wi="654" he="166" />量测更新:<img file="FDA0000508589600000021.GIF" wi="707" he="241" />式中,<img file="FDA0000508589600000022.GIF" wi="627" he="83" /><img file="FDA0000508589600000023.GIF" wi="62" he="82" />为滤波输出的状态估计,<img file="FDA0000508589600000024.GIF" wi="280" he="83" />为状态和及其方差一步预测,K<sub>k</sub>为滤波增益,滤波和曲线拟合可以同步并行运行,通过曲线拟合覆盖滤波估计值提高定位精度;(3)拟合数据保存在5秒时间内进行滤波,并保存每一个时刻的状态,将保存的滤波数据用多项式拟合方法拟合4个数据点,得到曲线的方程,并由曲线方程预测出5秒后从艇的位置,并将拟合数据保存;(4)选取拟合向量当每5秒进行量测更新时,从艇k时刻的位置坐标x<sub>k</sub>、y<sub>k</sub>都是时刻k的函数,把x<sub>k</sub>、y<sub>k</sub>分别分成p组,每组包含5个时刻的x<sub>k</sub>、y<sub>k</sub>,当组数p为偶数且k是5的整数倍时,拟合x<sub>k</sub>和y<sub>k</sub>关于时间k的多项式,设位置坐标x<sub>k</sub>、y<sub>k</sub>对应的拟合向量为u<sub>l</sub>、v<sub>l</sub>,令u<sub>l</sub>=x<sub>k‑5+l</sub>,v<sub>l</sub>=y<sub>k‑5+l</sub>,l=1,2,3,4;(5)曲线拟合当从艇每5秒测量和主艇的距离进行量测更新,把向量u<sub>l</sub>,v<sub>l</sub>关于l=1,2,3,4按照3阶多项式拟合方法拟合成多项式a(l),b(l),以下具体说明三阶多项式a(l)的拟合,设a(l)=a+bl+cl<sup>2</sup>+dl<sup>3</sup>,l=1,2,3,4,用最小二乘原理求三次拟合曲线的系数,拟合方差为:<img file="FDA0000508589600000025.GIF" wi="416" he="148" />求a(l)的系数使拟合方差最小,令<img file="FDA0000508589600000026.GIF" wi="594" he="138" />代入u<sub>l</sub>(l=1,2,3,4)四个点,得到四个方程分别解出a(l)的系数a,b,c,d;(6)利用拟合曲线预测下一个4秒从艇的位置每5秒进行量测更新,当时间k%5≠0时,令x<sub>k</sub>=a(5+k%5),y<sub>k</sub>=b(5+k%5),利用拟合多项式a(l),b(l)预测5秒后的4秒内从艇的位置,并将预测值赋给x<sub>k</sub>,y<sub>k</sub>,即覆盖滤波估计值,再进行新的量测更新。
地址 150001 黑龙江省哈尔滨市南岗区南通大街145号哈尔滨工程大学科技处知识产权办公室