主要内容gydF4y2Ba

iscolagydF4y2Ba

确定window-overlap组合是可乐兼容gydF4y2Ba

所有的页面崩溃gydF4y2Ba

语法gydF4y2Ba

noverlap tf = iscola(窗口)gydF4y2Ba
tf = iscola(窗口、noverlap方法)gydF4y2Ba
(tf, m) = iscola (gydF4y2Ba___gydF4y2Ba)gydF4y2Ba
[tf, m, maxDeviation] = iscola (gydF4y2Ba___gydF4y2Ba)gydF4y2Ba

描述gydF4y2Ba

例子gydF4y2Ba

特遣部队gydF4y2Ba= iscola (gydF4y2Ba窗口gydF4y2Ba,gydF4y2BanoverlapgydF4y2Ba)gydF4y2Ba检查满足指定的窗口和重叠gydF4y2Ba不断交叠相加(可乐)的约束gydF4y2Ba以确保gydF4y2Ba短时傅里叶反变换gydF4y2Ba结果完美重建nonmodified光谱。gydF4y2Ba

特遣部队gydF4y2Ba= iscola (gydF4y2Ba窗口gydF4y2Ba,gydF4y2BanoverlapgydF4y2Ba,gydF4y2Ba方法gydF4y2Ba)gydF4y2Ba指定要使用的反演方法。gydF4y2Ba

(gydF4y2Ba特遣部队gydF4y2Ba,gydF4y2Ba米gydF4y2Ba)= iscola (gydF4y2Ba___gydF4y2Ba)gydF4y2Ba还返回可乐求和的值。您可以使用这些输入输出参数与任何以前的语法。gydF4y2Ba

例子gydF4y2Ba

(gydF4y2Ba特遣部队gydF4y2Ba,gydF4y2Ba米gydF4y2Ba,gydF4y2BamaxDeviationgydF4y2Ba)= iscola (gydF4y2Ba___gydF4y2Ba)gydF4y2Ba返回的最大偏差值gydF4y2Ba米gydF4y2Ba。gydF4y2Ba

例子gydF4y2Ba

全部折叠gydF4y2Ba

打开生活的脚本gydF4y2Ba

创建一个周期性root-Hann窗口长度为120。测试窗口是否可乐符合50%的重叠。gydF4y2Ba

赢得=√损害(120gydF4y2Ba“周期”gydF4y2Ba));noverlap = 60;gydF4y2Ba

检查窗口是否可乐符合50%的重叠。gydF4y2Ba

tf = iscola(赢,noverlap)gydF4y2Ba
tf =gydF4y2Ba逻辑gydF4y2Ba1gydF4y2Ba
打开生活的脚本gydF4y2Ba

创建一个周期256年汉明窗的长度。设置交叠相加的方法gydF4y2Ba“ola”gydF4y2Ba。gydF4y2Ba

窗口=汉明(256gydF4y2Ba“周期”gydF4y2Ba);方法=gydF4y2Ba“ola”gydF4y2Ba;noverlap = 128;gydF4y2Ba

测试窗口是否可乐符合50%的重叠。也算可乐求和的值和最大偏差求和。gydF4y2Ba

[tf, m, maxDeviation] = iscola(窗口、noverlap方法)gydF4y2Ba
tf =gydF4y2Ba逻辑gydF4y2Ba1gydF4y2Ba
m = 1.0800gydF4y2Ba
maxDeviation = 2.2204 e-16gydF4y2Ba

输入参数gydF4y2Ba

全部折叠gydF4y2Ba

分析窗口中,指定为一个向量。gydF4y2Ba

例子:gydF4y2Ba赢得=巴特利特(120)gydF4y2Ba是120年Bartlett窗口长度。gydF4y2Ba

数据类型:gydF4y2Ba双gydF4y2Ba|gydF4y2Ba单gydF4y2Ba

重叠的样本数量,指定为一个正整数的长度小于gydF4y2Ba窗口gydF4y2Ba。gydF4y2Ba

数据类型:gydF4y2Ba双gydF4y2Ba|gydF4y2Ba单gydF4y2Ba

