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）;
A=[CSmatrix-CSmatrix];
b=measurement;
lb=zeros（2*m1）;
options=optimset（‘LargeScale‘‘off‘‘Simplex‘‘on‘）;
[zfvalexitflagoutput]=linprog（c[][]Ablb[][]options）;
%?output
u=z（1:m1）;
v=z（m+1:2*m1）;
dct_coef1（ii:）=u-v;
dcterror（ii）=norm（dct_coef0（ii:）-dct_coef1（ii:）1）;
x1（ii:）=idctbasis*dct_coef1（ii:）‘;
x_error（ii）=norm（x0（ii:）-x1（ii:）1）;
signal1=[signal1x1（ii:）];
end
signal0=signal0（1:P*320）;
snr=signal0‘*signal0/（（signal0‘-signal1）*（signal0-signal1‘））;
snr=10*log10（snr）;
signal_error=norm（signal0‘-signal11）;
ratio=sum（N）/（L*P）;
wavwrite（signal116000‘M1_263_dctcs_csmr05‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1263??2009-03-12?11:35??dct_cs.m

-----------?---------??----------?-----??----

?????????????????1263????????????????????1