matlib中的xcorr函数的c语言版本

以下是基于C语言实现的xcorr函数，该函数计算两个向量之间的互相关。

#include <stdio.h>
#include <stdlib.h>
#include <math.h>

double* xcorr(double* x, int nx, double* y, int ny, int maxlag) {
    int len = 2 * maxlag + 1;
    double* r = (double*) calloc(len, sizeof(double));

    for (int lag = -maxlag; lag <= maxlag; lag++) {
        double sum = 0.0;
        int i, j;

        if (lag >= 0) {
            i = lag;
            j = 0;
        } else {
            i = 0;
            j = abs(lag);
        }

        for (; i < nx && j < ny; i++, j++) {
            sum += x[i] * y[j];
        }
        r[lag + maxlag] = sum;
    }

    return r;
}

int main() {
    double x[] = {1, 2, 3, 4, 5};
    double y[] = {5, 4, 3, 2, 1};
    int nx = sizeof(x) / sizeof(double);
    int ny = sizeof(y) / sizeof(double);
    int maxlag = 2;

    double* r = xcorr(x, nx, y, ny, maxlag);

    for (int i = 0; i < 2 * maxlag + 1; i++) {
        printf("%f ", r[i]);
    }
    printf("\n");

    free(r);
    return 0;
}

该函数接受两个向量x和y，以及它们的长度nx和ny。maxlag参数指定要计算的最大滞后。函数返回一个长度为2*maxlag+1的数组，其中包含从-lag到+lag的所有滞后值的相关性。

该算法使用两个指针i和j，分别在x和y中移动，以便计算两个向量之间的点积。循环变量lag从-maxlag到+maxlag循环，每次计算一个相关性值。最后，将所有相关性值存储在r数组中并返回它。