CN104219188A - 利用压缩感知的双端时域波束搜索方法 - Google Patents

Abstract

本发明属于无线通信技术领域，具体涉及在多天线正交频分复用(Orthogonal Frequency Division Multiplexing，OFDM)通信系统中的采用压缩感知来搜索最优波束矢量的方法。本发明提供了一种在多天线OFDM通信系统中的一种利用压缩感知的双端时域波束来搜索最优波束矢量的方法。该方法利用离开角、到达角的稀疏性将波束搜索的问题转化为压缩感知的问题，结合信道的对称性，通过反复迭代的方法确定最优的发射和接收波束矢量。

Description

利用压缩感知的双端时域波束搜索方法

技术领域

本发明属于无线通信技术领域，具体涉及在多天线正交频分复用(OrthogonalFrequency Division Multiplexing，OFDM)通信系统中的采用压缩感知来搜索最优波束矢量的方法。

背景技术

UWB系统和60GHz系统主要用于短距离高速传输，应用范围广泛，包括无线个域网(WPAN，Wireless Personal Area Network)，无线高清多媒体接口，医疗成像，车载雷达等等。为了适应高数据率和高系统容量等方面的需要，UWB系统和60GHz系统往往利用多天线多载波技术用于传输数据。

多天线技术包括多输入多输出(Multiple Input Multiple Output，MIMO)，多输入单输出(Multiple Input Single Output，MISO)和单输入多输出(Single Input MultipleOutput，SIMO)。基于阵列天线的波束成形技术利用传输信号的方向性提高信噪比(Signal to Noise Ratio，SNR)，抑制干扰，改善系统性能。

阵列天线在空间的分布情况影响了信道空间的相关性，智能天线中的波束成形技术利用了这种相关性对信号进行处理，在期望方向上产生方向性强的辐射波束增强有用信号，零瓣方向对准干扰源达到抑制作用，由此提高信噪比和增加传输距离。在收/发端应用天线阵列波束成形具有以下优势：首先，降低对功率放大器的要求。发射端如果使用单个天线时，对PA增益要求很高。如果发射端使用天线阵列发送信号，每个天线阵元前面增加一个功放，这样通过使用多个较低功率增益的PA就能够满足发射功率要求。其次，天线阵列波束成形便于定向传输。在发射功率不变情况下，等效增加接收机接收信号的功率，同时还可以有效降低多径时延扩展。这样可以简化收发机的基带设计，降低模拟数字转换器的分辨率指标。最后，天线阵列系统动态地调整波束的方向，以使期望方向获得最大的功率并减小其他方向的功率。这样不仅改善了信号干扰比，还提高了系统的容量，扩大了系统通信覆盖范围，降低了发射功率要求。

OFDM是多载波调制的一种。其主要思想是：将信道分成若干正交子信道，将高速数据信号转换成并行的低速子数据流，调制到在每个子信道上进行传输。正交信号可以通过在接收端采用相关技术来分开，这样可以减少子信道间相互干扰ISI。每个子信道上的信号带宽小于信道的相干带宽，因此每个子信道上的可以看成平坦性衰落，从而可以消除符号间干扰。而且由于每个子信道的带宽仅仅是原信道带宽的一小部分，信道均衡变得相对容易。基于OFDM的波束成形需要在发射端天线前做快速傅里叶逆变换变换，接收端做快速傅里叶变换变换来解调。

波束切换是一种波束搜索规则，它在发射机和接收机两端都预先设置好波束控制矢量码本，使用时只需要从中选取。因此，切换波束形成也称为基于码本的波束成形，使用开关天线阵列，在发送数据包前，发射机要多次发送携带不同波束控制矢量的信息。

基于信道状态信息的波束成形技术，发射机和接收机都可以找到一个最优的波束成形控制矢量。其详细方法可参考：Yoon S,Jeon T,Lee W.Hybrid beam-formingand beam-switching for OFDM based wireless personal area networks[J].SelectedAreas in Communications,IEEE Journal on,2009,27(8):1425-1432.物理层(PHY)解决方案能够提供最优的系统性能，波束成形操作往往考虑在物理层进行，但获取完整的信道状态信息要很高的时间成本和开销。基于码本的波束成形技术有助于降低复杂度和开销，而且码本既可以完全根据基带信号处理而设计，也可以结合控制层(MAC)实现。

