没想到1kHz正弦波做时频分析竟然能出这种花里胡哨pattern(每1ms出两个,正好为信号频率的2倍)。 仔细一看,原来是分析参数设置不正确,主要是短时傅里叶变换的每段长度及频率分辨率没设置对。时频分析的方法之一是采用短时傅里叶变换,就是将采集到的整段数据分为多段,每段之间可部分重叠,然后对每段做FFT(可加窗)得到频谱,最后得到频谱随时间的变化图,应用包括语音识别等。
LZ图中的采样频率为800kHz, 短时傅里叶变换的每段长度为64点,才80us,还远远不到1KHZ信号的一个周期长度。频率分辨率为800000/64=12.5kHz,也远远不能分辨1kHz。
下面采用相近的参数复制一下LZ的这个pattern,采样频率为784kHz,先按错误的参数设置,每段38个点,约50us,无重叠,平头窗, FFT点数1024。
声谱图:
瀑布图:
下面按比较正确的参数设置,每段7680个点,10ms,无重叠,平头窗,FFT点数8192。
瀑布图:
|