分布式压缩感知中充分稀疏源信号的交替迭代估计方法
【专利摘要】一种分布式压缩感知中充分稀疏源信号的交替迭代估计方法，具体涉及由混合信号的压缩观测值，重构源信号的方法，解决现有技术中，由于没有考虑混合信号压缩观测值的特点，采用先恢复混合信号再分离源信号的算法重构源信号而导致源信号精度低的问题。本发明所述的方法：在分布式压缩感知框架下，利用充分稀疏源信号的特点，采用交替估计的方法，从混合信号的观测值中估计源信号，每次迭代过程中，在不对混合信号进行重构的前提下，直接恢复得到源信号的估计值。本发明方法在信号压缩域从混合信号的观测值中分离源信号，降低了分离过程的复杂度。
【专利说明】分布式压缩感知中充分稀疏源信号的交替迭代估计方法
【技术领域】
[0001]本发明涉及分布式压缩感知领域以及盲源分离领域，具体涉及在分布式压缩感知信号处理框架下，从混合信号的观测值恢复充分稀疏源信号的方法。
【背景技术】
[0002] 传统的信号采样一般都基于奈奎斯特采样定理，即采样速率至少需要大于信号最高频率的2倍，才能利用采样得到的离散数据不失真的分离得到源信号。但是随着信息技术的发展，以奈奎斯特采样定理为基础的信号处理框架对前端模数转换器(ADC)的采样速率和处理速度提出了更高的要求，也对后端信息的传输、存储环节带来了巨大的挑战。
[0003]压缩感知(CompressedSensing,CS)理论是在 2004 年由 CandSs 和 Donoho 等人提出的。该理论指出，当信号满足稀疏性时，能以远低于奈奎斯特采样频率的速度对源信号进行全局观测，然后采用合适的重构算法精确重构出源信号。CS理论将采样和压缩合并成一步进行，极大的降低了信号的采样频率以及数据的传输、存储代价。若考虑稀疏度为K且长度为N的离散实信号X，即X满足:
[0004]X e Rn, I IXI 10 ^ K < < N (1)
[0005]其中I I.I Itl表示信号的Itl范数，即信号值不为零的个数。
[0006]令y为信号X的M次观测，则具体的压缩观测模型如式(2)所示，
[0007]y=Φχ (2)
[0008]其中，y是长度为M的观测信号，Φ为MXN维的观测矩阵。当观测矩阵满足RIP(Restricted Isometry Property)条件时，即可用I。范数意义下的优化问题来重构或者逼近源信号。即
[0009]arg min | | x | 10, s.t.y=Φx (3)
[0010]RIP 条件:对于任意的 qe Rl11 Jc{l，2，...,iV}且 |I| ≤ K,以及 K<m,0< δ ( I,如果有式(4)成立，则称观测矩阵满足RIP条件。其中O1为观测矩阵Φ中由索引I所指示的相关列组成大小为mX 111的子矩阵空间。
[0011](l--<5)|q[u <i#,q|2 <(l + 5)||q||j (4)
[0012]近年来，随着网络通信、多媒体技术、存储技术的发展，网络规模的越来越大，人们对分布式信息处理技术等先进技术的需求越来越迫切。压缩感知理论一般是应用于单信号场景，D.Baron等人充分发掘了多通道信号的信号内和信号间的相关性结构，提出了分布式压缩感知(DistributedCompressive Sensing, DCS)理论。在DCS信号处理框架下，先对各通道信号进行独立的CS观测，然后利用观测值采用联合重构算法恢复出源信号。在实际的分布式的场景下，传感器往往接收到的是多个信号的混合信号，并且源信号以及混合信号均是未知的。如果利用混合信号之间的相关性，采用现有的分布式压缩感知重构算法求解，我们仅仅能得到混合信号的估计值。但是，人们往往对源信号更感兴趣，而非混合信号。因此，分布式压缩感知框架下，研究从混合信号中估计源信号的方法是十分有必要的。[0013]目前，通用的在分布式压缩感知框架下，利用混合信号的观测值恢复源信号的方法如说明书附图图1所示。图中的“盲源分离”过程，也是现在的一个研究热点，主要研究如何从混合信号中恢复源信号及混合参数。利用CS采样，在不降低源信号的恢复精度的前提下，可以很大程度的缓解对硬件的采样速度以及传输带宽的压力。
[0014]图1中的S=[Si S2...S/表示J个长度为N的源信号，混合系统可以用PXJ维矩阵A表示，CS观测系统可以用观测矩阵Φ表示。则混合信号X可以表示为式(5)，观测信号Y可以表不为式(6)。
[0015]X=AS (5)
[0016]Y= Φ Xt (6)
[0017]由图1可知，通用的基于压缩感知的源信号估计方法主要包括两个步骤:
[0018]步骤一、采用现有的CS重构方法对混合信号进行重构，得到混合信号的估计值x。[0019]步骤二、结合现有的盲源分离算法，对混合信号i分离，得到源信号的估计值
[0020]由于压缩感知理论应用的前提是信号满足稀疏性或者可压缩性。根据源信号的稀疏程度，稀疏信号又可以分为充分稀疏源信号和非充分稀疏源信号。充分稀疏源信号是指多个源信号在同一时刻只有一个源信号非零，其他源信号为零或者趋近于零。非充分稀疏源信号是指源信号在同一时刻不止一个源信号的取值为零。
[0021]一般情况下，我们仅仅对源信号感兴趣，对混合信号并不感兴趣。上述方法没有充分考虑混合信号的压缩观测值的特性以及所携带的源信号的信息，算法分离源信号的效率低，并且分离源信号的精度在一定程度上受混合信号重构精度的影响。
【发明内容】
[0022]本发明主要针对稀疏源信号中的充分稀疏源信号的估计，为了解决现有技术中，由于没有考虑混合信号压缩观测值的特点，而采用通用的重构算法重构源信号而导致源信号精度低的问题，而提出了一种分布式压缩感知中充分稀疏源信号的交替迭代估计方法。
[0023]分布式压缩感知中稀疏源信号的交替迭代估计方法，所述交替迭代估计方法包括如下步骤:
[0024]步骤一、获取观测信号Y并转换为向量形式:通过采集混合信号Xi得到观测信号yi，其中i e {1，2，..，P}，P表示混合信号的通道数，对各通道进行独立观测，观测矩阵为Φ，观测过程71表示为式(7)，
[0025]Yi=Oxi (7)
[0026]其中，Xi的长度为N，Yi的长度为Μ，Φ是MXN维的矩阵，且M〈〈N，即Xi e Rnxi,Φ e Rmxn, Yi e Rmxi ；
[0027]将充分稀疏源信号S、混合信号X以及观测信号Y转变为矩阵形式，如式(8)所示，S = [S1 S2...SjJeRnxj
X = [X1 X2 …xf]eR'VxP (8)
Y-- fm/ mr?/Ir O ^
一 LyI y，.** y j 01?[0028]X=SAt
[0029]则Υ=ΦΧ (9)
[0030]其中，A表示PX J维的混合系统，式(9)可以表示为，
[0031]Y= Φ SAt (10)
[0032]将观测信号Y和源信号S变形为列向量的形式，如式(11)所示，
【权利要求】
1.分布式压缩感知中充分稀疏源信号的交替迭代估计方法，其特征在于:所述交替迭代估计方法包括如下步骤: 步骤一、获取观测信号Y并转换为向量形式:通过采集混合信号Xi得到观测信号yi，其中i e {1，2，..，P}，P表示混合信号的通道数，对各通道进行独立观测，观测矩阵为Φ，观测过程Yi表示为式(7)，
Υ?=Φχ? (7) 其中，Xi的长度为N，Yi的长度为Μ，Φ是MXN维的矩阵，且Μ〈〈Ν，即Xi e Rnxi,Φ e Rmxn, Yi e Rmxi ； 将充分稀疏源信号S、混合信号X以及观测信号Y转变为矩阵形式，如式(8)所示，
2.根据权利要求1所述的分布式压缩感知中充分稀疏源信号的交替迭代估计方法，其特征在于:步骤一中所述的观测矩阵Φ服从高斯分布。
3.根据权利要求1所述的分布式压缩感知中充分稀疏源信号的交替迭代估计方法，其特征在于:在迭代次数小于20时，步骤三中所述的迭代步长λ取0.01 ;在迭代次数大于20时，步骤三中所述的迭代步长λ取0.001。
【文档编号】H03M7/30GK103929186SQ201410155538
【公开日】2014年7月16日 申请日期:2014年4月17日 优先权日:2014年4月17日
【发明者】徐红伟, 付宁, 殷聪如, 张毅刚, 彭喜元 申请人:哈尔滨工业大学