波束搜索时的搜索策略是至关重要的，高效的波束搜索策略能够有效降低搜索时间，假设发射端有N个发射波束矢量，M个接收波束矢量，则最多需要N×M次搜索，802.15.3c中采用了两级的码本结构：一个扇形码本和一个波束码本，波束码本的每个列向量表示一个波束，每个波束图案都表示一个精确的方向，每个扇区都是几个波束的集合，在空间中表示较宽的方向，所有的扇区加起来覆盖整个空间。搜索过程也分为两阶段：第一阶段在根据信噪比找到最优的扇区，第二阶段在最优的扇区中找到最优的波束。其详细方法可参考：Wang J,Lan Z,Pyo C W,et al.Beamcodebook based beamforming protocol for multi-Gbps millimeter-wave WPANsystems[J].Selected Areas in Communications,IEEE Journal on,2009,27(8):1390-1399.。

分阶段的波束搜索策略可以大幅减低搜索次数，但是当天线阵列很大时，需要的搜索次数仍然是巨大的。因此，研究一种快速有效的波束搜索算法是一项有创新性和重要实际意义且具挑战性的任务。

发明内容

本发明提供了一种在多天线OFDM通信系统中的一种利用压缩感知的双端频域波束来搜索最优波束矢量的方法。该方法利用发射角、到达角的稀疏性将波束搜索的问题转化为压缩感知的问题，通过发射端和接收端使用不同的发射和接收矢量，由接收端单独确定最优的发射/接收波束矢量。

本发明的目的是通过如下步骤来实现的：

S1、令设备1的收发天线数为Nt，所述设备1的码本中的波束数目为Ct，所述设备1使用全向天线向设备2发射，发射波束矢量为 ${\stackrel{\→}{w}}_1=\left[\begin{array}{c}1\\ 0\\ .\\ .\\ .\\ 0\end{array}\right],$ 所述发射波束长度为Nt，所述设备1使用OFDM技术在频域发射的序列为[1,1,...,1]，所述发射的序列长度为N

令设备2的收发天线数为Nr，所述设备2的码本中的波束数目为Cr，所述设备2的天线接收到的第m个子载波信号向量为所述设备2的接收端使用P_r个接收矢量来接收信号，任意一个接收矢量都是长度为Nr的向量，所述接收矢量中的每个元素的值从集合[1,i,-1,-i]中随机选择，组成一个测量矩阵所述测量矩阵Φ_r每一行都对应一次接收，子载波信号为测量信号向量为其中，是长度为Nr的噪声向量，H_m为第m个频点的阶数为Nr×Nt信道矩阵，所述信道矩阵中第x行第y列的元素表示从发射端第y根天线到接收端第x根天线间的频域信道冲击响应，其中，m＝1,2,...,N，y＝1,2,...,Nt，x＝1,2,...,Nr，d＝1,2,...P_r，i为虚数单位，是噪声向量，中的每个元素对应一个测量值，()^T是矩阵的转置运算，P_r为大于1的整数，N、Nt、Nr、Ct和Cr均为大于1的整数；

S2、根据S1所述构建字典矩阵为D，D的每一列对应[-90°，90°]中的一个角度，S1所述信号可以在D下展开，并且是稀疏的，展开系数为复数，是在D下的展开系数；

S3、、使用N个任务的正交匹配追踪算法(Orthogonal Matching Pursuit，OMP)联合所有S1所述子载波信号共同恢复出所有的具体为：

S31、 ${\stackrel{\→}{y}}_{m,r}={\mathrm{\Φ}}_r{\stackrel{\→}{g}}_{m,r}={\mathrm{\Φ}}_{r}D{\stackrel{\→}{\mathrm{\ω}}}_{m,r}+{\stackrel{\→}{n}}_{m,r}^{\′}=V_r{\stackrel{\→}{\mathrm{\ω}}}_{m,r}+{\stackrel{\→}{n}}_{m,r}^{\′},$ V_r＝Φ_rD，所述可以在V_r下展开，是在V_r下的展开系数；

