www.pudn.com > bfunction.rar > wave.m

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) 
  t=real(ifftshift(h)); 
  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('小波函数');