|
楼主 |
发表于 2008-2-3 20:55
|
显示全部楼层
本帖最后由 牛小贱 于 2014-7-3 20:46 编辑
我构造了一个信号,M文件如下:
- clc;
- clear;
- fs=5000;
- nfft=1024;
- n=1:nfft/4; %(1,2,3,…,256)
- f=0:fs/nfft:fs/2-fs/nfft; %频率向量-横坐标
- i=0;
- for k=1:6
- for t=0:1/fs:1
- i=i+1;
- x(1,i)=4*sin(2*pi*25*k*t)*sin(2*pi*350*t)+2.5*sin(2*pi*350*t);
- end
- end
- subplot(3,1,1)
- y=abs(fft(x,nfft)); %FFT
- plot(f(n),y(n))
- xlabel('Frequency');
- ylabel('|FFT)|')
- grid on
- subplot(3,1,2)
- w=hanning(nfft);
- z=psd(x,nfft,fs,w,nfft/2); %PSD
- plot(f(n),abs(z(n)));
- xlabel('Frequency(Hz)');
- ylabel('|psd()|')
- grid on
- subplot(3,1,3)
- t=1000*(0:1/fs:(nfft-1)/fs); %倒频谱的时间向量-横坐标ms
- zw=rceps(x); %rceps
- plot(t(n),abs(zw(n)));
- xlabel('ms')
- ylabel('|Cepstrum|');
- grid on;
复制代码 得到如下图:图三中1000/40=25Hz
|
|