logfft.m

小波阈值去噪 VC编译通过,效果很好,有很大参考价值!

%%%傅立叶变换对数形式%%
clear
clc
load signal
load syn
load noise
detat = 0.005;                                                    %采样间隔
n = length(signal);
offt=fft(signal,n);                                                 %输入数据进行傅立叶变换
offt=abs(offt)/n/detat;
sfft=fft(syn,n);
sfft=abs(sfft)/n/detat;
nofft=fft(noise,n);                                             %输入数据进行傅立叶变换
nofft=abs(nofft)/n/detat;
f = (0:n-1)/(n*detat);                                           %频率序列
figure
loglog(f(1:n/2),offt(1:n/2),'r');
hold on
loglog(f(1:n/2),sfft(1:n/2),'g')
hold on
loglog(f(1:n/2),nofft(1:n/2),'k')
legend('原始信号','合成信号','噪声','2');
title('重建小波以及原信号的FFT对比')
ylabel('傅里叶振幅(cm/s)');xlabel('频率(Hz)');
figure
subplot(211);
hold on
plot(signal);
title('原始信号');
subplot(212);
hold on;
plot(syn);
title('重建信号');