多次采样求均值的所谓滤波算法到底用的对不对？

在采样信号时，干扰几乎很难避免的，通过软件滤波是非常常见的手段。而最常用的是连续多个采样并ADC，然后将多个采样值相加除以数据个数求得平均值；

也有的同样多个采样，去掉一个最高值，去掉一个最低值，其他中间数值求平均值。

最近在几家单位做设计审核时，是几家，都存在这个问题，可见这不是冷门的现象。看到设计时使用了这个方法，我直接问出了扎心几连击的问题。

这种方法其实是一个数字低通滤波器，但现实情况是不是真的需要的是一个低通滤波，即有用信号的频率低于干扰噪声的频率。只有符合这个前提的抗干扰类型才适用于求均值方式。然后就是以下的几个问题。

1、准备处理的有用信号的频率是多少？

2、存在的干扰的频段是多少？

3、干扰的分布规律是服从什么分布的？

4、采样并ADC一个数据的周期是多长时间？

5、连续采样多个求均值，具体多个数据的个数是多少个？为什么？

如果各位能准确地回答出这几个问题，下面的文字就不用看了。这是给回答不出问题来的朋友看的。

1、准备处理的有用信号的频率是多少？

2、存在的干扰的频段是多少？

以上两个问题的答复：能自己分析的就自己分析出来，不能分析出来的就用测试方式，启动示波器的FFT功能，由时域波形转换成频谱图，提炼出的频谱毛刺除了基频就是基频的整数倍频率。这个道理请参阅2021-01出版的《电子产品世界》中的《电路信号线间串扰机理及措施研究》一文（知网www.cnki.net可以检索得到）。

3、干扰的分布规律是服从什么分布的？

这个问题暂时我也不知道怎么回答，凭常识和经验能知道，但用讲道理的方式讲出来暂时没想好。欲深入探究建议去研究下统计学里的几种分布规律，看看均匀分布/正态分布/二项分布/指数分布/威布尔分布/...分别会在哪类事件特征下存在。一般的EMI干扰服从正态分布、脉冲式工况导致的瞬变干扰服从指数分布、温漂导致的数据变化规律属于非线性漂移... ...

4、采样并ADC一个数据的周期是多长时间？

采样定理我们大家是都学过的了，众所周知，采样频率必须大于欲采集信号的最高频率的至少2倍以上，才能保证被采集的信号能恢复出其原始形状来，但这是不够的，见（图1）。

           

  图1

A为原始被采样信号，1234为采样点，采样频率刚好是被采样信号的2倍，B为根据采样点的电压恢复出来的信号，由示意图可以看出，如果采样频率fsa=2*fs信号频率的话，恢复后的信号的幅值会与原信号有出入。也就是说，2fs的采样频率，充其量只能保证信号的基本形状能被保留下来，但幅值数据并不很准确。工程上，一个周期中的采样点越多，根据采样数据拟合出来的波形才越逼近于真实波形。所以，一般至少要fsa > 5fs（推荐）。示波器常标注的5GSa指标就是这个道理。

由以上道理，如果有用信号的最高频率为fsmax，则ADC速率fsamin > 5fsmax。如1kHz的fs，fsa至少5kHz，采样周期为1/fsa =Tsa <0.2ms，这里采样周期假设取20us。

5、连续采样多个数值求均值，具体这里的多个应该是多少个？为什么？

求均值这种计算方式，针对就是正态分布的这种干扰。正态分布的分布特点，正偏差负偏差对称，大部分的随机干扰都是服从这种分布形式。按照这种形式，取一个周期内的好多个数据，正偏差和负偏差的个数和偏差程度应该是对称相加求均值是能正负抵消掉的。以图2为例，正常信号是一个1kHz稳定信号，干扰是5kHz，则多点采样的总采样时间必须至少涵盖一个完整的噪声周期，如果采样的点太少，数据个数和时间间隔不足一个周期，相加后正负不能完全抵消，噪声就消除不了。5kHz的周期时间为200us。按第4条，采样时间周期Tsa=20us；按第5条，需要采集的数据总时间n*Tsa > 200us。因此采样点须取至少n > 10个。若低于10个，则滤波达不到期望滤除效果。

图2

不知大家是否会问，我的干扰信号不是正弦波，这么用还可以吗？当然可以，君莫忘，任何规律性波形是不是都可以傅立叶展开分解成一系列正弦波相加的形式，反正是个低通滤波，我们取其最低正弦波作为计算便可以了，最低的能滤掉，高频的自然更能滤掉了。

以上几个问题就回答完了，总结起来主要是四点：

a）、需要低通滤波的数据才适用此方法；

b）、服从随机正态分布、随机均匀分布的干扰，才适用此方法；

c）、fsamin > 5fsmax

fsamin：可选择的最低采样频率；

fsmax：信号的最高频率

d）、n*Tsa > TNmax

n：累加求均值的采样数据个数；

Tsa：采样一个数据的周期；

TNmax：噪声的最大周期，也即最低频噪声的周期。

最后一句：GOOD  LUCK。