MP算法实现

%基于MP算法
clc;clear
%观测向量y的长度M=80，即采样率M/N=0.3
N=256;
K=15;             %信号稀疏度为15
M=80;             %
x = zeros（N1）;
q = randperm（N）;
x（q（1:K）） =randn（K1）;    %原始信号

%构造高斯测量矩阵，用以随机采样
Phi = randn（MN）*sqrt（1/M）; 
for i = 1:N
    Phi（:i） = Phi（:i）/norm（Phi（:i））;
end
y=Phi*x ;       %获得线性测量

%用MP算法开始迭代重构
m=2*K;              %总的迭代次数
r_n=y;              %  残差值初始值
x_find=zeros（N1）;      %x_find为MP算法恢复的信号

for times