发明名称 一种雷达跟踪状态下预测目标RCS的方法
摘要 本发明公开了一种雷达跟踪状态下预测目标RCS的方法,包括:步骤1,设定目标雷达截面积RCS的预测滤波器的阶数为M,步骤2,雷达接收目标自n-M+1时刻至n时刻的回波s<sub>i</sub>;并且记录自n-M+1时刻至n时刻目标的距离,第n个跟踪时刻目标的速度,以及自n-M+1时刻至n时刻目标速度与目标距离的夹角,步骤3,求得在自n-M+1时刻至n时刻的目标复幅度的估计值,步骤4,得到目标在第n+1时刻的目标速度与目标距离的夹角预测值,步骤5,求得目标RCS值的自相关矩阵和目标RCS值的互相关列向量,步骤6,求得目标RCS的预测滤波器系数,步骤7,得到在第n+1时刻的目标RCS的预测值。
申请公布号 CN104076342B 申请公布日期 2016.05.18
申请号 CN201410289649.X 申请日期 2014.06.25
申请人 西安电子科技大学 发明人 刘宏伟;吴梦;臧会凯;周生华;王英华
分类号 G01S7/41(2006.01)I 主分类号 G01S7/41(2006.01)I
代理机构 西安睿通知识产权代理事务所(特殊普通合伙) 61218 代理人 惠文轩
主权项 一种雷达跟踪状态下预测目标RCS的方法,其特征在于,包括以下步骤:步骤1,设定目标的雷达截面积RCS的预测滤波器的阶数M的取值范围为1至tan(2arcsin(c/(4fd)))r/(v<sub>max</sub>T<sub>r</sub>)之间的整数,其中,tan(·)表示正切函数,arcsin(·)表示反正弦函数,c表示光速,f表示雷达的载频,d表示目标长度的平均值,r表示目标距离的平均值,v<sub>max</sub>表示目标的最大速度,T<sub>r</sub>表示脉冲重复间隔;步骤2,雷达接收目标自n‑M+1时刻至n时刻的回波s<sub>i</sub>;并且记录自n‑M+1时刻至n时刻目标的距离r<sub>i</sub>,第n个跟踪时刻目标的速度v<sub>n</sub>,以及自n‑M+1时刻至n时刻目标速度与目标距离的夹角α<sub>i</sub>,i=n‑M+1,n‑M+2,…,n,M为预测滤波器的阶数,n为自然数;步骤3,根据目标自n‑M+1时刻至n时刻的目标回波s<sub>i</sub>和目标距离r<sub>i</sub>,求得在自n‑M+1时刻至n时刻的目标复幅度的估计值η<sub>i</sub>=s<sub>i</sub>×(r<sub>i</sub>/r<sub>1</sub>)<sup>2</sup>,r<sub>i</sub>表示第i时刻目标的距离,r<sub>1</sub>为第1个时刻目标的距离,i=n‑M+1,n‑M+2,…,n,M为预测滤波器的阶数;×表示乘积;步骤4,根据目标在第n时刻的目标距离r<sub>n</sub>、第n时刻的目标的速度v<sub>n</sub>、第n时刻的目标速度与目标距离的夹角α<sub>n</sub>,所述目标速度与目标距离的夹角为目标速度方向与目标到雷达的矢量方向所成的夹角,按照下式得到目标在第n+1时刻的目标速度与目标距离夹角的预测值<img file="FDA0000871627030000011.GIF" wi="111" he="70" /><maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mover><mi>&alpha;</mi><mo>^</mo></mover><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><mi>&pi;</mi><mo>-</mo><mi>a</mi><mi>r</mi><mi>c</mi><mi>c</mi><mi>o</mi><mi>s</mi><mrow><mo>(</mo><mfrac><mrow><msub><mi>T</mi><mi>r</mi></msub><msub><mi>v</mi><mi>n</mi></msub><mo>-</mo><msub><mi>r</mi><mi>n</mi></msub><mi>c</mi><mi>o</mi><mi>s</mi><mrow><mo>(</mo><msub><mi>&alpha;</mi><mi>n</mi></msub><mo>)</mo></mrow></mrow><msqrt><mrow><msubsup><mi>r</mi><mi>n</mi><mn>2</mn></msubsup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>T</mi><mi>r</mi></msub><msub><mi>v</mi><mi>n</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><mn>2</mn><msub><mi>T</mi><mi>r</mi></msub><msub><mi>v</mi><mi>n</mi></msub><msub><mi>r</mi><mi>n</mi></msub><mi>c</mi><mi>o</mi><mi>s</mi><mrow><mo>(</mo><msub><mi>&alpha;</mi><mi>n</mi></msub><mo>)</mo></mrow></mrow></msqrt></mfrac><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000871627030000012.GIF" wi="1022" he="207" /></maths>其中,arccos(·)表示反余弦函数,cos(·)表示余弦函数,T<sub>r</sub>表示脉冲重复间隔;步骤5,根据自n‑M+1时刻至n时刻的目标速度与目标距离的夹角α<sub>i</sub>,i=n‑M+1,n‑M+2,…,n,求得目标RCS值的自相关矩阵R,自相关矩阵R是一个M行M列矩阵,自相关矩阵R形式为:<img file="FDA0000871627030000021.GIF" wi="1550" he="359" />其中,ρ(·)表示相关系数函数,相关系数函数的表达式为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>&rho;</mi><mrow><mo>(</mo><msub><mi>u</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>=</mo><mi>j</mi><mi>i</mi><mi>n</mi><mi>c</mi><mrow><mo>(</mo><mfrac><mrow><mn>4</mn><mi>f</mi><mi>d</mi><mi>s</mi><mi>i</mi><mi>n</mi><mrow><mo>(</mo><msub><mi>u</mi><mn>1</mn></msub><mo>/</mo><mn>2</mn><mo>)</mo></mrow></mrow><mi>c</mi></mfrac><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000871627030000022.GIF" wi="605" he="159" /></maths>式中,u<sub>1</sub>表示相关系数函数的输入变量,f表示雷达的载频,d表示目标长度的平均值,c表示光速,jinc(·)表示jinc函数,jinc函数定义为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>j</mi><mi>i</mi><mi>n</mi><mi>c</mi><mrow><mo>(</mo><msub><mi>u</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mn>2</mn><mo>&times;</mo><msub><mi>J</mi><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>&pi;u</mi><mn>2</mn></msub><mo>)</mo></mrow></mrow><mrow><msub><mi>&pi;u</mi><mn>2</mn></msub></mrow></mfrac></mrow>]]></math><img file="FDA0000871627030000023.GIF" wi="465" he="149" /></maths>其中,J<sub>1</sub>(·)表示第一类一阶贝塞尔函数,u<sub>2</sub>为jinc函数的输入变量;根据自n‑M+1时刻至n时刻目标速度与目标距离的夹角α<sub>i</sub>和第n+1时刻目标速度与目标距离夹角的预测值<img file="FDA0000871627030000026.GIF" wi="111" he="70" />i=n‑M+1,n‑M+2,…,n,求得目标RCS值的互相关列向量P,互相关向量P是一个M维的列向量,互相关向量P具体形式为:P=[‑ρ(α<sub>n+1</sub>‑α<sub>n</sub>) ‑ρ(α<sub>n+1</sub>‑α<sub>n‑1</sub>) … ‑ρ(α<sub>n+1</sub>‑α<sub>n‑M+1</sub>)]<sup>T</sup>其中,[·]<sup>T</sup>表示转置,ρ(·)表示相关系数函数;步骤6,根据自相关矩阵R和互相关列向量P,求得目标RCS的预测滤波器系数W=R<sup>‑1</sup>P,预测滤波器系数W的表达式为:W=[w<sub>1</sub>,…,w<sub>k</sub>,…,w<sub>M</sub>]<sup>T</sup>其中,w<sub>k</sub>为预测滤波器系数W的第k个元素,k=1,2,…,M;M表示预测滤波器的阶数;(·)<sup>‑1</sup>表示矩阵求逆运算;步骤7,根据预测滤波器系数W以及自n‑M+1时刻至n时刻目标复幅度的估计值η<sub>i</sub>,i=n‑M+1,n‑M+2,…,n,得到在第n+1时刻的目标复幅度的预测值<img file="FDA0000871627030000024.GIF" wi="454" he="135" />求取第n+1时刻目标RCS的预测值<img file="FDA0000871627030000025.GIF" wi="357" he="79" />其中,Σ表示求和运算,k=1,2,…,M,w<sub>k</sub>为预测滤波器系数W的第k个元素,(·)<sup>*</sup>表示取共轭。
地址 710071 陕西省西安市太白南路2号