跳至正文

Select your region & language

Global

Region

技术报告关于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