发明名称 基于非凸优化的MIMO雷达动目标检测方法
摘要 本发明公开了一种基于非凸优化的MIMO雷达动目标检测方法,根据雷达参数对观测区域进行距离-角度-多普勒网格划分;根据雷达发射、接收阵列的位置、雷达观测区域参数计算发射、接收引导矢量;根据发射信号波形、引导矢量以及距离单元时延变换矩阵构造一个稀疏字典矩阵;将接收天线阵列接收到的回波信号进行序列化;根据以上模型将MIMO雷达动目标检测构造成一个<img file="dest_path_image002.GIF" wi="26" he="25" />范数最小化约束的优化问题;利用重赋权<img file="dest_path_image004.GIF" wi="17" he="25" />范数正则化方法求解上述优化问题,得到MIMO雷达动目标的距离-角度-多普勒成像,并检测出观测区域中的动目标。本发明能够获得比凸优化算法更加准确的动目标检测结果,检测结果具有更高的分辨率。
申请公布号 CN103744076B 申请公布日期 2016.07.06
申请号 CN201310726840.1 申请日期 2013.12.25
申请人 河海大学 发明人 鹿浩;陈亮;王佳希;胡晓雯;曹宁
分类号 G01S13/50(2006.01)I;G01S13/89(2006.01)I;G01S7/41(2006.01)I 主分类号 G01S13/50(2006.01)I
代理机构 南京经纬专利商标代理有限公司 32200 代理人 朱小兵
主权项 一种基于非凸优化的MIMO雷达动目标检测方法,所述MIMO雷达是一个收发天线阵列均为均匀线性阵的窄带相干MIMO雷达系统,且收发阵列中第一个天线位置相同,其特征在于,检测方法包括如下步骤:步骤1,根据雷达参数对观测区域进行距离‑角度‑多普勒网格划分,将远场雷达观测场景划分成N<sub>r</sub>个距离分辨单元,N<sub>a</sub>个角度分辨单元,N<sub>d</sub>个关心的多普勒频移分辨单元,用度数来描述目标多普勒频移;步骤2,根据雷达发射、接收阵列的位置和雷达观测区域参数计算发射、接收引导矢量,具体如下:设定雷达系统具有M<sub>t</sub>个发射天线和M<sub>r</sub>个接收天线,则第a个角度分辨单元的发射和接收引导矢量分别表示如下:<maths num="0001"><math><![CDATA[<mrow><msub><mi>a</mi><mi>a</mi></msub><mo>=</mo><msup><mfenced open = "[" close = "]"><mtable><mtr><mtd><mi>1</mi></mtd><mtd><msup><mi>e</mi><mrow><mo>-</mo><mfrac><mrow><mi>j</mi><mn>2</mn><msub><mi>&pi;&Delta;</mi><mi>t</mi></msub><mi>s</mi><mi>i</mi><mi>n</mi><mrow><mo>(</mo><msub><mi>&theta;</mi><mi>a</mi></msub><mo>)</mo></mrow></mrow><msub><mi>&lambda;</mi><mn>0</mn></msub></mfrac></mrow></msup></mtd><mtd><mo>...</mo></mtd><mtd><msup><mi>e</mi><mrow><mo>-</mo><mfrac><mrow><mi>j</mi><mn>2</mn><mi>&pi;</mi><mrow><mo>(</mo><msub><mi>M</mi><mi>t</mi></msub><mo>-</mo><mn>1</mn><mo>)</mo></mrow><msub><mi>&Delta;</mi><mi>t</mi></msub><mi>s</mi><mi>i</mi><mi>n</mi><mrow><mo>(</mo><msub><mi>&theta;</mi><mi>a</mi></msub><mo>)</mo></mrow></mrow><msub><mi>&lambda;</mi><mn>0</mn></msub></mfrac></mrow></msup></mtd></mtr></mtable></mfenced><mi>T</mi></msup></mrow>]]></math><img file="FDA0000908402900000011.GIF" wi="907" he="184" /></maths><maths num="0002"><math><![CDATA[<mrow><msub><mi>b</mi><mi>a</mi></msub><mo>=</mo><msup><mfenced open = "[" close = "]"><mtable><mtr><mtd><mi>1</mi></mtd><mtd><msup><mi>e</mi><mrow><mo>-</mo><mfrac><mrow><mi>j</mi><mn>2</mn><msub><mi>&pi;&Delta;</mi><mi>r</mi></msub><mi>s</mi><mi>i</mi><mi>n</mi><mrow><mo>(</mo><msub><mi>&theta;</mi><mi>a</mi></msub><mo>)</mo></mrow></mrow><msub><mi>&lambda;</mi><mn>0</mn></msub></mfrac></mrow></msup></mtd><mtd><mo>...</mo></mtd><mtd><msup><mi>e</mi><mrow><mo>-</mo><mfrac><mrow><mi>j</mi><mn>2</mn><mi>&pi;</mi><mrow><mo>(</mo><msub><mi>M</mi><mi>r</mi></msub><mo>-</mo><mn>1</mn><mo>)</mo></mrow><msub><mi>&Delta;</mi><mi>r</mi></msub><mi>s</mi><mi>i</mi><mi>n</mi><mrow><mo>(</mo><msub><mi>&theta;</mi><mi>a</mi></msub><mo>)</mo></mrow></mrow><msub><mi>&lambda;</mi><mn>0</mn></msub></mfrac></mrow></msup></mtd></mtr></mtable></mfenced><mi>T</mi></msup></mrow>]]></math><img file="FDA0000908402900000012.GIF" wi="918" he="187" /></maths>式中,j表示虚数单位;Δ<sub>t</sub>和Δ<sub>r</sub>分别为MIMO雷达发射天线阵元间距和接收天线阵元间距;λ<sub>0</sub>为雷达发射信号工作波长;θ<sub>a</sub>为第a个角度分辨单元相对于阵列标准方向的方位角,a=1,2,…N<sub>a</sub>;步骤3,将接收天线阵列接收到的回波信号进行序列化,矢量化后的回波信号为:y=vec(Y)式中,Y为接收天线阵列接收到的回波信号;<img file="FDA0000908402900000013.GIF" wi="346" he="81" />即Y为一个M<sub>r</sub>×(L+N<sub>r</sub>‑1)维的复数矩阵;C为复数矩阵;Y的定义如下:<maths num="0003"><math><![CDATA[<mrow><mi>Y</mi><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>r</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>r</mi></msub></munderover><munderover><mo>&Sigma;</mo><mrow><mi>a</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>a</mi></msub></munderover><munderover><mo>&Sigma;</mo><mrow><mi>d</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>d</mi></msub></munderover><msub><mi>&alpha;</mi><mrow><mi>r</mi><mo>,</mo><mi>a</mi><mo>,</mo><mi>d</mi></mrow></msub><msub><mi>b</mi><mi>a</mi></msub><msubsup><mi>a</mi><mi>a</mi><mi>T</mi></msubsup><msub><mover><mi>S</mi><mo>~</mo></mover><mi>d</mi></msub><msub><mi>J</mi><mi>r</mi></msub><mo>+</mo><mi>E</mi></mrow>]]></math><img file="FDA0000908402900000015.GIF" wi="651" he="147" /></maths>式中,L为发射信号的长度;{α<sub>r,a,d</sub>}为关心区域内的分辨单元复反射系数;<img file="FDA0000908402900000016.GIF" wi="330" he="83" />为包含观测区域内回波延时的补零发射波形矩阵;E为一个包括杂波、干扰和加性噪声的矩阵;J<sub>r</sub>为一个变换矩阵,通过J<sub>r</sub>描述包含不同距离分辨单元时间延迟的反射波形;r=1,2,…,N<sub>r</sub>;d=1,2,…,N<sub>d</sub>;J<sub>r</sub>的定义如下:<img file="FDA0000908402900000021.GIF" wi="712" he="368" />步骤4,根据发射信号波形、引导矢量以及距离分辨单元时延变换矩阵构造一个稀疏字典矩阵:<maths num="0004"><math><![CDATA[<mrow><mi>A</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>v</mi><mrow><mn>1</mn><mo>,</mo><mn>1</mn><mo>,</mo><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>v</mi><mrow><mn>1</mn><mo>,</mo><mn>1</mn><mo>,</mo><mn>2</mn></mrow></msub></mtd><mtd><mo>...</mo></mtd><mtd><msub><mi>v</mi><mrow><msub><mi>N</mi><mi>r</mi></msub><mo>,</mo><msub><mi>N</mi><mi>a</mi></msub><mo>,</mo><msub><mi>N</mi><mi>d</mi></msub></mrow></msub></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000908402900000022.GIF" wi="667" he="98" /></maths>式中,<maths num="0005"><math><![CDATA[<mrow><msub><mi>v</mi><mrow><mi>r</mi><mo>,</mo><mi>a</mi><mo>,</mo><mi>d</mi></mrow></msub><mo>=</mo><mi>v</mi><mi>e</mi><mi>c</mi><mo>&lsqb;</mo><msub><mi>b</mi><mi>a</mi></msub><msubsup><mi>a</mi><mi>a</mi><mi>T</mi></msubsup><msub><mover><mi>S</mi><mo>~</mo></mover><mi>d</mi></msub><msub><mi>J</mi><mi>r</mi></msub><mo>&rsqb;</mo><mo>&Element;</mo><msup><mi>C</mi><mrow><mo>(</mo><msub><mi>M</mi><mi>r</mi></msub><mo>(</mo><mi>L</mi><mo>+</mo><msub><mi>N</mi><mi>r</mi></msub><mo>-</mo><mn>1</mn><mo>)</mo><mo>)</mo><mo>&times;</mo><mn>1</mn></mrow></msup><mo>;</mo></mrow>]]></math><img file="FDA0000908402900000023.GIF" wi="792" he="97" /></maths>矢量化的回波信号表示为:y=Ax+e式中,<maths num="0006"><math><![CDATA[<mrow><msub><mi>a</mi><mi>a</mi></msub><mo>=</mo><msup><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>&alpha;</mi><mrow><mn>1</mn><mo>,</mo><mn>1</mn><mo>,</mo><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>&alpha;</mi><mrow><mn>1</mn><mo>,</mo><mn>1</mn><mo>,</mo><mn>2</mn></mrow></msub></mtd><mtd><mo>...</mo></mtd><mtd><msub><mi>&alpha;</mi><mrow><msub><mi>N</mi><mi>r</mi></msub><mo>,</mo><msub><mi>N</mi><mi>a</mi></msub><mo>,</mo><msub><mi>N</mi><mi>d</mi></msub></mrow></msub></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>;</mo></mrow>]]></math><img file="FDA0000908402900000024.GIF" wi="699" he="105" /></maths>e=vec(E);步骤5,根据以上模型将MIMO雷达动目标检测构造成一个l<sub>1/2</sub>范数最小化约束的优化问题:<maths num="0007"><math><![CDATA[<mrow><munder><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow><mi>x</mi></munder><mo>{</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>|</mo><mo>|</mo><mi>y</mi><mo>-</mo><mi>A</mi><mi>x</mi><mo>|</mo><msubsup><mo>|</mo><mn>2</mn><mn>2</mn></msubsup><mo>+</mo><mi>&lambda;</mi><mo>|</mo><mo>|</mo><mi>x</mi><mo>|</mo><msubsup><mo>|</mo><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msubsup><mo>}</mo></mrow>]]></math><img file="FDA0000908402900000025.GIF" wi="563" he="147" /></maths>式中,λ为正则化参数;步骤6,利用重赋权l<sub>1</sub>范数正则化方法求解步骤5中的优化问题,得到MIMO雷达动目标的距离‑角度‑多普勒成像,并检测出观测区域中的动目标。
地址 210098 江苏省南京市鼓楼区西康路1号