其实关于二氧化碳流量计信号的特点及稀疏傅里叶变换是有理论在的

其实当测量使用的二氧化碳流量计信号的特点为什么说是有理论的，小编总结了一下在这篇文章中的二氧化碳流量计信号的特点及稀疏傅里叶变换的详细介绍，接下来就跟随着小编的步伐让我们来看看究竟都是如何的吧？

传统的涡街信号处理方法主要是用傅里叶变换，由于传统的傅里叶变换的时间复杂度与其分析的信号长度成正比，采样点数越多其频谱谱线越接近理想状态，但需要较大的运算时间。针对这个问题，提出了采用稀疏傅里叶变换分析涡街流量信号的方法，利用实验数据进行 Matlab 仿真，验证了该方法的性能优于传统的方法，不仅提高了对强噪声的抗干扰能力还加快了计算速度。 近半个世纪以来， 二氧化碳计量表因其测量精度高、 无可动部件、测量精度高等优点得到了迅猛的发展。二氧化碳流量计主要测量部件为压电传感器，其易受到噪声的干扰，如管道振动、电磁干扰、流体的低频摆动等。在含有噪声的信号的中提取出有用的涡街信号，国内外众多研究学者对涡街信号的处理方式主要有FFT 的周期图法、互相关法、自适应陷波滤波法、小波分析法和数字跟踪滤波方法等。但是这些方法对于含有强噪声的信号测量精度不高或错误，即噪声频率在涡街信号频率范围内，而噪声的幅值高于涡街信号的幅值。 对此提出一种基于稀疏傅里叶变换的涡街信号分析方法， 该方法不仅具有很高的实时性而且对含有强噪声的信号也能够保证测量的准确性。 

1、稀疏傅里叶变换的理论分析 

快速傅里叶变换（ Fast Fourier Transform ， FFT ）的时间复杂度为 O （ nlogn ），与离散傅里叶变换（ Discrete Fourier Trans-form ， DFT ）的复杂度 O （ N 2 ）相比，运算速度发生了质的飞跃，尤其是随着采样点数 N 的增加这种优势就越加明显。但是随着时代的发展，需要实时处理的信号越来越多，即便是 FFT 对于这样的需求也难以满足。传统的 FFT 只考虑到了信号的长度 N 需要为2 的整数次幂，并未考虑到信号的自身的特性，如稀疏性。 在实际生活中常见的信号的傅里叶系数只有小部分是我们感兴趣的，其大部分都是可以忽略的，如图像和语音信号。针对这样的信号能否找到一种更加快速的算法来计算其傅里叶变换，MIT的团队给出了答案。该团队提出了稀疏傅里叶变换（ Sparse Fourier Transform ， SFT ）， 该算法利用了信号频域的稀疏性，先对信号进行分“桶”，将长的 DFT 运算变为较短的运算，再根据一定的规则重构了信号的频谱，其运算速度为 FFT 的十倍甚至百倍。稀疏傅里叶变换使用的先决条件就是分析的信号具有稀疏性，设 x （ n ）是长度为 N 点的有限长序列，则该序列的 N 点离散傅里叶变换逆变换为：

其中 Ω N 表示集合邀0 ， 1 ，……， N－1妖 。 只有 K （ K垲N ）个非零的傅里叶系数，只通过信号 x （ n ）的部分 采样值来确定这 K 非零傅里叶系数与位置。 

2、频域降采样 参数 B 整除 N ，若想要以等间隔 N／B 对信号频域进行降采样，即：

混叠后频域谱线由 N 减少到 B ，信号点数成倍较少，这正是SFT算法复杂度为亚线性的关键原因之一。 

3、稀疏傅里叶变换运算步骤 

稀疏傅里叶变换包括频谱重排、加窗函数、频域降采样、定位、估值与迭代等运算过程。

