基于B样条小波变换的matlab仿真【包括程序，注释，参考文献，操作步骤】

clc; clear; close all; warning off; addpath(genpath(pwd)); m=6; %求样条函数在整数点的值 N=func_nm(2*m,2*m+1) ; w=0:0.1*pi:10*pi %求尺度函数; mea_w=func_w(m,w,N); figure(1); subplot(1,3,1),plot(w/pi,abs(mea_w)); title('尺度函数的幅频图'); %计算相位角 mea_ang=atan(imag(mea_w)./(real(mea_w)+eps)); subplot(1,3,2),plot(w/pi,mea_ang); title('尺度函数的相频图'); %求反傅立叶变换，求出尺度函数的时域函数。 mea_x=ifft(mea_w); %对时域函数进行做一些搬移 mea_x=ifftshift(mea_x); subplot(1,3,3),plot(w/pi,real(mea_x)); title('尺度函数的时域图'); mea_2w=func_w(m,2*w,N); %求H(w), Hw=mea_2w./(mea_w+eps); %求h（k） H_k=ifft(Hw(1:21)) ; H_kreal=real(H_k) wdiv2=w/2; %求H(w+pi) Hwpidiv2=func_w(m,2*(wdiv2+pi),N)./func_w(m,wdiv2+pi,N); %求H(w/2+pi)的共轭 Hwpidiv2_=real(Hwpidiv2)-i*imag(Hwpidiv2); %求G(w/2) G_wdiv2=exp(-j*wdiv2).*Hwpidiv2_; %求尺度函数 mea_wdiv2=func_w(m,wdiv2,N); %求小波函数 bwav_w=G_wdiv2.*mea_wdiv2; figure(2); subplot(1,3,1),plot(w/pi,abs(bwav_w)); title('小波函数的幅频图'); %计算相位角 bwav_ang=atan(imag(bwav_w)./(real(bwav_w)+eps)); subplot(1,3,2),plot(w/pi,bwav_ang); title('小波函数的相频图'); %求反傅立叶变换，求出小波函数的时域函数。 bwav_x=ifft(bwav_w); bwav_x=ifftshift(bwav_x); subplot(1,3,3),plot(w/pi,real(bwav_x)); title('小波函数的时域图');