|
回复8楼kleeblatt 的贴子
本帖最后由 wdhd 于 2016-9-19 13:38 编辑
你的问题是:若1秒内得128个采样,如何作apFFT?
选定fft阶数N后,当采样频率fs=N赫/秒,作apFFT需2秒采样, 当采样频率fs=N/2赫/秒,作apFFT只需1秒采样,
所以若得1秒内128个采样(即采样频率为128赫/秒),相当采样频率为64赫/秒2秒内得128个采样. 计算频率时注意采样频率为64赫/秒.
若只得128个采样,你用127采样作N=64的apFFT即可,反正作N阶FFT/apFFT须2*N-1个采样.如何作apFFT 5楼程序中都写出了
类同的有一种FFT时移相位差法,它须二组相邻序列作FFT才能校正振幅等,所以也要采样2秒
由信号频谱求振幅的一个长处是一次FFT可得多种频率式成份的振幅值. 用fft/apfft校正法得的振幅校正谱可用下aaaa.doc附图示
N=256;
t=-N+1:N-1;
A=1.0;f=29.2;phi=40;
y=A*cos(2*pi*t*f/N+phi*pi/180)+2*A*cos(2*pi*t*2*f/N+phi*pi/180)+3*A*cos(2*pi*t*3*f/N+phi*pi/180);
附图是3个成分信号y的振幅谱校正图;(a)为fft振幅谱a1,(b)为apfft振幅谱a2,(c)校正振幅谱.
3个成分信号振幅校正为1,2,3.而且在峰值附近形成一个平台,很醒目. 有噪时峰值处最准.
振幅校正值(111db)
0.9999 1.9999 3.0000
振幅校正值(10db)
1.0536 1.9113 3.0705
从图也可见,未经校正的(a)fft振幅谱a1和(b)apfft振幅谱a2峰值都不是1,2,3. 只有校正后(C)的峰值才是1,2,3.
aaaa.doc
(30.5 KB, 下载次数: 5)
[ 本帖最后由 zhwang554 于 2009-11-17 23:33 编辑 ]
|
|