[发明专利]一种测量值缺失的压缩感知重构方法

权利要求书

1.一种测量值缺失的压缩感知重构方法，其特征在于，包括如下步骤：

S1，获取稀疏度为k的N维源信号x；

S2，使用M×N维测量矩阵Φ对所述源信号x进行线性投影，获得x的低维测量向量y＝[y(1),y(2),...,y(M)]^T，令f＝[f(1),f(2),...,f(M)]^T表示测量值缺失部分，测量向量为：y＝Φx+f，

其中，M、N为正整数，f为缺失值向量，其非零元素位置为缺失测量位置，其非零元素值为缺失测量值的相反数；

S3，利用测量矩阵构造恢复矩阵[Φ,I]，设定迭代次数初始值i＝0，重构信号初始值u_i＝0，残差初始值d_i＝y，本次迭代残差值：d_i＝y-[Φ,I]u_i；

S4，计算恢复矩阵的转置[Φ,I]^T，按u_i支撑集抽取对应行，得到其组成的子矩阵计算本次迭代步长：

μi=||[Φ,I]ΓiT(y-[Φ,I]ui)||22||[Φ,I]Γi[Φ,I]T(y-[Φ,I]ui)||22]]>

S5，计算b_i＝u_i+μ_i[Φ,I]^Td_i；

S6，保留b_i中幅值最大的k个元素，将其他元素置零，令u_i+1＝b_i；

S7，计算残差值d_i+1＝y-[Φ,I]u_i+1，并判断是否满足终止条件d_i+1≥d_i，若满足，执行步骤S8；若不满足，令i＝i+1，返回步骤S2；

S8，依次选取输出u_i中前N个元素，组成的向量赋值给信号重构值依次选取输出u_i中第N+1个到第N+M个元素，组成的向量赋值给恢复缺失测量值

2.如权利要求1所述的测量值缺失的压缩感知重构方法，其特征在于，所述稀疏度k小于p，测量矩阵行数大于q，所述p、q为正整数。