4、频谱重排频谱重排的目的是使各大值点均匀分布， 分桶时大频点不要分到同一个桶中，当两个或两个以上大值点在同一桶中时， 无法求解取大值点的频率和位置。p （ n ） ＝x ｛ mod ［ σ · n ， N ］｝， nε ［ 1 ， N ］ （ 4 ）式中 σ 为一个随 机数，且为奇数，并满足 mod ［ σ×σ －1 ， N ］ ＝1 ，这就保证了 σ 与 N 互为质数， σ －1 为 σ 的模逆算子。 根 据傅里叶变换可知上式中的 p （ n ）， x （ n ）满足：P （ k ） ＝X ｛ mod ［ σ－1 · k ， N ］｝， σ ， kε ［ 1 ， N ］ （ 5 ）通过式（ 4 ）、（ 5 ）知道信号时域上的重排也会导致频谱信号位置上发生变换。

5、窗函数滤波器 

为了保证算法的效率且防止频谱泄漏，需要设计一个在时域和频域能量都集中的滤波器， 根据文献该滤波器的为 sinc 窗函数与高斯窗函数的卷积，该窗函数具有过渡带陡峭、通带平滑等特点。 

6、哈希映射定义一个映射区间 Ω N →Ω B 的哈希函数： h σ （ k ） ＝round （ σ ·k · N／B ）， round 表示四舍五入， 将 Ω N 中每一个点都映射到 Ω B中。 定义偏移量： o σ （ k ） ＝σ · k－h σ （ k ）·（ N／B ）；定义集合 J ，集合 J包含了 Z （ k 中 K 个较大幅值的坐标 k ；通过哈希反映射得到 I r ，即 I r ＝邀kε ［ 0 ， N－1 ］ ｜h σ （ k ） εJ妖 ，后 从中取出 K 个大值点对原信号的频率估计。

7、循环投票对于每一个 kεI ， X‘（ k ） ＝Z （ h σ （ k ） W Nτk／G （ o σ （ k ））频率估计值。 每一次定位循环得到一个坐 标集合 I r ，在 L＝O （ log 2 N ）次循环中，对任意坐标 kεI＝I 1 U …… υI r ，若出现次数大于 L／2 ，则将其归入集 合 I‘ 中，并认为集合 I‘ 包含所有目标频点坐标。 对每一个kεI‘ ，取 L 次循环得到 X （ k ）的中值作为终的频率 值，即：X （ k ） ＝median （ 邀X r （ k ） ｜rε邀1 ，……， L) 2 涡街信号的特点 在一定范围内，流体流速 V 与涡街频率 f 有以下关系： f＝πK 1 VD2／4 （ 6 ） 其中 K 1 为仪表系数， D 为管道直径。在管道口径 D 不变，流体密度不变的情况下，涡街传感器的输出幅值与 f 2 成正比，具体表达形式可以根据实验测出。本文以 50mm 口径气体实验为例，数据如表 1 所示：

从表 1 的气体流量的实际幅值和拟合幅值的误差可以看出，涡街信号的幅值在理论值附近波动，且波动的范围一定，则幅频关系更一般的形式表达如下： 其中 c 为系数， δ 为相对误差限，其示意图如图 2 所示，图中实线为幅频关系的理论拟合曲线，而虚线为幅值波动的阈值曲线。根据实验的数据，涡街信号幅值波动的相对误差为 ±10％ 。

本文提出了一种噪声情形下新颖的室内定位与环境构建算法，该算法不需要无线传感器网络中的节点配备激光传感模块，只需要通过节点间的射频信号即可实现对周围环境的感知，为室内定位与环境构建技术提供了一种经济且可靠的解决方案。

好啦，以上内容的描述总结就是这些啦，对此不知道大家都认为是如何的呢？小编在此可是欢迎大家随时来电或在线与小编分享，了解更多流量计知识尽在威海坤科流量仪表，在下期的内容中大家还有哪些其他需要补充得地方呢？关注小编，下期的内容一定会更加精彩，点击右下角弹窗与小编进一步详询吧。