发明名称 利用压缩感知的单端时域波束搜索方法
摘要 本发明属于无线通信技术领域,具体涉及在多天线波束成形通信系统中的采用压缩感知来搜索最优波束矢量的方法。本发明提供了一种在多天线波束成形通信系统中的一种利用压缩感知的单端时域波束来搜索最优波束矢量的方法。该方法利用离开角、到达角的稀疏性将波束搜索的问题转化为压缩感知的问题,通过发射端和接收端使用不同的发射和接收矢量,由接收端单独确定最优的发射/接收波束矢量。本发明使用范围极广,可用于所有的慢衰落视距或者非视距信道。
申请公布号 CN104168047A 申请公布日期 2014.11.26
申请号 CN201410396381.X 申请日期 2014.08.13
申请人 电子科技大学 发明人 王梦瑶;成先涛
分类号 H04B7/04(2006.01)I 主分类号 H04B7/04(2006.01)I
代理机构 成都宏顺专利代理事务所(普通合伙) 51227 代理人 李玉兴
主权项 利用压缩感知的单端时域波束搜索方法,其特征在于,包括以下步骤:S1、令设备1的收发天线数为Nt,所述设备1的码本中的波束数目为Ct,所述设备1采用P<sub>t</sub>种发射矢量进行发射,任意一个发射矢量<img file="FDA0000553086370000011.GIF" wi="80" he="89" />都是长度为Nt的向量,所述发射矢量中每个位置的元素的值从集合[1,i,‑1,‑i]中随机选择,组成一个测量矩阵<img file="FDA0000553086370000012.GIF" wi="456" he="92" />所述测量矩阵Φ<sub>t</sub>每一行都对应一次发射,所述设备1在时域发射时间序列为[1,0,...,0],所述时间序列长度为N,其中,d=1,2,...P<sub>t</sub>,i为虚数单位,令设备2的收发天线数为Nr,所述设备2的码本中的波束数目为Cr,对于设备1的每个发射矢量设备2都有P<sub>r</sub>个接收矢量来接收,任意一个接收矢量<img file="FDA0000553086370000013.GIF" wi="90" he="88" />都是长度为Nr的向量,所述接收矢量中每个位置的元素的值从集合[1,i,‑1,‑i]中随机选择,组成一个测量矩阵<img file="FDA0000553086370000014.GIF" wi="473" he="92" />所述测量矩阵Φ<sub>r</sub>每一行都对应一次接收,接收端能够得到第n个时间点测量信号矩阵为Y<sub>n</sub>=Φ<sub>r</sub>h<sub>n</sub>Φ<sup>T</sup><sub>t</sub>+n<sub>n</sub>,其中,d′=1,2,...P<sub>r</sub>,n=1,2,...,N,矩阵Y<sub>n</sub>的阶数为P<sub>r</sub>×P<sub>t</sub>,n<sub>n</sub>是噪声矩阵,h<sub>n</sub>为第n个时间点的阶数为Nr×Nt的信道矩阵,矩阵中第x行第y列的元素表示从发射端第y根天线到接收端第x根天线间的频域信道冲击响应,x=1,2,...,Nr,y=1,2,...,Nt,()<sup>T</sup>是矩阵的转置运算,N、Nt、Nr、Ct、Cr、P<sub>r</sub>和P<sub>t</sub>为大于1的整数;S2、根据S1所述构建字典矩阵为D,D的每一列对应[‑90°,90°]中的一个角度;S3、根据S1所述Φ<sub>r</sub>和Φ<sub>t</sub>恢复出h<sub>n</sub>,即已知Y<sub>n</sub>、Φ<sub>r</sub>和Φ<sub>t</sub>,算出矩阵h<sub>n</sub>,具体为:S31、根据每个时间点信号Y<sub>n</sub>恢复出Y′<sub>n</sub>,所述Y′<sub>n</sub>的每一列都可以在S2所述字典矩阵D下展开,即Y′<sub>n</sub>的每一列都可以表示为字典矩阵中少数列与对应不为0的展开系数相乘后的线性加和,展开系数为复数,具体如下:对每个时间点信号Y<sub>n</sub>使用P<sub>t</sub>个任务的正交匹配追踪算法(Orthogonal Matching Pursuit,OMP)联合所述时间点信号Y<sub>n</sub>的每一列共同恢复出Y′<sub>n</sub>,所述时间点信号Y<sub>n</sub>的任意第l列为<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mover><mi>y</mi><mo>&RightArrow;</mo></mover><mrow><mi>n</mi><mo>,</mo><mi>l</mi></mrow></msub><mo>=</mo><msub><mi>&Phi;</mi><mi>r</mi></msub><mi>D</mi><msub><mover><mi>&omega;</mi><mo>&RightArrow;</mo></mover><mrow><mi>n</mi><mo>,</mo><mi>l</mi></mrow></msub><mo>+</mo><msub><mover><mi>n</mi><mo>&RightArrow;</mo></mover><mrow><mi>n</mi><mo>,</mo><mi>l</mi></mrow></msub><mo>=</mo><msub><mi>V</mi><mi>r</mi></msub><msub><mover><mi>&omega;</mi><mo>&RightArrow;</mo></mover><mrow><mi>n</mi><mo>,</mo><mi>l</mi></mrow></msub><mo>+</mo><msub><mover><mi>n</mi><mo>&RightArrow;</mo></mover><mrow><mi>n</mi><mo>,</mo><mi>l</mi></mrow></msub><mo>,</mo></mrow>]]></math><img file="FDA0000553086370000015.GIF" wi="697" he="77" /></maths>其中,V<sub>r</sub>=Φ<sub>r</sub>D,<img file="FDA0000553086370000016.GIF" wi="70" he="71" />是nn的第l列,<img file="FDA0000553086370000017.GIF" wi="75" he="75" />可以在V<sub>r</sub>下展开,<img file="FDA0000553086370000021.GIF" wi="81" he="76" />就是<img file="FDA0000553086370000022.GIF" wi="66" he="75" />在V<sub>r</sub>下的展开系数,l=1,2,...,P<sub>t</sub>;S32、对每个时间点n使用Nr个任务的OMP联合S31所述Y′<sub>m</sub>的每一列共同恢复出h<sub>n</sub><sup>T</sup>,所述h<sub>n</sub><sup>T</sup>的每一列都可以在S2所述字典矩阵D下展开,其中,Y′<sub>n</sub>=h<sub>n</sub>Φ<sup>T</sup><sub>t</sub>;S33、根据S32所述h<sub>n</sub><sup>T</sup>恢复出h<sub>n</sub>;S4、恢复出所有时间点的的h<sub>n</sub>,做N点的离散傅里叶变换得到每个频点的频域信道矩阵<img file="FDA0000553086370000023.GIF" wi="130" he="84" />从码本中找到一个最优的<img file="FDA0000553086370000024.GIF" wi="45" he="57" />和<img file="FDA0000553086370000025.GIF" wi="51" he="57" />使得频谱效率最大,即<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mrow><mo>(</mo><mover><mi>c</mi><mo>&RightArrow;</mo></mover><mo>,</mo><mover><mi>w</mi><mo>&RightArrow;</mo></mover><mo>)</mo></mrow><mo>=</mo><munder><mi>arg</mi><mrow><mover><mi>c</mi><mo>&RightArrow;</mo></mover><mo>,</mo><mover><mi>w</mi><mo>&RightArrow;</mo></mover></mrow></munder><mi>max</mi><mrow><mo>(</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>log</mi><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msub><mi>&gamma;</mi><mrow><mi>r</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>)</mo></mrow><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000553086370000026.GIF" wi="776" he="154" /></maths>其中<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>&gamma;</mi><mrow><mi>r</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>=</mo><mfrac><msup><mrow><mo>|</mo><msup><mover><mi>w</mi><mo>&RightArrow;</mo></mover><mi>T</mi></msup><msub><mover><mi>H</mi><mo>^</mo></mover><mi>n</mi></msub><mi></mi><mover><mi>c</mi><mo>&RightArrow;</mo></mover><mo>|</mo></mrow><mn>2</mn></msup><mrow><msub><mi>N</mi><mi>r</mi></msub><msub><mi>N</mi><mi>t</mi></msub><msup><mi>&sigma;</mi><mn>2</mn></msup></mrow></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000553086370000027.GIF" wi="360" he="173" /></maths>σ<sup>2</sup>是噪声的功率,<img file="FDA0000553086370000028.GIF" wi="53" he="60" />和<img file="FDA0000553086370000029.GIF" wi="46" he="55" />是长度为Nr的复向量。
地址 611731 四川省成都市高新区(西区)西源大道2006号
您可能感兴趣的专利