主要内容GydF4y2Ba

iscola.GydF4y2Ba

确定窗口重叠组合是否符合ColaGydF4y2Ba

描述GydF4y2Ba

例子GydF4y2Ba

TF.GydF4y2Ba= iscola (GydF4y2Ba窗户GydF4y2Ba那GydF4y2BanoverlapGydF4y2Ba的)GydF4y2Ba检查指定的窗口和重叠是否满足GydF4y2Ba恒定重叠 - 添加(可乐)约束GydF4y2Ba确保GydF4y2Ba逆短时傅里叶变换GydF4y2Ba结果对未修饰的光谱进行了完美的重建。GydF4y2Ba

TF.GydF4y2Ba= iscola (GydF4y2Ba窗户GydF4y2Ba那GydF4y2BanoverlapGydF4y2Ba那GydF4y2Ba方法GydF4y2Ba的)GydF4y2Ba指定要使用的反转方法。GydF4y2Ba

[GydF4y2BaTF.GydF4y2Ba那GydF4y2BaM.GydF4y2Ba] = iscola(GydF4y2Ba___GydF4y2Ba的)GydF4y2Ba也返回COLA求和的中位数。可以将这些输出参数与前面的任何输入语法一起使用。GydF4y2Ba

例子GydF4y2Ba

[GydF4y2BaTF.GydF4y2Ba那GydF4y2BaM.GydF4y2Ba那GydF4y2BamaxDeviationGydF4y2Ba] = iscola(GydF4y2Ba___GydF4y2Ba的)GydF4y2Ba返回与中值的最大偏差GydF4y2BaM.GydF4y2Ba.GydF4y2Ba

例子GydF4y2Ba

全部折叠GydF4y2Ba

创建长度为120的周期性根- hann窗口。测试窗口是否符合COLA,重叠50%。GydF4y2Ba

WIN = SQRT(HANN(120,GydF4y2Ba'定期'GydF4y2Ba));noverlap = 60;GydF4y2Ba

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

tf = iscola(win,noverlap)GydF4y2Ba
tf =GydF4y2Ba逻辑GydF4y2Ba1GydF4y2Ba

创建长度为256的周期性汉明窗口。设置重叠添加的方法为GydF4y2Ba“ola”GydF4y2Ba.GydF4y2Ba

窗口=汉明(256,GydF4y2Ba'定期'GydF4y2Ba);方法=GydF4y2Ba“ola”GydF4y2Ba;noverlap = 128;GydF4y2Ba

测试窗口是否符合COLA,重叠50%。还要计算COLA求和的中位数和与该求和的最大偏差。GydF4y2Ba

[tf,m,maxdeviation] = iscra(窗口,noverlap,方法)GydF4y2Ba
tf =GydF4y2Ba逻辑GydF4y2Ba1GydF4y2Ba
m = 1.0800.GydF4y2Ba
maxDeviation = 2.2204 e-16GydF4y2Ba

输入参数GydF4y2Ba

全部折叠GydF4y2Ba

分析窗口,指定为矢量。GydF4y2Ba

例子:GydF4y2BaWin = Bartlett(120)GydF4y2Ba是一个长度120的巴特利特窗口。GydF4y2Ba

数据类型:GydF4y2Ba双GydF4y2Ba|GydF4y2Ba单身的GydF4y2Ba

重叠样本的数目,指定为小于的长度的正整数GydF4y2Ba窗户GydF4y2Ba.GydF4y2Ba

数据类型:GydF4y2Ba双GydF4y2Ba|GydF4y2Ba单身的GydF4y2Ba

重叠的方法 - 添加,指定为:GydF4y2Ba

  • 'wola'GydF4y2Ba- 加权重叠 - 添加。GydF4y2Ba

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

输出参数GydF4y2Ba

全部折叠GydF4y2Ba

COLA遵从性,作为逻辑标量返回。如果函数返回aGydF4y2Ba1GydF4y2Ba(true),然后窗口和重叠长度满足Cola约束。GydF4y2Ba

中位数的可乐求和,作为一个真正的标量返回。如果输入是符合可乐的,那么GydF4y2BaM.GydF4y2Ba等于可乐求和常数。GydF4y2Ba

