%裝載語音信號
N=1024;
s=wavread('m?l.wav',N);
figure(1);
plot(1:N,s,'LineWidth',2);
xlabel('時間 n');
ylabel('幅值 A');

%選用小波函數
wavelet='db5';
%分解級數
level=5;
%壓縮參數
alpha=1.5;
%閾值類型
sorh='h';
%小波分解
[c,l]=wavedec(s,level,wavelet);
%計算閾值
[thr,nkeep]=wdcbm(c,l,alpha);
%使用硬閾值壓縮信號
[xc,cxc,lxc,perf0,perfl2]=wdencmp('lvd',c,l,wavelet,level,thr,sorh);
%對系數編碼
cC=encode(cxc);
%傳輸的系數
tC=cC;
%傳輸的系數向量長度
tL=length(tC);
%零系數的百分比
pzeros=perf0
%保留的能量大小
pnormen=perfl2
%編碼的壓縮比
compratio=length(s)/length(tC)