CN103944578A - 一种多信号的重构方法 - Google Patents

Abstract

本发明属于无线通信技术领域，具体涉及一种在超宽带通信系统中采用压缩感知来联合恢复多个信号的低复杂度的改进贝叶斯算法。本发明提供了一种在无线通信系统中，利用低复杂度的改进贝叶斯压缩感知算法联合重建多个信号的方法。该方法首先对多个接收信号分段、滤波，再使用不同的测量矩阵对每个滤波之后的信号重新线性组合，在一系列利用了这多个信号之间相关性的低复杂度迭代运算后，可以测量出每个原始信号在同一特征基下的展开系数，从而实现对每个原始信号更加精确的重建。

Description

一种多信号的重构方法

技术领域

本发明属于无线通信技术领域，具体涉及在超宽带通信系统中的一种采用低复杂度的改进贝叶斯压缩感知算法来重构多个信号的方法。

背景技术

由于超带宽（UWB，Ultra Wideband）信号多为时域上的窄脉冲信号，导致其频域带宽过宽，按照Nyquist采样定律的要求，其采样速率将高达数GHz，这样就给模数转换器的设计带来巨大的挑战，并且大量的数据处理给数字信号处理器带来很大的负担，目前的硬件水平难以实现以上要求。压缩感知（CS,Compressed sensing）作为一种新颖的信号低速采样理论，突破了Nyquist采样定律的限制，已经有文献对其在UWB信道估计方面的应用展开了理论研究，并指出在合作通信的前提下，使用压缩感知理论可以将UWB信号的采样速率从GHz级降到目前硬件能够实现的范围，大大降低了UWB信号的采样难度。

近年来很多研究表明，无线多径信道的冲激响应往往呈现为数量相对较少的几簇重要路径集，尤其是在传输带宽较大或采用多天线系统时，这种稀疏特性会更加明显。随着无线通信系统的宽带化和高速化发展，对稀疏信道估计的研究已成为新的焦点。近年来新兴的压缩感知理论指出，可以根据少量的观测值高效率重构稀疏信号。与传统信道估计方法相比，基于压缩感知的信道估计方法可以大大降低导频开销，提高频谱利用率和估计性能。

压缩感知作为一个新的采样理论，通过开发信号的稀疏特性，不必直接测量信号，而是通过一个随机的测量矩阵测量原来信号的线性组合，压缩感知需要的测量样本数要远小于Nyquist采样率条件下的样本数，通过对这些测量样本使用各种各样的非线性重建算法，可以完美地重建信号。让代表原始信号，Ψ是它的一个基，在Ψ中展开的展开系数大部分趋近于0，我们不用直接测量而是测量样本点较少的 Ψ^T是Ψ的转置，Φ是K×N阶测量矩阵（K＜＜N），通过各种的压缩感知算法可实现对的重建，所以简单来说，压缩感知的目的就是实现

在利用压缩感知前要考虑的一个重要的问题是：需要找到一个合适的完备或非完备展开基让信号基于此展开基是稀疏的，因此设计展开基就成为了运用压缩感知的关键。常用于超宽带的压缩感知的展开基有时间稀疏基和多路基,后者的性能要比前者好，其详细方法可参考文献：Jose L.Paredes,Gonzalo R.Arce and Zhongmin Wang.Ultra-widebandcompressed sensing:Channel estimation[J].IEEE J.Select.Topics SignalProcess,2007,1(3):383-395.接收信号近似为多路成分的叠加，使用多路字典，信道估计问题可以直接转换为压缩感知的问题，但是不同的多路成分会呈现出不同程度的由散射和折射引起的失真。为了更充分利用UWB信号内部的统计特性，我们考虑使用基于信号统计特性的特征展开基。

常用的压缩感知的方法主要有：基追踪算法（BP，Basis Pursuit）、匹配追踪算法（MP，Matching Pursuit）、正交匹配追踪算法（OMP，Orthogonal Matching Pursuit）以及贝叶斯压缩感知算法（BCS，Bayesian Compressive Sensing）。其中，贝叶斯压缩感知理论将统计学习理论（SLT，Statistical Learning Theory）和相关向量机（RVM，Relevance VectorMachine）思想引入到压缩感知理论中来，给需要重构向量中的每个值设置由超参数（Hyperparameters）控制的后验概率密度函数，在超参数的更新过程中，趋近于零的非重要多径所对应的超参数将变得无限大，导致该多径的后验概率趋近于零，而真正的重要多径所对应的矩阵中的列向量起到了相关向量的作用，其详细方法可参考文献：Shihao Ji,YaXue and Lawrence Carin.Bayesian compressive sensing.IEEE Transactions on SignalProcessing[J],2008,56(6):2346-2356.