与中位数的最大偏差GydF4y2BaM.GydF4y2Ba.如果GydF4y2Ba窗户GydF4y2Ba和GydF4y2BanoverlapGydF4y2Ba可乐是合规的吗GydF4y2BamaxDeviationGydF4y2Ba接近可乐求和的预期数值误差。GydF4y2Ba

笔记GydF4y2Ba

你可以得出结论,严格遵守可乐,如果GydF4y2Bam = 1GydF4y2Ba和GydF4y2BamaxDeviationGydF4y2Ba接近数值精度错误。GydF4y2Ba

更多关于GydF4y2Ba

全部折叠GydF4y2Ba

逆短时傅里叶变换GydF4y2Ba

通过拍摄STFT的每个DFT向量的IFFT和重叠添加反相信号来计算逆短时傅里叶变换。ISTFT计算如下:GydF4y2Ba

XGydF4y2Ba (GydF4y2Ba N.GydF4y2Ba 的)GydF4y2Ba =GydF4y2Ba ∫GydF4y2Ba -GydF4y2Ba 1GydF4y2Ba /GydF4y2Ba 2GydF4y2Ba 1GydF4y2Ba /GydF4y2Ba 2GydF4y2Ba σ.GydF4y2Ba M.GydF4y2Ba =GydF4y2Ba -GydF4y2Ba ∞GydF4y2Ba ∞GydF4y2Ba XGydF4y2Ba M.GydF4y2Ba (GydF4y2Ba FGydF4y2Ba 的)GydF4y2Ba E.GydF4y2Ba jGydF4y2Ba 2GydF4y2Ba πGydF4y2Ba FGydF4y2Ba N.GydF4y2Ba D.GydF4y2Ba FGydF4y2Ba =GydF4y2Ba σ.GydF4y2Ba M.GydF4y2Ba =GydF4y2Ba -GydF4y2Ba ∞GydF4y2Ba ∞GydF4y2Ba ∫GydF4y2Ba -GydF4y2Ba 1GydF4y2Ba /GydF4y2Ba 2GydF4y2Ba 1GydF4y2Ba /GydF4y2Ba 2GydF4y2Ba XGydF4y2Ba M.GydF4y2Ba (GydF4y2Ba FGydF4y2Ba 的)GydF4y2Ba E.GydF4y2Ba jGydF4y2Ba 2GydF4y2Ba πGydF4y2Ba FGydF4y2Ba N.GydF4y2Ba D.GydF4y2Ba FGydF4y2Ba =GydF4y2Ba σ.GydF4y2Ba M.GydF4y2Ba =GydF4y2Ba -GydF4y2Ba ∞GydF4y2Ba ∞GydF4y2Ba XGydF4y2Ba M.GydF4y2Ba (GydF4y2Ba N.GydF4y2Ba 的)GydF4y2Ba

在哪里GydF4y2Ba R.GydF4y2Ba 是连续DFT之间的跳尺寸,GydF4y2Ba XGydF4y2Ba M.GydF4y2Ba 是窗口数据的DFT以时间为中心GydF4y2Ba M.GydF4y2Ba R.GydF4y2Ba 和GydF4y2Ba XGydF4y2Ba M.GydF4y2Ba (GydF4y2Ba N.GydF4y2Ba 的)GydF4y2Ba =GydF4y2Ba XGydF4y2Ba (GydF4y2Ba N.GydF4y2Ba 的)GydF4y2Ba GGydF4y2Ba (GydF4y2Ba N.GydF4y2Ba -GydF4y2Ba M.GydF4y2Ba R.GydF4y2Ba 的)GydF4y2Ba .逆STFT是对原始信号的完美再现,只要GydF4y2Ba σ.GydF4y2Ba M.GydF4y2Ba =GydF4y2Ba -GydF4y2Ba ∞GydF4y2Ba ∞GydF4y2Ba GGydF4y2Ba 一种GydF4y2Ba +GydF4y2Ba 1GydF4y2Ba (GydF4y2Ba N.GydF4y2Ba -GydF4y2Ba M.GydF4y2Ba R.GydF4y2Ba 的)GydF4y2Ba =GydF4y2Ba CGydF4y2Ba ∀GydF4y2Ba N.GydF4y2Ba ∈GydF4y2Ba ℤGydF4y2Ba 在哪里GydF4y2Ba分析窗口GydF4y2Ba GGydF4y2Ba (GydF4y2Ba N.GydF4y2Ba 的)GydF4y2Ba 是用来屏蔽原始信号的GydF4y2Ba CGydF4y2Ba 是一个常数。下图描绘了重建原始信号之后的步骤。GydF4y2Ba

