发明名称 一种基于矩阵恢复的降采样阵列信号处理方法
摘要 本发明属于雷达阵列信号处理技术领域,特别涉及一种基于矩阵恢复的降采样阵列信号处理方法。该基于矩阵恢复的降采样阵列信号处理方法包括以下步骤:S1:利用雷达发射信号,雷达的接收阵列接收回波信号,所述雷达的接收阵列为均匀线阵,构造回波信号的模型;S2:分别恢复出x(t<sub>1</sub>)至x(t<sub>N</sub>),得出整个回波信号X(t);恢复雷达在t<sub>i</sub>时刻的接收信号x(t<sub>i</sub>)包括以下步骤:构造低秩矩阵,通过矩阵恢复方法恢复出x(t<sub>i</sub>)。
申请公布号 CN103885035B 申请公布日期 2015.11.11
申请号 CN201410106241.4 申请日期 2014.03.20
申请人 西安电子科技大学 发明人 朱圣棋;胡海洋;杨东;廖桂生;陶海红;高永婵;张学攀
分类号 G01S7/02(2006.01)I;G01S7/40(2006.01)I 主分类号 G01S7/02(2006.01)I
代理机构 西安睿通知识产权代理事务所(特殊普通合伙) 61218 代理人 惠文轩
主权项 一种基于矩阵恢复的降采样阵列信号处理方法,其特征在于,包括以下步骤:S1:利用雷达发射信号,雷达的接收阵列接收回波信号,所述雷达的接收阵列为均匀线阵,由M个阵元组成;在所述M个阵元中,存在部分失效阵元;所述回波信号的模型为:X(t)=a(θ)s(t)+N(t),其中,N(t)表示设定的观测噪声,X(t)表示回波信号,X(t)=[x(t<sub>1</sub>),...,x(t<sub>i</sub>),...,x(t<sub>N</sub>)],t<sub>1</sub>至t<sub>N</sub>表示长度为N的时间序列,x(t<sub>i</sub>)表示雷达在t<sub>i</sub>时刻的接收信号,i取1至N;s(t)=[s(t<sub>1</sub>),...,s(t<sub>i</sub>),...,s(t<sub>N</sub>)],s(t<sub>i</sub>)为雷达在t<sub>i</sub>时刻的接收信号的幅度;a(θ)为雷达的接收阵列的导向矢量,并且有a(θ)=[a<sub>1</sub>(θ),...,a<sub>m</sub>(θ),...,a<sub>M</sub>(θ)]<sup>H</sup>,H表示矩阵的共轭转置,a<sub>m</sub>(θ)表示雷达的接收阵列中第m个阵元的导向矢量,m取1至M,θ为波达方向;S2:分别恢复出x(t<sub>1</sub>)至x(t<sub>N</sub>),得出整个回波信号X(t);恢复雷达在t<sub>i</sub>时刻的接收信号x(t<sub>i</sub>)包括以下步骤:S21:x(t<sub>i</sub>)表示为x(t<sub>i</sub>)=[x<sub>1</sub>(t<sub>i</sub>),...,x<sub>m</sub>(t<sub>i</sub>),...,x<sub>M</sub>(t<sub>i</sub>)]<sup>T</sup>,其中,T表示矩阵的转置,x<sub>m</sub>(t<sub>i</sub>)为雷达的接收阵列中第m个阵元在t<sub>i</sub>时刻的接收信号;当第m个阵元为失效阵元时,x<sub>m</sub>(t<sub>i</sub>)为未知量,反之,x<sub>m</sub>(t<sub>i</sub>)为已知量;根据x(t<sub>i</sub>),构造对应的矩阵<img file="FDA0000479813310000011.GIF" wi="186" he="84" /><maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mover><mi>x</mi><mo>~</mo></mover><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>=</mo><mo>[</mo><msub><mover><mi>x</mi><mo>~</mo></mover><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mover><mi>x</mi><mo>~</mo></mover><mi>p</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mover><mi>x</mi><mo>~</mo></mover><mi>l</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>]</mo><mo>,</mo></mrow>]]></math><img file="FDA0000479813310000012.GIF" wi="804" he="86" /></maths>其中,<img file="FDA0000479813310000013.GIF" wi="698" he="94" />p取1至l,l为小于M的自然数,k为小于M的自然数,且l和k满足以下关系式:l=M‑k+1;S22:求解以下优化问题:<maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open='' close=''><mtable><mtr><mtd><mi>min</mi></mtd><mtd><msub><mrow><mo>|</mo><mo>|</mo><msub><mover><mi>x</mi><mo>~</mo></mover><mi>est</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>|</mo><mo>|</mo></mrow><mo>*</mo></msub></mtd></mtr><mtr><mtd><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo></mtd><mtd><msub><mrow><mo>|</mo><mo>|</mo><msub><mi>P</mi><mi>&Omega;</mi></msub><mrow><mo>(</mo><msub><mover><mi>x</mi><mo>~</mo></mover><mi>est</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>-</mo><mover><mi>x</mi><mo>~</mo></mover><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mo>|</mo><mo>|</mo></mrow><mi>F</mi></msub><mo>&lt;</mo><mi>&delta;</mi></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000479813310000014.GIF" wi="704" he="201" /></maths>其中,<img file="FDA0000479813310000015.GIF" wi="161" he="84" />表示x(t<sub>i</sub>)的估计值,<img file="FDA0000479813310000016.GIF" wi="207" he="93" />表示求矩阵<img file="FDA0000479813310000017.GIF" wi="148" he="80" />的核范数,||·||<sub>F</sub>表示求矩阵的F范数;Ω是与<img file="FDA0000479813310000018.GIF" wi="118" he="84" />具有相同维度的矩阵,Ω中的元素只有0和1,每个元素1的位置与<img file="FDA0000479813310000021.GIF" wi="112" he="84" />中已知元素的位置相同,每个元素0的位置与<img file="FDA0000479813310000022.GIF" wi="114" he="74" />中未知元素的位置相同;<img file="FDA0000479813310000023.GIF" wi="414" he="78" />表示矩阵<img file="FDA0000479813310000024.GIF" wi="304" he="68" />在Ω上的投影;δ为设置的常数;然后,通过求解上述优化问题恢复雷达在t<sub>i</sub>时刻的接收信号x(t<sub>i</sub>)。
地址 710071 陕西省西安市太白南路2号