qt学习--自相关函数代码

最新推荐文章于 2024-11-13 08:26:19 发布

流心蛋黄

最新推荐文章于 2024-11-13 08:26:19 发布

阅读量606

点赞数

CC 4.0 BY-SA版权

分类专栏： qt学习文章标签： matlab c++

本文链接：https://blog.csdn.net/weixin_43139923/article/details/120327761

qt学习专栏收录该内容

1 篇文章

订阅专栏

这篇博客介绍了如何在C++中实现自相关函数，代码基于MATLAB的xcorr()函数，用于信号处理中的序列相关性分析。博主分享了完整的C++代码，并指出处理大量数据时欢迎探讨更高效的实现方式。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

自相关函数c++代码

最近在做信号处理，将matlab中的xcorr()自相关函数用c++写出来，mark一下，自相关的概念有空再补充

void  xcorr_func(vector<double> &in , vector<double> &out)
{
    int i,j;
    for (i=0;i<=lenth;i++)
    {
        for (j=0;j<=i;j++)
        {
            out[i]+=in[j]*in[lenth-1-i+j];
        }

        //如果对应的是自相关的无偏估计的话，还应加上此步骤，对应的是xcorr(x,'unbiased')
        out[i]=out[i]/j;
    }

    //自相关具有对称性，直接用前一半的数值来对应后一半的数值
    for (int i=lenth;i<(lenth<<1)-1;i++)
    {
        out[i]=out[(lenth<<1)-2-i];
    }
}

完整代码如下，因为处理数据很大，如果有更快的更省内存的计算方法欢迎来交流

//自相关函数，是对信号自身的互相关，表示同一序列不同时刻的相关程度
//一维数据的自相关函数计算，对应matlab中的"xcorr"函数
#include <iostream>
#include<vector>
#include<math.h>

using namespace std;
int lenth;

//自相关函数
void  xcorr_func(vector<double> &in , vector<double> &out)
{
    int i,j;
    for (i=0;i<=lenth;i++)
    {
        for (j=0;j<=i;j++)
        {
            out[i]+=in[j]*in[lenth-1-i+j];
        }

        //如果对应的是自相关的无偏估计的话，还应加上此步骤，对应的是xcorr(x,'unbiased')
        out[i]=out[i]/j;
    }

    //自相关具有对称性，直接用前一半的数值来对应后一半的数值
    for (int i=lenth;i<(lenth<<1)-1;i++)
    {
        out[i]=out[(lenth<<1)-2-i];
    }
}

int main()
{
    //自己定义的数组
    double a[]={1,2,3,4,5,6,7,8,9,10,11,12,13,14,15};
    //vector向量o从数组中赋初值
    vector<double> in(a,a+15);
    lenth = in.size();                            //数组长度
    int out_len = lenth*2 - 1;                //自相关结果长度，也可以不用这一步，直接在计算中使用push_back添加
    vector<double> out(out_len);  //定义vector向量存储结果

        xcorr_func(in,out);
    for (int	i=0;i<=28;i++)
    {
        std::cout<<out[i]<<  std::endl;
    }
    return 0;
}