% fig27
clear;clc;close all;
N=512;
z=zeros(1,N);
z(1,[129:256])=sin((abs([128:255]-128).^1.7)/128);
z(1,[385:448])=sin((abs([384:447]-128).^2)/128);
plot(z);axis tight

N=4;
[c,l] = wavedec(z,N,'db6');
ca4 = appcoef(c,l,'db6',4);
cd4 = detcoef(c,l,4);
cd3 = detcoef(c,l,3);
cd2 = detcoef(c,l,2);
cd1 = detcoef(c,l,1);

figure;plot(cd1,'x');axis tight;
figure;plot(cd2,'x');axis tight;
figure;plot(cd3,'x');axis tight;
figure;plot(cd4,'x');axis tight;
figure;plot(ca4,'x');axis tight;


%figure;
%plot(C(1:L(1)),'x');

%for k=1:N,
%   figure;
%   x=plot(C(L(k)+1:L(k)+L(k+1)),'x');axis tight;
%end;