技术报告关于FFT分析器5
3.傅里叶级数和傅里叶变换
3.6傅立叶展开计算示例
本节将计算实际时间波形中一个周期的傅里叶级数。
下图3-9所示为根据采样数N=12的1周期波形数据计算k=2的2次谐波的傅里叶系数的步骤,以及进行同样的计算至k=0~6 *所得的该波形的傅里叶级公式3-20。(※:根据抽样定理,在样本数N=12的情况下,可分析的谐波阶f m =6)
下面的图3-10显示了当采样数据是周期T=Nh=0.02秒 (ω 0 =50 Hz) 时FFT分析器显示的数据。
通过这样的手动计算,即使最多12个采样数据也需要花费大量时间。FFT分析器采用将N取2 n的FFT方法(例如1024=2 10),可在短时间内完成上述计算,并将数据保存到内存中,通过显示时间轴波形和频谱等波形,将输入波形的各种特征传达给我们。
-
图3-9
傅立叶系数a 0、a k、b k为:
二次谐波计算表
| n | Xi (角度) | Yi (数据序列) | Cos (2xi) | YiCos (2xi) |
| 0 | 0 | 0 | - | - |
| 1 | 30 | 16 | Cos(60)=0.5 | 8 |
| 2 | 60 | 26 | Cos(120)=-0.5 | -13 |
| 3 | 90 | 31 | Cos(180)=-1 | -31 |
| 4 | 120 | 36 | Cos(240)=-0.5 | -18 |
| 5 | 150 | 46 | Cos(300)=0.5 | 23 |
| 6 | 180 | 21 | Cos(360)=1 | 21 |
| 7 | 210 | 10 | Cos(420)=0.5 | 5 |
| 8 | 240 | -4 | Cos(480)=-0.5 | 2 |
| 9 | 270 | -19 | Cos(540)=-1 | 19 |
| 10 | 300 | -16 | Cos(600)=-0.5 | 8 |
| 11 | 330 | -14 | Cos(660)=0.5 | -7 |
| 12 | 360 | 0 | Cos(720)=1 | 0 |
| Σ | - | 133 | - | 17 |
根据上表,k=2时的傅立叶系数;
a 0 =133/12=11.08, a 2 =2/N∑x (n) cos (2n) =2.83
) 中被调用,将出现故障。在k=0~6之前进行同样的计算,可以求出下式。
(公式3-20)
| 频率 (Hz) |
Ak (Cos项) |
Bk (Sin项) |
r (振幅) |
φ (相位) |
r (振幅对数) |
|
| A0 | 0 | 11.08 | 0.00 | 11.08 | 0.00 | 20.89 |
| ω0 | 50 | -13.13 | 25.67 | 28.83 | -62.91 | 29.20 |
| 2ω0 | 100 | 2.83 | -0.58 | 2.89 | -11.58 | 9.21 |
| 3ω0 | 150 | 0.17 | 2.67 | 2.68 | 86.36 | 8.55 |
| 4ω0 | 200 | -2.83 | -1.15 | 3.05 | 22.11 | 9.70 |
| 5ω0 | 250 | 2.46 | 2.00 | 3.17 | 39.11 | 10.02 |
| 6ω0 | 300 | -0.59 | 0.00 | 0.59 | 0.00 | -4.66 |
图3-10
功率谱
相位
振幅
傅里叶光谱
4.查看FFT分析器屏幕
4.1时间轴分解能和频率分解能
在使用FFT分析器实际分析简单波形的同时,让我们来看看与前面解释的FFT和FFT特征的关系。
下面的数据1是分析1个原始波形的结果,其特征在10 kHz的频率范围内未知。相对于原波形,采样时间较短,因此用时间轴表示时,可以细致观察到原波形,而频率轴的频率分辨率不够,因此出现频谱相互连接的状态 (包络线) 。
-
<数据1>
此时频域的X轴与频率范围的关系可计算如下。
<时间轴 (顶部数据) >
以时间轴分辨率ΔT=约0.039 ms为间隔对原波形进行采样的2048个点的数据显示80 ms,其为基频T。
<频率轴 (底部数据) >
DFT时间波形显示从0 Hz到10 kHz的800光谱数据每12.5 Hz。
12.5 Hz是基频80 ms的倒数 (1/T) 。下表总结了样本数,数据长度和频率分辨率之间的关系,供您参考。
| 抽样点数 | 数据长度T (秒) | 频率分辨率1/T |
| 256 | ||
| 512 | ||
| 1024 | ||
| 2048 | ||
| 4096 |
<f:频率范围>
接下来请看数据2。这是在1kHz频率范围内分析与数据1相同波形的结果。虽然时间轴分辨率变差,但是*,相反频率分辨率提高到1.25 Hz,因此能够详细读取基波10 Hz的奇数倍的谐波成分。结果显示原始波形是10 Hz的矩形波。(*注:与信号频率相比,采样频率足够高,因此时间轴波形没有明显恶化。)
-
<数据2>
让我们看一下数据3,在相同的数据波形中,采样点数为1/4=512个点,频率范围为1 kHz,类似于数据2。对于数据2,数据长度从800 ms减少到200 ms,因此频率分辨率变为5 Hz,即1 kHz的1/200,并且出现不确定性。
-
<数据3>
虽然我们在数据1~3的每个频谱波形 (底部) 中显示了50 Hz的频谱值,但您可以看到它们几乎相同。这表明DFT频谱的包络线与原始波形的频谱一致,如前一期图3-10中的G所述。
-
图4-1