发明名称 一种SAR聚束工作模式的瞄准点的获取方法
摘要 本发明公开了一种SAR聚束工作模式的瞄准点的获取方法,包括以下几个步骤:(1)获取雷达飞行轨迹数据及接收信号功率数据,初始化瞄准点中心坐标及距离范围,构建二维搜索点矩阵。(2)获取接收信号功率误差平方和。(3)获取瞄准点坐标。本发明提供的方法可以在瞄准点位置未知情况下,通过测量接收雷达发射信号功率来得出SAR聚束模式瞄准点坐标位置。本发明提供的方法可以通过迭代,从而获取较为精确的SAR聚束模式瞄准点坐标。
申请公布号 CN102331576A 申请公布日期 2012.01.25
申请号 CN201110213627.1 申请日期 2011.07.28
申请人 北京航空航天大学 发明人 孙兵;王鹏;唐小明;李春升
分类号 G01S13/90(2006.01)I 主分类号 G01S13/90(2006.01)I
代理机构 北京永创新实专利事务所 11121 代理人 赵文利
主权项 1.一种SAR聚束工作模式的瞄准点的获取方法,其特征在于,包括以下几个步骤:步骤一、获取雷达飞行轨迹数据及接收信号功率数据,初始化瞄准点中心坐标及距离范围,构建二维搜索点矩阵;具体为:(a)设定地面接收站R点坐标为(Xr,Yr,0),其中,Xr表示R点X轴坐标,Yr表示R点Y轴坐标,Xr>0,Yr>0,根据采样时间间隔划分得到采样时间序列t<sub>k</sub>,k为正整数,k=1、2…N;获取雷达飞行轨迹数据,得到对应于采样时间序列t<sub>k</sub>的雷达P点坐标(Xp<sub>k</sub>,Yp<sub>k</sub>,Zp<sub>k</sub>),其中,Xp<sub>k</sub>表示P点在t<sub>k</sub>时刻X轴坐标,Yp<sub>k</sub>表示P点在t<sub>k</sub>时刻Y轴坐标,Zp<sub>k</sub>表示P点在t<sub>k</sub>时刻Z轴坐标,获取对应于采样时间序列t<sub>k</sub>的接收信号功率数组Pm<sub>k</sub>;(b)设定迭代次数,并初始化迭代计数器;(c)初始化瞄准点中心坐标和距离范围,将地面接收站R点坐标值设为瞄准点中心坐标初始值(Xc,Yc,0),其中,Xc表示瞄准点中心X轴坐标,Yc表示瞄准点中心Y轴坐标,并且,设定距离范围初始值为d;(d)在XOY平面中,将X轴方向上点(Xc-d/2)至点(Xc+d/2)之间的距离等分为(M-1)份,得到M个X轴上坐标点,将Y轴方向上点(Yc-d/2)至点(Yc+d/2)之间的距离等分为(M-1)份,得到M个Y轴上坐标点,M个X轴上坐标点和M个Y轴上坐标点组成M×M个搜索坐标点T(Xt<sub>i</sub>,Yt<sub>j</sub>,0),i,j为正整数,i,j=1、2…M,Xt<sub>i</sub>为T点X轴坐标,Yt<sub>j</sub>为T点Y轴坐标;步骤二、获取接收信号功率误差平方和;具体为:(a)获取地面接收站的接收距离,地面接收站的接收距离序列R<sub>k</sub>为:<maths num="0001"><![CDATA[<math><mrow><msub><mi>R</mi><mi>k</mi></msub><mo>=</mo><msqrt><msup><mrow><mo>(</mo><mi>Xr</mi><mo>-</mo><msub><mi>Xp</mi><mi>k</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>Yr</mi><mo>-</mo><msub><mi>Yp</mi><mi>k</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><msub><mi>Zp</mi><mi>k</mi></msub><mn>2</mn></msup></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>(b)设辅助点E坐标为(Xr,Yt<sub>j</sub>,0),辅助点F坐标为(Xt<sub>i</sub>,Yr,0),在每一个设定的瞄准点T坐标(Xt<sub>i</sub>,Yt<sub>j</sub>,0)下,获取对应雷达天线方位向离轴角序列θ<sub>k</sub>(i,j),距离向离轴角序列<img file="FDA0000079235610000012.GIF" wi="168" he="49" />具体如下所示:<maths num="0002"><![CDATA[<math><mrow><msub><mrow><mi>a</mi><mn>1</mn></mrow><mi>k</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><msqrt><msup><mrow><msub><mrow><mo>(</mo><mi>Xt</mi></mrow><mi>i</mi></msub><mo>-</mo><msub><mi>Xp</mi><mi>k</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>Yt</mi><mi>j</mi></msub><mo>-</mo><msub><mi>Yp</mi><mi>k</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><msub><mi>Zp</mi><mi>k</mi></msub><mn>2</mn></msup></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><msub><mrow><mi>b</mi><mn>1</mn></mrow><mi>k</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><msqrt><mrow><mo>(</mo><mi>Xr</mi></mrow><msup><mrow><mo>-</mo><msub><mi>Xp</mi><mi>k</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>Yt</mi><mi>j</mi></msub><mo>-</mo><msub><mi>Yp</mi><mi>k</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><msub><mi>Zp</mi><mi>k</mi></msub><mn>2</mn></msup></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths>c1<sub>k</sub>(i,j)=|Xt<sub>i</sub>-Xr|                      (4)式中:a1<sub>k</sub>(i,j)为PT两点间距离,b1<sub>k</sub>(i,j)为PE两点间距离,c1<sub>k</sub>(i,j)为ET两点间距离;<maths num="0004"><![CDATA[<math><mrow><msub><mrow><mi>a</mi><mn>2</mn></mrow><mi>k</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><msqrt><msup><mrow><msub><mrow><mo>(</mo><mi>Xt</mi></mrow><mi>i</mi></msub><mo>-</mo><msub><mi>Xp</mi><mi>k</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>Yt</mi><mi>j</mi></msub><mo>-</mo><msub><mi>Yp</mi><mi>k</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><msub><mi>Zp</mi><mi>k</mi></msub><mn>2</mn></msup></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><msub><mrow><mi>b</mi><mn>2</mn></mrow><mi>k</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><msqrt><msup><mrow><msub><mrow><mo>(</mo><mi>Xt</mi></mrow><mi>i</mi></msub><mo>-</mo><msub><mi>Xp</mi><mi>k</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>Yr</mi><mo>-</mo><msub><mi>Yp</mi><mi>k</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><msub><mi>Zp</mi><mi>k</mi></msub><mn>2</mn></msup></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths>c2<sub>k</sub>(i,j)=|Yt<sub>j</sub>-Yr|                    (7)式中:a2<sub>k</sub>(i,j)为PT两点间距离,b2<sub>k</sub>(i,j)为PF两点间距离,c2<sub>k</sub>(i,j)为FT两点间距离;<maths num="0006"><![CDATA[<math><mrow><msub><mi>&theta;</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><msup><mi>cos</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mfrac><mrow><mo>(</mo><msub><mrow><mi>a</mi><mn>1</mn></mrow><mi>k</mi></msub><msup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msub><mrow><mi>b</mi><mn>1</mn></mrow><mi>k</mi></msub><msup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><msub><mrow><mi>c</mi><mn>1</mn></mrow><mi>k</mi></msub><msup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>)</mo></mrow><mrow><mn>2</mn><mo>&CenterDot;</mo><msub><mrow><mi>a</mi><mn>1</mn></mrow><mi>k</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mrow><mi>b</mi><mn>1</mn></mrow><mi>k</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow></math>]]></maths><img file="FDA0000079235610000024.GIF" wi="1520" he="97" />(c)根据功率方向增益图及离轴角序列θ<sub>k</sub>(i,j),<img file="FDA0000079235610000025.GIF" wi="167" he="49" />得到雷达在PR方向上的增益G<sub>k</sub>(i,j),,通过式(10)得到接收信号功率数组Pr<sub>k</sub>(i,j):<maths num="0007"><![CDATA[<math><mrow><msub><mi>Pr</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mi>Pt</mi><mo>&CenterDot;</mo><msub><mi>G</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mfrac><mn>1</mn><mrow><mn>4</mn><mi>&pi;</mi><msup><msub><mi>R</mi><mi>k</mi></msub><mn>2</mn></msup></mrow></mfrac><mo>&CenterDot;</mo><mfrac><mrow><msup><mi>&lambda;</mi><mn>2</mn></msup><mo>&CenterDot;</mo><mi>Gr</mi></mrow><mrow><mn>4</mn><mi>&pi;</mi></mrow></mfrac><mo>&CenterDot;</mo><mfrac><mn>1</mn><mi>L</mi></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中:Pt为雷达发射功率,Gr接收天线增益,λ为雷达波长,L为系统损耗;(d)将M×M个搜索坐标点中每个点的接收信号功率数组和地面接收信号功率数组比较,获取每个搜索坐标点接收信号功率误差平方和δ(i,j):<maths num="0008"><![CDATA[<math><mrow><mi>&delta;</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>&Sigma;</mi><mi>k</mi></munder><msup><mrow><mo>|</mo><msub><mi>Pr</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>Pm</mi><mi>k</mi></msub><mo>|</mo></mrow><mn>2</mn></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow></math>]]></maths>步骤三、获取瞄准点坐标;逐个搜索M×M个搜索坐标点中,得到误差平方和最小值对应的坐标点即为雷达聚束模式瞄准点,更新迭代参数,获取新的瞄准点坐标,进行迭代,当达到迭代次数后,所得到的瞄准点坐标即为最终的瞄准点坐标。
地址 100191 北京市海淀区学院路37号