S32、从S31所述V_r中找到一列使得最大，构造矩阵算出所有在V_c下的展开系数 $[{\stackrel{\→}{\mathrm{\ω}}}_{1,r,c},{\stackrel{\→}{\mathrm{\ω}}}_{2,r,c},...,{\stackrel{\→}{\mathrm{\ω}}}_{N,r,c}]={\left({V_c}^H{V}_c\right)}^{-1}{V_c}^H[{\stackrel{\→}{y}}_{1,r},{\stackrel{\→}{y}}_{2,r},...,{\stackrel{\→}{y}}_{N,r}],$ 表示当前恢复程度的剩余量 $e_r=[{\stackrel{\→}{y}}_{1,r},{\stackrel{\→}{y}}_{2,r},...,{\stackrel{\→}{y}}_{N,r}]-V_c[{\stackrel{\→}{\mathrm{\ω}}}_{1,r,c},{\stackrel{\→}{\mathrm{\ω}}}_{2,r.c}...,{\stackrel{\→}{\mathrm{\ω}}}_{N,r,c}],$ 其中，()^-1是矩阵的求逆运算，()^H是矩阵的共轭转置运算，|·|表示取复数的幅度，||·||₂表示向量的二范数运算；

S33、从V_r中找到使得最大，其中是矩阵e_r中的第m列，将加到S32所述V_c中得到更新后的算出在更新后的V_c下的展开系数；

S34、循环S34到S33，直到e_r的F范数小于的F范数的a倍时停止，结合S2所述字典矩阵D中的列向量与对应位置的系数线性组合恢复出其中，a为门限值，0<a<1,且a为实数；

S4、恢复出所有频点的记作从码本中找到一个适合的波束矢量使得频谱效率最大， ${\stackrel{\&RightArrow;}{c}}_1=\underset{\stackrel{\&RightArrow;}{c}}{\arg }\max \left(\frac{1}{N}\underset{m=1}{\overset{N}{\mathrm{\&Sigma;}}}{\log }_2(1+{\mathrm{\&gamma;}}_{r,m})\right),$ 其中 ${\mathrm{\&gamma;}}_{r,m}=\frac{{\left|{\stackrel{\&RightArrow;}{c}}^T{\widehat{\stackrel{\&RightArrow;}{g}}}_{r,m}\right|}^2}{N_r{N}_t{\mathrm{\&sigma;}}^2},$ σ²是噪声的功率，是长度为Nr的复向量；

S5、使用N个任务的OMP联合所有S1所述测量信号共同恢复出所有的具体如下：