利用特征基结合贝叶斯压缩感知的算法需要对矩阵求逆，特别的，特征展开基需要利用信号的统计特性来提高对信号的恢复能力，信号的长度比较长，此时矩阵的阶数会比较大，求逆运算会大大增加运算的复杂度，使得应用推广变得十分困难，因此，研究一种低复杂度的改进贝叶斯压缩感知算法是一项有创新性和重要实际意义且具挑战性的任务。

发明内容

本发明的目的在于一种在无线通信系统中，采用低复杂度的改进贝叶斯压缩感知算法来重构多个信号的方法。

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

S1、假设有D个信号，包括D≥1，且D为整数，将均匀分为L段，每段信号的长度相等且为大于1的整数，使用完备基Ψ使得这D个信号在Ψ中是稀疏的，所述每段信号都可以在同一个完备基中展开，每段信号对应不同的展开系数其中，Ψ是由特征基向量构成的正交方阵，L＞1,且L为自然数，完备基是一种特殊的矩阵，矩阵的列向量之间是线性无关的，任意一个信号都可以用这个矩阵中的列向量和对应的展开系数线性加和来表示，特征基向量是指对矩阵进行特征值分解后得到的特征向量，这些向量是线性无关的，可用于组成完备基；

S2、对S1所述的分为L段的信号进行滤波处理，得到信号其中，信号的长度和的长度相同，每段信号的长度相等且为大于1的整数；

S3、使用测量矩阵对S2所述信号进行线性组合，其中，测量矩阵记作Φ₁,Φ₂,...,Φ_D，Φ_D是K×N阶的矩阵，N为不使用压缩感知时需要的测量样本数，K为使用压缩感知时需要的测量样本数，Φ_D中的每个元素独立服从均值为0的正态分布，K＞0,且K为自然数，N＞0,且N为自然数，L＞0，且L为自然数，N＝N_e×L，K根据实际需要设置具体数值，且K＜＜N；

S4、初始化，包括：

S41、初始化D个信号共同的参数其中，D个信号是同分布的，所述的D个信号展开系数向量有相同的方差向量，方差的倒数向量记为参数 $\stackrel{\→}{\mathrm{\β}}=[{\mathrm{\β}}_1,{\mathrm{\β}}_2,...,{\mathrm{\β}}_N],\stackrel{\→}{\mathrm{\β}}>0,$ 且为实数；

S42、初始化Σ_i和i＝1,2,...,D，其中，第i个信号的展开系数的后验概率密度函数是均值为协方差为Σ_i的高斯分布，是长度为N的向量，Σ_i是N×N阶的方阵，经运算，Σ_i＝(α₀Φ_i ^TΦ_i+A)^-1，A是N×N阶的对角矩阵，主对角线位置上的元素是中元素按顺序排列，其余位置上的元素都是0，（）^-1是矩阵的求逆运算；对D个信号的均值和协方差都进行初始化，σ₀ ²为噪声方差，σ₀ ²＝测量样本的方差/100；

S43、初始化l_i，i＝1,2,...,D，其中，l_i为的展开系数的后验概率密度函数的对数似然函数，l_i≤0,且l_i为实数，D个信号的展开系数的联合后验概率密度函数的对数似然函数就是将他们单独的似然函数加和，记为l；

S44、初始化i＝1,2,...,D，其中，是长度为N的向量，是长度为N的向量；

S5、迭代更新Σ_i、得到

S6、重构信号：就是信号展开系数的近似估计，用作为展开系数，计算得到的就是使用低复杂度的改进贝叶斯压缩感知算法对原信号的重构，由此，D个信号可以同时得到恢复。

进一步地，S1所述展开系数范围为实数，其中，所述系数范围由在Ψ_S上的展开决定，取Ψ的N_e个列向量，构成Ψ的子空间，记作Ψ'，由此构造矩阵其中，U'＞N_e≥1，且N_e为自然数，U为原始信号长度。

