[发明专利]一种浅海匀速运动目标初始状态估计方法有效
申请号: | 202010134094.7 | 申请日: | 2020-03-02 |
公开(公告)号: | CN111273302B | 公开(公告)日: | 2022-08-23 |
发明(设计)人: | 李辉;徐哲臻;杨坤德 | 申请(专利权)人: | 西北工业大学 |
主分类号: | G01S15/88 | 分类号: | G01S15/88 |
代理公司: | 西安凯多思知识产权代理事务所(普通合伙) 61290 | 代理人: | 王鲜凯 |
地址: | 710072 *** | 国省代码: | 陕西;61 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 一种 浅海 匀速运动 目标 初始 状态 估计 方法 | ||
1.一种浅海匀速运动目标初始状态估计方法,其特征在于步骤如下:
步骤一:基于观测时间段内的n个观测值n为观测次数,k为观测时刻序号,k=1,2,3,…,n,为第k观测时刻目标水平方位角观测值,为第k观测时刻阵列不变量观测值,按下式构造测量矩阵和状态矩阵
其中,Xok=[xok,yok,vxok,vyok]T表示拖曳线列阵在第k时刻的状态,xok表示拖曳线列阵在第k时刻的横坐标,yok表示拖曳线列阵在第k时刻的纵坐标,vxok表示拖曳线列阵在第k时刻的速度沿x轴的分量,vyok表示拖曳线列阵在第k时刻速度沿y轴分量;表示第k时刻的测量向量1,表示第k时刻的测量向量2,Fk为第k时刻的匀速运动模型的状态转移矩阵,表示第k时刻的状态向量1,表示第k时刻的状态向量2,分别按照下列公式计算:
其中,c为声速,I2表示2阶单位矩阵,02表示2阶零矩阵,02×1表示大小为2×1的零矩阵,01×2表示大小为1×2的零矩阵,T为观测总时长;是拖曳线列阵在第k时刻的正横方向和目标方位之间的夹角的观测值;
步骤二:根据线性关系:
解出目标初始状态X0的估计值为:
其中,为一个5×1的向量,其前2个维度为目标初始时刻位置[x0,y0]T的估计值,第3、4个维度为初始时刻目标运动速度[vx0,vy0]T的估计值,最后一个维度为初始时刻当前环境波导不变量β0的估计值;x0表示初始时刻目标所在位置横坐标值,y0表示初始时刻目标所在位置纵坐标值,vx0表示初始时刻目标速度沿x轴的分量,vy0表示初始时刻目标速度沿y轴的分量;
步骤三:将目标初始状态的估计值带入匀速目标运动模型,得到第k观测时刻目标状态的估计值
其中,表示第k时刻目标所在位置横坐标的估计值,表示第k时刻目标所在位置纵坐标的估计值,表示第k时刻目标速度沿x轴分量的估计值,表示第k时刻目标速度沿y轴分量的估计值,表示第k时刻当前环境波导不变量的估计值;
分别按以下公式计算得到第k时刻目标水平方位角和阵列不变量的估计值:
其中,[xok,yok]T表示第k时刻拖曳线列阵的位置,xok表示拖曳线列阵在第k时刻的横坐标,yok表示拖曳线列阵在第k时刻的纵坐标,表示拖曳线列阵在第k时刻的正横方向和目标方位之间夹角的估计值,表示第k时刻目标与拖曳线列阵之间距离的估计值,按下式进行计算:
然后,按下式构造辅助变量矩阵G:
其中,表示第k时刻的测量向量1的估计值,表示第k时刻的测量向量2的估计值,分别按照下列公式计算:
步骤四:按下式计算得到最终的目标初始状态估计值
2.如权利要求1所述的一种浅海匀速运动目标初始状态估计方法,其特征在于:所述的拖曳线列阵在观测时间内处于运动状态且其运动复杂度高于目标运动的复杂度。
3.如权利要求1所述的一种浅海匀速运动目标初始状态估计方法,其特征在于:所述的目标水平方位角观测值利用拖曳线列阵的波束形成方法观测得到,所述的阵列不变量利用时间-到达角平面的曲线拟合的方法观测得到。
4.如权利要求1所述的一种浅海匀速运动目标初始状态估计方法,其特征在于:假设目标所处浅海环境的波导不变量在观测时间段内保持恒定。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于西北工业大学,未经西北工业大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/202010134094.7/1.html,转载请声明来源钻瓜专利网。