clc;
clear;
m=3;
w=0.001:0.3:10*pi;
x=w./5;
f=scale(m,w); 
figure(1); subplot(2,2,1);
 plot(x,abs(f));
 axis([0 10  0  1]);
  title('尺度函數幅頻圖');
 xlabel('w');
 ylabel('幅度');
 
 f_seq=ifft(f);
 y=real(ifftshift(f_seq));
 axis([0 10  0  1]);
 subplot(2,2,2);
 
 plot(x,y);
 title('尺度函數');
 xlabel('x');
 ylabel('尺度函數');
  s1=scale(m,w./2);
  s2=scale(m,w);
  H=s1./s2;
  h=ifft(H)
  h=real
  G=exp(-j*w./2).*H;
  lw=G.*scale(m,w./2);
   subplot(2,2,3);
 plot(x,abs(lw));
 title('小波函數幅頻圖');
 xlabel('w');
 ylabel('幅度');
  lw_seq=ifft(lw);
 y=real(ifftshift(lw_seq));
 subplot(2,2,4);
 plot(x,y);
 title('小波函數');
 xlabel('x');
 ylabel('小波函數');