[发明专利]分布式压缩感知中充分稀疏源信号的交替迭代估计方法

摘要

一种分布式压缩感知中充分稀疏源信号的交替迭代估计方法，具体涉及由混合信号的压缩观测值，重构源信号的方法，解决现有技术中，由于没有考虑混合信号压缩观测值的特点，采用先恢复混合信号再分离源信号的算法重构源信号而导致源信号精度低的问题。本发明所述的方法在分布式压缩感知框架下，利用充分稀疏源信号的特点，采用交替估计的方法，从混合信号的观测值中估计源信号，每次迭代过程中，在不对混合信号进行重构的前提下，直接恢复得到源信号的估计值。本发明方法在信号压缩域从混合信号的观测值中分离源信号，降低了分离过程的复杂度。

主权项

分布式压缩感知中充分稀疏源信号的交替迭代估计方法，其特征在于：所述交替迭代估计方法包括如下步骤：步骤一、获取观测信号Y并转换为向量形式：通过采集混合信号xi得到观测信号yi，其中i∈{1,2,..,P}，P表示混合信号的通道数，对各通道进行独立观测，观测矩阵为Φ，观测信号yi表示为式(7)，yi＝Φxi                     (7)其中，xi的长度为N，yi的长度为M，Φ是M×N维的矩阵，且M＜＜N，即xi∈RN×1，Φ∈RM×N，yi∈RM×1；将充分稀疏源信号S、混合信号X以及观测信号Y转变为矩阵形式，如式(8)所示，S=S1S2...SJ∈RN×JX=x1x2...xP∈RN×PY=y1y2...yP∈RM×P---(8)]]>则X＝SATY＝ΦX                      (9)其中，A表示P×J维的混合系统，式(9)表示为：Y＝ΦSAT                    (10)将观测信号Y和充分稀疏源信号S变形为列向量的形式，如式(11)所示，S‾=S1TS2T...SJTT=s1,1s1,2...s1,N.......sJ,1sJ,2...sJ,NT∈RNJ×1Y‾=y1Ty2T...yPTT=y1,1y1,2...y1,M.......yP,1yP,2...yP,MT∈RMP×1---(11)]]>令ai,j,i∈{1,2,...,P},j∈{1,2,...,J}表示中的元素，为混合系统的估计值，则式(10)变形为：Y‾=a1,1Φa1,2Φ...a1,J-1Φa1,JΦa2,1Φa2,2Φ...a2,J-1Φa2,JΦ...............aP-1,1ΦaP-1,2Φ...aP-1,J-1ΦaP-1,JΦaP,1ΦaP,2Φ...aP,J-1ΦaP,JΦS‾---(12)]]>步骤二、迭代估计条件选择：随机选取一个混合矩阵作为混合系统的估计值令源信号的估计值令迭代次数变量t＝1；步骤三、利用得到的混合系统的估计值估计充分稀疏源信号S：式(12)采用现有的压缩感知重构算法求解，令Φ1=a1,1Φa1,2Φ...a1,J-1Φa1,JΦa2,1Φa2,2Φ...a2,J-1Φa2,JΦ...............aP-1,1ΦaP-1,2Φ...aP-1,J-1ΦaP-1,JΦaP,1ΦaP,2Φ...aP,J-1ΦaP,JΦ---(13)]]>则，采用线性规划的方法求解，如式(14)所示，argmin|S‾|1s.t.Y‾=Φ1S‾---(14)]]>其中，Φ1为变形后的观测矩阵Φ，将变形为N×J维的矩阵，则充分稀疏源信号的估计值即为变形后的矩阵；步骤四、更新混合系统的估计值ΔA=-A1(BS^T+I)]]>A^=A1+λΔA---(15)]]>其中，A1为上一次迭代得到的混合系统的估计值，λ为迭代步长，I为单位矩阵，为第j个源信号的先验概率分布；充分稀疏源信号服从超高斯分布，则其概率密度函数为P(Sm)∝αexp(α|Sm|)；步骤五、归一化混合系统的估计值的列向量；步骤六、如果相邻两次的混合系统的估计值相差很小，则停止迭代，进行步骤七；否则，t＝t+1，且返回步骤二；步骤七、输出混合系统的估计值和充分稀疏源信号的估计值