wavelet transform program.rar_four level wavelet_wavelet transfo

clear; x=load('B.dat'); x=x(6000:7000); %N=length(x); N=1024; x=x(1:N); x=x'; [ld,hd]=wfilters('db4','d'); %%%%%%%%%%%%%%% di yi ji fen jie %%%%%%%%%%%%%%%%%%% ld1=[ld,zeros(1,N-length(ld))]; hd1=[hd,zeros(1,N-length(hd))]; t10=ifft(fft(x).*fft(ld1)); t11=ifft(fft(x).*fft(hd1)); t10=dyaddown(t10); t11=dyaddown(t11); %%%%%%%%%%%%%%%%%%% di er ji fen jie %%%%%%%%%%%%%%%%%%%%% ld2=[ld,zeros(1,N/2-length(ld))]; hd2=[hd,zeros(1,N/2-length(hd))]; t20=ifft(fft(t10).*fft(ld2)); t21=ifft(fft(t10).*fft(hd2)); t22=ifft(fft(t11).*fft(ld2)); t23=ifft(fft(t11).*fft(hd2)); t20=dyaddown(t20); t21=dyaddown(t21); t22=dyaddown(t22); t23=dyaddown(t23); %%%%%%%%%%%%%%%%%%%% di san ji fen jie %%%%%%%%%%%%%% ld3=[ld,zeros(1,N/4-length(ld))]; hd3=[hd,zeros(1,N/4-length(hd))]; t30=ifft(fft(t20).*fft(ld3)); t31=ifft(fft(t20).*fft(hd3)); t32=ifft(fft(t21).*fft(ld3)); t33=ifft(fft(t21).*fft(hd3)); t34=ifft(fft(t22).*fft(ld3)); t35=ifft(fft(t22).*fft(hd3)); t36=ifft(fft(t23).*fft(ld3)); t37=ifft(fft(t23).*fft(hd3)); t30=dyaddown(t30); t31=dyaddown(t31); t32=dyaddown(t32); t33=dyaddown(t33); t34=dyaddown(t34); t35=dyaddown(t35); t36=dyaddown(t36); t37=dyaddown(t37); figure(1); subplot(511) plot(real(t30)) subplot(512) plot(real(t31)) subplot(513) plot(real(t32)) subplot(514) plot(real(t33)) subplot(515) plot(real(t34)) %%%%%%%%%%%%%%%%%%%%%%%%% di si ji fen jie %%%%%%%%%%%%%%%%%%% ld4=[ld,zeros(1,N/8-length(ld))]; hd4=[hd,zeros(1,N/8-length(hd))]; t40=ifft(fft(t30).*fft(ld4)); t41=ifft(fft(t30).*fft(hd4)); t42=ifft(fft(t31).*fft(ld4)); t43=ifft(fft(t31).*fft(hd4)); t44=ifft(fft(t32).*fft(ld4)); t45=ifft(fft(t32).*fft(hd4)); t46=ifft(fft(t33).*fft(ld4)); t47=ifft(fft(t33).*fft(hd4)); t48=ifft(fft(t34).*fft(ld4)); t49=ifft(fft(t34).*fft(hd4)); t410=ifft(fft(t35).*fft(ld4)); t411=ifft(fft(t35).*fft(hd4)); t412=ifft(fft(t36).*fft(ld4)); t413=ifft(fft(t36).*fft(hd4)); t414=ifft(fft(t37).*fft(ld4)); t415=ifft(fft(t37).*fft(hd4)); t40=dyaddown(t40); t41=dyaddown(t41); t42=dyaddown(t42); t43=dyaddown(t43); t44=dyaddown(t44); t45=dyaddown(t45); t46=dyaddown(t46); t47=dyaddown(t47); t48=dyaddown(t48); t49=dyaddown(t49); t410=dyaddown(t410); t411=dyaddown(t411); t412=dyaddown(t412); t413=dyaddown(t413); t414=dyaddown(t414); t415=dyaddown(t415); %%%%%%%%%%%%%%%%%%%%%%%%%%%% di si ji chong gou %%%%%%%%%%%%%% t40=dyadup(t40); t41=dyadup(t41); t42=dyadup(t42); t43=dyadup(t43); t44=dyadup(t44); t45=dyadup(t45); t46=dyadup(t46); t47=dyadup(t47); t48=dyadup(t48); t49=dyadup(t49); t410=dyadup(t410); t411=dyadup(t411); t412=dyadup(t412); t413=dyadup(t413); t414=dyadup(t414); t415=dyadup(t415); t40=t40(1:length(t40)-1); t41=t41(1:length(t41)-1); t42=t42(1:length(t42)-1); t43=t43(1:length(t43)-1); t44=t44(1:length(t44)-1); t45=t45(1:length(t45)-1); t46=t46(1:length(t46)-1); t47=t47(1:length(t47)-1); t48=t48(1:length(t48)-1); t49=t49(1:length(t49)-1); t410=t410(1:length(t410)-1); t411=t411(1:length(t411)-1); t412=t412(1:length(t412)-1); t413=t413(1:length(t413)-1); t414=t414(1:length(t414)-1); t415=t415(1:length(t415)-1); lr4=ld4(end:-1:1);%为什么lr4是ld4的逆？？？？？ hr4=hd4(end:-1:1); lr4=circshift(lr4',1)';%为什么lr4还要做此变换？？？？ hr4=circshift(hr4',1)'; t40=ifft(fft(lr4).*fft(t40)); t41=ifft(fft(hr4).*fft(t41)); t30=real(t40+t41);%为什么t30只取t40和t41的实部？？？ t42=ifft(fft(lr4).*fft(t42)); t43=ifft(fft(hr4).*fft(t43)); t31=real(t42+t43); t44=ifft(fft(lr4).*fft(t44)); t45=ifft(fft(hr4).*fft(t45)); t32=real(t44+t45); t46=ifft(fft(lr4).*fft(t46)); t47=ifft(fft(hr4).*fft(t47)); t33=real(t46+t47); t48=ifft(fft(lr4).*fft(t48)); t49=ifft(fft(hr4).*fft(t49)); t34=real(t48+t49); t410=ifft(fft(lr4).*fft(t410)); t411=ifft(fft(hr4).*fft(t411)); t35=real(t410+t411); t412=ifft(fft(lr4).*fft(t412)); t413=ifft(fft(hr4).*fft(t413)); t36=real(t412+t413); t414=ifft(fft(lr4).*fft(t414)); t415=ifft(fft(hr4).*fft(t415)); t37=real(t414+t415); figure(2); subplot(511) plot(t30) subplot(512) plot(t31) subplot(513) plot(t32) subplot(514) plot(t33) subplot(515) plot(t34) %%%%%%%%%%%%%%%%%%% di san ji chong gou %%%%%%%%%%% t30=dyadup(t30); t31=dyadup(t31); t32=dyadup(t32); t33=dyadup(t33); t34=dyadup(t34); t35=dyadup(t35); t36=dyadup(t36); t37=dyadup(t37); t30=t30(1:length(t30)-1); t31=t31(1:length(t31)-1); t32=t32(1:length(t32)-1); t33=t33(1:length(t33)-1); t34=t34(1:length(t34)-1); t35=t35(1:length(t35)-1); t36=t36(1:length(t36)-1); t37=t37(1:length(t37)-1); lr3=ld3(end:-1:1); hr3=hd3(end:-1:1); lr3=circshift(lr3',1)'; hr3=circshift(hr3',1)'; t30=ifft(fft(lr3).*fft(t30)); t31=ifft(fft(hr3).*fft(t31)); t20=real(t30+t31); t32=ifft(fft(lr3).*fft(t32)); t33=ifft(fft(hr3).*fft(t33)); t21=real(t32+t33); t34=ifft(fft(lr3).*fft(t34)); t35=ifft(fft(hr3).*fft(t35)); t22=real(t34+t35); t36=ifft(fft(lr3).*fft(t36)); t37=ifft(fft(hr3).*fft(t37)); t23=real(t36+t37); %%%%%%%%%%%%%%%%%%%%%%% di er ji chong gou %%%%%%%%%%%%%%%%%% t20=dyadup(t20); t21=dyadup(t21); t22=dyadup(t22); t23=dyadup(t23); t20=t20(1:length(t20)-1); t21=t21(1:length(t21)-1); t22=t22(1:length(t22)-1); t23=t23(1:length(t23)-1); lr2=ld2(end:-1:1); hr2=hd2(end:-1:1); lr2=circshift(lr2',1)'; hr2=circshift(hr2',1)'; t20=ifft(fft(lr2).*fft(t20)); t21=ifft(fft(hr2).*fft(t21)); t10=real(t20+t21); t22=ifft(fft(lr2).*fft(t22)); t23=ifft(fft(hr2).*fft(t23)); t11=real(t22+t23); %%%%%%%%%%%%%%%%% di yi ji chong gou %%%%%%%%%%%%%% t10=dyadup(t10); t11=dyadup(t11); t10=t10(1:length(t10)-1); t11=t11(1:length(t11)-1); lr1=ld1(end:-1:1); hr1=hd1(end:-1:1); lr1=circshift(lr1',1)'; hr1=circshift(hr1',1)'; t10=ifft(fft(lr1).*fft(t10)); t11=ifft(fft(hr1).*fft(t11)); sig=t10+t11; sig=real(sig); figure(3); subplot(211) plot(x) subplot(212) plot(sig,'r') legend('original signal','synthesis signal')