如何使用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本身并不能直接改善信号质量，但它能显示信号的频率成分结构，并提供了如何移除不需要的频谱成分的信息。
fft产生的频谱有一个离散的时间轴，就像时域信号有离散的采样时间一样。频谱中的谱线，通常称为“bins”或“cells”，每个谱线之间由分辨率带宽（f）隔开，而分辨率带宽与采集到的原始信号长度成反比，因此，要增加fft频谱的频率分辨率，必须增加采集的原始信号长度。而频谱显示的总频率范围，或者说频率跨度，是原始信号采样率的一半，因此，要增加总频率跨度，也必须增加采样率。
在德思特sbench 6软件中，fft结果的纵轴缩放可以设为线性单位伏特，也可以是对数单位分贝（db）。对数刻度可以参考数字化仪的满偏刻度（dbfs）、1 毫瓦（dbm）、1 微伏（dbμv）或设为针对假定载波频谱中的最大峰值（dbc）。
（1）权重函数理论上的傅立叶变换假设输入记录的长度为无穷大。有限的采样长度则会在频域边缘引入不连续性，为频域引入伪频率，并一定程度扭曲实际的频谱。例如，当信号的开始和结束相位不同，或者信号频率落在两个相邻的谱线之间时，使频谱变宽。
频谱展宽，以至于扩展到许多相邻的频带，称为泄漏。对于这种问题的理论应对方法是确保在显示网格内包含整数个周期，或在边缘处不出现不连续性。两者都需要信号波形频率和数字化采样率之间非常精确的同步，并且需要准确地设置采集长度，这通常只可能在实验室中实现，而不能作用于真实世界的信号。另一种方法是使用窗函数（加权）来平滑信号的边缘。
为了尽量减少这些负面影响，对获取的信号应用加权函数，使记录的端点为零。德思特sbench 6软件中的fft功能为用户提供八种加权函数的选择，而这些加权函数则可以改变谱线的形状。对此没什么概念的读者，可以把它想象为组合到一起的一些列并行的理想带通滤波器，这些滤波器之间的频率间距为分辨率带宽，而加权函数则会影响滤波器频率响应的形状。图1比较了四种最常用的加权函数对应的频率响应。
图1：四种最常用的加权函数的谱形状比较
理想情况下，主瓣应尽可能窄，且两侧平坦，以代表真实的频谱成分，而所有旁瓣应无限衰减。窗函数类型定义了在fft处理中使用的等效滤波器的带宽和形状。表1中展示了频谱响应的最大旁瓣幅度。较大的旁瓣电平将有助于区分密集的频谱成分。
表1：不同加权函数的关键特性
如前所述，fft结果的频率轴是离散的，具有以分辨率带宽的倍数间隔的频点。如果输入信号频率落在两个相邻的谱线之间，能量将被分配到两个谱线中，且峰值幅度将降低。这被称为栅栏效应或波浪状效应，而扩宽频谱响应可以减少幅度变化。表1中的栅栏损失列给出了不同加权函数该效应幅值。
同时，加权函数还会影响频谱响应的带宽。等效噪声带宽（enbw）说明了相对于矩形窗加权带宽的相对变化。将功率谱归一化到测量带宽（功率谱密度）需要将功率谱除以enbw与分辨率带宽的积（f✖️enbw）。
相干增益描述了给定加权函数相对于矩形窗加权的频谱幅度变化。这是一个应用在所有频率上的固定增益，可以很容易地归一化。
矩形窗加权函数是采集信号没有任何加权的响应。它具有最窄的带宽，但旁瓣幅值则相当高。由于采集时域记录中的所有点的幅度响应都是均匀的，它常用于具有瞬态性质的信号（或相对记录总长短得多的信号）。当需要以最佳频率精度进行分析时，也会使用它。
hanning和hamming加权函数具有良好的通用的频率响应，能提供较好的频率分辨率以及合理的旁瓣响应。blackman-harris则旨在获得最佳幅度精度和优秀的旁瓣抑制。
（2）fft应用实例图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成分之外的其它成分。现在我们将带着这个目标，来进行滤波器的设置。
（3）滤波在德思特sbench 6专业版软件中，我们提供了低通、带通或高通选项的有限冲激响应(fir)数字滤波器。通过输入所需的滤波器类型、截止频率或频率，以及滤波器阶数，用户可以直接在图形界面中创建这些滤波器。sbench 6软件会在滤波器无法实现时给出提示，并提出解决建议。或者，您也可以输入从其他来源获取的滤波器系数。我们将这些滤波器应用于采集到的信号，然后将滤波结果与原始采样结果或平均采样结果进行比较。在图3中，应用的带通fir滤波器截止频率为30和50 khz，以帮助提取目标信号。
图3：原始波形和滤波后信号及其fft结果的对比
左上角的窗格显示了原始波形，下面对应我们之前看到的原始信号fft结果。右上角的窗格则为通过带通滤波后的波形，滤波信号的fft结果在右下角的窗格中。可以发现，带通滤波器消除了低频拾取噪声和80 khz位置的二次谐波。滤波后的信号在时间域视图现在有了一个相对平坦的基线，从而能使反射信号更清楚地分辨出来，而这就是我们滤波处理的目标。由此我们也可以看出fft对深入了解信号提供的帮助。
结论使用德思特sbench 6软件提供的信号处理工具，如模拟运算、平均值、fft、滤波和直方图等，将有助于加深对采集信号的见解，此外，还能生成一系列有利于进一步分析的二级信号波形。
CES:博通将发新WiFi芯片 比现有技术快三倍
中国联通同中兴和腾讯开展网络切片项目
解锁你的 Tesla 完全自动驾驶能力
iPhone12Pro全面支持杜比视界视频，能力之强让许多电影摄像机都羡慕
数字化有什么作用？数字化转型是好还是坏？
如何使用SBench 6对数字化仪采集信号进行处理？（三）——快速傅立叶变换（FFT）
联发科高端旗舰市场再进一步，可能是最好的移动游戏平台发布，开发生态更强
关于模拟弧长反馈新技术分析探讨介绍
LDO的基础知识介绍
CD4051模拟开关芯片实现IO口扩展
压力变送器工作原理及分解图如何选型？
首届广州国际电子及电器博览会圆满落幕
临界流文丘里喷嘴气体流量计检定系统
波士顿动力或将再度被卖，买家或是现代汽车
汽车AFS的电机驱动方案及应用设计要点
应用在触控一体机中的电容式触摸芯片
GP328GP88S系对讲机维修实例
3日内更新两次 苹果向开发者推出iOS10.3 beta7！
冰箱响的原因及解决办法浅析
8英寸晶圆到底有多缺？
原文链接：https://mp.weixin.qq.com/s/mys1idxfnvfrecggtf78mw
欢迎关注虹科，为您提供最新资讯！
#信号采集 #信号处理 #测量测试
上一篇文章介绍了德思特sbench 6的平均运算功能。本章将继续为大家介绍sbench 6的快速傅立叶变换（fft）。
前文回顾：https://mp.weixin.qq.com/s/j-in_2jrn9zhgmaaaysdjg
快速傅立叶变换（ fft ）快速傅立叶变换（fft）将获取的波形从时域（幅度对时间）映射到频谱域（幅度对频率）。这使得用户可以观察构成信号的频率成分。fft本身并不能直接改善信号质量，但它能显示信号的频率成分结构，并提供了如何移除不需要的频谱成分的信息。
fft产生的频谱有一个离散的时间轴，就像时域信号有离散的采样时间一样。频谱中的谱线，通常称为“bins”或“cells”，每个谱线之间由分辨率带宽（f）隔开，而分辨率带宽与采集到的原始信号长度成反比，因此，要增加fft频谱的频率分辨率，必须增加采集的原始信号长度。而频谱显示的总频率范围，或者说频率跨度，是原始信号采样率的一半，因此，要增加总频率跨度，也必须增加采样率。
在德思特sbench 6软件中，fft结果的纵轴缩放可以设为线性单位伏特，也可以是对数单位分贝（db）。对数刻度可以参考数字化仪的满偏刻度（dbfs）、1 毫瓦（dbm）、1 微伏（dbμv）或设为针对假定载波频谱中的最大峰值（dbc）。
（1）权重函数理论上的傅立叶变换假设输入记录的长度为无穷大。有限的采样长度则会在频域边缘引入不连续性，为频域引入伪频率，并一定程度扭曲实际的频谱。例如，当信号的开始和结束相位不同，或者信号频率落在两个相邻的谱线之间时，使频谱变宽。
频谱展宽，以至于扩展到许多相邻的频带，称为泄漏。对于这种问题的理论应对方法是确保在显示网格内包含整数个周期，或在边缘处不出现不连续性。两者都需要信号波形频率和数字化采样率之间非常精确的同步，并且需要准确地设置采集长度，这通常只可能在实验室中实现，而不能作用于真实世界的信号。另一种方法是使用窗函数（加权）来平滑信号的边缘。
为了尽量减少这些负面影响，对获取的信号应用加权函数，使记录的端点为零。德思特sbench 6软件中的fft功能为用户提供八种加权函数的选择，而这些加权函数则可以改变谱线的形状。对此没什么概念的读者，可以把它想象为组合到一起的一些列并行的理想带通滤波器，这些滤波器之间的频率间距为分辨率带宽，而加权函数则会影响滤波器频率响应的形状。图1比较了四种最常用的加权函数对应的频率响应。
图1：四种最常用的加权函数的谱形状比较
理想情况下，主瓣应尽可能窄，且两侧平坦，以代表真实的频谱成分，而所有旁瓣应无限衰减。窗函数类型定义了在fft处理中使用的等效滤波器的带宽和形状。表1中展示了频谱响应的最大旁瓣幅度。较大的旁瓣电平将有助于区分密集的频谱成分。
表1：不同加权函数的关键特性
如前所述，fft结果的频率轴是离散的，具有以分辨率带宽的倍数间隔的频点。如果输入信号频率落在两个相邻的谱线之间，能量将被分配到两个谱线中，且峰值幅度将降低。这被称为栅栏效应或波浪状效应，而扩宽频谱响应可以减少幅度变化。表1中的栅栏损失列给出了不同加权函数该效应幅值。
同时，加权函数还会影响频谱响应的带宽。等效噪声带宽（enbw）说明了相对于矩形窗加权带宽的相对变化。将功率谱归一化到测量带宽（功率谱密度）需要将功率谱除以enbw与分辨率带宽的积（f✖️enbw）。
相干增益描述了给定加权函数相对于矩形窗加权的频谱幅度变化。这是一个应用在所有频率上的固定增益，可以很容易地归一化。
矩形窗加权函数是采集信号没有任何加权的响应。它具有最窄的带宽，但旁瓣幅值则相当高。由于采集时域记录中的所有点的幅度响应都是均匀的，它常用于具有瞬态性质的信号（或相对记录总长短得多的信号）。当需要以最佳频率精度进行分析时，也会使用它。
hanning和hamming加权函数具有良好的通用的频率响应，能提供较好的频率分辨率以及合理的旁瓣响应。blackman-harris则旨在获得最佳幅度精度和优秀的旁瓣抑制。
（2）fft应用实例图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成分之外的其它成分。现在我们将带着这个目标，来进行滤波器的设置。
（3）滤波在德思特sbench 6专业版软件中，我们提供了低通、带通或高通选项的有限冲激响应(fir)数字滤波器。通过输入所需的滤波器类型、截止频率或频率，以及滤波器阶数，用户可以直接在图形界面中创建这些滤波器。sbench 6软件会在滤波器无法实现时给出提示，并提出解决建议。或者，您也可以输入从其他来源获取的滤波器系数。我们将这些滤波器应用于采集到的信号，然后将滤波结果与原始采样结果或平均采样结果进行比较。在图3中，应用的带通fir滤波器截止频率为30和50 khz，以帮助提取目标信号。
图3：原始波形和滤波后信号及其fft结果的对比
左上角的窗格显示了原始波形，下面对应我们之前看到的原始信号fft结果。右上角的窗格则为通过带通滤波后的波形，滤波信号的fft结果在右下角的窗格中。可以发现，带通滤波器消除了低频拾取噪声和80 khz位置的二次谐波。滤波后的信号在时间域视图现在有了一个相对平坦的基线，从而能使反射信号更清楚地分辨出来，而这就是我们滤波处理的目标。由此我们也可以看出fft对深入了解信号提供的帮助。
结论使用德思特sbench 6软件提供的信号处理工具，如模拟运算、平均值、fft、滤波和直方图等，将有助于加深对采集信号的见解，此外，还能生成一系列有利于进一步分析的二级信号波形。
CES:博通将发新WiFi芯片 比现有技术快三倍
中国联通同中兴和腾讯开展网络切片项目
解锁你的 Tesla 完全自动驾驶能力
iPhone12Pro全面支持杜比视界视频，能力之强让许多电影摄像机都羡慕
数字化有什么作用？数字化转型是好还是坏？
如何使用SBench 6对数字化仪采集信号进行处理？（三）——快速傅立叶变换（FFT）
联发科高端旗舰市场再进一步，可能是最好的移动游戏平台发布，开发生态更强
关于模拟弧长反馈新技术分析探讨介绍
LDO的基础知识介绍
CD4051模拟开关芯片实现IO口扩展
压力变送器工作原理及分解图如何选型？
首届广州国际电子及电器博览会圆满落幕
临界流文丘里喷嘴气体流量计检定系统
波士顿动力或将再度被卖，买家或是现代汽车
汽车AFS的电机驱动方案及应用设计要点
应用在触控一体机中的电容式触摸芯片
GP328GP88S系对讲机维修实例
3日内更新两次 苹果向开发者推出iOS10.3 beta7！
冰箱响的原因及解决办法浅析
8英寸晶圆到底有多缺？