%基于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