matlab中hamming窗的 c/c++ 版本的实现

matlab中的hamming窗调用：

>> hamming(5)
ans =
    0.0800
    0.5400
    1.0000
    0.5400
    0.0800

原理: hammin 窗的公式：

当a0 = 0.54时，即为hamming窗；
当a0 = 0.5 时，即为hanning窗

故其对应的 c/c++ 代码为:

//hamming窗 : a = 0.54;   hann窗 : a = 0.5
void hamming(float* win, int len)
{
	float a = 0.54; // 0.504是hamming   0.5是hann
	for (size_t i = 0; i < len; i++)
	{
		win[i] = a - (1.0f - a) * cos(2 * PI * i / (len - 1));
	}
}

原理介绍：

当 a0=0.53836 ，称作 Hamming窗；当 a0=0.5 则叫作 Hann窗。

• Hann窗又称升余弦窗。Hann窗可以看作是3个矩形时间窗的频谱之和，或者说是 3个 sinc(t) 型函数之和，而括号中的两项相对于第一个谱窗向左、右各移动了π/T，从而使旁瓣互相抵消，消去高频干扰和漏能。
  从减小泄漏观点出发，Hann窗优于矩形窗。但Hann窗主瓣加宽，相当于分析带宽加宽，频率分辨力下降。

• 如果我们将a0设为接近 0.53836 的数值，或是更精确来说是 25/46，便会得到Hamming窗，而设定这个数值的用意，是在频率为 5π/(N − 1) 处产生零交会处(zero-crossing)，使原先Hann窗的第一个旁瓣(sidelobe)可以被大幅消除，产生只有Hann窗 1/5 高度的旁瓣。

• 一般来说，为了减少运算量，可以把公式精简为：
  

本文参考： (如有侵权，请联系作者)

1. 什么是汉明窗？加Hanmming窗的作用？：https://www.cnblogs.com/lxy2017/p/4049124.html
2. matlab矩形窗代码_Hamming(汉明)窗的原理介绍及实例解析：https://blog.csdn.net/weixin_36260304/article/details/113489274