www.gusucode.com > 局部均值分解源代码 难得的matlab程序代码源码 > lmd/hua_fft1.m

    function hua_fft1(y,fs)
%该FFT是为了和hua_xihua.m这个画细化频率的函数的幅度值相吻合，来使用的，
%如果不画细化普，还是用hua_fft.m这个画FFT的结果
%当style=1,画幅值谱；当style=2,画功率谱;当style=其他的，那么花幅值谱和功率谱
nfft= 2^nextpow2(length(y));%找出大于y的个数的最大的2的指数值
y_ft=fft(y,nfft);%对y信号进行DFT，得到频率的幅值分布
y_p=y_ft.*conj(y_ft)/nfft;%conj()函数是求y函数的共轭复数，实数的共轭复数是他本身。
y_f=fs*(0:nfft/2-1)/nfft;%DFT变换后对应的频率的序列
y_p=y_ft.*conj(y_ft)/nfft;%conj()函数是求y函数的共轭复数，实数的共轭复数是他本身。
plot(y_f,abs(y_ft(1:nfft/2)));%论坛上画FFT的方法
%plot(y_f,2*abs(y_ft(1:nfft/2))/length(y));%matlab的帮助里画FFT的方法
end