恒定重叠 - 添加(可乐)约束GydF4y2Ba

为了保证非修饰谱的成功重构,分析窗口必须满足COLA约束。一般情况下,如果分析窗口满足条件GydF4y2Ba σ.GydF4y2Ba M.GydF4y2Ba =GydF4y2Ba -GydF4y2Ba ∞GydF4y2Ba ∞GydF4y2Ba GGydF4y2Ba 一种GydF4y2Ba +GydF4y2Ba 1GydF4y2Ba (GydF4y2Ba N.GydF4y2Ba -GydF4y2Ba M.GydF4y2Ba R.GydF4y2Ba 的)GydF4y2Ba =GydF4y2Ba CGydF4y2Ba ∀GydF4y2Ba N.GydF4y2Ba ∈GydF4y2Ba ℤGydF4y2Ba ,窗户被认为是符合可乐。另外,可乐依从性可以分为弱遵从性和强遵从性。GydF4y2Ba

  • 弱COLA顺应性意味着分析窗口的傅立叶变换在帧率谐波下为零,因此GydF4y2Ba

    GGydF4y2Ba (GydF4y2Ba FGydF4y2Ba K.GydF4y2Ba 的)GydF4y2Ba =GydF4y2Ba 0.GydF4y2Ba 那GydF4y2Ba K.GydF4y2Ba =GydF4y2Ba 1GydF4y2Ba 那GydF4y2Ba 2GydF4y2Ba 那GydF4y2Ba ......GydF4y2Ba 那GydF4y2Ba R.GydF4y2Ba -GydF4y2Ba 1GydF4y2Ba 那GydF4y2Ba FGydF4y2Ba K.GydF4y2Ba ≜GydF4y2Ba K.GydF4y2Ba R.GydF4y2Ba .GydF4y2Ba

    别名取消受光谱修改的干扰。弱可乐依赖于频域中的别名消除。因此,只要信号没有经历任何频谱修改,可以使用弱Cola标准的窗口来实现完美的重建。GydF4y2Ba

  • 对于强的COLA遵从性,窗口的傅里叶变换必须与帧速率下采样的带宽限制一致GydF4y2Ba

    GGydF4y2Ba (GydF4y2Ba FGydF4y2Ba 的)GydF4y2Ba =GydF4y2Ba 0.GydF4y2Ba 那GydF4y2Ba FGydF4y2Ba ≥GydF4y2Ba 1GydF4y2Ba 2GydF4y2Ba R.GydF4y2Ba .GydF4y2Ba

    该等式表明,强大的可乐约束不允许别名。另外,对于强大的可乐符合性,常数的值GydF4y2Ba CGydF4y2Ba 必须等于1.一般情况下,如果以任何方式修改短时间频谱,则优选符合更强的可乐窗口。GydF4y2Ba

你可以使用GydF4y2Baiscola.GydF4y2Ba检查可乐合规性的功能。用于检查可乐合规性的总和的数量由窗口长度和跳尺寸决定。一般来说,它是常用的GydF4y2Ba 一种GydF4y2Ba =GydF4y2Ba 1GydF4y2Ba 在GydF4y2Ba σ.GydF4y2Ba M.GydF4y2Ba =GydF4y2Ba -GydF4y2Ba ∞GydF4y2Ba ∞GydF4y2Ba GGydF4y2Ba 一种GydF4y2Ba +GydF4y2Ba 1GydF4y2Ba (GydF4y2Ba N.GydF4y2Ba -GydF4y2Ba M.GydF4y2Ba R.GydF4y2Ba 的)GydF4y2Ba =GydF4y2Ba CGydF4y2Ba ∀GydF4y2Ba N.GydF4y2Ba ∈GydF4y2Ba ℤGydF4y2Ba 对于加权重叠 - 添加(WOLA),和GydF4y2Ba 一种GydF4y2Ba =GydF4y2Ba 0.GydF4y2Ba 交叠相加(OLA)。默认情况下,GydF4y2BaistftGydF4y2Ba使用wola方法,通过应用aGydF4y2Ba综合窗口GydF4y2Ba在执行重叠添加方法之前。GydF4y2Ba

