如何使用vDSP_conv来模仿MATLAB的xcorr函数?

时间:2012-12-10 15:30:15

标签: ios matlab fft accelerate-framework cross-correlation

我目前正在将MATLAB算法转换为C,以便在iOS应用程序中使用它。

我一直在努力学习MATLAB的xcorr功能。这是相关的MATLAB代码。

xcr = xcorr(A,A,maxlags);

根据MATLAB文档

  

返回滞后范围[-maxlags:maxlags]上的互相关序列。输出c的长度为2 * maxlags + 1.

Apple Accelerate.Framework提供了一个名为vDSP_conv的卷积/相关函数,但我没有看到如何使用它产生与xcorr相同的输出。这可能吗 ?如果有,任何人都可以帮助我。

最诚挚的问候,

Acacio

1 个答案:

答案 0 :(得分:4)

要复制MATLAB xcorr的结果,需要在前后填充带零的向量:

#include <stdio.h>

#include <Accelerate/Accelerate.h>


int main(void)
{
    #define NF  3
    #define NC  (2*NF+1)

    float A[3*NF] = {0, 0, 0, 1, 2, 3, 0, 0, 0};
    float C[NC];

    vDSP_conv(A, 1, A+NF, 1, C, 1, NC, NF);

    for (vDSP_Length i = 0; i < NC; ++i)
        printf("C[%u] = %g.\n", (unsigned) i, C[i]);

    return 0;
}