我编写这段代码是为了平滑曲线。 它在一个点旁边取 5 个点并将它们相加并取平均值。

/* Smoothing */
void smoothing(vector<Point2D> &a)
{
    //How many neighbours to smooth
    int NO_OF_NEIGHBOURS=10;
    vector<Point2D> tmp=a;
    for(int i=0;i<a.size();i++)
    {

        if(i+NO_OF_NEIGHBOURS+1<a.size())
        {
            for(int j=1;j<NO_OF_NEIGHBOURS;j++)
            {
                a.at(i).x+=a.at(i+j).x;
                a.at(i).y+=a.at(i+j).y;
            }
            a.at(i).x/=NO_OF_NEIGHBOURS;
            a.at(i).y/=NO_OF_NEIGHBOURS;

        }
        else
        {
            for(int j=1;j<NO_OF_NEIGHBOURS;j++)
            {
                a.at(i).x+=tmp.at(i-j).x;
                a.at(i).y+=tmp.at(i-j).y;
            }
            a.at(i).x/=NO_OF_NEIGHBOURS;
            a.at(i).y/=NO_OF_NEIGHBOURS;
        }

    }

}

但是我得到的每个点的值都非常高，而不是与前一个点相似的值。形状最大化了很多，这个算法哪里出了问题？

最佳答案

您在这里看到的是一个实现了 boxcar window function 的有限脉冲响应 (FIR) 滤波器的低音响应实现。 .考虑 DSP 方面的问题，您需要使用 NO_OF_NEIGHBOURS 相等的 FIR 系数过滤传入的 vector，每个系数的值为 1/NO_OF_NEIGHBOURS。通常最好使用既定算法而不是重新发明轮子。

这是一个非常简陋的实现，我很快就敲定了它，它使过滤器加倍。您可以轻松修改它以过滤您的数据类型。该演示仅出于演示目的显示了对上升锯函数 (0,.25,.5,1) 的几个周期进行过滤。它会编译，因此您可以使用它。

#include <iostream>
#include <vector>

using namespace std;

class boxFIR
{
    int numCoeffs; //MUST be > 0
    vector<double> b; //Filter coefficients
    vector<double> m; //Filter memories

public:
    boxFIR(int _numCoeffs) :
    numCoeffs(_numCoeffs)
    {
        if (numCoeffs<1)
            numCoeffs = 1; //Must be > 0 or bad stuff happens

        double val = 1./numCoeffs;
        for (int ii=0; ii<numCoeffs; ++ii) {
            b.push_back(val);
            m.push_back(0.);
        }
    }    

    void filter(vector<double> &a)
    {
        double output;

        for (int nn=0; nn<a.size(); ++nn)
        {
            //Apply smoothing filter to signal
            output = 0;
            m[0] = a[nn];
            for (int ii=0; ii<numCoeffs; ++ii) {
                output+=b[ii]*m[ii];
            }

            //Reshuffle memories
            for (int ii = numCoeffs-1; ii!=0; --ii) {
                m[ii] = m[ii-1];
            }                        
            a[nn] = output;
        }
    }


};

int main(int argc, const char * argv[])
{
    boxFIR box(1); //If this is 1, then no filtering happens, use bigger ints for more smoothing

    //Make a rising saw function for demo
    vector<double> a;
    a.push_back(0.); a.push_back(0.25); a.push_back(0.5); a.push_back(0.75); a.push_back(1.);
    a.push_back(0.); a.push_back(0.25); a.push_back(0.5); a.push_back(0.75); a.push_back(1.);
    a.push_back(0.); a.push_back(0.25); a.push_back(0.5); a.push_back(0.75); a.push_back(1.);
    a.push_back(0.); a.push_back(0.25); a.push_back(0.5); a.push_back(0.75); a.push_back(1.);

    box.filter(a);

    for (int nn=0; nn<a.size(); ++nn)
    {
        cout << a[nn] << endl;
    }
}

使用这条线增加滤波器系数的数量，以查看逐渐更加平滑的输出。只有 1 个滤波器系数，没有平滑。

boxFIR box(1);

代码非常灵活，您甚至可以根据需要更改窗口形状。通过修改构造函数中定义的系数来做到这一点。

注意:这将为您的实现提供略微不同的输出，因为这是一个因果过滤器(仅取决于当前样本和先前样本)。您的实现不是因果关系，因为它会及时预测 future 的样本以求平均值，这就是为什么您需要条件语句来处理接近 vector 末尾的情况。如果您希望输出像您尝试使用此算法对过滤器执行的操作，请通过此算法反向运行您的 vector (只要窗口函数是对称的，这就可以正常工作)。这样您就可以获得类似的输出，而无需算法的讨厌的条件部分。

关于c++ - 平滑算法，我们在Stack Overflow上找到一个类似的问题： https://stackoverflow.com/questions/12973255/