发明名称 一种圆环形阵列恒定束宽波束形成器设计方法
摘要 本发明涉及一种圆环形阵列恒定束宽波束形成器设计方法,针对均匀圆环形阵列提出一种简单而精确的恒定束宽波束形成器设计方法。该方法利用循环矩阵的特性,给出了圆环阵最小平方误差波束图综合问题的精确解,并将最终合成的波束以及最小平方误差简单地表示为子分量叠加之和。当期望波束转换为适当形式时,可得到恒定束宽波束形成器权值向量的闭式表达式。本发明克服了现有技术操作复杂和不够精确的不足:最优解精确地表示为子分量的叠加和,且有闭式的表达式。
申请公布号 CN103903609B 申请公布日期 2016.11.30
申请号 CN201410146364.0 申请日期 2014.04.11
申请人 西北工业大学 发明人 杨益新;汪勇;马远良;何正耀
分类号 G10K11/34(2006.01)I 主分类号 G10K11/34(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 王鲜凯
主权项 一种圆环形阵列恒定束宽波束形成器设计方法,其特征在于步骤如下:步骤1:确定期望波束为水平范围内的二维波束,由圆环形阵列在参考频率理论计算得到,其表达式为:<img file="FDA0001022581950000011.GIF" wi="765" he="133" />其中M为阵元个数,(·)<sup>*</sup>表示求共轭,p<sub>n</sub>(k<sub>r</sub>a<sub>r</sub>,φ)为第n个阵元接收到的从方向φ入射的单位幅度平面波信号,ω<sub>m</sub>(k<sub>r</sub>a<sub>r</sub>)为给定的期望权值向量元素,且需满足<img file="FDA0001022581950000012.GIF" wi="395" he="63" />m=1,2,…,M/2‑1;所述<img file="FDA0001022581950000013.GIF" wi="534" he="71" />式中:<img file="FDA0001022581950000014.GIF" wi="198" he="62" />a<sub>r</sub>为参考圆环形阵列的半径,k<sub>r</sub>=2π/λ<sub>r</sub>,λ<sub>r</sub>表示参考入射平面波的波长,φ<sub>n</sub>=nβ,β=2π/M;所述<img file="FDA0001022581950000015.GIF" wi="813" he="135" />步骤2:计算期望波束与合成波束的最小均方误差:<img file="FDA0001022581950000016.GIF" wi="998" he="135" />其中:<img file="FDA0001022581950000017.GIF" wi="52" he="70" />和λ<sub>m</sub>分别为循环矩阵ρ在参考频率k<sub>r</sub>a<sub>r</sub>和所要计算频率ka时的特征值,<img file="FDA0001022581950000018.GIF" wi="57" he="71" />为循环矩阵<img file="FDA0001022581950000019.GIF" wi="43" he="54" />的特征值,满足关系<img file="FDA00010225819500000110.GIF" wi="230" he="78" />λ<sub>m</sub>=λ<sub>M‑m</sub>和<img file="FDA00010225819500000111.GIF" wi="230" he="71" />所述<img file="FDA00010225819500000112.GIF" wi="574" he="127" />ρ<sub>s</sub>(k<sub>r</sub>a<sub>r</sub>)=J<sub>0</sub>(k<sub>r</sub>·Δr<sub>sr</sub>);所述<img file="FDA00010225819500000113.GIF" wi="502" he="135" />ρ<sub>s</sub>(ka)=J<sub>0</sub>(k·Δr<sub>s</sub>);所述<img file="FDA00010225819500000114.GIF" wi="1324" he="127" />式中:J<sub>0</sub>(·)是第0阶柱形Bessel函数,Δr<sub>sr</sub>=2a<sub>r</sub> sin(sβ/2)和Δr<sub>s</sub>=2a sin(sβ/2)分别是参考圆环阵和所要研究的圆环阵各自第m和m'个传感器之间的距离,s=|m‑m′|;步骤3:计算恒定束宽波束形成器在确定的频带范围内的权值向量,权值向量的第m个元素表达式为:<maths num="0001"><math><![CDATA[<mrow><msub><mi>&omega;</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>k</mi><mi>a</mi><mo>)</mo></mrow><mo>=</mo><mfrac><msub><mover><mi>&lambda;</mi><mo>&OverBar;</mo></mover><mi>m</mi></msub><msub><mi>&lambda;</mi><mi>m</mi></msub></mfrac><msub><mi>&omega;</mi><mi>m</mi></msub><mrow><mo>(</mo><msub><mi>k</mi><mi>r</mi></msub><msub><mi>a</mi><mi>r</mi></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA00010225819500000115.GIF" wi="443" he="142" /></maths>构成权值向量为ω=[ω<sub>0</sub>,ω<sub>1</sub>,…,ω<sub>M‑1</sub>]<sup>T</sup>,上标T表示转置;步骤4:根据权值向量的不同形式,采用以下两种方法合成最终波束:1、直接计算:<img file="FDA0001022581950000021.GIF" wi="805" he="132" />2、阵元域计算:<img file="FDA0001022581950000022.GIF" wi="797" he="134" />其中w=Vω,V=[v<sub>0</sub>,v<sub>1</sub>,…,v<sub>M‑1</sub>],P(ka,φ)=[p<sub>0</sub>(ka,φ),p<sub>1</sub>(ka,φ),…,p<sub>M‑1</sub>(ka,φ)]<sup>T</sup>为阵列流形向量,v<sub>m</sub>=M<sup>‑1/2</sup>[1 e<sup>imβ</sup> … e<sup>i(M‑1)mβ</sup>]<sup>T</sup>为循环矩阵特征向量且满足关系<img file="FDA0001022581950000023.GIF" wi="230" he="69" />ka在波束的频率范围内取值,该频率范围根据设计精度,由步骤2得到的最小均方误差曲线δ(ka)确定。
地址 710072 陕西省西安市友谊西路127号