交叠相加的方法,指定为:gydF4y2Ba

  • “呼啦”gydF4y2Ba——加权交叠相加。gydF4y2Ba

  • “ola”gydF4y2Ba——交叠相加。gydF4y2Ba

输出参数gydF4y2Ba

全部折叠gydF4y2Ba

可乐的依从性,作为一个逻辑标量返回。如果函数返回gydF4y2Ba1gydF4y2Ba(真正的),然后和重叠的窗口长度满足可乐约束。gydF4y2Ba

可乐求和的值,作为一个真正的标量返回。如果输入是可乐兼容,那么gydF4y2Ba米gydF4y2Ba等于可乐总和不变。gydF4y2Ba

最大偏差值gydF4y2Ba米gydF4y2Ba。如果gydF4y2Ba窗口gydF4y2Ba和gydF4y2BanoverlapgydF4y2Ba可乐兼容,是吗gydF4y2BamaxDeviationgydF4y2Ba接近预期的数值精度误差的可乐求和。gydF4y2Ba

请注意gydF4y2Ba

你可以强烈COLA-compliance如果结束gydF4y2Bam = 1gydF4y2Ba和gydF4y2BamaxDeviationgydF4y2Ba接近的数值精度误差。gydF4y2Ba

更多关于gydF4y2Ba

全部折叠gydF4y2Ba

短时傅里叶反变换gydF4y2Ba

逆短时傅里叶变换是通过传输线的DFT计算向量的STFT和反向信号重叠相加。ISTFT计算如下:gydF4y2Ba

xgydF4y2Ba (gydF4y2Ba ngydF4y2Ba )gydF4y2Ba =gydF4y2Ba ∫gydF4y2Ba −gydF4y2Ba 1gydF4y2Ba /gydF4y2Ba 2gydF4y2Ba 1gydF4y2Ba /gydF4y2Ba 2gydF4y2Ba ∑gydF4y2Ba 米gydF4y2Ba =gydF4y2Ba −gydF4y2Ba ∞gydF4y2Ba ∞gydF4y2Ba XgydF4y2Ba 米gydF4y2Ba (gydF4y2Ba fgydF4y2Ba )gydF4y2Ba egydF4y2Ba jgydF4y2Ba 2gydF4y2Ba πgydF4y2Ba fgydF4y2Ba ngydF4y2Ba dgydF4y2Ba fgydF4y2Ba =gydF4y2Ba ∑gydF4y2Ba 米gydF4y2Ba =gydF4y2Ba −gydF4y2Ba ∞gydF4y2Ba ∞gydF4y2Ba ∫gydF4y2Ba −gydF4y2Ba 1gydF4y2Ba /gydF4y2Ba 2gydF4y2Ba 1gydF4y2Ba /gydF4y2Ba 2gydF4y2Ba XgydF4y2Ba 米gydF4y2Ba (gydF4y2Ba fgydF4y2Ba )gydF4y2Ba egydF4y2Ba jgydF4y2Ba 2gydF4y2Ba πgydF4y2Ba fgydF4y2Ba ngydF4y2Ba dgydF4y2Ba fgydF4y2Ba =gydF4y2Ba ∑gydF4y2Ba 米gydF4y2Ba =gydF4y2Ba −gydF4y2Ba ∞gydF4y2Ba ∞gydF4y2Ba xgydF4y2Ba 米gydF4y2Ba (gydF4y2Ba ngydF4y2Ba )gydF4y2Ba

