原创 德思特干货｜如何使用SBench 6对数字化仪采集信号进行处理？（三）——快速傅立叶变换（FFT）

来源：德思特测量测试 德思特干货｜如何使用SBench 6对数字化仪采集信号进行处理？（三）——快速傅立叶变换（FFT）

原文链接：https://mp.weixin.qq.com/s/mYS1iDXFNVfReCGGtF78mw

欢迎关注虹科，为您提供最新资讯！

#信号采集 #信号处理 #测量测试

上一篇文章介绍了德思特SBench 6的平均运算功能。本章将继续为大家介绍SBench 6的快速傅立叶变换（FFT）。

前文回顾：https://mp.weixin.qq.com/s/j-iN_2Jrn9ZHGMaaAYsDJg

快速傅立叶变换（FFT）将获取的波形从时域（幅度对时间）映射到频谱域（幅度对频率）。这使得用户可以观察构成信号的频率成分。FFT本身并不能直接改善信号质量，但它能显示信号的频率成分结构，并提供了如何移除不需要的频谱成分的信息。

FFT产生的频谱有一个离散的时间轴，就像时域信号有离散的采样时间一样。频谱中的谱线，通常称为“bins”或“cells”，每个谱线之间由分辨率带宽（f）隔开，而分辨率带宽与采集到的原始信号长度成反比，因此，要增加FFT频谱的频率分辨率，必须增加采集的原始信号长度。而频谱显示的总频率范围，或者说频率跨度，是原始信号采样率的一半，因此，要增加总频率跨度，也必须增加采样率。

在德思特SBench 6软件中，FFT结果的纵轴缩放可以设为线性单位伏特，也可以是对数单位分贝（dB）。对数刻度可以参考数字化仪的满偏刻度（dBFS）、1 毫瓦（dBm）、1 微伏（dbμV）或设为针对假定载波频谱中的最大峰值（dBc）。

理论上的傅立叶变换假设输入记录的长度为无穷大。有限的采样长度则会在频域边缘引入不连续性，为频域引入伪频率，并一定程度扭曲实际的频谱。例如，当信号的开始和结束相位不同，或者信号频率落在两个相邻的谱线之间时，使频谱变宽。

频谱展宽，以至于扩展到许多相邻的频带，称为泄漏。对于这种问题的理论应对方法是确保在显示网格内包含整数个周期，或在边缘处不出现不连续性。两者都需要信号波形频率和数字化采样率之间非常精确的同步，并且需要准确地设置采集长度，这通常只可能在实验室中实现，而不能作用于真实世界的信号。另一种方法是使用窗函数（加权）来平滑信号的边缘。

为了尽量减少这些负面影响，对获取的信号应用加权函数，使记录的端点为零。德思特SBench 6软件中的FFT功能为用户提供八种加权函数的选择，而这些加权函数则可以改变谱线的形状。对此没什么概念的读者，可以把它想象为组合到一起的一些列并行的理想带通滤波器，这些滤波器之间的频率间距为分辨率带宽，而加权函数则会影响滤波器频率响应的形状。图1比较了四种最常用的加权函数对应的频率响应。

图1：四种最常用的加权函数的谱形状比较

理想情况下，主瓣应尽可能窄，且两侧平坦，以代表真实的频谱成分，而所有旁瓣应无限衰减。窗函数类型定义了在FFT处理中使用的等效滤波器的带宽和形状。表1中展示了频谱响应的最大旁瓣幅度。较大的旁瓣电平将有助于区分密集的频谱成分。

表1：不同加权函数的关键特性

如前所述，FFT结果的频率轴是离散的，具有以分辨率带宽的倍数间隔的频点。如果输入信号频率落在两个相邻的谱线之间，能量将被分配到两个谱线中，且峰值幅度将降低。这被称为栅栏效应或波浪状效应，而扩宽频谱响应可以减少幅度变化。表1中的栅栏损失列给出了不同加权函数该效应幅值。

