发明名称 基于雷达旋转发射阵列的微波关联成像方法
摘要 本发明公开了一种基于雷达旋转发射阵列的微波关联成像方法,涉及雷达信息获取与处理技术领域,其步骤为:步骤1,设定雷达旋转发射天线;步骤2,设定雷达的探测面,将该探测面划分成多个探测单元;步骤3,求取时间延迟;步骤4,利用时间延迟求取导向矢量,得到阵列流形矩阵;步骤5,设定回波信号向量;步骤6,构建辐射场向量;步骤7,建立雷达接收的M个脉冲下的原始回波向量与辐射场向量的关联形式;步骤8,在目标函数的约束条件求解目标反射系数向量的目标函数,得到目标反射系数向量;步骤9,得到重排的目标反射系数向量矩阵,即目标的成像矩阵。本发明提高辐射场空时随机性和空间自由度,实现对大场景目标的高分辨率成像。
申请公布号 CN104199028A 申请公布日期 2014.12.10
申请号 CN201410446365.7 申请日期 2014.09.03
申请人 西安电子科技大学 发明人 李军;李小敏;朱圣棋;邵自立;廖桂生;吕利;马玉芳
分类号 G01S13/89(2006.01)I;G06F19/00(2011.01)I;G01S7/41(2006.01)I 主分类号 G01S13/89(2006.01)I
代理机构 西安睿通知识产权代理事务所(特殊普通合伙) 61218 代理人 惠文轩
主权项 一种基于雷达旋转发射阵列的微波关联成像方法,其特征在于,包括以下步骤:步骤1,雷达旋转发射阵列是具有N个阵元的二维面阵,雷达旋转发射阵列的阵元分布在一个圆面上,雷达的接收机位于圆面的圆心,阵元所在圆面绕过圆心的轴线以角速度ω转动,所述圆心的轴线垂直于圆面;N为自然数;步骤2,雷达旋转发射阵列向探测面发射M个脉冲的微波信号,并且接收M个脉冲下的原始回波数据Y;M为自然数;将所述探测面划分成E×F个探测单元,探测单元总数Q=E×F,其中,E表示横向探测单元数,F表示纵向探测单元数,E、F分别为自然数,且Q&gt;M,M表示雷达旋转发射阵列发射脉冲数;步骤3,在第m个脉冲下第n个阵元向第q个探测单元发射微波信号表示为s<sub>mn</sub>(t),雷达接收机接收第q个探测单元的回波信号表示为s<sub>mn</sub>(t‑τ<sub>mnq</sub>),求取发射微波信号s<sub>mn</sub>(t)与回波信号s<sub>mn</sub>(t‑τ<sub>mnq</sub>)之间的时间延迟τ<sub>mnq</sub>,其中,t表示时间,m=1,2,...,M,n=1,2,...,N,q=1,2,...,Q,M表示雷达旋转发射阵列发射脉冲数,N表示雷达旋转发射阵列的发射阵元数,Q表示探测单元总数;步骤4,利用时间延迟τ<sub>mnq</sub>求取在第m个脉冲下N个阵元相对第q个探测单元的导向矢量V<sub>mq</sub>,V<sub>mq</sub>维数为N×1;进而得到在第m个脉冲下N个阵元相对Q个探测单元的阵列流形矩阵V<sub>m</sub>,V<sub>m</sub>维数为N×Q;步骤5,设定第m个脉冲下第n个阵元向第q个探测单元发射的微波信号s<sub>mn</sub>(t)等于接收机接收的第q个探测单元的回波信号s<sub>mn</sub>(t‑τ<sub>mnq</sub>),即s<sub>mn</sub>(t‑τ<sub>mnq</sub>)=s<sub>mn</sub>(t),则第m个脉冲下的回波信号向量:s<sub>m</sub>(t)=[s<sub>m1</sub>(t),s<sub>m2</sub>(t),...,s<sub>mn</sub>(t),...,s<sub>mN</sub>(t)]<sub>1×N</sub>;m=1,2,...,M,n=1,2,...,N,q=1,2,...,Q;步骤6,利用第m个脉冲下的N个阵元相对Q个探测单元的阵列流形矩阵V<sub>m</sub>和第m个脉冲下的回波信号向量s<sub>m</sub>(t),构建第m个脉冲下的辐射场向量A<sub>m</sub>=s<sub>m</sub>(t)·V<sub>m</sub>,维数为1×Q进而得到M个脉冲下的辐射场矩阵<img file="FDA0000565157160000011.GIF" wi="604" he="101" />A维数为M×Q,其中[·]<sup>T</sup>表示转置;步骤7,设定目标反射系数向量σ的向量化形式为[σ<sub>1</sub>,σ<sub>2</sub>,...,σ<sub>q</sub>,...,σ<sub>Q</sub>]<sup>T</sup>,其中[·]<sup>T</sup>表示转置,σ<sub>q</sub>表示第q个探测单元的目标反射系数,q=1,2,…,Q;利用M个脉冲下的辐射场矩阵A和目标反射系数向量σ将雷达接收的M个脉冲下的原始回波向量Y写成如下形式:Y=Aσ+w其中,A表示M个脉冲下的Q个探测单元的辐射场矩阵,w表示设定的维数为M×1的高斯白噪声矢量;步骤8,建立目标反射系数向量σ的目标函数,将Y=Aσ+w作为目标函数的约束条件,在目标函数的约束条件下求解目标反射系数向量σ的目标函数,得到目标反射系数向σ,σ=[σ<sub>1</sub>,σ<sub>2</sub>,...,σ<sub>q</sub>,...,σ<sub>Q</sub>]<sup>T</sup>=[σ<sub>1</sub>,σ<sub>2</sub>,...,σ<sub>E</sub>,σ<sub>E+1</sub>,σ<sub>E+2</sub>,...,σ<sub>2E</sub>,...,σ<sub>(F‑1)E+1</sub>,σ<sub>(F‑1xE+2</sub>,...,σ<sub>Q</sub>]<sup>T</sup>,其中[·]<sup>T</sup>表示转置,Q=E×F,E表示横向探测单元数,F表示纵向探测单元数;步骤9,将目标反射系数向量σ由Q×1维排列为E×F维矩阵,得到重排的目标反射系数矩阵<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mover><mi>&sigma;</mi><mo>&OverBar;</mo></mover><mo>=</mo><msub><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>&sigma;</mi><mn>1</mn></msub></mtd><mtd><msub><mi>&sigma;</mi><mrow><mi>E</mi><mo>+</mo><mn>1</mn></mrow></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mi>&sigma;</mi><mrow><mrow><mo>(</mo><mi>F</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>E</mi><mo>+</mo><mn>1</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>&sigma;</mi><mn>2</mn></msub></mtd><mtd><msub><mi>&sigma;</mi><mrow><mi>E</mi><mo>+</mo><mn>2</mn></mrow></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mi>&sigma;</mi><mrow><mrow><mo>(</mo><mi>F</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>E</mi><mo>+</mo><mn>2</mn></mrow></msub></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>&sigma;</mi><mi>E</mi></msub></mtd><mtd><msub><mi>&sigma;</mi><mrow><mn>2</mn><mi>E</mi></mrow></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mi>&sigma;</mi><mi>Q</mi></msub></mtd></mtr></mtable></mfenced><mrow><mi>E</mi><mo>&times;</mo><mi>F</mi></mrow></msub><mo>,</mo></mrow>]]></math><img file="FDA0000565157160000021.GIF" wi="728" he="316" /></maths>重排的目标反射系数矩阵<img file="FDA0000565157160000022.GIF" wi="49" he="75" />即目标的成像矩阵。
地址 710071 陕西省西安市太白南路2号