主权项 |
一种基于稀疏约束的波束展宽和旁瓣抑制方法,该方法的具体步骤为:步骤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>Π</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>‾</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>‾</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>λ</mi><mo>~</mo></mover><msup><mrow><mo>(</mo><msub><mi>R</mi><mi>xx</mi></msub><mo>+</mo><mover><mi>λ</mi><mo>~</mo></mover><mover><mi>A</mi><mo>‾</mo></mover><mi>Π</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>‾</mo></mover><mi>H</mi></msup><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mover><mi>A</mi><mo>‾</mo></mover><mi>Π</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>λ</mi><mo>=</mo><mi>γ</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>θ</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>θ</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(θ)|。 |