chgoudslow.m

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

function out=chgoudslow(a)
b=a;
%提取点蚀小波分解系数a的低频段
a(1:64)=0;
a(65:128)=fft(b(65:128));
a(129:256)=fft(b(129:256));
a(257:512)=fft(b(257:512));
a(513:2048)=0;
a(2049:4096)=0;
a(4097:8192)=0;
a(8193:16384)=0;
%对提取的频段做复共轭
c=zeros(1,16383);
c(64:127)=fft(conj(b(65:128)));
c(128:255)=fft(conj(b(129:256)));
c(256:511)=fft(conj(b(257:512)));

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

a1(32768:-1:16386)=c;
out=real(ifft(a1));
plot(out);
outputl=out;