进一步地，S1所述展开系数为服从零均值高斯分布的随机变量，是在Ψ_S下的展开系数，其中，i＝1,2,...,D。

进一步地，S3所述测量样本记为对于任意的一个 和是噪声向量，（）^T是矩阵的转置运算。

进一步地，S5所述迭代更新过程就是寻求最大l的过程，将l的变化量Δl作为判断迭代终止的条件，对应每个信号引入向量和和矩阵C_i来辅助运算，C_i是N×N阶的方阵，通过公式C_i＝σ₀ ²I+Φ_iA^-1Φ_i ^T来计算，I是N阶单位矩阵，主对角线上元素均为1，其余位置为0，中的第m个元素为S_i,m＝φ_i.m ^TC_i ^-1φ_i,m，中的第m个元素为 为Φ_i的第m个列向量，其中，m＝1,2,...,N。

进一步地，S4所述的初始化是通过逐个增加中一个元素β_i，D个信号的Σ_i、 初始化都分别计算，其中，i＝1,2,...,D。

进一步的，S5所述的更新分为两个部分，包括：对的更新和对的更新，分别对和中元素求导，中任意一个元素β_i的更新公式为 是更新后的β_i，M_i是中含有β_i成分的位置号的集合，则是位置号m所在的段落号，Σmm是Σ对角线上的第m个元素，μ_m是中第m个元素，中元素k_l(l＝2,3,...,L)的更新公式为k_l ^new是更新后的k_l，N_l是第l段位置号的集合，则是位置号n对应于第一段中的位置号，其中，为的第一段N_e个系数方差的倒数向量，为的第一段N_e个系数方差的倍数向量， 和k_L每个元素都是大于0的实数，可以重新写为 $\stackrel{\→}{\mathrm{\β}}=[{\stackrel{\→}{\mathrm{\β}}}_1,k_2{\stackrel{\→}{\mathrm{\β}}}_1,...,k_L{\stackrel{\→}{\mathrm{\β}}}_1].$

进一步地，S5所述代更新Σ_i、为：更新一个β_j，有L个元素同时变化，同理，更新一个中的元素都会导致中N_e个元素的同时变化,通过逐步改变中元素并循环更新Σ_i、的办法可以实现，每更新一个β_j，这种内部循环需要L次，每更新一个k_l这种内部循环需要N_e次，当更新完全部的β_j和k_l就完成了一次迭代，其中，i＝1,2,...,D，j＝1,2,...,N_e，l＝2,3,...,L。

进一步地，S5所述迭代更新终止于变化率＜Th时，即当变化率＜Th时认为收敛，终止迭代，得到其中，

$\mathrm{\Δ}l=\frac{1}{2}[\frac{{Q_i}^2{\mathrm{\Δ\α}}_i}{S_i\mathrm{\Δ}{\mathrm{\α}}_i-{\widetilde{\mathrm{\α}}}_i{\mathrm{\α}}_i}-\log ({\widetilde{\mathrm{\α}}}_i{\mathrm{\α}}_i-S_i\mathrm{\Δ}{\mathrm{\α}}_i)+\log {\widetilde{\mathrm{\α}}}_i{\mathrm{\α}}_i],$ Δα_i是α_i更新后的值减去更新前的值，其中，Th为经验值，1＜Th＜0,且Th为实数，i＝1,2,...,D。

本发明的有益效果是：本发明充分利用特征基的特性改进了贝叶斯压缩感知算法，提高了恢复性能，避开了贝叶斯算法中复杂的矩阵求逆过程，有效地减少信号恢复运算复杂度。特别的，利用多个信号内部的统计相关特性，同时测量出每个原始信号在同一展开基下的展开系数，从而实现对每个原始信号更加精确的重建。

附图说明

图1是本发明利用低复杂度的改进贝叶斯压缩感知算法联合重建多个信号的方法的结构图。

图2是本发明的流程图。

图3是本发明用于超宽带信道估计的相关系数的曲线图。

图4是本发明用于超宽带信道估计的误比特率（BER）性能曲线图。

具体实施方式

