clear?all;
clc;
signal0=wavread（‘M1_263‘）;
L=160;?%窗長
P=floor（length（signal0）/L）;?%幀數
%?dct_coef=zeros（NP）;
signal1=[];
n=L*0.5;
Measurement=[];
for?ii=1:P
%?ii=10;
x0（ii:）=signal0（（ii-1）*L+1:ii*L）;
dct_coef0（ii:）=dct（x0（ii:））;
%求dctbasis和idctbasis
for?kk=2:L
????for?nn=1:L
????????dctbasis（kknn）=（2/L）^0.5*cos（（2*（nn-1）+1）*（kk-1）*pi/2/L）;
????end
end
for?nn=1:L
????dctbasis（1nn）=（1/L）^0.5*cos（（2*（nn-1）+1）*（1-1）*pi/2/L）;
end
idctbasis=dctbasis‘;
%%CS重構
N（ii）=L*0.5;
n=N（ii）;
guancematrix=randn（nL）;??
measurement=guancematrix*x0（ii:）‘;
CSmatrix=guancematrix*idctbasis;
m=L;
c=ones（2*m1