同时，加权函数还会影响频谱响应的带宽。等效噪声带宽（ENBW）说明了相对于矩形窗加权带宽的相对变化。将功率谱归一化到测量带宽（功率谱密度）需要将功率谱除以ENBW与分辨率带宽的积（f✖️ENBW）。

相干增益描述了给定加权函数相对于矩形窗加权的频谱幅度变化。这是一个应用在所有频率上的固定增益，可以很容易地归一化。

矩形窗加权函数是采集信号没有任何加权的响应。它具有最窄的带宽，但旁瓣幅值则相当高。由于采集时域记录中的所有点的幅度响应都是均匀的，它常用于具有瞬态性质的信号（或相对记录总长短得多的信号）。当需要以最佳频率精度进行分析时，也会使用它。

Hanning和Hamming加权函数具有良好的通用的频率响应，能提供较好的频率分辨率以及合理的旁瓣响应。Blackman-Harris则旨在获得最佳幅度精度和优秀的旁瓣抑制。

图2展现了一个典型的FFT应用实例。在该实例中，我们使用宽带的仪器级麦克风和德思特TS-M4i系列14位数字化仪，获得了超声波测距仪的信号。

图2：40 kHz的超声波脉冲（左）及其FFT结果（右下全频段，右上局部放大）

采集到的时域信号在左窗格中。时域记录包含了在3.90625 MS/s采样率下采集到的16384个样本点，持续时间为4.2毫秒。得到的FFT结果（右窗格）共有8192条谱线，每条谱线之间间隔238 Hz分辨率带宽（记录时长的倒数），总频率跨度为1.95 MHz（采样率的一半）。其中，右下角的频谱是全频段总览图，右上角则是局部放大视图，只显示了前100 kHz频率范围的内容，以便更好地观察主要频谱成分。

FFT能帮助我们更好地理解构成这个信号的成分。首先我们看时域，可以知道这是一个持续时间小于采集记录长度的瞬态信号，在这种情况下我们选用了矩形窗进行加权。FFT结果中幅值最大的频率成分，明显就是我们的主要信号——40 kHz左右的脉冲。而在80 kHz附近频率的一个小峰，则是40 kHz信号成分的二次谐波。其幅度大约比40 kHz信号成分低45 dB。此外，在0~10 kHz之间也有很多低频干扰成分，其中最高的、接近DC的那部分，对应的是设备所在房间中的环境噪声。

在该实例中，我们的目标是能够测量发射脉冲和40 kHz回波之间的时延。为了实现这一点，我们需要改进这个测量过程。第一步，我们希望移除频谱中40 kHz成分之外的其它成分。现在我们将带着这个目标，来进行滤波器的设置。

在德思特SBench 6专业版软件中，我们提供了低通、带通或高通选项的有限冲激响应(FIR)数字滤波器。通过输入所需的滤波器类型、截止频率或频率，以及滤波器阶数，用户可以直接在图形界面中创建这些滤波器。SBench 6软件会在滤波器无法实现时给出提示，并提出解决建议。或者，您也可以输入从其他来源获取的滤波器系数。我们将这些滤波器应用于采集到的信号，然后将滤波结果与原始采样结果或平均采样结果进行比较。在图3中，应用的带通FIR滤波器截止频率为30和50 kHz，以帮助提取目标信号。

图3：原始波形和滤波后信号及其FFT结果的对比

左上角的窗格显示了原始波形，下面对应我们之前看到的原始信号FFT结果。右上角的窗格则为通过带通滤波后的波形，滤波信号的FFT结果在右下角的窗格中。可以发现，带通滤波器消除了低频拾取噪声和80 kHz位置的二次谐波。滤波后的信号在时间域视图现在有了一个相对平坦的基线，从而能使反射信号更清楚地分辨出来，而这就是我们滤波处理的目标。由此我们也可以看出FFT对深入了解信号提供的帮助。

使用德思特SBench 6软件提供的信号处理工具，如模拟运算、平均值、FFT、滤波和直方图等，将有助于加深对采集信号的见解，此外，还能生成一系列有利于进一步分析的二级信号波形。