一种多信号的重构方法

摘要

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

权利要求

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

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

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

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

S4、初始化，包括：

S41、初始化D个信号共同的参数 其中，D个信号 是同分布的，所述的D个信号展开系数向量 有相同的方差向量，方差的倒数向量记为参数且 为实数；

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

S43、初始化 其中，为 的展开系数 的后验概率密度函数的对数似

然函数， 且 为实数，D个信号的展开系数的联合后验概率密度函数的对数似然函数就是将他们单独的似然函数加和，记为S44、初始化 其中，是长度为N的向量， 是长度为N的向量；

S5、迭代更新 Σi、 得到

所述 的更新分为两个部分，包括：对 的更新和对 的更新，分别对 和 中元素求导，中任意一个元素βi的更新公式为 是更新后的βi，Mi是 中含有βi成分的位置号的集合，则是位置号m所在的段落号，Σmm是Σ对角线上的第m个元素，μm是中第m个元素，中元素kl，l＝2,3,...,L的更新公式为 klnew是更新后的kl，Nl是第l段位置号的集合，则是位置号n对应于第一段中的位置号，其中， 为的第一段Ne个系数方差的倒数向量，为 的第一段Ne个系数方差的倍数向量， 和kL每个元素都是大于0的实数，

可以重新写为

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

2,...,Ne，l＝2,3,...,L；

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

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

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

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

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

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

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

是变化后的αi，其中，Th为经验值，1＜Th＜0,且Th为实数，i＝1,2,...,D。

说明书

一种多信号的重构方法

技术领域

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

背景技术

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

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

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

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

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

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

发明内容

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

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

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

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

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

[0013] S4、初始化，包括：

[0014] S41、初始化D个信号共同的参数 其中，D个信号 是同分布的，所述的D个信号展开系数向量 有相同的方差向量，方差的倒数向量记为参数且 为实数；

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

[0016] S43、初始化 其中，为 的展开系数 的后验概率密度函数的对数似然函数， 且 为实数，D个信号的展开系数的联合后验概率密度函数的对数似然函数就是将他们单独的似然函数加和，记为

[0017] S44、初始化 其中， 是长度为N的向量， 是长度为N的向量；

[0018] S5、迭代更新 Σi、 得到

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

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

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

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

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

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

[0025] 进一步的，S5所述 的更新分为两个部分，包括：对 的更新和对 的更新，分别对 和 中元素求导，中任意一个元素βi的更新公式为 是更新后的βi，Mi是 中含有βi成分的位置号的集合，则是位置号m所在的段落号，Σmm是Σ对角线上的第m个元素，μm是 中第m个元素， 中元素kl(l＝2,3,...,L)的更新公式为klnew是更新后的kl，Nl是第l段位置号的集合，则是位置号n对应
于第一段中的位置号，其中， 为 的第一段Ne个系数方差的倒数向量，为
的第一段Ne个系数方差的倍数向量，
和kL每个元素都是大于0的实数，可以重新写为

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

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

[0028] Δαi是αi更新后的值减去更新前的值，其中，Th为经验值，1＜Th＜0,且Th为实数，i＝1,2,...,D。

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

附图说明

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

[0031] 图2是本发明的流程图。

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

[0033] 图4是本发明用于超宽带信道估计的误比特率(BER)性能曲线图。

具体实施方式

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

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

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

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

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

[0039] 初始化Σi、 通过逐个增加 中一个元素βi(i＝1,2,..,.N)，D个信号的Σi、 初始化都分别计算，为了表示一般性，下面叙述过程和公式中都把表示第i个信号的下标i去掉；从 中的第一个元素β1开始，此时，对于任意一个信号，A＝β1， 算出
之后逐次增加β2,β3,...,βN；代表增加一个元素之后的Σ，
每添加一个βi，Σ的阶数加1，对角线上加一个实数
元素Σii， 是Σ右上角增加的列向量， 是左下角增加的
行向量， 是与Σ同阶数的矩阵； 代表增加一个元素之后的
同样的 的长度加1,增加一个实数元素μi，μi＝ΣiiQi， 是与
同样长度的列向量。 代表增加一个元素之后的Sm， 代表增加
一个元素之后的Qm， 对于 和 长度不变，其中的每个元素都要更新，
是长度为N的列向量。

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

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

[0042]

[0043]

[0044]

[0045]

[0046]

[0047]

[0048] 记 是更新了α1,α2,...,αi后的Σ， 是更新了α1,α2,...,αi之后的 是更新了α1,α2,...,αi之后的 的第m个元素， 是更新了α1,α2,...,αi的 的第m个元素， 是更新了α1,α2,...,αi-1的Σ的第pi列， 是更新了α1,α2,...,αi-1的 的第pi个元素， 是更新了α1,α2,...,αi-1的 的第pi个元素，是更新了α1,α2,...,αi-1的 的第pi个元素。

[0049] 分别计算每个信号此次更新完 或 所带来 对于每个信号来说，中有多个元素同时变化时的总的 就是依次变化 中一个元素的 的叠加，当只有一个αj变化 时 的计算公式为：

[0050]

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

[0052] 每一次迭代的 和上一次迭代的 比较变化率，若变化率小于设定的值则认为收敛，终止迭代，得到最后的 进入下面的步骤；若不满足收敛条件则继续迭代更新。

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

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

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

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

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

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

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