function out=chgouds(a)
b=a;
%提取點蝕分解小波系數a的高頻段
a(1:64)=0;
a(65:128)=0;
a(129:256)=0;
a(257:512)=0;
a(513:1024)=0;
a(1025:2048)=0;
a(2049:4096)=fft(b(2049:4096));
a(4097:8192)=0;
a(8193:16384)=0;
c=zeros(1,16383);
%對提取出的頻段做復共軛
c(2048:4095)=fft(conj(b(2049:4096)));

a1=zeros(1,32768);
a1(1:16384)=a;
a1(32768:-1:16386)=c;
out=real(ifft(a1));
plot(out);%畫重構系數圖
output=out;