假设有D个（D是大于等于1的整数）信号：将这多个原始信号每个都均分为L段，每段的长度相等且为大于1的整数，任意一个信号的每一段都可以被同一个完备基和不同的展开系数完全表示，每一段对应不同的展开系数，使用特殊的完备基Ψ使得这D个信号在这个完备集中均是稀疏的，即展开系数大部分为0或接近0，Ψ是由特征基向量构成的正交方阵，我们只取Ψ的N_e个主要的列向量，构成Ψ的子空间，记为Ψ'，由此构造矩阵考虑在Ψ_S上的展开系数；对每个信号滤波，让的每一段在Ψ的子空间Ψ'上投影，滤波之后的信号记为（i＝1,2,...,D）是在Ψ_S下的展开系数；L为大于1的整数；N_e是大于等于1的整数，N_e小于中每一段的长度。

使用不同的测量矩阵对重新线性组合，记测量矩阵分别为Φ₁,Φ₂,...,Φ_D，均是K×N阶的矩阵，对于N长的信号来说，不使用压缩感知时需要N个测量样本数，K是使用压缩感知时需要的测量样本数。矩阵Φ₁,Φ₂,...,Φ_D中的每个元素独立服从均值为0的正态分布；对应的D个压缩感知测量样本记为对于任意的一个（i＝1,2,...,D）， 和是噪声向量，（）^T是矩阵的转置运算；K和N均为大于0的整数，N=N_e×L。

初始化D个信号共同的参数D个信号是同分布的，它们的展开系数向量有相同的方差向量，方差的倒数向量记为参数其中第一段的N_e个元素记为 元素可以任意取大于0的数，以后每一段中任意一个系数的方差都是第一段中相同位置系数方差的倍数，与所在段落有关，倍数向量记为其中第一个元素一直是1，其余元素都是大于0的实数，可以重新写为噪声方差σ₀ ²在整个算法过程中保持不变，且

初始化Σ_i和（i＝1,2,...,D）；与只有一个信号的情况类似，经过贝叶斯法则分析，第i个信号的展开系数的后验概率密度函数是均值为协方差为Σ_i的高斯分布，是长度为N的向量，Σ_i是N×N阶的方阵，经运算，Σ_i＝(α₀Φ_i ^TΦ_i+A)^-1，A是N×N阶的对角矩阵，主对角线位置上的元素是中元素按顺序排列，其余位置上的元素都是0，（）^-1是矩阵的求逆运算；对D个信号的均值和协方差都进行初始化。

的展开系数的后验概率密度函数的对数似然函数为l_i，D个信号的展开系数的联合后验概率密度函数的对数似然函数就是将他们单独的似然函数加和，记为l，和l_i一样是小于等于0的实数，迭代过程就是寻求最大l的过程，将l的变化量Δl作为判断迭代终止的条件；对应每个信号引入向量和和矩阵C_i（i＝1,2,...,D）来辅助运算，同样进行初始化，和都是长度为N的向量，C_i是N×N阶的方阵，通过公式C_i＝σ₀ ²I+Φ_iA^-1Φ_i ^T来计算，I是N阶单位矩阵，主对角线上元素均为1，其余位置为0，中的第m（m＝1,2,...,N）个元素为S_i,m＝φ_i.m ^TC_i ^-1φ_i,m，中的第m个元素为 为Φ_i的第m个列向量。

