在 Matlab 中创建长度为 64 且重叠率为 60% 的汉明窗答案

【问题标题】：To create Hamming window of length 64 with overlap 60% in Matlab在 Matlab 中创建长度为 64 且重叠率为 60% 的汉明窗
【发布时间】：2013-12-24 16:07:24
【问题描述】：

我正在尝试提供长度为 64 且重叠率为 60% 的 Wigner-Ville Distribution Hamming Window here。我可以通过

创建长度为 64 的汉明窗

h=hamming(64);

Here 是关于这个问题的一些理论建议。该窗口似乎是三个汉明波的某种convolution，卷积概率为 60%。

重叠似乎是三个函数的某种卷积。我尝试三个窗口及其重叠

conv(conv(hamming(64), hamming(64)), conv(hamming(64), hamming(64)))

我尝试两个窗口及其重叠

h = conv(hamming(64), hamming(64));

这两个结果似乎都没有给我更好的 Wigner-Ville 分布结果。许多云峰仍然可见。所以关键似乎在时间上分开了窗口，因为窗口函数的当前结果返回与 hamming(64) 窗口完全相同的图片。

考虑 60% 的重叠

hamming(64) 的尺寸为 64x1 双倍，而conv(hamming(64), hamming(64)) 的尺寸为 127x1 双倍。制作一个 60% 几率的概率算法并不简单，因为我们不能线性地迭代这两个函数。

如何创建具有 60% 重叠的汉明窗？

【问题讨论】：

根据您的“理论建议”链接，我认为重叠来自于选择跳跃长度（换句话说，如果连续窗口之间有一些重叠，则只能重叠；如果你的窗口是 64 长，你可能需要 36 跳才能在连续样本之间获得 60% 的重叠）。然后你可以不理会汉明窗本身。抱歉，我没有时间更深入地研究这个（并且现在没有方便的 matlab 来测试它）。
@Floris 请根据您的评论查看我的回答。你是这个意思吗？感谢您的精彩评论！

标签： matlab signal-processing time-frequency hamming-window

【解决方案1】：

提示buffer 的快速搜索可能值得一试。

h = hamming(64);
y = buffer(h, 1, floor(64 * 0.6));

但是我的Matlab版本不支持这个功能，所以没有尝试。

【讨论】：

我认为参数的顺序应该颠倒。重叠的数量排在最后，例如 y = buffer(h, floor(64 * 0.6), 1);
我不确定结果是否有效。传递到 Wigner-Ville 分布时，情况似乎是一样的。
缓冲区的尺寸是 38x2 双倍，而 hamming(64) 是 64x1 双倍。我不确定缓冲区是否通过遗漏 40% 的点并将列数乘以 2 来做正确的事情。大概，卷积现在应该连接到这个东西了。

【解决方案2】：

基于弗洛里斯的评论。运行 Wigner-Ville 分布，Hamming 窗口 64 有 60% 的重叠。

h = hamming(64);
h2 = hamming(38);
h = conv(h, h2);
[B,T,F] = tfrwv(data, 1:length(data), length(data), h);

图片似乎与 hamming-64 窗口完全相同。图片不应该是一样的，因为向量 hamming-64 和这个窗口函数的值不同。所以可能应该研究规范来估计这个东西。

log(abs(data)) 在左侧给出，而原始在右侧

这里现在将对数函数应用于使用 Hamming 64 的分布和使用 Hamming 64 但重叠 60% 的分布

对数函数后的图片看起来也一样。

【讨论】：

不明白为什么和这里的图片相比变化这么小stackoverflow.com/questions/20750660/…图片好像一模一样。
你比较了两个频谱图的大小吗？也许尝试 log(abs(B)) 来观察更多细节。
@lennon310 取对数时图片变黄。应该是这样吗？
玛西，它们还是一模一样吗？我的意思是 tfrwv 和频谱图的幅度对数？
似乎 h 参数根本不起作用（确保每次都清除所有内容，然后再使用另一个配置重新运行程序..）。您可能需要向他们询问这个问题......或者如果您可以在 tfrwv.m 上获得原始代码，您可以调查 h 的使用方式。