技术报告关于FFT分析器4
3.傅里叶级数和傅里叶变换
3.4离散傅立叶变换
FFT分析器将时间x (t) 波形视为采样时间间隔h内的离散数据。设x (n) (n=0, ..... N-1)为在采样时间间隔h采样的有限数据的N个采样值序列,则此数据同样适用傅里叶变换。傅里叶展开由于“乘以时间波形x (n) 中的波形e-j2πft,对其1个周期进行积分,就可以求出某个波形的振幅。”,所以参考公式3-9b将该离散数据应用到公式中。;
因此,第k次谐波的傅立叶变换方程为;
(公式3-16)
,公式3-16称为离散傅里叶变换 (DFT) 。此外,与公式3-7b相同,公式3-11的关系公式也成立。类似地,离散傅里叶逆变换 (IDFT) :
(公式3-17)
的双曲余切值。DFT的公式3-16在概念上如下图3-5所示,请参考。
3.5抽样定理
那么,谐波k可以计算到什么程度?
要获取原始波形,必须以至少两倍于波形频率的频率进行采样。这称为抽样定理。现在考虑周期T的波形中采样数为N个的情况。此时,采样时间h及其频率f s;
(公式3-18)
因此,采样定理可以分析的频率f m;
(公式3-19)
即,直到第N/2次谐波。
用f s对限制为0到f m的波形进行采样并计算DFT,得到频率为-∞到+∞的频谱。图3-6给出了频率范围从-f m到4f m的频谱示意图。
-
图3-6
正如您在此图中看到的那样,得到的频谱重复相同的值,就像信件是f m的整数倍一样。另外,即使对从f m〜2f m受到频带限制的波形使用相同的f s进行采样并进行DFT,也同样如图3-6所示,在计算上求得实际上不存在的0〜f m的分量。这称为区域现象 (折返失真),此时的f m称为奈奎斯特频率。只观察得到的光谱的0〜f m时,无法判断原波形是0〜f m还是f m〜2f m。因此,在FFT分析仪中,通过预先限制为f m的低通滤波器 (称为区域滤波器),确定小于f m的信号,以满足采样定理。此外,由于需要考虑f m +α的频率以获得该低通滤波器充分衰减的频率,因此FFT分析器仅显示高达f m的1/1.28的光谱 (图3-7) 。
例如,如果FFT频率范围设置为1 kHz,则FFT分析器将以频率fs=2.56 kHz的频率采样,该频率是频率范围的2.56倍。
-
图3-7
图3-8示出了从时间波形考虑的采样定理,图3-9a和b示出了基于区域滤波器的ON/OFF的FFT分析器中的显示波形。
-
图3-8
● (红圈) :采样点
--:采样点再现的波形
如果原信号的频率小于采样频率fs的1/2 (=fm),则如 (a) 所示能够再现原信号,但如果原信号的频率大,则不能再现,如 (b) 的虚线波形那样产生乱频。另外,当原信号的频率与采样频率fm一致时,波形如虚线 (c) 所示。
-
图3-9a
由于过热过滤效果,未发生过热现象。然而,由于高频分量被面积滤波器切断,因此时间波形的矩形波受到波形的影响,并且不是完整的矩形波。
-
图3-9b
由于没有施加区域化滤波器,所以时间波形是原来的矩形波,而FFT的结果是产生区域化滤波器,出现了原来没有的频率成分。
这个故事变得非常复杂和数学,但你能理解到目前为止的内容吗?
在这里,我想从FFT分析器的角度组织傅里叶级数,转换,DFT和采样,并加入一些形象的概念。我将引用Kenichi Shito的“数字信号处理简介”并解释它。
请参阅图3-10。“A”波形是无限重复的时间波形,其被限制在奈奎斯特频率fm以下,并且对于“A”的频谱,获得“B”,其以原点为中心分布在±fm处,并且在该范围之外为零。其次,“B”的光谱X (f) 在频率轴上以2fx的周期无限排列,设其为X (f) 。;
-
-
图3-10
fx=fm时,“C”的光谱波形
当 fx < fm 时,频谱波形为“D”。
当 fx > fm 时,频谱波形为“E”。
) 中被调用,将出现故障。这里,光谱波形“E”是光谱交迭和区域发生。此外,在光谱波形“C”和“D”中,在-f m到+f m的范围内,光谱X (f) 与x (t) 相同。从另一个角度考虑采样列x n,其中x (t) 以2f x采样。此时的x n光谱;
(公式3-20)
它是C D E的表达式。即,用2f x对受f x频带限制的连续波形x (t) 进行采样,根据其无限的采样序列{x n}的数据求频谱X (f),观察其-f x〜+f x的频率范围时,就是x (t) 的频谱X (f) 本身 (“C”“D”) 。
将X (f) 进行傅里叶逆变换得到x (t),因此每1/ (2f x) 的样本值序列有足够的信息来表示连续波形x (t) 。如果原始波形的频率高于采样频率,则由于区域划分,产生重叠的光谱,如“E”,这与X (f) 不同,因此x (t) 不能再现。
到目前为止,我们讨论的是傅里叶变换。接下来,如果“A”的波形在周期T中无限重复,例如“F”,则通过傅里叶扩展波形获得的频谱是在频率轴上以1/T的间隔排列的频谱。该波形的频率成分被限制在f x以下,形成如“G”的光谱,该线谱的包络线为“B”的连续光谱。将“G”的光谱与“C”“D”同样考虑以2fx的周期无限重复的光谱“H”,将“H”进行逆傅立叶变换后的时间波形如“I”所示,成为周期T无限重复的样本值系列。换句话说,时间波形和光谱都是无限重复的周期性采样序列。此关系由公式3-16和3-17中的DFT和IDFT表示。公式3-16、3-17中的;
与Sin和Cos波形一样,是在某个周期内重复取相同值的周期函数,因此您可以理解DFT和IDFT在图中表示“H”和“I”。
另外,如果正确求出0~f x的光谱,我想您也可以理解作为信息是充分的。为了便于参考,这里给出了从光谱返回到连续波形的表达式。示例列{x n};
(公式3-21)
,其光谱X (f) 的-n项的傅里叶系数为x n。另外,根据样本值系列{x n},连续的原波形x (t) 可通过下式求出 (参照图3-11) 。
(公式3-22)
-
图3-11
那么,最后总结一下DFT,就完成本期的傅里叶级数、傅里叶变换的说明。
总结
如果在周期T中无限重复的信号x (t) 被限制在某一频率以下,则通过从根据采样定理以一定间隔h对周期T中的信号进行采样的N个离散采样列{x n}以每1/T=1/Nh执行DFT,可以获得离散频谱 (线频谱) X (f),并且可以通过IDFT再现离散采样列{x n}。
FFT分析器用时域的0到N-1个采样序列表示时间波形,用频域的0到f m /1.28显示频谱 (线谱) 。