初始化Σ_i、（i＝1,2,...,D），通过逐个增加中一个元素β_i（i＝1,2,..,.N），D个信号的Σ_i、初始化都分别计算，为了表示一般性，下面叙述过程和公式中都把表示第i个信号的下标i去掉；从中的第一个元素β₁开始，此时，对于任意一个信号，A＝β₁ 算出 $\mathrm{\Σ}=\frac{1}{({\mathrm{\α}}_0{{\stackrel{\→}{\mathrm{\φ}}}_1}^T{\stackrel{\→}{\mathrm{\φ}}}_1+{\mathrm{\β}}_1)},\stackrel{\→}{\mathrm{\μ}}={\mathrm{\α}}_0\mathrm{\Σ}{{\stackrel{\→}{\mathrm{\φ}}}_1}^T\stackrel{\→}{y},S_m={\mathrm{\α}}_0{{\stackrel{\→}{\mathrm{\φ}}}_m}^T{\stackrel{\→}{\mathrm{\φ}}}_m-{{\mathrm{\α}}_0}^2{{\stackrel{\→}{\mathrm{\φ}}}_m}^T{\stackrel{\→}{\mathrm{\φ}}}_m\mathrm{\Σ}{{\stackrel{\→}{\mathrm{\φ}}}_1}^T{\stackrel{\→}{\mathrm{\φ}}}_m,$ 之后逐次增加β₂,β₃,...,β_N；代表增加一个元素之后的Σ， $\widetilde{\mathrm{\Σ}}=\left[\begin{array}{cc}\mathrm{\Σ}+{{\mathrm{\α}}_0}^2{\mathrm{\Σ}}_{\mathrm{ii}}\mathrm{\Σ}{\mathrm{\Φ}}^T{\stackrel{\→}{\mathrm{\φ}}}_i{{\stackrel{\→}{\mathrm{\φ}}}_i}^T\mathrm{\Φ\Σ}& -{\mathrm{\α}}_0{\mathrm{\Σ}}_{\mathrm{ii}}\mathrm{\Σ}{\mathrm{\Φ}}^T{\stackrel{\→}{\mathrm{\φ}}}_i\\ -{\mathrm{\α}}_0{\mathrm{\Σ}}_{\mathrm{ii}}{\left(\mathrm{\Σ}{\mathrm{\Φ}}^T{\stackrel{\→}{\mathrm{\φ}}}_i\right)}^T& {\mathrm{\Σ}}_{\mathrm{ii}}\end{array}\right],$ 每添加一个β_i，Σ的阶数加1，对角线上加一个实数元素Σ_ii， 是Σ右上角增加的列向量，是左下角增加的行向量，是与Σ同阶数的矩阵；代表增加一个元素之后的 $\stackrel{\→}{\widetilde{\mathrm{\μ}}}=\left[\begin{array}{c}\stackrel{\→}{\mathrm{\μ}}-{\mathrm{\α}}_0{\mathrm{\μ}}_i\mathrm{\Σ}{\mathrm{\Φ}}^T{\stackrel{\→}{\mathrm{\φ}}}_i\\ {\mathrm{\μ}}_i\end{array}\right],$ 同样的的长度加1,增加一个实数元素μ_i，μ_i＝Σ_iiQ_i，是与同样长度的列向量。代表增加一个元素之后的S_m， 代表增加一个元素之后的Q_m，对于和长度不变，其中的每个元素都要更新，是长度为N的列向量。

的更新可以分为两个部分：对的更新和对的更新。更新的原则是使得这多个信号的联合对数似然值l的值最大，分别对和中元素求导；中任意一个元素的更新公式为： 是更新后的β_j，M_j是中含有β_j成分的位置号的集合，则是位置号m所在的段落号，Σ_i,mm是第i个信号的协方差Σ_i对角线上的第m个元素，μ_i,m是第i个信号的中第m个元素；中元素k_l(l＝2,3,...,L)的更新公式为 是更新后的k_l，N_l是第l段位置号的集合，则是位置号n对应于第一段中的位置号。

随着和的变化，用低复杂度的方法更新Σ_i、（i＝1,2,...,D），若更新β_i(i＝1,2,...,N_e)，有L个元素同时变化，同理，更新一个中的元素都会导致中N_e个元素的同时变化。通过逐步改变中元素并循环更新Σ_i、的办法可以实现，每更新一个β_i这种内部循环需要L次，每更新一个k_l这种内部循环需要N_e次，当更新完全部的β_i和k_l就完成了一次迭代；假设有任意n个中的元素需要更新，记为[α₁,α₂,...,α_n]，更新之后记为它们在中相应的位置为P＝[p₁,p₂,...,p_n]，Σ_i、的更新公式为：（特别的，D个信号的Σ_i、在更新时互不干扰，每个信号都单独计算，为了表示一般性，下面的公式都把表示第i个信号的下标i去掉，如将Σ_i表示为Σ）：

${\widetilde{\mathrm{\κ}}}_{(\mathrm{1,2},...i-1),n}=\frac{\mathrm{\Δ}{\mathrm{\α}}_i}{{\widetilde{\mathrm{\Σ}}}_{(\mathrm{1,2},...,i-1),p_i{p}_i}\mathrm{\Δ}{\mathrm{\α}}_i+1}$

$\mathrm{\Δ}{\mathrm{\α}}_i={\widetilde{\mathrm{\α}}}_i-{\mathrm{\α}}_i$

