马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
本帖最后由 牛小贱 于 2014-9-12 13:16 编辑
以下是我用matlab求的FFT频谱和HHT频谱,为何纵轴差这么远啊,这数量级差的太远了吧?
一下是我求边际谱FFT谱的程序:求傅立叶谱程序: - y2=fft(y);
- Fourier_Amp=sqrt(real(y2).^2+imag(y2).^2);
- fori=1:length(Fourier_Amp)/2
- Real_Fourier_Amp(i)=Fourier_Amp(i)+Fourier_Amp(length(Fourier_Amp)-i);%折半
- end
- Ts=2e-6;
- fs=1/Ts;
- f_Fourier=(0:length(Fourier_Amp)/2-1)/length(Fourier_Amp)*fs;
- plot(f_Fourier,Real_Fourier_Amp)
- title('FourierSpectrum');
- xlabel('Frequency/Hz');
- ylabel('amplitude')
复制代码求边际谱程序: >>imf=emd(y); >>[A,fa,tt]=hhspectrum(imf); >>[E,tt1]=toimage(A,fa,tt,length(tt)); >>for k=1:size(E,1) bip(k)=sum(E(k,:))*(2e-6); end >>Ts=2e-6; >>fs=1/Ts; >>f=(0:(length(y)-3))/(length(y))*(fs/2); >>plot(f,bip) >>titile('Marginal Spectrum of Ultrasonic Signal') >>xlabel('Frequency/Hz') >>ylabel('Amplitude')
真切期望得到大家的解答,谢谢!!!
|