用C++实现matlab中的finddelay函数

Matlab中的finddelay函数用于计算两个时间序列之间的延迟，可以通过以下C++代码实现：

#include <iostream>
#include <vector>
#include <algorithm>
#include <cmath>

using namespace std;

int finddelay(const vector<double>& x, const vector<double>& y)
{
    int len1 = x.size();
    int len2 = y.size();
    
    // 计算互相关函数
    vector<double> corr(len1 + len2 - 1);
    for(int i = 0; i < len1 + len2 - 1; i++)
    {
        double sum = 0;
        for(int j = max(0, i - len2 + 1); j <= min(i, len1 - 1); j++)
        {
            sum += x[j] * y[i - j];
        }
        corr[i] = sum;
    }
    
    // 找到最大值的位置
    int max_index = 0;
    double max_value = corr[0];
    for(int i = 1; i < corr.size(); i++)
    {
        if(corr[i] > max_value)
        {
            max_value = corr[i];
            max_index = i;
        }
    }
    
    // 计算延迟
    int delay = max_index - len2 + 1;
    return delay;
}

int main()
{
    vector<double> x = {1, 2, 3, 4, 5};
    vector<double> y = {4, 5, 1, 2, 3};
    
    int delay = finddelay(x, y);
    cout << "Delay: " << delay << endl;
    
    return 0;
}

该代码定义了一个名为`finddelay`的函数，输入为两个`vector<double>`类型的时间序列，输出为它们之间的延迟。其中，计算互相关函数部分使用了嵌套循环，计算延迟部分使用了最大值的位置来计算。

在主函数中，我们定义了两个时间序列，并调用了`finddelay`函数来计算它们之间的延迟。最后输出结果。