记是更新了α₁,α₂,...,α_i后的Σ，是更新了α₁,α₂,...,α_i之后的是更新了α₁,α₂,...,α_i之后的的第m个元素，是更新了α₁,α₂,...,α_i的的第m个元素，是更新了α₁,α₂,...,α_i-1的Σ的第p_i列，是更新了α₁,α₂,...,α_i-1的的第pi个元素，是更新了α₁,α₂,...,α_i-1的的第pi个元素，是更新了α₁,α₂,...,α_i-1的的第p_i个元素。

分别计算每个信号此次更新完或所带来Δl_i（i＝1,2,...,D）；对于每个信号来说，中有多个元素同时变化时的总的Δl_i就是依次变化中一个元素的Δl_i的叠加，当只有一个α_j变化时Δl_i的计算公式为：

$\mathrm{\Δ}{l}_i=\frac{1}{2}[\frac{{Q_{i,j}}^2{\mathrm{\Δ\α}}_j}{S_{i,j}\mathrm{\Δ}{\mathrm{\α}}_j-{\widetilde{\mathrm{\α}}}_j{\mathrm{\α}}_j}-\log ({\widetilde{\mathrm{\α}}}_j{\mathrm{\α}}_j-S_{i,j}\mathrm{\Δ}{\mathrm{\α}}_j)+\log {\widetilde{\mathrm{\α}}}_j{\mathrm{\α}}_j].$

