压缩感知离散余弦变换基-小波基MATLAB程序

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