在哪里gydF4y2Ba RgydF4y2Ba 是跳连续阶之间的大小,gydF4y2Ba XgydF4y2Ba 米gydF4y2Ba 的DFT窗口的数据集中的时间呢gydF4y2Ba 米gydF4y2Ba RgydF4y2Ba 和gydF4y2Ba xgydF4y2Ba 米gydF4y2Ba (gydF4y2Ba ngydF4y2Ba )gydF4y2Ba =gydF4y2Ba xgydF4y2Ba (gydF4y2Ba ngydF4y2Ba )gydF4y2Ba ggydF4y2Ba (gydF4y2Ba ngydF4y2Ba −gydF4y2Ba 米gydF4y2Ba RgydF4y2Ba )gydF4y2Ba 。逆STFT是一个完美的,只要重建原始信号gydF4y2Ba ∑gydF4y2Ba 米gydF4y2Ba =gydF4y2Ba −gydF4y2Ba ∞gydF4y2Ba ∞gydF4y2Ba ggydF4y2Ba 一个gydF4y2Ba +gydF4y2Ba 1gydF4y2Ba (gydF4y2Ba ngydF4y2Ba −gydF4y2Ba 米gydF4y2Ba RgydF4y2Ba )gydF4y2Ba =gydF4y2Ba cgydF4y2Ba ∀gydF4y2Ba ngydF4y2Ba ∈gydF4y2Ba ℤgydF4y2Ba 在哪里gydF4y2Ba分析窗口gydF4y2Ba ggydF4y2Ba (gydF4y2Ba ngydF4y2Ba )gydF4y2Ba 用于窗口原始信号和gydF4y2Ba cgydF4y2Ba 是一个常数。下图描述了重构原始信号之后的步骤。gydF4y2Ba

不断交叠相加(可乐)的约束gydF4y2Ba

确保成功的重建nonmodified光谱,分析窗口必须满足可乐约束。一般来说,如果满足条件的分析窗口gydF4y2Ba ∑gydF4y2Ba 米gydF4y2Ba =gydF4y2Ba −gydF4y2Ba ∞gydF4y2Ba ∞gydF4y2Ba ggydF4y2Ba 一个gydF4y2Ba +gydF4y2Ba 1gydF4y2Ba (gydF4y2Ba ngydF4y2Ba −gydF4y2Ba 米gydF4y2Ba RgydF4y2Ba )gydF4y2Ba =gydF4y2Ba cgydF4y2Ba ∀gydF4y2Ba ngydF4y2Ba ∈gydF4y2Ba ℤgydF4y2Ba ,被认为是COLA-compliant的窗口。另外,可乐合规可以被描述为或弱或强。gydF4y2Ba

  • 弱可乐遵从性暗示的傅里叶变换分析窗口有0的帧频谐波等gydF4y2Ba

    GgydF4y2Ba (gydF4y2Ba fgydF4y2Ba kgydF4y2Ba )gydF4y2Ba =gydF4y2Ba 0gydF4y2Ba ,gydF4y2Ba kgydF4y2Ba =gydF4y2Ba 1gydF4y2Ba ,gydF4y2Ba 2gydF4y2Ba ,gydF4y2Ba …gydF4y2Ba ,gydF4y2Ba RgydF4y2Ba −gydF4y2Ba 1gydF4y2Ba ,gydF4y2Ba fgydF4y2Ba kgydF4y2Ba ≜gydF4y2Ba kgydF4y2Ba RgydF4y2Ba 。gydF4y2Ba

    别名取消由光谱干扰的修改。弱可乐依靠别名取消在频域。因此,完美的重建是可能使用弱COLA-compliant windows只要信号没有经历任何光谱的修改。gydF4y2Ba

  • 强劲可乐合规,窗口的傅里叶变换必须与将采样带宽受限的一贯的帧率等gydF4y2Ba

    GgydF4y2Ba (gydF4y2Ba fgydF4y2Ba )gydF4y2Ba =gydF4y2Ba 0gydF4y2Ba ,gydF4y2Ba fgydF4y2Ba ≥gydF4y2Ba 1gydF4y2Ba 2gydF4y2Ba RgydF4y2Ba 。gydF4y2Ba

    这个方程表明,强劲的可乐不允许混叠的约束。此外,强劲可乐合规,常量的值gydF4y2Ba cgydF4y2Ba 必须等于1。一般来说,如果短时谱以任何方式修改,更可乐的窗口是首选。gydF4y2Ba

