chgoulwlow.m

本源代码主要实现谐波小波的分解及重构的matlab程序。

function out=chgoulwlow(a)
b=a;
%提取裂纹小波分解系数a的低频段
a(1:32)=0;
a(33:64)=fft(b(33:64));
a(65:128)=fft(b(65:128));
a(129:256)=fft(b(129:256));
a(257:512)=0;
a(513:1024)=0;
a(1025:2048)=0;
a(2049:32768)=0;
%对提取的频段做复共轭
c=zeros(1,32767);
c(32:63)=fft(conj(b(33:64)));
c(64:127)=fft(conj(b(65:128)));
c(128:255)=fft(conj(b(129:256)));

a1=zeros(1,65536);
a1(1:32768)=a;

a1(65536:-1:32770)=c;
out=real(ifft(a1));
plot(out);
outputl=out;