chgoulw.m

来自「本源代码主要实现谐波小波的分解及重构的matlab程序。」· M 代码 · 共 18 行

function out=chgoulw(a)
b=a;
%提取裂纹分解小波系数a的高频段
a(1:512)=0;
a(513:1024)=0;
a(1025:2048)=fft(b(1025:2048));
a(2049:32768)=0;
%对提取出的频段做复共轭
c=zeros(1,32767);
c(1024:2047)=fft(conj(b(1025:2048)));

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

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