一般情况下,合成窗口与分析窗口是相同的。通过对强OLA窗口取平方根,可以构建有用的WOLA窗口。您可以对所有非负OLA窗口使用此方法。例如,根- hann窗口就是一个很好的WOLA窗口示例。GydF4y2Ba

完美的重建GydF4y2Ba

一般来说,计算输入信号的短时傅立叶变换并对其进行反变换并不能得到完美的重构。如果想要ISTFT的输出与原始输入信号尽可能的匹配,信号和窗口必须满足以下条件:GydF4y2Ba

  • 输入大小 - 如果逆变输出GydF4y2BastGydF4y2Ba使用GydF4y2BaistftGydF4y2Ba并且希望结果与输入信号相同GydF4y2BaXGydF4y2Ba, 的价值GydF4y2Ba K.GydF4y2Ba =GydF4y2Ba (GydF4y2Ba L.GydF4y2Ba E.GydF4y2Ba N.GydF4y2Ba GGydF4y2Ba T.GydF4y2Ba HGydF4y2Ba (GydF4y2Ba XGydF4y2Ba 的)GydF4y2Ba -GydF4y2Ba N.GydF4y2Ba O.GydF4y2Ba V.GydF4y2Ba E.GydF4y2Ba R.GydF4y2Ba L.GydF4y2Ba 一种GydF4y2Ba P.GydF4y2Ba 的)GydF4y2Ba (GydF4y2Ba L.GydF4y2Ba E.GydF4y2Ba N.GydF4y2Ba GGydF4y2Ba T.GydF4y2Ba HGydF4y2Ba (GydF4y2Ba W.GydF4y2Ba 一世GydF4y2Ba N.GydF4y2Ba D.GydF4y2Ba O.GydF4y2Ba W.GydF4y2Ba 的)GydF4y2Ba -GydF4y2Ba N.GydF4y2Ba O.GydF4y2Ba V.GydF4y2Ba E.GydF4y2Ba R.GydF4y2Ba L.GydF4y2Ba 一种GydF4y2Ba P.GydF4y2Ba 的)GydF4y2Ba 必须是整数。GydF4y2Ba

  • 可乐合规则 - 使用可调整的窗口,假设您没有修改信号的短时傅立叶变换。GydF4y2Ba

  • 填充-如果输入信号的长度等于GydF4y2BaK.GydF4y2Ba在计算短时傅里叶变换之前,不是整数,零键盘信号。在反转信号后取出额外的零。GydF4y2Ba

参考文献GydF4y2Ba

[1] j·B·艾伦用离散傅里叶变换进行短期光谱分析、合成和修正。GydF4y2Ba声学,语音和信号处理的IEEE交易。GydF4y2Ba1977年6月,第25卷第3期,第235-238页。GydF4y2Ba

[2]格里芬,丹尼尔W.和Jae S. LIM。“从修改的短时傅里叶变换的信号估计。”GydF4y2Ba声学,语音和信号处理的IEEE交易。GydF4y2Ba第32卷第2期,1984年4月,236-243页。GydF4y2Ba

[3]夏普,布鲁斯。GydF4y2Ba重叠叠加处理的可逆性GydF4y2Ba.GydF4y2Bahttps://gauss256.github.io/blog/cola.htmlGydF4y2Ba已于2019年7月通过。GydF4y2Ba

史密斯,我是朱利叶斯·奥里翁。GydF4y2Ba光谱音频信号处理GydF4y2Ba.GydF4y2Bahttps://ccrma.stanford.edu/~jos/sasp/GydF4y2Ba,在线图书,2011版,访问2018年11月。GydF4y2Ba

扩展能力GydF4y2Ba

C / c++代码生成GydF4y2Ba
使用MATLAB®Coder™生成C和c++代码。GydF4y2Ba

也可以看看GydF4y2Ba

职能GydF4y2Ba

在R2019A引入GydF4y2Ba