您可以使用gydF4y2BaiscolagydF4y2Ba函数来检查弱可乐合规。合计的数量用于检查可乐合规是由窗口长度和大小。一般来说,使用是很常见的gydF4y2Ba 一个gydF4y2Ba =gydF4y2Ba 1gydF4y2Ba 在gydF4y2Ba ∑gydF4y2Ba 米gydF4y2Ba =gydF4y2Ba −gydF4y2Ba ∞gydF4y2Ba ∞gydF4y2Ba ggydF4y2Ba 一个gydF4y2Ba +gydF4y2Ba 1gydF4y2Ba (gydF4y2Ba ngydF4y2Ba −gydF4y2Ba 米gydF4y2Ba RgydF4y2Ba )gydF4y2Ba =gydF4y2Ba cgydF4y2Ba ∀gydF4y2Ba ngydF4y2Ba ∈gydF4y2Ba ℤgydF4y2Ba 加权交叠相加(呼啦)gydF4y2Ba 一个gydF4y2Ba =gydF4y2Ba 0gydF4y2Ba 交叠相加(OLA)。默认情况下,gydF4y2BaistftgydF4y2Ba使用呼啦方法,通过应用gydF4y2Ba合成窗口gydF4y2Ba在执行之前交叠相加的方法。gydF4y2Ba

一般来说,合成窗口是一样的分析窗口。你可以构造有用呼啦windows以√OLA窗口。您可以使用这种方法对所有非负OLA窗口。例如,root-Hann呼啦窗口的窗口是一个很好的例子。gydF4y2Ba

完美的重建gydF4y2Ba

一般来说,计算输入信号的STFT和反相这不会导致完美的重建。如果你想让ISTFT的输出尽可能匹配原始输入信号,信号和窗口必须满足下列条件:gydF4y2Ba

  • 输入大小——如果你反转的输出gydF4y2BastftgydF4y2Ba使用gydF4y2BaistftgydF4y2Ba和想要的结果是相同的长度作为输入信号gydF4y2BaxgydF4y2Ba的价值,gydF4y2Ba

    kgydF4y2Ba =gydF4y2Ba NgydF4y2Ba xgydF4y2Ba −gydF4y2Ba lgydF4y2Ba 米gydF4y2Ba −gydF4y2Ba lgydF4y2Ba

    必须是一个整数。的方程,gydF4y2BaNgydF4y2BaxgydF4y2Ba是信号的长度,gydF4y2Ba米gydF4y2Ba窗口的长度,gydF4y2BalgydF4y2Ba重叠长度。gydF4y2Ba

  • 可乐合规使用COLA-compliant windows,假设你没有修改了短时傅里叶变换的信号。gydF4y2Ba

  • 填充,如果输入信号的长度,这样的价值gydF4y2BakgydF4y2Ba不是一个整数,计算短时傅里叶变换在信号之前。反相的信号后去除多余的零。gydF4y2Ba

引用gydF4y2Ba

[1]艾伦,j·B。“短期内光谱分析,合成,通过离散傅里叶变换和修改。”gydF4y2BaIEEE声学,演讲和信号处理。gydF4y2Ba第三卷。25日,1977年6月,第235 - 238页。gydF4y2Ba

[2]格里芬,丹尼尔·W。,Jae Lim。“信号短时傅里叶变换估计从修改。”gydF4y2BaIEEE声学,演讲和信号处理。gydF4y2Ba32卷,2号,1984年4月,第236 - 243页。gydF4y2Ba

[3]夏普,布鲁斯。gydF4y2Ba可逆性的交叠相加处理gydF4y2Ba。gydF4y2Bahttps://gauss256.github.io/blog/cola.htmlgydF4y2Ba2019年7月,访问。gydF4y2Ba

史密斯[4],朱利叶斯猎户座。gydF4y2Ba光谱音频信号处理gydF4y2Ba。gydF4y2Bahttps://ccrma.stanford.edu/乔斯/ sasp /gydF4y2Ba在线图书,2011年版,2018年11月访问。gydF4y2Ba

扩展功能gydF4y2Ba

C / c++代码生成gydF4y2Ba
生成C和c++代码使用MATLAB®编码器™。gydF4y2Ba

版本历史gydF4y2Ba

介绍了R2019agydF4y2Ba

另请参阅gydF4y2Ba

功能gydF4y2Ba