发明名称 一种基于稀疏约束的波束展宽和旁瓣抑制方法
摘要 该发明公开了一种基于稀疏约束的波束展宽和旁瓣抑制方法,属于接收无线传输信号的天线阵列信号处理范畴,特别是一种利用稀疏重构来实现波束展宽和旁瓣抑制从而使天线阵优化的方法。利用稀疏约束算法实现波束展宽的同时,能够得到较低的旁瓣电平,进而增强了对干扰信号的抑制能力。通过建立天线阵模型,确定导向矢量参数,再利用超完备字典表示观测和干扰方向信息的矩阵,随后初始化波束形成器的权值向量,并对其做稀疏迭代,已得到满足条件的权值向量,从而使输出峰值在原有基础上得到有更低的旁瓣。相对于其他的算法,本算法在实现主瓣窄的同时做到超分辨,旁瓣低的同时又提高稳定性。
申请公布号 CN104615854A 申请公布日期 2015.05.13
申请号 CN201510003601.2 申请日期 2015.01.05
申请人 电子科技大学 发明人 张瑛;赵丹旎;王婷静;陈垠江;康宁;赵华鹏
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 电子科技大学专利中心 51203 代理人 张杨
主权项 一种基于稀疏约束的波束展宽和旁瓣抑制方法,该方法的具体步骤为:步骤1、参数设置:步骤1‑1、根据所给阵列结构设置导向矢量的参数γ、λ、p,其中γ和λ为后续步骤中需要的拉格朗日乘子,λ为确定稀疏约束的有效性,γ为调节无失真响应约束及稀疏约束的影响,p调节波束的稀疏度;步骤1‑2、根据观测方向和旁瓣区域生成<img file="FDA0000651839200000011.GIF" wi="85" he="69" /><img file="FDA0000651839200000012.GIF" wi="394" he="80" />其中A表示L×N矩阵,包括了导向矢量角度范围,包含了可能存在的干扰方向;L是表示传感器数目,N代表了超过干扰范围的空间抽样数;而α=γ/λ为拉格朗日乘数,a(θ<sub>0</sub>)为信号源的方向向量;步骤2、初始化:令迭代指数i=0,初始权值为w(0),其中w表示空域滤波器的权向量;w(i)表示w在第i次迭代时的当前权值;步骤3、迭代:步骤3‑1、计算Π(w(i))的值;其中<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>&Pi;</mi><mrow><mo>(</mo><mi>w</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>=</mo><mi>diag</mi><mo>{</mo><msup><mrow><mo>|</mo><msub><mrow><mo>(</mo><msup><mover><mi>A</mi><mo>&OverBar;</mo></mover><mi>H</mi></msup><mi>w</mi><mo>-</mo><msup><mi>d</mi><mo>*</mo></msup><mo>)</mo></mrow><mn>1</mn></msub><mo>|</mo></mrow><mrow><mi>p</mi><mo>-</mo><mn>2</mn></mrow></msup><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msup><mrow><mo>|</mo><msub><mrow><mo>(</mo><msup><mover><mi>A</mi><mo>&OverBar;</mo></mover><mi>H</mi></msup><mi>w</mi><mo>-</mo><msup><mi>d</mi><mo>*</mo></msup><mo>)</mo></mrow><mi>N</mi></msub><mo>|</mo></mrow><mrow><mi>p</mi><mo>-</mo><mn>2</mn></mrow></msup><mo>}</mo><mo>,</mo></mrow>]]></math><img file="FDA0000651839200000014.GIF" wi="1129" he="116" /></maths>而d=[0|α·1],d*表示d的共轭复数,<img file="FDA0000651839200000015.GIF" wi="80" he="70" />是矩阵<img file="FDA0000651839200000016.GIF" wi="60" he="73" />的共轭转置;步骤3‑2、通过式子<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>w</mi><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><mover><mi>&lambda;</mi><mo>~</mo></mover><msup><mrow><mo>(</mo><msub><mi>R</mi><mi>xx</mi></msub><mo>+</mo><mover><mi>&lambda;</mi><mo>~</mo></mover><mover><mi>A</mi><mo>&OverBar;</mo></mover><mi>&Pi;</mi><mrow><mo>(</mo><mi>w</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>)</mo></mrow><msup><mover><mi>A</mi><mo>&OverBar;</mo></mover><mi>H</mi></msup><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mover><mi>A</mi><mo>&OverBar;</mo></mover><mi>&Pi;</mi><mrow><mo>(</mo><mi>w</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>)</mo></mrow><msup><mi>d</mi><mo>*</mo></msup></mrow>]]></math><img file="FDA0000651839200000017.GIF" wi="1069" he="102" /></maths>来更新的w的,其中<img file="FDA0000651839200000018.GIF" wi="179" he="84" />R<sub>xx</sub>是接收信号x(k)的自协方差矩阵;步骤3‑3、根据L‑曲线确定新的λ值,即<img file="FDA0000651839200000019.GIF" wi="201" he="152" />其中||A||代表的是矩阵A的范数,<img file="FDA00006518392000000110.GIF" wi="420" he="161" />a<sub>ij</sub>为矩阵A中的第i行第j列元素;步骤3‑4、在通过式子<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>&lambda;</mi><mo>=</mo><mi>&gamma;</mi><msup><mrow><mo>(</mo><mfrac><mrow><mi>w</mi><msup><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mi>H</mi></msup><msub><mi>R</mi><mi>xx</mi></msub><mi>w</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></mrow><mrow><mrow><mo>|</mo><mi>w</mi><msup><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mi>H</mi></msup></mrow><mi>a</mi><msup><mrow><mrow><mo>(</mo><msub><mi>&theta;</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>-</mo><mn>1</mn><mo>|</mo></mrow><mrow><mi>p</mi><mo>-</mo><mn>2</mn></mrow></msup><mrow><mo>(</mo><mi>w</mi><msup><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mi>H</mi></msup><mi>a</mi><mrow><mo>(</mo><msub><mi>&theta;</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></mfrac><mo>)</mo></mrow><mrow><mn>1</mn><mo>-</mo><mi>p</mi></mrow></msup></mrow>]]></math><img file="FDA00006518392000000111.GIF" wi="938" he="245" /></maths>获得新的γ值;步骤3‑5、通过i=i+1使w的值改变;步骤4、如果结果符合要求,算法终止,若不满足,则返回步骤3利用新获得的λ和γ继续计算;算法停止的准则可以为迭代次数达到根据经验预设的次数;也可以将<img file="FDA0000651839200000021.GIF" wi="502" he="104" />与根据经验预设的阈值比较,若其小于阈值,则停止计算。步骤5、根据得到的迭代后的空域滤波器权向量w,获得阵列天线接收信号方向图增益20lg|w<sup>H</sup>a(θ)|。
地址 611731 四川省成都市高新区(西区)西源大道2006号