CN107390199A - 一种雷达机动目标跟踪波形设计方法 - Google Patents
一种雷达机动目标跟踪波形设计方法 Download PDFInfo
- Publication number
- CN107390199A CN107390199A CN201710855364.1A CN201710855364A CN107390199A CN 107390199 A CN107390199 A CN 107390199A CN 201710855364 A CN201710855364 A CN 201710855364A CN 107390199 A CN107390199 A CN 107390199A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mover
- mtd
- msup
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000013461 design Methods 0.000 title claims abstract description 17
- 239000011159 matrix material Substances 0.000 claims abstract description 98
- 230000004927 fusion Effects 0.000 claims abstract description 29
- 230000007704 transition Effects 0.000 claims abstract description 22
- 238000005259 measurement Methods 0.000 claims abstract description 19
- 238000001914 filtration Methods 0.000 claims abstract description 9
- 239000013598 vector Substances 0.000 claims description 20
- 239000002131 composite material Substances 0.000 claims description 13
- 230000008569 process Effects 0.000 claims description 9
- 230000005540 biological transmission Effects 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 6
- 239000007787 solid Substances 0.000 claims description 4
- 238000012546 transfer Methods 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 abstract description 4
- 238000004891 communication Methods 0.000 abstract description 3
- 150000001875 compounds Chemical class 0.000 abstract 1
- 238000000342 Monte Carlo simulation Methods 0.000 description 4
- 230000001149 cognitive effect Effects 0.000 description 4
- 230000003993 interaction Effects 0.000 description 4
- 230000001133 acceleration Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000007711 solidification Methods 0.000 description 2
- 230000008023 solidification Effects 0.000 description 2
- 241000375392 Tana Species 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000008685 targeting Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
Abstract
一种雷达机动目标跟踪波形设计方法,属于雷达通信技术领域,具体涉及雷达机动目标跟踪波形设计方法。本发明首先对机动目标构造运动模型,计算滤波器更新矩阵权值、运动模型的混合输入状态及对应的估计误差协方差矩阵,然后利用线性或非线性滤波算法获得各子模型的局部无偏滤波器估值和目标状态估计误差协方差矩阵,进行复合矩阵更新以获得最优融合状态估计及融合估计误差协方差矩阵,在此基础上获得发射波形旋转参数,利用分数阶傅里叶变换旋转用户设定波形得到新的量测误差椭圆及发射波形,最后进行马尔科夫转移概率矩阵的更新,以达到更好的跟踪精度。本发明解决了机动目标跟踪稳健性不强、准确性低的问题。本发明可运用于雷达通信技术。
Description
技术领域
本发明属于雷达通信技术领域,具体涉及雷达机动目标跟踪波形设计方法。
背景技术
战场中机动目标状态往往呈现随机性、多样性特点,使得传统雷达难以有效追踪,成为当前研究的难点。多数研究从接收端数据处理出发,侧重目标状态建模及滤波算法改进(见文献:New interacting multiple model algorithms for the tracking of themaneuvering target,FU X,JIA Y,DU J,et al.;IET control theory&applications,2010,4(10):2184-2194;Dynamic waveform selection for maneuvering targettracking in clutter,WANG Jiantao,QIN Yuliang,WANG Hongqiang,et al;IET Radar,Sonar&Navigation,2013,7(7):815-825);忽视目标跟踪精度不仅与数据处理方式有关,且受发射波形影响,导致跟踪误差大、鲁棒性差等问题。
新近,波形设计与多传感器接收数据联合处理成为提升机动目标跟踪性能的热点(见文献:Cognitive waveform and receiver selection mechanism for multistaticradar,KILANI M B,NIJSURE Y,GAGNON G,et al;IET Radar,Sonar&Navigation,2016,10(2):417-425;基于波形捷变的多传感器机动目标跟踪,盛丹,王国宏,张翔宇;系统工程与电子技术,2015,(03):485-491)。但多传感器配置耗资过大,数据实时性也难以兼顾;另外,目标模型失配或先验知识缺失下目标跟踪性能也难以保证。为解决目标模型失配问题,常见思路为构造交互多模型(IMM)框架(见文献:跟踪机动目标的雷达波形选择新方法,檀甲甲,张建秋;系统工程与电子技术,2011,(03):515-522+54;Waveform selection formaneuvering targets within an IMM framework,SAVAGE C O,MORAN B;IEEETransactions on Aerospace and Electronic Systems,2007,43(3);针对机动目标跟踪的雷达发射波形选择,靳标,纠博,苏涛等;电子与信息学报,2014,(08):1912-1918);但IMM中采用概率加权来联接不同子模型的状态输出,忽视目标状态向量及其误差协方差矩阵中不同维度元素间的相关性,进而影响目标状态匹配性能,使得跟踪误差随时间积累而逐步增大。另外,IMM中固定马尔科夫转移矩阵亦使模型转换速度变慢、精度降低(见文献:马尔可夫矩阵修正IMM跟踪算法,封普文,黄长强,曹林平等;系统工程与电子技术,2013,11:2269-2274.时变转移概率IMM-SRCKF机动目标跟踪算法,郭志,董春云,蔡远利等;系统工程与电子技术,2015,01:24-30)。
以上提出的IMM、MIMM、AIMM等技术跟踪精度差、稳定效果欠佳,不适合工程应用。针对机动目标跟踪模型失配、多模型融合中概率权重难以兼顾速度与位置估计等棘手问题,急需一种有望获得更低的跟踪误差和更强的稳健性的跟踪波形设计方法。
发明内容
本发明为解决现有技术忽视各模型位置误差与速度误差相关性,导致跟踪稳健性不强、准确性低的问题,提供了一种雷达机动目标跟踪波形设计方法。
本发明所述雷达机动目标跟踪波形设计方法,通过以下技术方案实现:
步骤一、对机动目标构造运动模型;
步骤二、计算与运动模型相对应的滤波器更新矩阵权值;
步骤三、计算运动模型的混合输入状态及对应的估计误差协方差矩阵;
步骤四、利用线性或非线性滤波算法获得各子模型的局部无偏滤波器估值和目标状态估计误差协方差矩阵;
步骤五、进行复合矩阵更新以获得最优融合状态估计及融合估计误差协方差矩阵;
步骤六、基于融合估计误差协方差矩阵获得发射波形旋转参数,利用分数阶傅里叶变换旋转用户设定波形得到新的量测误差椭圆及发射波形;
步骤七、进行马尔科夫转移概率矩阵的更新。
本发明与现有技术相比较,最为突出的特点和显著的有益效果是:
本发明将多模型融合与雷达波形设计相结合应用于机动目标跟踪,采用多模型复合矩阵加权的认知跟踪雷达波形,相比采用现有的跟踪波形在线设计方法,本发明具有更好的自由性、使用性,且效率高、可有效避免机动目标失配、跟踪丢失问题,跟踪稳健性较高。应用本发明的机动目标多模型复合矩阵加权跟踪波形设计方法,可使雷达机动目标跟踪准确性提高。
本发明所提设计方法可兼顾目标状态向量及误差协方差矩阵中元素间关系,均衡各模型位置误差与速度误差,克服采用近似概率密度函数作为融合权值时出现的概率函数失效问题,所设计的波形增强跟踪稳健性、准确性,更有利于目标的检测与跟踪,其优点如下:
(1)就其采用多模型融合中复合矩阵加权产生方式而言,本发明提出的认知雷达跟踪波形设计思路具有精确的表达式及具体的步骤程序,避免了现有技术中涉及的传统概率加权或常值加权难以兼顾位置误差与速度误差的问题。
(2)就多模型中转移概率矩阵设计问题,本发明提出的指数形式模型概率变化率构造方法,考虑相邻时刻模型概率之差反映模型与实际运动模式间的匹配程度,具有较高的灵敏度和转换速度,相比现有技术解决了传统IMM中马尔科夫转移概率矩阵固化问题。
(3)就跟踪波形参数优化而言,本发明提出的基于融合后估计误差协方差矩阵知识的设计方法,通过旋转特定波形(用户设定波形)模糊函数使其量测误差椭圆与融合估计误差椭圆正交,来获得最优波形,避免信息论方法中线性调频参数优选只保留最大或最小调频率的贫化问题,增加了波形参数设计的自由度。其他性能指标(跟踪误差RMSE)完全优于已有技术IMM、MIMM、AIMM等的相应指标。
以1000次独立蒙特卡洛实验来评估各算法统计性能,量测误差的强度为250,并以目标状态估计的均方根误差RMSE作为评价指标,IMM算法、AIMM算法、MIMM算法以及本发明的平均跟踪误差分别为:0.58、0.40、0.37、0.31,平均速度跟踪误差为:0.15、0.13、0.12、0.09;本发明明显优于其他算法。
附图说明
图1是本发明实施例中,目标位置轨迹图;
图2是本发明实施例中,目标的速度随时间变化的线性图;
图3是本发明实施例中,噪声强度r=250下的X方向位置跟踪误差图;
图4是本发明实施例中,噪声强度r=250下的Y方向位置跟踪误差图;
图5是本发明实施例中,噪声强度r=250下的X方向速度跟踪误差图;
图6是本发明实施例中,噪声强度r=250下的Y方向速度跟踪误差图。
具体实施方式
具体实施方式一:本实施方式给出的一种雷达机动目标跟踪波形设计方法,具体是按照以下步骤进行的:
步骤一、对机动目标构造运动模型;
步骤二、计算与运动模型相对应的滤波器更新矩阵权值;
步骤三、计算运动模型的混合输入状态及对应的估计误差协方差矩阵;
步骤四、利用线性或非线性滤波算法获得各子模型的局部无偏滤波器估值和目标状态估计误差协方差矩阵;
步骤五、进行复合矩阵更新以获得最优融合状态估计及融合估计误差协方差矩阵;
步骤六、基于融合估计误差协方差矩阵获得发射波形旋转参数,利用分数阶傅里叶变换旋转用户设定波形得到新的量测误差椭圆及发射波形,用户设定波形指的是可由用户自行决定的波形,体现了本发明的普适性;
步骤七、进行马尔科夫转移概率矩阵的更新,利用量测数据及模型概率变化率自适应地修正转移概率矩阵,以达到更好的跟踪精度。
具体实施方式二:本实施方式与具体实施方式一不同的是:所述步骤一中对机动目标构造运动模型具体为:
式中包括目标状态方程x(k+1)=Fj(x(k))+wj(k)和量测方程z(k)=Hj(x(k))+v(k);x(k)表示k时刻的目标状态向量,维度为n×1,含X、Y方向位置和速度状态,z(k)为量测向量;j∈{1,...,s}表示模型库中的模型序号,s为模型数目;当上式表示线性运动模型时,Fj(·)和Hj(·)为线性转移矩阵,非线性运动模型时Fj(·)和Hj(·)表示非线性函数;wj(k)表示均值为零、协方差矩阵为Qj的高斯过程噪声,v(k)表示均值为零、协方差矩阵为R的量测噪声;表示k时刻机动目标的第j个运动模型,则所述运动模型从k-1时刻跳变到k时刻的转移过程,能够用具有转移概率矩阵的一阶马尔科夫过程表示,其中转移概率矩阵中第i行第j列的元素为:
式中:Prob(·)表示概率密度函数,表示条件下的转移概率。
具体实施方式三:本实施方式与具体实施方式二不同的是:所述步骤二计算与运动模型相对应的滤波器更新矩阵权值具体包括以下步骤:
对状态向量x而言,设s个子模型相应的局部无偏滤波器估值为表示局部估值误差,表示第i个运动模型的目标状态估计误差协方差矩阵,表示运动模型和的估计误差互协方差矩阵,且i≠j时有其中E(·)为求取期望值函数;局部无偏滤波器估值可视为第i个运动模型相对应的滤波器对x的量测,即有:
进而可定义:
其中:
e=[In ... In]T
由无偏性知即In表示n×n维单位矩阵,e=[In ... In]T为列满秩矩阵;的协方差矩阵可构造为:
因为目标状态向量为n×1维度,则局部无偏滤波器估值也为n×1维,则Pii与Pij为n×n维,所以P为ns×ns维;
由线性无偏最小方差估计定理,n×ns维无偏性约束加权矩阵为Ω,即Ωe=In,用分块矩阵表示为:
Ω=[Ω1 ... Ωs]
各子模型滤波误差互不相关,即Pij=0(i≠j),则在线性无偏最小方差意义下最优加权矩阵为:
另构造矩阵其元素能够通过以下计算求得:
通过a=1,…,n,可以求得
而中的元素属于
因此可以求得矩阵中的元素;
其中,为中间过程参数,为Pii的第a个对角元素;
进而获得复合矩阵更新权值:
式中diag(·)表示取矩阵对角元素构成行向量,Diag(·)表示将行向量构成对角阵;
与运动模型相对应的滤波器更新矩阵如下:
其中,Zk表示1至k时刻的量测向量序列。
具体实施方式四:本实施方式与具体实施方式三不同的是:所述步骤三中运动模型的混合输入状态及对应的估计误差协方差矩阵Poj计算方法为:
其它步骤及参数与具体实施方式三相同。
具体实施方式五:本实施方式与具体实施方式四不同的是:所述步骤五进行复合矩阵更新以获得最优融合状态估计及融合估计误差协方差矩阵具体为:
获得复合矩阵更新权值:
最优融合状态估值为:
且融合估计误差协方差矩阵为:
具体实施方式六:
本实施方式与具体实施方式五不同的是:所述步骤六中基于融合估计误差协方差矩阵获得发射波形旋转参数,利用分数阶傅里叶变换旋转用户设定波形得到新的量测误差椭圆及发射波形包括以下具体步骤:
对融合估计误差协方差矩阵PM(k)进行特征值分解,令λi,vi分别代表特征值及其特征向量,利用最大特征值对应的特征向量vmax中元素间的关系进行误差椭圆角度计算:
β=arctan(vmax(2)/vmax(1))
式中:β∈[-π/4,π/4],vmax(1)指的是特征向量vmax的第一个元素,vmax(2)指的是特征向量vmax的第二个元素;令R0表示用户设定波形估计位置与速度的克拉美罗界(CRLB),同理可得R0的误差椭圆角度参数ψ;令θk+1为发射波形旋转参数,则量测误差椭圆旋转关系可表示为:
通过旋转用户设定波形模糊函数以使其量测误差椭圆R(θk+1)与估计误差椭圆正交,得到旋转参数为:
θk+1=-π/2-ψ+β
进而利用分数阶傅里叶变换获得相应波形。
具体实施方式七:
本实施方式与具体实施方式六不同的是:所述步骤七中进行马尔科夫转移概率矩阵的更新具体为:
为解决传统IMM中马尔科夫转移概率矩阵固化问题,考虑相邻时刻模型概率之差反映模型与实际运动模式间的匹配程度,构造指数形式模型概率变化率为:
式中:ρ表示转换因子;实时修正模型转移概率,即:
当运动模型概率随时间增大时,有κj(k)>1,矩阵的第j列元素也随之增大,滤波前模型交互中概率大的子模型状态估计输出在交互过程中所占比重更大;反之,当运动模型概率随时间减小时,有κj(k)<1,矩阵的第j列元素也随之减小,滤波前模型交互中概率小的子模型状态估计输出在交互过程中所占比重减小。利用量测数据及模型概率变化率自适应地修正转移概率矩阵,放大匹配模型作用,抑制非匹配模型作用,从而提高误差收敛速度。
具体实施方式八:
本实施方式与具体实施方式七不同的是:所述步骤七中转换因子的范围是:ρ∈[1,2]。
实施例
如图3~图6所示,结合步骤一~步骤七给出的基于本发明的技术方案,从而可设计针对单个或多个目标跟踪的多模型复合矩阵加权的认知跟踪波形的实施例:
目标初始状态为:x(0)=[2100 0 10000 -15];先匀速运动200s,然后在X方向进行机动转弯,加速度为ax=ay=0.1m/s2,并在t=270s时完成转弯,加速度为0;第二个状态机动在t=410s处以ax=0.05m/s2,ay=0.03m/s2至t=610s;再以ax=-0.3m/s2,ay=0.2m/s2持续至t=640s。其X-Y平面的位置、速度轨迹如图1、图2所示。采用两个具有不同过程噪声强度的恒速度模型来构造多模型框架,如下:
式中:Fi表示第i个模型的状态转移矩阵,Hi表示量测矩阵;Qi表示过程噪声协方差矩阵,取q1=0.5,q2=50;R表示量测误差协方差矩阵,r表示量测误差的强度;采样间隔T=10s;另定义目标初始状态误差协方差矩阵为:
模型转移概率矩阵为:
模型初始概率为以1000次独立蒙特卡洛实验来评估各算法统计性能,并以目标状态估计的均方根误差RMSE作为评价指标,即:
其中,xn(k)和分别表示第k时刻第n次蒙特卡洛实验的真实目标状态及其估计值;N表示蒙特卡罗实验次数。位置跟踪误差比较(r=250)如图3、图4所示,速度跟踪误差比较(r=250)如图5、图6所示;各算法平均跟踪误差及耗时比较(r=250)见表1。为便于比较,特将本专利所提方法命名为EAMIMM。
表1各算法平均跟踪误差及耗时比较(r=250)
本发明还可有其它多种实施例,在不背离本发明精神及其实质的情况下,本领域技术人员当可根据本发明作出各种相应的改变和变形,但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。
Claims (8)
1.一种雷达机动目标跟踪波形设计方法,其特征在于,所述方法包括以下步骤:
步骤一、对机动目标构造运动模型;
步骤二、计算与运动模型相对应的滤波器更新矩阵权值;
步骤三、计算运动模型的混合输入状态及对应的估计误差协方差矩阵;
步骤四、利用线性或非线性滤波算法获得各子模型的局部无偏滤波器估值和目标状态估计误差协方差矩阵;
步骤五、进行复合矩阵更新以获得最优融合状态估计及融合估计误差协方差矩阵;
步骤六、基于融合估计误差协方差矩阵获得发射波形旋转参数,利用分数阶傅里叶变换旋转用户设定波形得到新的量测误差椭圆及发射波形;
步骤七、进行马尔科夫转移概率矩阵的更新。
2.根据权利要求1所述的一种雷达机动目标跟踪波形设计方法,其特征在于,所述步骤一中对机动目标构造运动模型具体为:
<mrow>
<msubsup>
<mi>M</mi>
<mi>j</mi>
<mi>k</mi>
</msubsup>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>F</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>w</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>z</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>H</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
式中:x(k)表示k时刻的目标状态向量,维度为n×1,z(k)为量测向量;j∈{1,...,s}表示模型库中的模型序号,s为模型数目;Fj(·)、Hj(·)分别为线性转移矩阵,或者Fj(·)、Hj(·)分别为非线性函数;wj(k)表示均值为零、协方差矩阵为Qj的高斯过程噪声,v(k)表示均值为零、协方差矩阵为R的量测噪声;表示k时刻机动目标的第j个运动模型,则所述运动模型从k-1时刻跳变到k时刻的转移过程,能够用具有转移概率矩阵的一阶马尔科夫过程表示,其中转移概率矩阵中第i行第j列的元素为:
<mrow>
<msub>
<mover>
<mi>&pi;</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>Pr</mi>
<mi>o</mi>
<mi>b</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>M</mi>
<mi>j</mi>
<mi>k</mi>
</msubsup>
<mo>|</mo>
<msubsup>
<mi>M</mi>
<mi>i</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>0</mn>
<mo>&le;</mo>
<msub>
<mover>
<mi>&pi;</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>&le;</mo>
<mn>1</mn>
<mo>,</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>s</mi>
</munderover>
<msub>
<mover>
<mi>&pi;</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mn>1</mn>
<mo>;</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>&Element;</mo>
<mo>{</mo>
<mn>1</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>s</mi>
<mo>}</mo>
</mrow>
式中:Prob(·)表示概率密度函数,表示条件下的转移概率。
3.根据权利要求2所述的一种雷达机动目标跟踪波形设计方法,其特征在于,所述步骤二计算与运动模型相对应的滤波器更新矩阵权值具体包括以下步骤:
对状态向量x而言,设s个子模型相应的局部无偏滤波器估值为表示局部估值误差,表示第i个运动模型的目标状态估计误差协方差矩阵,表示运动模型和的估计误差互协方差矩阵,且i≠j时有其中E(·)为求取期望值函数;局部无偏滤波器估值能够视为第i个运动模型相对应的滤波器对x的量测,即有:
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mi>i</mi>
</msub>
<mo>=</mo>
<mi>x</mi>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>s</mi>
</mrow>
设:
<mrow>
<mi>y</mi>
<mo>=</mo>
<mi>e</mi>
<mi>x</mi>
<mo>+</mo>
<mover>
<mi>v</mi>
<mo>^</mo>
</mover>
<mo>,</mo>
</mrow>
其中:
<mrow>
<mi>y</mi>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mi>i</mi>
</msub>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mi>s</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
</mrow>
e=[In ... In]T
<mrow>
<mover>
<mi>v</mi>
<mo>^</mo>
</mover>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mi>i</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mi>s</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
</mrow>
由无偏性知即In表示n×n维单位矩阵,e=[In ... In]T为列满秩矩阵;的协方差矩阵能够构造为:
由线性无偏最小方差估计定理,n×ns维无偏性约束加权矩阵为Ω,即Ωe=In,用分块矩阵表示为:
Ω=[Ω1 ... Ωs]
各子模型滤波误差互不相关,即Pij=0且i≠j,则在线性无偏最小方差意义下最优加权矩阵为:
<mrow>
<msub>
<mi>&Omega;</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>s</mi>
</munderover>
<msubsup>
<mi>P</mi>
<mrow>
<mi>j</mi>
<mi>j</mi>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mi>i</mi>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>,</mo>
<mi>i</mi>
<mo>&Element;</mo>
<mo>{</mo>
<mn>1</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>s</mi>
<mo>}</mo>
</mrow>
另构造矩阵其元素能够通过以下计算求得:
<mrow>
<msub>
<mover>
<mi>w</mi>
<mo>&OverBar;</mo>
</mover>
<mi>a</mi>
</msub>
<mo>=</mo>
<mo>&lsqb;</mo>
<mtable>
<mtr>
<mtd>
<msub>
<mi>w</mi>
<mrow>
<mn>1</mn>
<mi>a</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>w</mi>
<mrow>
<mn>2</mn>
<mi>a</mi>
</mrow>
</msub>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msub>
<mi>w</mi>
<mrow>
<mi>s</mi>
<mi>a</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<msub>
<mover>
<mi>w</mi>
<mo>&OverBar;</mo>
</mover>
<mi>a</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<mi>e</mi>
<mi>T</mi>
</msup>
<msup>
<mrow>
<mo>(</mo>
<msup>
<mi>P</mi>
<mrow>
<mi>a</mi>
<mi>a</mi>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
</mrow>
<mrow>
<msup>
<mi>e</mi>
<mi>T</mi>
</msup>
<msup>
<mrow>
<mo>(</mo>
<msup>
<mi>P</mi>
<mrow>
<mi>a</mi>
<mi>a</mi>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>e</mi>
</mrow>
</mfrac>
<mo>,</mo>
<mi>a</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>n</mi>
</mrow>
其中,为Pii的第a个对角元素;
进而获得复合矩阵更新权值:
<mrow>
<msub>
<mover>
<mi>&Omega;</mi>
<mo>~</mo>
</mover>
<mi>i</mi>
</msub>
<mo>=</mo>
<msub>
<mi>&Omega;</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mi>D</mi>
<mi>i</mi>
<mi>a</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>d</mi>
<mi>i</mi>
<mi>a</mi>
<mi>g</mi>
<mo>(</mo>
<msub>
<mi>&Omega;</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mover>
<mi>&Omega;</mi>
<mo>^</mo>
</mover>
<mi>i</mi>
</msub>
</mrow>
式中diag(·)表示取矩阵对角元素构成行向量,Diag(·)表示将行向量构成对角阵;
与运动模型相对应的滤波器更新矩阵如下:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>&Omega;</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>i</mi>
<mo>|</mo>
<mi>j</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>P</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>M</mi>
<mi>i</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>|</mo>
<msubsup>
<mi>M</mi>
<mi>j</mi>
<mi>k</mi>
</msubsup>
<mo>,</mo>
<msup>
<mi>Z</mi>
<mi>k</mi>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mover>
<mi>&pi;</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mover>
<mi>&Omega;</mi>
<mo>~</mo>
</mover>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>s</mi>
</munderover>
<msub>
<mover>
<mi>&pi;</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mover>
<mi>&Omega;</mi>
<mo>~</mo>
</mover>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>s</mi>
</munderover>
<msub>
<mover>
<mi>&pi;</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mover>
<mi>&Omega;</mi>
<mo>~</mo>
</mover>
<mi>i</mi>
</msub>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msub>
<mover>
<mi>&pi;</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mover>
<mi>&Omega;</mi>
<mo>~</mo>
</mover>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
2
其中,Zk表示1至k时刻的量测向量序列。
4.根据权利要求3所述的一种雷达机动目标跟踪波形设计方法,其特征在于,所述步骤三中运动模型的混合输入状态及对应的估计误差协方差矩阵Poj计算方法为:
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>o</mi>
<mi>j</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>s</mi>
</munderover>
<msub>
<mover>
<mi>&Omega;</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>i</mi>
<mo>|</mo>
<mi>j</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>o</mi>
<mi>j</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>s</mi>
</munderover>
<msub>
<mover>
<mi>&Omega;</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>i</mi>
<mo>|</mo>
<mi>j</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>{</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mo>&lsqb;</mo>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>o</mi>
<mi>j</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>&times;</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>o</mi>
<mi>j</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>}</mo>
<mo>.</mo>
</mrow>
5.根据权利要求4所述的一种雷达机动目标跟踪波形设计方法,其特征在于,所述步骤五进行复合矩阵更新以获得最优融合状态估计及融合估计误差协方差矩阵具体为:
<mrow>
<msub>
<mi>&Omega;</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>s</mi>
</munderover>
<msubsup>
<mi>P</mi>
<mrow>
<mi>j</mi>
<mi>j</mi>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mi>i</mi>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>,</mo>
<mi>i</mi>
<mo>&Element;</mo>
<mo>{</mo>
<mn>1</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>s</mi>
<mo>}</mo>
</mrow>
获得复合矩阵更新权值:
<mrow>
<msub>
<mover>
<mi>&Omega;</mi>
<mo>~</mo>
</mover>
<mi>i</mi>
</msub>
<mo>=</mo>
<msub>
<mi>&Omega;</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mi>D</mi>
<mi>i</mi>
<mi>a</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>d</mi>
<mi>i</mi>
<mi>a</mi>
<mi>g</mi>
<mo>(</mo>
<msub>
<mi>&Omega;</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mover>
<mi>&Omega;</mi>
<mo>^</mo>
</mover>
<mi>i</mi>
</msub>
</mrow>
最优融合状态估值为:
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mi>M</mi>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>s</mi>
</munderover>
<msub>
<mover>
<mi>&Omega;</mi>
<mo>~</mo>
</mover>
<mi>i</mi>
</msub>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mi>i</mi>
</msub>
</mrow>
且融合估计误差协方差矩阵为:
<mrow>
<msub>
<mi>P</mi>
<mi>M</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>s</mi>
</munderover>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mi>i</mi>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>.</mo>
</mrow>
6.根据权利要求5所述的一种雷达机动目标跟踪波形设计方法,其特征在于,所述步骤六中基于融合估计误差协方差矩阵获得发射波形旋转参数,利用分数阶傅里叶变换旋转用户设定波形得到新的量测误差椭圆及发射波形包括以下具体步骤:
对融合估计误差协方差矩阵PM(k)进行特征值分解,令λi,vi分别代表特征值及其特征向量,利用最大特征值对应的特征向量vmax中元素间的关系进行误差椭圆角度计算:
β=arctan(vmax(2)/vmax(1))
式中:β∈[-π/4,π/4],vmax(1)指的是特征向量vmax的第一个元素,vmax(2)指的是特征向量vmax的第二个元素;令R0表示用户设定波形估计位置与速度的克拉美罗界,并得到R0的误差椭圆角度参数ψ;令θk+1为发射波形旋转参数,则量测误差椭圆旋转关系表示为:
<mrow>
<mi>R</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>cos&theta;</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>sin&theta;</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>sin&theta;</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>cos&theta;</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msub>
<mi>R</mi>
<mn>0</mn>
</msub>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>cos&theta;</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>sin&theta;</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>sin&theta;</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>cos&theta;</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<mo>-</mo>
<mi>T</mi>
</mrow>
</msup>
</mrow>
通过旋转用户设定波形模糊函数以使其量测误差椭圆R(θk+1)与估计误差椭圆正交,得到旋转参数为:
θk+1=-π/2-ψ+β
进而利用分数阶傅里叶变换获得相应波形。
7.根据权利要求6所述的一种雷达机动目标跟踪波形设计方法,其特征在于,所述步骤七中进行马尔科夫转移概率矩阵的更新具体为:
构造指数形式模型概率变化率为:
<mrow>
<msub>
<mi>&kappa;</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mi>&rho;</mi>
<mo>&CenterDot;</mo>
<mo>(</mo>
<mrow>
<msubsup>
<mi>u</mi>
<mi>j</mi>
<mi>k</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>u</mi>
<mi>j</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
式中:ρ表示转换因子;实时修正模型转移概率,即:
<mrow>
<msub>
<mover>
<mi>&pi;</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>&kappa;</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>&pi;</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>,</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>s</mi>
</mrow>
利用量测数据及模型概率变化率自适应地修正转移概率矩阵。
8.根据权利要求7所述的一种雷达机动目标跟踪波形设计方法,其特征在于,所述步骤七中转换因子的范围是:ρ∈[1,2]。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710855364.1A CN107390199B (zh) | 2017-09-20 | 2017-09-20 | 一种雷达机动目标跟踪波形设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710855364.1A CN107390199B (zh) | 2017-09-20 | 2017-09-20 | 一种雷达机动目标跟踪波形设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107390199A true CN107390199A (zh) | 2017-11-24 |
CN107390199B CN107390199B (zh) | 2019-06-18 |
Family
ID=60349958
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710855364.1A Active CN107390199B (zh) | 2017-09-20 | 2017-09-20 | 一种雷达机动目标跟踪波形设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107390199B (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109975765A (zh) * | 2019-03-28 | 2019-07-05 | 中国人民解放军国防科技大学 | 一种雷达嵌入式通信的正交波形设计方法 |
CN110471034A (zh) * | 2019-09-19 | 2019-11-19 | 上海无线电设备研究所 | 一种超宽带雷达波形设计方法 |
CN110706265A (zh) * | 2019-11-05 | 2020-01-17 | 中国人民解放军国防科技大学 | 一种改进srckf强跟踪滤波的机动目标跟踪方法 |
CN111175738A (zh) * | 2020-01-08 | 2020-05-19 | 中国船舶重工集团公司第七二四研究所 | 一种基于多模型隶属控制的相控阵雷达目标快速建航方法 |
CN112800889A (zh) * | 2021-01-18 | 2021-05-14 | 浙江工业大学 | 一种基于分布式矩阵加权融合高斯滤波的目标跟踪方法 |
CN113190960A (zh) * | 2021-03-26 | 2021-07-30 | 大连海事大学 | 一种基于非等维状态混合估计的并行imm机动目标跟踪方法 |
CN113238219A (zh) * | 2021-05-18 | 2021-08-10 | 电子科技大学 | 一种适用于rsn在目标跟踪任务中的波形优化方法 |
CN113468750A (zh) * | 2021-07-06 | 2021-10-01 | 东南大学 | 一种基于信息融合的天线罩误差斜率估计方法 |
CN113566821A (zh) * | 2021-06-28 | 2021-10-29 | 江南造船(集团)有限责任公司 | 基于交互式滤波的无人机航向估计方法、系统及电子设备 |
CN117169818A (zh) * | 2023-10-30 | 2023-12-05 | 哈尔滨工业大学(威海) | 一种面向海面机动目标跟踪的雷达波形设计方法 |
CN117648521A (zh) * | 2023-11-30 | 2024-03-05 | 常熟理工学院 | 基于角度量测序贯结构状态估计方法、系统及存储介质 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102721951A (zh) * | 2012-05-04 | 2012-10-10 | 西安电子科技大学 | 一种高机动目标跟踪方法 |
CN103853908A (zh) * | 2012-12-04 | 2014-06-11 | 中国科学院沈阳自动化研究所 | 一种自适应交互式多模型的机动目标跟踪方法 |
CN104020466A (zh) * | 2014-06-17 | 2014-09-03 | 西安电子科技大学 | 基于变结构多模型的机动目标跟踪方法 |
CN104880707A (zh) * | 2014-11-30 | 2015-09-02 | 中国科学院沈阳自动化研究所 | 一种基于自适应转移概率矩阵的交互多模型跟踪方法 |
CN105652250A (zh) * | 2016-01-15 | 2016-06-08 | 西北工业大学 | 一种基于双层期望最大化的机动目标跟踪技术 |
CN105954743A (zh) * | 2016-05-31 | 2016-09-21 | 西安电子科技大学 | 一种改进权值的变结构多模型机动目标跟踪方法 |
-
2017
- 2017-09-20 CN CN201710855364.1A patent/CN107390199B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102721951A (zh) * | 2012-05-04 | 2012-10-10 | 西安电子科技大学 | 一种高机动目标跟踪方法 |
CN103853908A (zh) * | 2012-12-04 | 2014-06-11 | 中国科学院沈阳自动化研究所 | 一种自适应交互式多模型的机动目标跟踪方法 |
CN104020466A (zh) * | 2014-06-17 | 2014-09-03 | 西安电子科技大学 | 基于变结构多模型的机动目标跟踪方法 |
CN104880707A (zh) * | 2014-11-30 | 2015-09-02 | 中国科学院沈阳自动化研究所 | 一种基于自适应转移概率矩阵的交互多模型跟踪方法 |
CN105652250A (zh) * | 2016-01-15 | 2016-06-08 | 西北工业大学 | 一种基于双层期望最大化的机动目标跟踪技术 |
CN105954743A (zh) * | 2016-05-31 | 2016-09-21 | 西安电子科技大学 | 一种改进权值的变结构多模型机动目标跟踪方法 |
Non-Patent Citations (3)
Title |
---|
檀甲甲 等: "跟踪机动目标的雷达波形选择新方法", 《系统工程与电子技术》 * |
郭志 等: "时变转移概率IMM-SRCKF机动目标跟踪算法", 《系统工程与电子技术》 * |
靳标 等: "针对机动目标跟踪的雷达发射波形选择", 《电子与信息学报》 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109975765A (zh) * | 2019-03-28 | 2019-07-05 | 中国人民解放军国防科技大学 | 一种雷达嵌入式通信的正交波形设计方法 |
CN110471034A (zh) * | 2019-09-19 | 2019-11-19 | 上海无线电设备研究所 | 一种超宽带雷达波形设计方法 |
CN110471034B (zh) * | 2019-09-19 | 2021-08-31 | 上海无线电设备研究所 | 一种超宽带雷达波形设计方法 |
CN110706265A (zh) * | 2019-11-05 | 2020-01-17 | 中国人民解放军国防科技大学 | 一种改进srckf强跟踪滤波的机动目标跟踪方法 |
CN110706265B (zh) * | 2019-11-05 | 2022-02-01 | 中国人民解放军国防科技大学 | 一种改进srckf强跟踪滤波的机动目标跟踪方法 |
CN111175738A (zh) * | 2020-01-08 | 2020-05-19 | 中国船舶重工集团公司第七二四研究所 | 一种基于多模型隶属控制的相控阵雷达目标快速建航方法 |
CN112800889B (zh) * | 2021-01-18 | 2024-02-23 | 浙江工业大学 | 一种基于分布式矩阵加权融合高斯滤波的目标跟踪方法 |
CN112800889A (zh) * | 2021-01-18 | 2021-05-14 | 浙江工业大学 | 一种基于分布式矩阵加权融合高斯滤波的目标跟踪方法 |
CN113190960A (zh) * | 2021-03-26 | 2021-07-30 | 大连海事大学 | 一种基于非等维状态混合估计的并行imm机动目标跟踪方法 |
CN113190960B (zh) * | 2021-03-26 | 2024-03-12 | 大连海事大学 | 一种基于非等维状态混合估计的并行imm机动目标跟踪方法 |
CN113238219A (zh) * | 2021-05-18 | 2021-08-10 | 电子科技大学 | 一种适用于rsn在目标跟踪任务中的波形优化方法 |
CN113566821A (zh) * | 2021-06-28 | 2021-10-29 | 江南造船(集团)有限责任公司 | 基于交互式滤波的无人机航向估计方法、系统及电子设备 |
CN113468750A (zh) * | 2021-07-06 | 2021-10-01 | 东南大学 | 一种基于信息融合的天线罩误差斜率估计方法 |
CN113468750B (zh) * | 2021-07-06 | 2022-11-18 | 东南大学 | 一种基于信息融合的天线罩误差斜率估计方法 |
CN117169818B (zh) * | 2023-10-30 | 2024-02-13 | 哈尔滨工业大学(威海) | 一种面向海面机动目标跟踪的雷达波形设计方法 |
CN117169818A (zh) * | 2023-10-30 | 2023-12-05 | 哈尔滨工业大学(威海) | 一种面向海面机动目标跟踪的雷达波形设计方法 |
CN117648521A (zh) * | 2023-11-30 | 2024-03-05 | 常熟理工学院 | 基于角度量测序贯结构状态估计方法、系统及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN107390199B (zh) | 2019-06-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107390199B (zh) | 一种雷达机动目标跟踪波形设计方法 | |
Li et al. | Reduced-dimension MUSIC for angle and array gain-phase error estimation in bistatic MIMO radar | |
CN104809326A (zh) | 一种异步传感器空间配准算法 | |
CN111985093A (zh) | 一种带噪声估计器的自适应无迹卡尔曼滤波状态估计方法 | |
CN101221238B (zh) | 基于高斯均值移动配准的动态偏差估计方法 | |
CN107561503A (zh) | 一种基于多重渐消因子的自适应目标跟踪滤波方法 | |
CN106443661A (zh) | 基于无迹卡尔曼滤波的机动扩展目标跟踪方法 | |
Niedfeldt et al. | Recursive RANSAC: Multiple signal estimation with outliers | |
CN104267375B (zh) | 一种外辐射源雷达网误差自配准方法 | |
Wan et al. | Tracking maneuvering target with angle-only measurements using IMM algorithm based on CKF | |
CN105954712A (zh) | 联合无线电信号复包络和载波相位信息的多目标直接定位方法 | |
CN111046591B (zh) | 传感器幅相误差与目标到达角度的联合估计方法 | |
CN106525042A (zh) | 一种基于蚁群与扩展卡尔曼滤波相结合的多auv协同定位方法 | |
CN109581281A (zh) | 基于到达时间差和到达频率差的移动目标定位方法 | |
CN107994885A (zh) | 一种同时估计未知输入和状态的分布式融合滤波方法 | |
CN109752710A (zh) | 一种基于稀疏贝叶斯学习的快速目标角度估计方法 | |
CN103323815A (zh) | 一种基于等效声速的水下声学定位方法 | |
CN104699965A (zh) | 基于旋转干涉仪的近场源参数估计方法 | |
CN102706345A (zh) | 一种基于衰减记忆序贯检测器的机动目标跟踪方法 | |
CN103487800B (zh) | 基于残差反馈的多模型高速高机动目标跟踪方法 | |
CN106383342B (zh) | 一种基于存在测量误差的阵列流形先验知识的稳健stap方法 | |
CN104181513A (zh) | 一种雷达天线阵元位置的校正方法 | |
CN107727097A (zh) | 基于机载分布式位置姿态测量系统的信息融合方法和装置 | |
CN111291471A (zh) | 一种基于l1正则无迹变换的约束多模型滤波方法 | |
CN104021285A (zh) | 一种具有最优运动模式切换参数的交互式多模型目标跟踪方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |