taylorwin(n,nbar,sll)窗函数 c/c++实现

这是一个实现 taylorwin(n,nbar,sll) 窗函数的 C/C++ 代码：

#include <cmath>

void taylorwin(int n, int nbar, double sll, double *w)
{
    double alpha = nbar * sll;
    double beta = std::sqrt(1 - std::pow(sll, 2));

    for (int i = 0; i < n; i++)
    {
        double x = (2 * i - n + 1) / (double)n;
        double ak = 0;
        for (int k = 0; k < nbar; k++)
        {
            ak += std::pow(-1, k) * std::pow(beta, 2 * k) / (std::tgamma(k + 1) * std::tgamma(nbar + k + 1)) * std::cos((2 * k + 1) * std::acos(x));
        }
        w[i] = ak / std::pow(std::pow(x, 2) - 1, nbar / 2.0);
    }
}

该函数接受四个参数：n、nbar、sll 和一个指向 double 类型数组的指针 w。n 表示窗口长度，nbar 是 taylor 系数的数量，sll 是过渡带衰减的 dB 数量，w 是输出的窗口数组。

该函数的实现基于 taylorwin.m 文件中的算法。在实现中，我们先计算了 alpha 和 beta，然后对于每个输出点，计算了对应的 taylor 系数 ak，最后计算了窗口值。具体的计算细节可以参考注释。