发明名称 一种大动态期望信号下的波束合成方法
摘要 一种大动态期望信号下的波束合成方法,其实现步骤为:1)阵列第k时刻接收的信号为x(k);2)接收实信号进行正交变换得解析信号<img file="DDA00003294071800011.GIF" wi="137" he="82" />3)设权矢量初值为参考角度的方向矢量v<sub>0</sub>,按照公式<img file="DDA00003294071800012.GIF" wi="888" he="106" />进行迭代计算,收敛后的权矢量为w<sub>f</sub>(K);4)根据得到的权值矢量w<sub>f</sub>(K)计算方向图<img file="DDA00003294071800013.GIF" wi="444" he="88" />5)方向图上参考角度<img file="DDA00003294071800014.GIF" wi="60" he="93" />对应的导数<img file="DDA00003294071800015.GIF" wi="614" he="157" />6)判断:若α的绝对值大于门限α<sub>0</sub>,选择相控阵的权矢量w<sub>q</sub>作为波束合成权矢量,即w<sub>e</sub>=w<sub>q</sub>;否则,使用权矢量w<sub>f</sub>(K)作为波束合成权矢量,即w<sub>e</sub>=w<sub>f</sub>(K);7)计算波束合成的输出<img file="DDA00003294071800016.GIF" wi="351" he="81" />本发明不仅克服了“自零陷”问题,相比于各种复杂稳健波束合成算法,实现简单,满足工程应用的需求。
申请公布号 CN103245934A 申请公布日期 2013.08.14
申请号 CN201310216613.4 申请日期 2013.06.03
申请人 重庆大学 发明人 曾浩;凤林锋;刘陆军;幸倩;刘玲
分类号 G01S7/292(2006.01)I;G01S7/527(2006.01)I;H04B7/08(2006.01)I 主分类号 G01S7/292(2006.01)I
代理机构 代理人
主权项 1.一种大动态期望信号下的波束合成方法,其实现步骤是:第一步,对N个阵元构成的均匀线阵,阵元间距d为半个波长λ/2,第k时刻阵列的接收信号可以表示为<maths num="0001"><![CDATA[<math><mrow><mi>x</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>x</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><msub><mi>x</mi><mi>N</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths>N为自然数,k=1,......;第二步,对第k时刻接收到的信号x(k)进行正交变换,将实信号变成相应的解析信号<maths num="0002"><![CDATA[<math><mrow><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mover><mi>x</mi><mo>^</mo></mover><mn>1</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><msub><mover><mi>x</mi><mo>^</mo></mover><mi>N</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>;</mo></mrow></math>]]></maths>第三步,根据解析信号<img file="FDA00003294071500013.GIF" wi="94" he="87" />和已知期望信号波达方向的参考值<img file="FDA00003294071500014.GIF" wi="91" he="90" />按Frost算法计算权矢量w<sub>f</sub>:(1)计算参考角度<img file="FDA00003294071500015.GIF" wi="62" he="93" />的方向矢量<maths num="0003"><![CDATA[<math><mrow><msub><mi>v</mi><mn>0</mn></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd></mtr><mtr><mtd><msup><mi>e</mi><mrow><mi>j&pi;</mi><mi>sin</mi><msub><mover><mi>&theta;</mi><mo>^</mo></mover><mn>0</mn></msub></mrow></msup></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><msup><mi>e</mi><mrow><mi>j&pi;</mi><mrow><mo>(</mo><mi>N</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>sin</mi><msub><mover><mi>&theta;</mi><mo>^</mo></mover><mn>0</mn></msub></mrow></msup></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths>并确定步长μ;(2)计算N×N维矩阵P=I-v<sub>0</sub>(v<sub>0</sub><sup>H</sup>v<sub>0</sub>)<sup>-1</sup>v<sub>0</sub><sup>H</sup>,I是N×N维单位矩阵;(3)计算N维列向量<img file="FDA00003294071500017.GIF" wi="205" he="118" />(4)取变量m=1,权矢量初值w<sub>f</sub>(1)=w<sub>q</sub>;(5)计算第m次输出<img file="FDA00003294071500018.GIF" wi="424" he="92" />符号H表示共轭转置;(6)更新权矢量<maths num="0004"><![CDATA[<math><mrow><msub><mi>w</mi><mi>f</mi></msub><mrow><mo>(</mo><mi>m</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><mi>P</mi><mrow><mo>(</mo><msub><mi>w</mi><mi>f</mi></msub><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>-</mo><mi>&mu;</mi><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><msup><mi>y</mi><mo>*</mo></msup><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>+</mo><msub><mi>w</mi><mi>q</mi></msub><mo>;</mo></mrow></math>]]></maths>(7)令m=m+1,重复执行上述(5)—(7)步,直到权矢量w<sub>f</sub>收敛为止,此时m=K;第四步,计算此时阵列方向图为<img file="FDA000032940715000110.GIF" wi="419" he="90" />第五步,设置角度偏移常数△θ,估计方向图上参考角度<img file="FDA000032940715000111.GIF" wi="60" he="88" />对应的导数<maths num="0005"><![CDATA[<math><mrow><mi>&alpha;</mi><mo>=</mo><mfrac><mrow><mi>B</mi><mrow><mo>(</mo><msub><mover><mi>&theta;</mi><mo>^</mo></mover><mn>0</mn></msub><mo>+</mo><mi>&Delta;&theta;</mi><mo>)</mo></mrow><mo>-</mo><mi>B</mi><mrow><mo>(</mo><msub><mover><mi>&theta;</mi><mo>^</mo></mover><mn>0</mn></msub><mo>-</mo><mi>&Delta;&theta;</mi><mo>)</mo></mrow></mrow><mrow><mn>2</mn><mi>&Delta;&theta;</mi></mrow></mfrac><mo>;</mo></mrow></math>]]></maths>第六步,设置门限值α<sub>0</sub>,若α的绝对值大于门限α<sub>0</sub>,则表示在期望方向附近的方向图有零陷,选择矢量w<sub>q</sub>作为最后的波束合成权矢量,即w<sub>e</sub>=w<sub>q</sub>;否则,使用Frost算法得到的权矢量w<sub>f</sub>(K)作为最后的波束合成权矢量,即w<sub>e</sub>=w<sub>f</sub>(K);第七步,阵列波束合成输出<img file="FDA00003294071500021.GIF" wi="338" he="86" />此时k=K+1,K+2,......。
地址 400044 重庆市沙坪坝区沙正街174号