信号采样中的滤波处理

以前曾经在别人的帖子里回复过这个算法, 今天又看到有人问类似的问题, 就特别提一下, 以备将来有人需要时参考.

在各种信号采集后处理的过程中,对原始信号中夹杂的噪音去除是一个比较广泛的需求.

经典的高斯滤波或中值滤波都可以消除高频噪音, 但是前者对原始波形的损伤较大,后者算法略为复杂,

本文将介绍一个我曾经采用过的高频噪音滤波算法. 复杂度和处理速度介于两者之间, 效果接近中值滤波.

缺点在于内存占用略多.

首先想到这个算法的时候是一次为一个朋友做的示波器数据采样分析.

因为信号源是一个高频冲击波,在波峰附近有很多高频谐波,并且谐波的峰值不小于主波的20%,如果不加以去除的话,后面的数据分析就会很离谱.(因为是谐波,所以其峰值和主波峰值有关联)

一般来说针对这样的高频杂波,用均值滤波器也就可以了,但是实际操作时发现问题很大,因为采样得到的波源频率已经接近示波器的采样频率极限,所以在关键的波峰波谷阶段的有效采样点并不多,使用均值滤波器的话,就会发生明显的波幅失真.

实际要作的工作就类似于拿起一支铅笔,把一张对焦问题严重的照片,把照片中间模糊的人物轮廓给"描"出来. (这实际上已经不是一般的去噪问题了)

我当时想出的一个办法就是:先插值,再均值. (后来才知道,3D里面的抗锯齿就是这么干的)

插值就不赘述了,与图片缩放不同的是,这里的采样数据只需要一维的线形插值即可,所以更简便了(线形插值可以参考我之前的文章,有详细的算法说明和历程).

值得一提的就是插值和均值的倍率取舍, 经过大致的估算, 如果杂波波幅占主波的10%的话,采用5倍插值(即两个采样数据中间插4个点)和均值(每个点取自身和前后各2个点,一共5点的均值)就可以得到一个比较完美波形描线.