每次迭代都更新一遍和中的值，将这些Δl_i加和(包括由的更新导致的Δl_i和由的更新导致的Δl_i），记为第i个信号本次迭代的总的Δl_i，将D个信号的Δl_i加和就是本次迭代的联合对数似然函数的总变化，记为Δl。

每一次迭代的Δl和上一次迭代的Δl比较变化率，若变化率小于设定的值则认为收敛，终止迭代，得到最后的（i＝1,2,...,D），进入下面的步骤；若不满足收敛条件则继续迭代更新。

重构信号，就是第i个信号展开系数的近似估计，用作为展开系数，计算得到的就是使用低复杂度的改进贝叶斯压缩感知算法对原信号的重建，由此，D个信号可以同时得到恢复。

以下将结合附图详细介绍本发明的具体方法步骤：

图1是本发明低复杂度的改进贝叶斯压缩感知算法联合重建多个信号的方法的结构图。中元素对应同一方差的倒数向量是分为L段，第一段记为 之后每一段都是第一段的倍数，倍数向量如：第二段第L段测量矩阵分别为Φ₁,Φ₂,...,Φ_D，均是K×N阶的矩阵，矩阵Φ₁,Φ₂,...,Φ_D中的元素服从均值为0的正态分布，压缩感知D个测量样本记为其中第i个测量样本 和是噪声向量。α₀是噪声方差的倒数，噪声方差一般可以取经验值σ₀ ²＝D个测量样本的方差/100，α₀在迭代更新的过程中始终保持不变。

图2是本发明的算法流程图。初始化 中元素和k₂,k₃,...,k_L的值可以随机取大于0的值，但是，最后组成的必须符合的对应关系，在迭代更新的过程中也要保持这种对应关系。为避免N×N阶矩阵的求逆运算，使用逐一添加中元素的方式来初始化Σ_i、（i＝1,2,...,D），如同单信号，D个信号要分开计算，每个任务要循环使用N次添加公式。设置迭代的最大次数，如：1000次。迭代开始后，首先更新中的元素的值，可以从β₁开始，更新β₁实际上意味着中有L个和β₁有关的元素需要同时更新：D个信号要分开计算，使用更新公式更新Σ_i、 （i＝1,2,...,D），使用同样的方法更新中剩下的元素，算出每个任务由的更新引起的总的Δl_i。同样的中的元素也需要更新，与更新β₁相类似，每更新一个元素，都有N_e个元素变化，使用更新公式更新Σ_i、（i＝1,2,...,D），类似的方法算出每个任务由的更新引起的总的Δl_i，最后可以将两个阶段加起来就是本次迭代的Δl。当相邻两次迭代的Δl变化不大时就可以判断为收敛，终止迭代，同时恢复D个信号，就是第i个信号展开系数的近似估计，用作为展开系数，计算得到的就是使用低复杂度的改进贝叶斯压缩感知算法对原信号的重建。

图3是使用图1的利用低复杂度的改进贝叶斯压缩感知算法联合重建多个信号的方法的结构图、图2的算法流程图，应用到具体的通信系统中，这里用来对超宽带通信进行信道估计，仿真得到本发明用于信道估计的相关系数的曲线图。本例的仿真系统是属于室内环境的超宽带通信系统，它主要仿真参数是：采样率为20.48GHz，2ns的clean pulse（CP），将信号截成50段，选择每一段的长度等于CP的长度41，并基于UWB信号的特征结构来构建本征字典，滤波时选择N_e＝5。每次传输5000个2PAM信号，前50个用于信道估计，K=150，选择单信号，2个信号和4个信号的情况做对比。图3呈现的是压缩感知算法的重建性能，相关系数的公式是：相关系数越大说明重建性能越好。从图中可以看出，在使用特征展开基时，多信号联合重建比单个信号重建有更好的重建性能，4个信号的联合重建要优于2个信号的联合重建，且使用低复杂度多信号改进贝叶斯算法的重建性能远好于普通多信号的BCS算法（不考虑内部和的对应特性）的性能。

图4是本发明用于超宽带信道估计的误比特率（BER）性能曲线图。仿真参数和图3中一样，除去信道估计50个符号的4950个符号，信噪比从8dB到18dB，每个点仿真200次。可以看出使用低复杂度的4信号改进贝叶斯算法（使用特征展开基）的性能仍然是最好的。

以上实例仅为本发明的优选例子而已，本发明的使用并不局限于该实例，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (9)

1.一种多信号的重构方法，其特征在于，包括如下步骤：

S1、假设有D个信号，包括D≥1，且D为整数，将均匀分为L段，每段信号的长度相等且为大于1的整数，使用完备基Ψ使得这D个信号在Ψ中是稀疏的，所述每段信号都可以在同一个完备基中展开，每段信号对应不同的展开系数其中，Ψ是由特征基向量构成的正交方阵，L＞1,且L为自然数，完备基是一种特殊的矩阵，矩阵的列向量之间是线性无关的，任意一个信号都可以用这个矩阵中的列向量和对应的展开系数线性加和来表示，特征基向量是指对矩阵进行特征值分解后得到的特征向量，这些向量是线性无关的，可用于组成完备基；

S2、对S1所述的分为L段的信号进行滤波处理，得到信号其中，信号的长度和的长度相同，每段信号的长度相等且为大于1的整数；

S3、使用测量矩阵对S2所述信号进行线性组合，其中，测量矩阵记作Φ₁,Φ₂,...,Φ_D，Φ_D是K×N阶的矩阵，N为不使用压缩感知时需要的测量样本数，K为使用压缩感知时需要的测量样本数，Φ_D中的每个元素独立服从均值为0的正态分布，K＞0,且K为自然数，N＞0,且N为自然数，L＞0，且L为自然数，N＝N_e×L，K根据实际需要设置具体数值，且K＜＜N；

S4、初始化，包括：

S41、初始化D个信号共同的参数其中，D个信号是同分布的，所述的D个信号展开系数向量有相同的方差向量，方差的倒数向量记为参数 $\stackrel{\→}{\mathrm{\β}}=[{\mathrm{\β}}_1,{\mathrm{\β}}_2,...,{\mathrm{\β}}_N],\stackrel{\→}{\mathrm{\β}}>0,$ 且为实数；

S42、初始化Σ_i和i＝1,2,...,D，其中，第i个信号的展开系数的后验概率密度函数是均值为协方差为Σ_i的高斯分布，是长度为N的向量，Σ_i是N×N阶的方阵，经运算，Σ_i＝(α₀Φ_i ^TΦ_i+A)^-1，A是N×N阶的对角矩阵，主对角线位置上的元素是中元素按顺序排列，其余位置上的元素都是0，（）^-1是矩阵的求逆运算；对D个信号的均值和协方差都进行初始化，σ₀ ²为噪声方差，σ₀ ²＝测量样本的方差/100；

S43、初始化l_i，i＝1,2,...,D，其中，l_i为的展开系数的后验概率密度函数的对数似然函数，l_i≤0,且l_i为实数，D个信号的展开系数的联合后验概率密度函数的对数似然函数就是将他们单独的似然函数加和，记为l；

S44、初始化i＝1,2,...,D，其中，是长度为N的向量，是长度为N的向量；

S5、迭代更新Σ_i、 得到

S6、重构信号：就是信号展开系数的近似估计，用作为展开系数，计算得到的就是使用低复杂度的改进贝叶斯压缩感知算法对原信号的重构，由此，D个信号可以同时得到恢复。

2.根据权利要求1所述的一种多信号的重构方法，其特征在于：S1所述展开系数范围为实数，其中，所述系数范围由在Ψ_S上的展开决定，取Ψ的N_e个列向量，构成Ψ的子空间，记作Ψ'，由此构造矩阵其中，U'＞N_e≥1，且N_e为自然数，U为原始信号长度。

3.根据权利要求1所述的一种多信号的重构方法，其特征在于：S1所述展开系数为服从零均值高斯分布的随机变量，是在Ψ_S下的展开系数，其中，i＝1,2,...,D。

4.根据权利要求1所述的一种多信号的重构方法，其特征在于：S3所述测量样本记为对于任意的一个 和是噪声向量，（）^T是矩阵的转置运算。

5.根据权利要求1所述的一种多信号的重构方法，其特征在于：S5所述迭代更新过程就是寻求最大l的过程，将l的变化量Δl作为判断迭代终止的条件，对应每个信号引入向量和和矩阵C_i来辅助运算，C_i是N×N阶的方阵，通过公式C_i＝σ₀ ²I+Φ_iA^-1Φ_i ^T来计算，I是N阶单位矩阵，主对角线上元素均为1，其余位置为0，中的第m个元素为S_i,m＝φ_i.m ^TC_i ^-1φ_i,m，中的第m个元素为 为Φ_i的第m个列向量，其中，m＝1,2,...,N。

6.根据权利要求1所述的一种多信号的重构方法，其特征在于：S4所述的初始化是通过逐个增加中一个元素β_i，D个信号的Σ_i、初始化都分别计算，其中，i＝1,2,...,D。

7.根据权利要求1所述的一种多信号的重构方法，其特征在于：S5所述的更新分为两个部分，包括：对的更新和对的更新，分别对和中元素求导，中任意一个元素β_i的更新公式为 是更新后的β_i，M_i是中含有β_i成分的位置号的集合，则是位置号m所在的段落号，Σ_mm是Σ对角线上的第m个元素，μ_m是中第m个元素，中元素k_l(l＝2,3,...,L)的更新公式为klnew是更新后的k_l，N_l是第l段位置号的集合，则是位置号n对应于第一段中的位置号，其中，为的第一段N_e个系数方差的倒数向量，为的第一段N_e个系数方差的倍数向量， ${\stackrel{\→}{\mathrm{\β}}}_1=[{\mathrm{\β}}_1,{\mathrm{\β}}_2,...,{\mathrm{\β}}_{N_e}],\stackrel{\→}{k}=[1,k_2,k_3,...,k_L],$ 和k_L每个元素都是大于0的实数，可以重新写为

8.根据权利要求1所述的一种多信号的重构方法，其特征在于：S5所述代更新Σ_i、 为：更新一个β_j，有L个元素同时变化，同理，更新一个中的元素都会导致中N_e个元素的同时变化,通过逐步改变中元素并循环更新Σ_i、的办法可以实现，每更新一个β_j，这种内部循环需要L次，每更新一个k_l这种内部循环需要N_e次，当更新完全部的β_j和k_l就完成了一次迭代，其中，i＝1,2,...,D，j＝1,2,...,N_e，l＝2,3,...,L。

9.根据权利要求1所述的一种多信号的重构方法，其特征在于：S5所述迭代更新终止于变化率＜Th时，即当变化率＜Th时认为收敛，终止迭代，得到其中，

$\mathrm{\Δ}l=\frac{1}{2}[\frac{{Q_i}^2{\mathrm{\Δ\α}}_i}{S_i\mathrm{\Δ}{\mathrm{\α}}_i-{\widetilde{\mathrm{\α}}}_i{\mathrm{\α}}_i}-\log ({\widetilde{\mathrm{\α}}}_i{\mathrm{\α}}_i-S_i\mathrm{\Δ}{\mathrm{\α}}_i)+\log {\widetilde{\mathrm{\α}}}_i{\mathrm{\α}}_i],$ Δα_i是α_i更新后的值减去更新前的值，其中，Th为经验值，1＜Th＜0,且Th为实数，i＝1,2,...,D。