%----------------------------------------------------------------------
% exa100701.m,  例10.7.1 和圖 10.7.5
%  用于信號的多分辨率分解和重建。調用函數子程序wavelet_decwavelet_rec
%  使用  A'trous 算法；在分解的過程中，同時給出各級的細節和概貌，它們
%  和原數據有著同樣的長度。
%----------------------------------------------------------------------
clf

f1=1;f2=20;f3=40;fs=200;N=400;
omiga=2*pi/fs;
t=0:N-1;
s1=sin(f1*omiga*t);
s2=sin(f2*omiga*t);
s3=sin(f3*omiga*t);
s=s1+s2+s3;
figure(1)
subplot(521);plot(s);Ylabel('x(t)')
subplot(522);plot(s);Ylabel('x(t)')

wname='sym8';
deccoef=wavelet_dec(s,4,wname);

sig = wavelet_rec(deccoef,4,wname)
figure(2)
subplot(211)
plot(s)
subplot(212)
plot(sig*sqrt(2))