S51、 ${\stackrel{\&RightArrow;}{y}}_{m,t}={\mathrm{\&Phi;}}_t{\stackrel{\&RightArrow;}{g}}_{m,t}={\mathrm{\&Phi;}}_{t}D{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{m,t}+{\stackrel{\&RightArrow;}{n}}_{m,t}^{\&prime;}=V_t{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{m,t}+{\stackrel{\&RightArrow;}{n}}_{m,t}^{\&prime;},$ 所述可以在可以在V_t下展开，是在V_t下的展开系数，其中，V_t＝Φ_tD；

S52、从S51所述V_t中找到一列使得最大，构造矩阵算出所有在V_t下的展开系数 $[{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{1,r,c},{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{2,r,c},...,{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{N,r,c}]={\left({V_c}^H{V}_c\right)}^{-1}{V_c}^H[{\stackrel{\&RightArrow;}{y}}_{1,r},{\stackrel{\&RightArrow;}{y}}_{2,r},...,{\stackrel{\&RightArrow;}{y}}_{N,r}],$ 表示当前恢复程度的剩余量 $e_r=[{\stackrel{\&RightArrow;}{y}}_{1,r},{\stackrel{\&RightArrow;}{y}}_{2,r},...,{\stackrel{\&RightArrow;}{y}}_{N,r}]-V_c[{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{1,r,c},{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{2,r.c}...,{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{N,r,c}],$ 其中，()^-1是矩阵的求逆运算，()^H是矩阵的共轭转置运算，|·|表示取复数的幅度，||·||₂表示向量的二范数运算；

S53、从V_r中找到使得最大，其中是矩阵e_r中的第m列，将加到S32所述V_c中得到更新后的算出在更新后的V_c下的展开系数；

S54、循环S51到S53，直到e_r的F范数小于的F范数的a倍时停止，结合S2所述字典矩阵D中的列向量与对应位置的系数线性组合恢复出记为从码本中找到一个最适合的波束矢量使得频谱效率最大，其中σ²是噪声的功率，是长度为Nr的复向量其中，a为门限值，0<a<1,且a为实数，σ²是噪声的功率，是长度为Nr的复向量；

S6、设备2的接收端使用P_t个接收矢量来接收信号，与S1.3的步骤一样，结合所有频点的信号运用多任务正交匹配追踪算法，测量信号为 ${\stackrel{\&RightArrow;}{y}}_{m,t}={\mathrm{\&Phi;}}_t{\stackrel{\&RightArrow;}{g}}_{m,t}={\mathrm{\&Phi;}}_{t}D{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{m,t}+{\stackrel{\&RightArrow;}{n}}_{m,t}^{\&prime;}=V_t{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{m,t}+{\stackrel{\&RightArrow;}{n}}_{m,t}^{\&prime;},$ 其中V_t＝Φ_tD， ${\mathrm{\&Phi;}}_t={\left[\begin{array}{ccc}{\stackrel{\&RightArrow;}{\mathrm{\&phi;}}}_{t,1}& ...& {\stackrel{\&RightArrow;}{\mathrm{\&phi;}}}_{t,P_t}\end{array}\right]}^T$ 是测量矩阵，每一行接收矢量对应一次测量，任意一个接收矢量(d＝1,2,...,P_r)都是长度为Nt的向量，每个元素的值从集合[1,i,-1,-i]中随机选择，是在V_t下的展开系数，重复测量P_t次，根据恢复出记为当所有频点的(m＝1,2,...,N)都恢复出来后，从码本中找到一个最适合的使得频谱效率最大，其中P_t为大于1的整数；是长度为Nt的复向量；

S7、设备1以作为发射波束矢量向设备2发射，设备2通过重复S1-S6的步骤找到最优的接收矢量

S8、经过反复迭代，对于设备1和设备2来说，当相邻两次找到的波束矢量，即和相同时迭代终止，并将最终找到的和作为设备1和设备2最优的波束矢量。

进一步地，对于任意角度q，S2所述字典矩阵D中的对应列为 $\left[\begin{array}{c}1\\ e^{\mathrm{i\&pi;}\sin \left(\mathrm{\&theta;}\right)}\\ e^{i2\mathrm{\&pi;}\sin \left(\mathrm{\&theta;}\right)}\\ .\\ .\\ .\\ e^{i(\mathrm{Nt}-1)\mathrm{\&pi;}\sin \left(\mathrm{\&theta;}\right)}\end{array}\right].$

进一步地，S54所述a＝0.05。

本发明的有益效果是：波束搜索所需次数与总的路径数有关，搜索复杂度不会随着天线数目而增加。本发明适用范围极广，可用于所有的慢衰落视距或者非视距信道。

附图说明

图1是本发明利用压缩感知的双端时域波束搜索算法的结构图。

图2是本发明用于802.11.ad信道波束搜索的的成功概率性能曲线图。

具体实施方式

下面结合实施例和附图，详细说明本发明的技术方案。

如图1所示，本发明利用压缩感知的双端时域波束搜索算法的结构图。整个过程在频域完成，设备1以全向天线向设备2发射，设备2重复接收P_r次，每次用不同的接收矢量，设备2根据P_r个测量值使用压缩感知对原始接收信号还原，关于信号的处理过程都是在频域完成，根据还原出的原始信号从码本中找到一个最优的接收矢量使得频谱效率最大。由于对称性，最优的接收矢量就是最优的发射矢量，随后，设备2以找到的最优接收矢量作为发射矢量向设备1发射，同样的，设备1重复接收P_t次，每次使用不同的接收矢量，设备1使用压缩感知对信号还原，根据还原出的信号从码本中找到最优的接收矢量使得频谱效率最大。这样的过程反复进行，每个设备在计算出此次最优的接收矢量后都与上次计算出的最优接收矢量做对比，当相邻两次计算出的最优接收矢量相同时迭代就可以终止。

S1、令设备1的收发天线数为Nt，所述设备1的码本中的波束数目为Ct，所述设备1使用全向天线向设备2发射，发射波束矢量为 ${\stackrel{\&RightArrow;}{w}}_1=\left[\begin{array}{c}1\\ 0\\ .\\ .\\ .\\ 0\end{array}\right],$ 所述发射波束长度为Nt，所述设备1使用OFDM技术在频域发射的序列为[1,1,...,1]，所述发射的序列长度为N

令设备2的收发天线数为Nr，所述设备2的码本中的波束数目为Cr，所述设备2的天线接收到的第m个子载波信号向量为所述设备2的接收端使用P_r个接收矢量来接收信号，任意一个接收矢量都是长度为Nr的向量，所述接收矢量中的每个元素的值从集合[1,i,-1,-i]中随机选择，组成一个测量矩阵所述测量矩阵Φ_r每一行都对应一次接收，子载波信号为测量信号向量为其中，是长度为Nr的噪声向量，H_m为第m个频点的阶数为Nr×Nt信道矩阵，所述信道矩阵中第x行第y列的元素表示从发射端第y根天线到接收端第x根天线间的频域信道冲击响应，其中，m＝1,2,...,N，y＝1,2,...,Nt，x＝1,2,...,Nr，d＝1,2,...P_r，i为虚数单位，是噪声向量，中的每个元素对应一个测量值，()^T是矩阵的转置运算，P_r为大于1的整数，N、Nt、Nr、Ct和Cr均为大于1的整数；

S2、根据S1所述构建字典矩阵为D，D的每一列对应[-90°，90°]中的一个角度，S1所述信号可以在D下展开，并且是稀疏的，展开系数为复数，是在D下的展开系数；

S3、、使用N个任务的正交匹配追踪算法(Orthogonal Matching Pursuit，OMP)联合所有S1所述子载波信号共同恢复出所有的具体为：

S31、 ${\stackrel{\&RightArrow;}{y}}_{m,r}={\mathrm{\&Phi;}}_r{\stackrel{\&RightArrow;}{g}}_{m,r}={\mathrm{\&Phi;}}_{r}D{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{m,r}+{\stackrel{\&RightArrow;}{n}}_{m,r}^{\&prime;}=V_r{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{m,r}+{\stackrel{\&RightArrow;}{n}}_{m,r}^{\&prime;},$ V_r＝Φ_rD，所述可以在V_r下展开，是在V_r下的展开系数；

S32、从S31所述V_r中找到一列使得最大，构造矩阵算出所有在V_c下的展开系数 $[{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{1,r,c},{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{2,r,c},...,{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{N,r,c}]={\left({V_c}^H{V}_c\right)}^{-1}{V_c}^H[{\stackrel{\&RightArrow;}{y}}_{1,r},{\stackrel{\&RightArrow;}{y}}_{2,r},...,{\stackrel{\&RightArrow;}{y}}_{N,r}],$ 表示当前恢复程度的剩余量 $e_r=[{\stackrel{\&RightArrow;}{y}}_{1,r},{\stackrel{\&RightArrow;}{y}}_{2,r},...,{\stackrel{\&RightArrow;}{y}}_{N,r}]-V_c[{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{1,r,c},{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{2,r.c}...,{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{N,r,c}],$ 其中，()^-1是矩阵的求逆运算，()^H是矩阵的共轭转置运算，|·|表示取复数的幅度，||·||₂表示向量的二范数运算；

S33、从V_r中找到使得最大，其中是矩阵e_r中的第m列，将加到S32所述V_c中得到更新后的算出在更新后的Vc下的展开系数；

S34、循环S34到S33，直到e_r的F范数小于的F范数的a倍时停止，结合S2所述字典矩阵D中的列向量与对应位置的系数线性组合恢复出其中，a为门限值，0<a<1,且a为实数；

S4、恢复出所有频点的记作从码本中找到一个适合的波束矢量使得频谱效率最大， ${\stackrel{\&RightArrow;}{c}}_1=\underset{\stackrel{\&RightArrow;}{c}}{\arg }\max \left(\frac{1}{N}\underset{m=1}{\overset{N}{\mathrm{\&Sigma;}}}{\log }_2(1+{\mathrm{\&gamma;}}_{r,m})\right),$ 其中 ${\mathrm{\&gamma;}}_{r,m}=\frac{{\left|{\stackrel{\&RightArrow;}{c}}^T{\widehat{\stackrel{\&RightArrow;}{g}}}_{r,m}\right|}^2}{N_r{N}_t{\mathrm{\&sigma;}}^2},$ σ²是噪声的功率，是长度为Nr的复向量；

S5、使用N个任务的OMP联合所有S1所述测量信号共同恢复出所有的具体如下：

S51、 ${\stackrel{\&RightArrow;}{y}}_{m,t}={\mathrm{\&Phi;}}_t{\stackrel{\&RightArrow;}{g}}_{m,t}={\mathrm{\&Phi;}}_{t}D{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{m,t}+{\stackrel{\&RightArrow;}{n}}_{m,t}^{\&prime;}=V_t{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{m,t}+{\stackrel{\&RightArrow;}{n}}_{m,t}^{\&prime;},$ 所述可以在可以在V_t下展开，是在V_t下的展开系数，其中，V_t＝Φ_tD；

S52、从S51所述V_t中找到一列使得最大，构造矩阵算出所有在V_t下的展开系数 $[{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{1,r,c},{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{2,r,c},...,{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{N,r,c}]={\left({V_c}^H{V}_c\right)}^{-1}{V_c}^H[{\stackrel{\&RightArrow;}{y}}_{1,r},{\stackrel{\&RightArrow;}{y}}_{2,r},...,{\stackrel{\&RightArrow;}{y}}_{N,r}],$ 表示当前恢复程度的剩余量 $e_r=[{\stackrel{\&RightArrow;}{y}}_{1,r},{\stackrel{\&RightArrow;}{y}}_{2,r},...,{\stackrel{\&RightArrow;}{y}}_{N,r}]-V_c[{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{1,r,c},{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{2,r.c}...,{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{N,r,c}],$ 其中，()^-1是矩阵的求逆运算，()^H是矩阵的共轭转置运算，|·|表示取复数的幅度，||·||₂表示向量的二范数运算；

S53、从V_r中找到使得最大，其中是矩阵e_r中的第m列，将加到S32所述V_c中得到更新后的算出在更新后的V_c下的展开系数；

S54、循环S51到S53，直到e_r的F范数小于的F范数的a倍时停止，结合S2所述字典矩阵D中的列向量与对应位置的系数线性组合恢复出记为从码本中找到一个最适合的波束矢量使得频谱效率最大，其中σ²是噪声的功率，是长度为Nr的复向量其中，a为门限值，0<a<1,且a为实数，σ²是噪声的功率，是长度为Nr的复向量；

S6、设备2的接收端使用P_t个接收矢量来接收信号，与S1.3的步骤一样，结合所有频点的信号运用多任务正交匹配追踪算法，测量信号为 ${\stackrel{\&RightArrow;}{y}}_{m,t}={\mathrm{\&Phi;}}_t{\stackrel{\&RightArrow;}{g}}_{m,t}={\mathrm{\&Phi;}}_{t}D{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{m,t}+{\stackrel{\&RightArrow;}{n}}_{m,t}^{\&prime;}=V_t{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{m,t}+{\stackrel{\&RightArrow;}{n}}_{m,t}^{\&prime;},$ 其中V_t＝Φ_tD， ${\mathrm{\&Phi;}}_t={\left[\begin{array}{ccc}{\stackrel{\&RightArrow;}{\mathrm{\&phi;}}}_{t,1}& ...& {\stackrel{\&RightArrow;}{\mathrm{\&phi;}}}_{t,P_t}\end{array}\right]}^T$ 是测量矩阵，每一行接收矢量对应一次测量，任意一个接收矢量(d＝1,2,...,P_r)都是长度为Nt的向量，每个元素的值从集合[1,i,-1,-i]中随机选择，是在V_t下的展开系数，重复测量P_t次，根据恢复出记为当所有频点的(m＝1,2,...,N)都恢复出来后，从码本中找到一个最适合的使得频谱效率最大，其中P_t为大于1的整数；是长度为Nt的复向量；

S7、设备1以作为发射波束矢量向设备2发射，设备2通过重复S1-S6的步骤找到最优的接收矢量

S8、经过反复迭代，对于设备1和设备2来说，当相邻两次找到的波束矢量，即和相同时迭代终止，并将最终找到的和作为设备1和设备2最优的波束矢量。

实施例1、

子载波总数为512，采样频率为1GHz，设备1和设备2都有20根天线，码本中的波束数目为40个，构造字典时以5度为一间隔，CM4是非视距信道，有多条多径。

如图2所示，802.11.ad信道波束搜索的的成功概率性能曲线图，图2中横坐标横坐标是每个设备每次接收时的重复测量次数，在信噪比为0dB的条件下，每个点都仿真1000次。

根据图2可以看出成功概率随着测量次数的增加而增大。

Claims (3)

1.利用压缩感知的双端时域波束搜索方法，其特征在于，包括如下步骤：

S1、令设备1的收发天线数为Nt，所述设备1的码本中的波束数目为Ct，所述设备1使用全向天线向设备2发射，发射波束矢量为 ${\stackrel{\&RightArrow;}{w}}_1=\left[\begin{array}{c}1\\ 0\\ .\\ .\\ .\\ 0\end{array}\right],$ 所述发射波束长度为Nt，所述设备1使用OFDM技术在频域发射的序列为[1,1,...,1]，所述发射的序列长度为N

令设备2的收发天线数为Nr，所述设备2的码本中的波束数目为Cr，所述设备2的天线接收到的第m个子载波信号向量为所述设备2的接收端使用P_r个接收矢量来接收信号，任意一个接收矢量都是长度为Nr的向量，所述接收矢量中的每个元素的值从集合[1,i,-1,-i]中随机选择，组成一个测量矩阵所述测量矩阵Φ_r每一行都对应一次接收，子载波信号为测量信号向量为其中，是长度为Nr的噪声向量，H_m为第m个频点的阶数为Nr×Nt信道矩阵，所述信道矩阵中第x行第y列的元素表示从发射端第y根天线到接收端第x根天线间的频域信道冲击响应，其中，m＝1,2,...,N，y＝1,2,...,Nt，x＝1,2,...,Nr，d＝1,2,...P_r，i为虚数单位，是噪声向量，中的每个元素对应一个测量值，()^T是矩阵的转置运算，P_r为大于1的整数，N、Nt、Nr、Ct和Cr均为大于1的整数；

S2、根据S1所述构建字典矩阵为D，D的每一列对应[-90°，90°]中的一个角度，S1所述信号可以在D下展开，并且是稀疏的，展开系数为复数，是在D下的展开系数；

S3、、使用N个任务的正交匹配追踪算法(Orthogonal Matching Pursuit，OMP)联合所有S1所述子载波信号共同恢复出所有的具体为：

S31、 ${\stackrel{\&RightArrow;}{y}}_{m,r}={\mathrm{\&Phi;}}_r{\stackrel{\&RightArrow;}{g}}_{m,r}={\mathrm{\&Phi;}}_{r}D{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{m,r}+{\stackrel{\&RightArrow;}{n}}_{m,r}^{\&prime;}=V_r{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{m,r}+{\stackrel{\&RightArrow;}{n}}_{m,r}^{\&prime;},$ V_r＝Φ_rD，所述可以在V_r下展开，是在V_r下的展开系数；

S32、从S31所述V_r中找到一列使得最大，构造矩阵算出所有在V_c下的展开系数 $[{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{1,r,c},{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{2,r,c},...,{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{N,r,c}]={\left({V_c}^H{V}_c\right)}^{-1}{V_c}^H[{\stackrel{\&RightArrow;}{y}}_{1,r},{\stackrel{\&RightArrow;}{y}}_{2,r},...,{\stackrel{\&RightArrow;}{y}}_{N,r}],$ 表示当前恢复程度的剩余量 $e_r=[{\stackrel{\&RightArrow;}{y}}_{1,r},{\stackrel{\&RightArrow;}{y}}_{2,r},...,{\stackrel{\&RightArrow;}{y}}_{N,r}]-V_c[{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{1,r,c},{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{2,r.c}...,{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{N,r,c}],$ 其中，()^-1是矩阵的求逆运算，()^H是矩阵的共轭转置运算，|·|表示取复数的幅度，||·||₂表示向量的二范数运算；

S33、从V_r中找到使得最大，其中是矩阵e_r中的第m列，将加到S32所述V_c中得到更新后的算出在更新后的V_c下的展开系数；

S34、循环S34到S33，直到e_r的F范数小于的F范数的a倍时停止，结合S2所述字典矩阵D中的列向量与对应位置的系数线性组合恢复出其中，a为门限值，0<a<1,且a为实数；

S4、恢复出所有频点的记作从码本中找到一个适合的波束矢量使得频谱效率最大， ${\stackrel{\&RightArrow;}{c}}_1=\underset{\stackrel{\&RightArrow;}{c}}{\arg }\max \left(\frac{1}{N}\underset{m=1}{\overset{N}{\mathrm{\&Sigma;}}}{\log }_2(1+{\mathrm{\&gamma;}}_{r,m})\right),$ 其中 ${\mathrm{\&gamma;}}_{r,m}=\frac{{\left|{\stackrel{\&RightArrow;}{c}}^T{\widehat{\stackrel{\&RightArrow;}{g}}}_{r,m}\right|}^2}{N_r{N}_t{\mathrm{\&sigma;}}^2},$ σ²是噪声的功率，是长度为Nr的复向量；

S5、使用N个任务的OMP联合所有S1所述测量信号共同恢复出所有的具体如下：

S51、 ${\stackrel{\&RightArrow;}{y}}_{m,t}={\mathrm{\&Phi;}}_t{\stackrel{\&RightArrow;}{g}}_{m,t}={\mathrm{\&Phi;}}_{t}D{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{m,t}+{\stackrel{\&RightArrow;}{n}}_{m,t}^{\&prime;}=V_t{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{m,t}+{\stackrel{\&RightArrow;}{n}}_{m,t}^{\&prime;},$ 所述可以在可以在V_t下展开，是在V_t下的展开系数，其中，V_t＝Φ_tD；

S52、从S51所述V_t中找到一列使得最大，构造矩阵算出所有在V_t下的展开系数 $[{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{1,r,c},{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{2,r,c},...,{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{N,r,c}]={\left({V_c}^H{V}_c\right)}^{-1}{V_c}^H[{\stackrel{\&RightArrow;}{y}}_{1,r},{\stackrel{\&RightArrow;}{y}}_{2,r},...,{\stackrel{\&RightArrow;}{y}}_{N,r}],$ 表示当前恢复程度的剩余量 $e_r=[{\stackrel{\&RightArrow;}{y}}_{1,r},{\stackrel{\&RightArrow;}{y}}_{2,r},...,{\stackrel{\&RightArrow;}{y}}_{N,r}]-V_c[{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{1,r,c},{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{2,r.c}...,{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{N,r,c}],$ 其中，()^-1是矩阵的求逆运算，()^H是矩阵的共轭转置运算，|·|表示取复数的幅度，||·||₂表示向量的二范数运算；

S53、从V_r中找到使得最大，其中是矩阵e_r中的第m列，将加到S32所述V_c中得到更新后的算出在更新后的V_c下的展开系数；

S54、循环S51到S53，直到e_r的F范数小于的F范数的a倍时停止，结合S2所述字典矩阵D中的列向量与对应位置的系数线性组合恢复出记为从码本中找到一个最适合的波束矢量使得频谱效率最大，其中σ²是噪声的功率，是长度为Nr的复向量其中，a为门限值，0<a<1,且a为实数，σ²是噪声的功率，是长度为Nr的复向量；

S6、设备2的接收端使用P_t个接收矢量来接收信号，与S1.3的步骤一样，结合所有频点的信号运用多任务正交匹配追踪算法，测量信号为 ${\stackrel{\&RightArrow;}{y}}_{m,t}={\mathrm{\&Phi;}}_t{\stackrel{\&RightArrow;}{g}}_{m,t}={\mathrm{\&Phi;}}_{t}D{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{m,t}+{\stackrel{\&RightArrow;}{n}}_{m,t}^{\&prime;}=V_t{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}_{m,t}+{\stackrel{\&RightArrow;}{n}}_{m,t}^{\&prime;},$ 其中V_t＝Φ_tD， ${\mathrm{\&Phi;}}_t={\left[\begin{array}{ccc}{\stackrel{\&RightArrow;}{\mathrm{\&phi;}}}_{t,1}& ...& {\stackrel{\&RightArrow;}{\mathrm{\&phi;}}}_{t,P_t}\end{array}\right]}^T$ 是测量矩阵，每一行接收矢量对应一次测量，任意一个接收矢量(d＝1,2,...,P_r)都是长度为Nt的向量，每个元素的值从集合[1,i,-1,-i]中随机选择，是在V_t下的展开系数，重复测量P_t次，根据恢复出记为当所有频点的(m＝1,2,...,N)都恢复出来后，从码本中找到一个最适合的使得频谱效率最大，其中P_t为大于1的整数；是长度为Nt的复向量；

S7、设备1以作为发射波束矢量向设备2发射，设备2通过重复S1-S6的步骤找到最优的接收矢量

S8、经过反复迭代，对于设备1和设备2来说，当相邻两次找到的波束矢量，即和相同时迭代终止，并将最终找到的和作为设备1和设备2最优的波束矢量。

2.根据权利要求1所述利用压缩感知的双端时域波束搜索方法，其特征在于：对于任意角度q，S2所述字典矩阵D中的对应列为 $\left[\begin{array}{c}1\\ e^{\mathrm{i\&pi;}\sin \left(\mathrm{\&theta;}\right)}\\ e^{i2\mathrm{\&pi;}\sin \left(\mathrm{\&theta;}\right)}\\ .\\ .\\ .\\ e^{i(\mathrm{Nt}-1)\mathrm{\&pi;}\sin \left(\mathrm{\&theta;}\right)}\end{array}\right].$

3.根据权利要求1所述利用压缩感知的双端时域波束搜索方法，其特征在于：S54所述a＝0.05。