几个简单的数据点平滑处理算法

河桐

已于 2022-08-05 20:22:12 修改

阅读量9.2k

点赞数 1

分类专栏： c++ 文章标签：点平滑 c

于 2018-04-19 18:13:54 首次发布

原文链接：https://blog.csdn.net/liyuanbhu/article/details/11119081?spm=1001.2101.3001.6661.1&utm_medium=distribute.pc_relevant_t0.none-task-blog-2%7Edefault%7ECTRLIST%7Edefault-1-11119081-blog-119898401.pc_relevant_multi_platform_featuressortv2removedup&depth_1-utm

版权

c++ 专栏收录该内容

33 篇文章 0 订阅

订阅专栏

最近在写一些数据处理的程序。经常需要对数据进行平滑处理。直接用FIR滤波器或IIR滤波器都有一个启动问题，滤波完成后总要对数据掐头去尾。因此去找了些简单的数据平滑处理的方法。

在一本老版本的《数学手册》中找到了几个基于最小二乘法的数据平滑算法。将其写成了C 代码，测试了一下，效果还可以。这里简单的记录一下，算是给自己做个笔记。

算法的原理很简单，以五点三次平滑为例。取相邻的5个数据点，可以拟合出一条3次曲线来，然后用3次曲线上相应的位置的数据值作为滤波后结果。简单的说就是 Savitzky-Golay 滤波器。只不过Savitzky-Golay 滤波器并不特殊考虑边界的几个数据点，而这个算法还特意把边上的几个点的数据拟合结果给推导了出来。

不多说了，下面贴代码。首先是线性拟合平滑处理的代码. 分别为三点线性平滑、五点线性平滑和七点线性平滑。

void linearSmooth5 ( double in[], double out[], int N )  
{  
    int i;  
    if ( N < 5 )  
    {  
        for ( i = 0; i <= N - 1; i++ )  
        {  
            out[i] = in[i];  
        }  
    }  
    else  
    {  
        out[0] = ( 3.0 * in[0] + 2.0 * in[1] + in[2] - in[4] ) / 5.0;  
        out[1] = ( 4.0 * in[0] + 3.0 * in[1] + 2 * in[2] + in[3] ) / 10.0;  
        for ( i = 2; i <= N - 3; i++ )  
        {  
            out[i] = ( in[i - 2] + in[i - 1] + in[i] + in[i + 1] + in[i + 2] ) / 5.0;  
        }  
        out[N - 2] = ( 4.0 * in[N - 1] + 3.0 * in[N - 2] + 2 * in[N - 3] + in[N - 4] ) / 10.0;  
        out[N - 1] = ( 3.0 * in[N - 1] + 2.0 * in[N - 2] + in[N - 3] - in[N - 5] ) / 5.0;  
    }  
}

几个简单的数据点平滑处理算法_liyuanbhu的博客-CSDN博客_如何平滑数据