N=1024;
r=4;
h=r/N;
x=0:h:r-h;

s1=bspline(x,1);
s2=bspline(x,2);
s3=bspline(x,3);
s4=bspline(x,4);

figure(1);
plot(x,s1,x,s2,x,s3,x,s4)

dwtmode('per');
[c1,d1]=dwt(s1,'haar');
[c2,d2]=dwt(s2,'haar');
[c3,d3]=dwt(s3,'haar');
[c4,d4]=dwt(s4,'haar');

x=x(1:2:end);

figure(2);
subplot(2,2,1);
plot(x,d1);
subplot(2,2,2);
plot(x,d2);
subplot(2,2,3);
plot(x,d3);
subplot(2,2,4);
plot(x,d4);

[c1,d1]=dwt(s1,'db3');
[c2,d2]=dwt(s2,'db3');
[c3,d3]=dwt(s3,'db3');
[c4,d4]=dwt(s4,'db3');

figure(3);
subplot(2,2,1);
plot(x,d1);
subplot(2,2,2);
plot(x,d2);
subplot(2,2,3);
plot(x,d3);
subplot(2,2,4);
plot(x,d4);


