clear
N?=?2^10;
M?=?N/8;
K?=?M/4;
t?=?1:1:N;
fs?=?1;
x?=?sin（t/fs）+2*cos（3*t/fs+2.4）+3*sin（7*t/fs+1）;??????????%若該行作為采樣信號，則所輸入的信號為非嚴格的稀疏信號
%x=0.3*sin（2*pi*50*t/800）+0.6*sin（2*pi*100*t/800）+0.1*sin（2*pi*200*t/800）+0.4*sin（2*pi*25*t/800）+0.7*sin（2*pi*150*t/800）+0.1*sin（2*pi*250*t/800）+0.7*sin（2*pi*125*t/800）;
%若該行作為采樣信號，則所輸入的信號為嚴格的稀疏信號，稀疏度為14
Phi?=?rand（MN）*2;?????%生成隨機矩陣，之后通過for循環，將隨機矩陣變為0/1二值獨立隨機矩陣
x_out?=?[];
for?i?=?1:M
????for?j?=?1:N
????????if?Phi（ij）>=1
????????????Phi（ij）?=?1;
????????else
????????????Phi（ij）?=?0;
????????end
????end
end

Psi=fft（eye（NN））/sqrt（N）;?
T?=?Phi;

%Kasi?=?randn（NM）;

for?i?=?0:0
????xi?=?x（i*N+1:i*N+N）;
????figure（2）??????????????????%該圖顯示了輸入信號通過FFT變換后，稀疏信號的頻域圖
????plot（tfft（xi）‘k.-‘）
????yi?=?T?*?fft（xi）.‘;
????r?=?yi;
????A?=?[];
????y_out?=?zeros（1N）;
????now_kasi?=?T;