stft_cola_check — ACOUSTICS transform op

• 数据种类:无 → table(仅由参数决定的算子 —— 不接受图像或数据输入)

• 调用:import acoustics; acoustics.stft_cola_check(window='hann', win=256, hop=None)(或 opsacoustics.get("stft_cola_check"))

用法

这组(窗, 跳距)是否满足 COLA,以及满足得有多精确。

> 以下的详细说明为原文 —— 摘要与标题已翻译。

COLA — the analysis windows summing to a constant over the hop lattice — is

what lets plain overlap-add work without a division. It is *not* required by

:func:istft, which is weighted, but it is required by anything that

overlap-adds modified frames without renormalising, and getting it wrong

produces a periodic amplitude ripple at `rate/hop` Hz that looks like

tremolo rather than like a bug.

Returns a dict: `cola (bool), constant` (the mean of the overlap sum),

`max_deviation (absolute), relative_deviation, nola` (bool),

`min_squared_sum`, plus the geometry.

Measured (periodic windows, `relative_deviation` of the plain sum):

======== ==== ==== ================== ======== ====

window win hop relative_deviation constant COLA

======== ==== ==== ================== ======== ====

hann 256 128 4.44e-16 1.00 yes

hann 256 64 2.22e-16 2.00 yes

hann 256 85 1.48e-03 1.506 no

hamming 256 128 2.06e-16 1.08 yes

blackman 256 128 1.91e-01 0.84 no

blackman 256 64 3.97e-16 1.68 yes

boxcar 256 128 0.00e+00 2.00 yes

======== ==== ==== ================== ======== ====

The two blackman rows are the useful ones: the same window is COLA at

hop = win/4 and 19 % off at hop = win/2, so "which window" is not the

question — the pair is. The boxcar row was worth measuring rather than

assuming: a rectangular window at 50 % overlap sums to exactly 2 and is

COLA, which is the opposite of the usual intuition about it. Note also that

the constant is not 1 in general — an overlap-add that does not divide by it

is off by a *gain*, which is the failure that looks like a working system.

Raises `ValueError: unknown / all-zero window, hop` outside

`[1, win], win outside [2, MAX_WINDOW]`.

详细使用指南

• acoustic_condition_monitoring 族使用指南

参考(示例数据・文献)

• 示例数据目录(下载 URL / 许可证) —— 2-D 用 skimage.data(BSD/公有领域)加合成图,3-D 给出真实数据源(Stanford/PDS 等)的下载 URL。

• 算子来历与参考文献 —— 该算子族所依据的研究/方法出处。

• 算法的正典(作者・年份)与用途见上面的族使用指南。

可运行的示例(实际调用该算子并已验证的样例)

• acoustic_condition_monitoring — py -3.11 examples/acoustic_condition_monitoring.py

类型可衔接的下一个算子(可接受 table 作为输入)

istft

同类别(transform)

stft · istft

*Provenance: acoustics.py — ACOUSTICS 算子登记表。本条目由 tools/opdocs.py md 自动生成(请勿手工编辑)。*

© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.