主要内容gydF4y2Ba

stftmag2siggydF4y2Ba

从STFT幅值重建信号gydF4y2Ba

折叠所有页面gydF4y2Ba

语法gydF4y2Ba

x = stftmag2sig(年代,nfft)gydF4y2Ba
x = stftmag2sig (s nfft fs)gydF4y2Ba
nfft x = stftmag2sig (s, t)gydF4y2Ba
x = stftmag2sig (gydF4y2Ba＿＿＿gydF4y2Ba、名称、值)gydF4y2Ba
[x, t,信息]= stftmag2sig (gydF4y2Ba＿＿＿gydF4y2Ba）gydF4y2Ba

描述gydF4y2Ba

例子gydF4y2Ba

xgydF4y2Ba= stftmag2sig (gydF4y2Ba年代gydF4y2Ba，gydF4y2BanfftgydF4y2Ba）gydF4y2Ba返回重构的时域实信号，gydF4y2BaxgydF4y2Ba，根据gydF4y2Ba短时傅里叶变换gydF4y2Ba(STFT)大小,gydF4y2Ba年代gydF4y2Ba，基于Griffin-Lim算法。函数假设gydF4y2Ba年代gydF4y2Ba采用离散傅里叶变换(DFT)的长度gydF4y2BanfftgydF4y2Ba．gydF4y2Ba

例子gydF4y2Ba

xgydF4y2Ba= stftmag2sig (gydF4y2Ba年代gydF4y2Ba，gydF4y2BanfftgydF4y2Ba，gydF4y2BafsgydF4y2Ba）gydF4y2Ba返回重构信号，假设gydF4y2Ba年代gydF4y2Ba抽样速率gydF4y2BafsgydF4y2Ba．gydF4y2Ba

xgydF4y2Ba= stftmag2sig (gydF4y2Ba年代gydF4y2Ba，gydF4y2BanfftgydF4y2Ba，gydF4y2BatsgydF4y2Ba）gydF4y2Ba返回重构信号，假设gydF4y2Ba年代gydF4y2Ba抽样时间是多少gydF4y2BatsgydF4y2Ba．gydF4y2Ba

xgydF4y2Ba= stftmag2sig (gydF4y2Ba＿＿＿gydF4y2Ba，gydF4y2Ba名称,值gydF4y2Ba）gydF4y2Ba使用名称-值对参数指定其他选项。选项包括FFT窗口和指定初始阶段的方法。这些参数可以添加到前面的任何输入语法中。例如,gydF4y2Ba“FrequencyRange”、“单向的”,“InitializePhaseMethod”、“随机”gydF4y2Ba指定信号由具有随机初始相位的单侧短时傅里叶变换重建而成。gydF4y2Ba

［gydF4y2BaxgydF4y2Ba，gydF4y2BatgydF4y2Ba，gydF4y2Ba信息gydF4y2Ba) = stftmag2sig (gydF4y2Ba＿＿＿gydF4y2Ba）gydF4y2Ba还返回重建信号的时刻以及包含重建过程信息的结构。gydF4y2Ba

例子gydF4y2Ba

全部折叠gydF4y2Ba

打开生活的脚本gydF4y2Ba

考虑一个正弦波的512个样本，其归一化频率为gydF4y2Ba πgydF4y2Ba /gydF4y2Ba 60gydF4y2Ba rad/样品，DC值为1。计算信号的短时傅立叶变换。gydF4y2Ba

n = 512;x = cos(π/ 60 * (0:n - 1)) + 1;S = stft (x);gydF4y2Ba

从STFT的大小重建正弦信号。绘制原始信号和重建信号。gydF4y2Ba

xr = stftmag2sig (abs (S)、大小(年代,1));情节(x)gydF4y2Ba在gydF4y2Ba情节(xr,gydF4y2Ba“——”gydF4y2Ba，gydF4y2Ba“线宽”gydF4y2Ba, 2)gydF4y2Ba从gydF4y2Ba传奇(gydF4y2Ba“原始”gydF4y2Ba，gydF4y2Ba“重建”gydF4y2Ba）gydF4y2Ba

图中包含一个Axis对象。Axis对象包含两个line类型的对象。这些对象表示原始的、重建的。gydF4y2Ba

重复计算，但现在用零填充信号以减少边缘效应。gydF4y2Ba

xz = circshift ([x;0 (n, 1)), n / 2);深圳= stft (xz);xr = stftmag2sig (abs(深圳),大小(深圳,1));xz = xz (n / 2 + (1: n));xr = xr (n / 2 + (1: n));情节(xz)gydF4y2Ba在gydF4y2Ba情节(xr,gydF4y2Ba“——”gydF4y2Ba，gydF4y2Ba“线宽”gydF4y2Ba, 2)gydF4y2Ba从gydF4y2Ba传奇(gydF4y2Ba“原始”gydF4y2Ba，gydF4y2Ba“重建”gydF4y2Ba）gydF4y2Ba

图中包含一个Axis对象。Axis对象包含两个line类型的对象。这些对象表示原始的、重建的。gydF4y2Ba

重复计算，但现在通过假设减少边缘效应gydF4y2BaxgydF4y2Ba是一段两倍长的信号。gydF4y2Ba

xx = cos(π/ 60 * (- n / 2: n / 2 + n - 1)”)+ 1;Sx = stft (xx);xr = stftmag2sig (abs (Sx),大小(Sx, 1));xx = xx (n / 2 + (1: n));xr = xr (n / 2 + (1: n));情节(xx)gydF4y2Ba在gydF4y2Ba情节(xr,gydF4y2Ba“——”gydF4y2Ba，gydF4y2Ba“线宽”gydF4y2Ba, 2)gydF4y2Ba从gydF4y2Ba传奇(gydF4y2Ba“原始”gydF4y2Ba，gydF4y2Ba“重建”gydF4y2Ba）gydF4y2Ba

图中包含一个Axis对象。Axis对象包含两个line类型的对象。这些对象表示原始的、重建的。gydF4y2Ba

打开生活的脚本gydF4y2Ba

加载包含两个渐减啁啾和一个宽带飞溅声的音频信号。信号采样频率为8192 Hz。绘制信号的STFT。将波形划分为128个采样段，并使用汉明窗口对这些采样段进行窗口处理。指定64个相邻段和1024个FFT点之间的重叠样本。gydF4y2Ba

负载gydF4y2Ba长条木板gydF4y2Ba泰=(0:长度(y) 1) / Fs;gydF4y2Ba要听，键入声音(y,Fs)gydF4y2Ba风=汉明(128);克拉= 64;nfft = 1024;stft (y, Fs,gydF4y2Ba“窗口”gydF4y2Ba、风能、gydF4y2Ba“OverlapLength”gydF4y2Ba克拉,gydF4y2Ba“FFTLength”gydF4y2Banfft)gydF4y2Ba

图中包含一个轴对象。标题为“短时间傅里叶变换”的axis对象包含一个类型为image的对象。gydF4y2Ba

计算STFT的幅度和相位。gydF4y2Ba

s=stft（y，Fs，gydF4y2Ba“窗口”gydF4y2Ba、风能、gydF4y2Ba“OverlapLength”gydF4y2Ba克拉,gydF4y2Ba“FFTLength”gydF4y2Ba, nfft);smag = abs (s);sph =角(s);gydF4y2Ba

根据STFT的大小重建信号。使用与计算STFT相同的参数。默认情况下,gydF4y2Bastftmag2siggydF4y2Ba将阶段初始化为零，并使用100次优化迭代。gydF4y2Ba

[x，tx，info]=stftmag2sig（smag，nfft，Fs，gydF4y2Ba“窗口”gydF4y2Ba、风能、gydF4y2Ba“OverlapLength”gydF4y2Ba克拉);gydF4y2Ba要听到，输入声音(x,Fs)gydF4y2Ba

绘制原始信号和重建信号。为了更好的比较，将重建信号向上和向右偏移。gydF4y2Ba

情节(泰,y, tx + 500 / Fs x + 1)传说(gydF4y2Ba“原始”gydF4y2Ba，gydF4y2Ba“重建”gydF4y2Ba，gydF4y2Ba“位置”gydF4y2Ba，gydF4y2Ba“最佳”gydF4y2Ba）gydF4y2Ba

图中包含一个Axis对象。Axis对象包含两个line类型的对象。这些对象表示原始的、重建的。gydF4y2Ba

输出最后两次迭代之间向收敛方向的相对改进。gydF4y2Ba

impr = info。在consistency
impr = 0.0424gydF4y2Ba

通过将优化迭代次数增加一倍并将初始阶段设置为STFT的实际阶段来改进重建。绘制原始信号和重构信号。为了更好地进行比较，绘制重建信号的负片，并将其向上向右偏移。gydF4y2Ba

[x，tx，info]=stftmag2sig（smag，nfft，Fs，gydF4y2Ba“窗口”gydF4y2Ba、风能、gydF4y2Ba“OverlapLength”gydF4y2Ba克拉,gydF4y2Ba．．．gydF4y2Ba“MaxIterations”gydF4y2Ba, 200,gydF4y2Ba“InitialPhase”gydF4y2Ba, sph);gydF4y2Ba要听到，输入声音(x,Fs)gydF4y2Ba情节(泰,y, tx + 500 / Fs - x + 1)传说(gydF4y2Ba“原始”gydF4y2Ba，gydF4y2Ba“重建”gydF4y2Ba，gydF4y2Ba“位置”gydF4y2Ba，gydF4y2Ba“最佳”gydF4y2Ba）gydF4y2Ba

图中包含一个Axis对象。Axis对象包含两个line类型的对象。这些对象表示原始的、重建的。gydF4y2Ba

输出最后两次迭代之间向收敛方向的相对改进。gydF4y2Ba

impr = info。在consistency
impr = 1.3874 e-16gydF4y2Ba

输入参数gydF4y2Ba

全部折叠gydF4y2Ba

STFT幅值，用矩阵表示。gydF4y2Ba年代gydF4y2Ba必须对应于一个单通道实值信号。gydF4y2Ba

例子:gydF4y2Ba腹肌gydF4y2Ba（gydF4y2BastftgydF4y2Ba（gydF4y2Ba罪gydF4y2Ba（π/2*（0:255）），'fftlelength'，128））gydF4y2Ba指定正弦信号的短时傅里叶变换幅度。gydF4y2Ba

例子:gydF4y2Ba腹肌gydF4y2Ba（gydF4y2BastftgydF4y2Ba（gydF4y2Ba尖声地说gydF4y2Ba1 (0:1/1e3:1, 25日,50)))gydF4y2Ba指定在1 kHz采样的啁啾的STFT幅度。gydF4y2Ba

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

DFT点的个数，指定为正整数标量。这个论点总是必需的。gydF4y2Ba

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

采样率，指定为正数值标量。gydF4y2Ba

采样时间，指定为gydF4y2Ba持续时间gydF4y2Ba标量。指定gydF4y2BatsgydF4y2Ba是否等同于设定抽样率gydF4y2BafgydF4y2Ba年代gydF4y2Ba= 1 /gydF4y2BatsgydF4y2Ba．gydF4y2Ba

例子:gydF4y2Ba秒(1)gydF4y2Ba是一个gydF4y2Ba持续时间gydF4y2Ba表示连续信号采样之间1秒时间差的标量。gydF4y2Ba

数据类型:gydF4y2Ba持续时间gydF4y2Ba

名称-值参数gydF4y2Ba

指定可选的逗号分隔的对gydF4y2Ba名称,值gydF4y2Ba参数。gydF4y2Ba的名字gydF4y2Ba参数名和gydF4y2Ba价值gydF4y2Ba为对应值。gydF4y2Ba的名字gydF4y2Ba必须出现在引号内。可以以任意顺序指定多个名称和值对参数gydF4y2Ba名称1，值1，…，名称，值gydF4y2Ba．gydF4y2Ba

例子:gydF4y2Ba“FrequencyRange”、“单向的”,“InitializePhaseMethod”、“随机”gydF4y2Ba指定信号由具有随机初始相位的单侧短时傅里叶变换重建而成。gydF4y2Ba

不一致显示选项，指定为逗号分隔对组成gydF4y2Ba“显示”gydF4y2Ba和一个逻辑值。如果该选项设置为gydF4y2Ba真正的gydF4y2Ba,然后gydF4y2Bastftmag2siggydF4y2Ba在每20次优化迭代之后显示规范化的不一致性，并在运行结束时显示停止信息。gydF4y2Ba

数据类型:gydF4y2Ba逻辑gydF4y2Ba

STFT幅值的频率范围，指定为逗号分隔对，由gydF4y2Ba“FrequencyRange”gydF4y2Ba和gydF4y2Ba“中心”gydF4y2Ba，gydF4y2Ba双侧的gydF4y2Ba,或gydF4y2Ba“单向的”gydF4y2Ba．gydF4y2Ba

  • “中心”gydF4y2Ba-治疗gydF4y2Ba年代gydF4y2Ba为双面居中STFT的大小。如果gydF4y2BanfftgydF4y2Ba是偶数,那么gydF4y2Ba年代gydF4y2Ba是在这个区间内计算出来的吗gydF4y2Ba(-gydF4y2BaπgydF4y2Ba，gydF4y2BaπgydF4y2Ba］gydF4y2Barad /样品。如果gydF4y2BanfftgydF4y2Ba是奇数,那么gydF4y2Ba年代gydF4y2Ba是在这个区间内计算出来的吗gydF4y2Ba(-gydF4y2BaπgydF4y2Ba，gydF4y2BaπgydF4y2Ba）gydF4y2Barad /样品。如果指定时间信息，则间隔为gydF4y2Ba(-gydF4y2BafgydF4y2Ba年代gydF4y2Ba，gydF4y2BafgydF4y2Ba年代gydF4y2Ba/ 2)gydF4y2Ba/单位时间和周期gydF4y2Ba(-gydF4y2BafgydF4y2Ba年代gydF4y2Ba，gydF4y2BafgydF4y2Ba年代gydF4y2Ba/ 2)gydF4y2Ba周期/单位时间，其中gydF4y2BafgydF4y2Ba年代gydF4y2Ba为抽样率。gydF4y2Ba

  • 双侧的gydF4y2Ba-治疗gydF4y2Ba年代gydF4y2Ba即在此区间内计算的双边短时傅里叶变换的幅度gydF4y2Ba[0, 2gydF4y2BaπgydF4y2Ba）gydF4y2Barad /样品。如果指定时间信息，则间隔为gydF4y2Ba[0,gydF4y2BafgydF4y2Ba年代gydF4y2Ba）gydF4y2Ba周期/单位时间。gydF4y2Ba

  • “单向的”gydF4y2Ba-治疗gydF4y2Ba年代gydF4y2Ba即单边短时傅里叶变换的大小。如果gydF4y2BanfftgydF4y2Ba是偶数,那么gydF4y2Ba年代gydF4y2Ba是在这个区间内计算出来的吗gydF4y2Ba[0,gydF4y2BaπgydF4y2Ba］gydF4y2Barad /样品。如果gydF4y2BanfftgydF4y2Ba是奇数,那么gydF4y2Ba年代gydF4y2Ba是在这个区间内计算出来的吗gydF4y2Ba[0,gydF4y2BaπgydF4y2Ba）gydF4y2Barad /样品。如果指定时间信息，则间隔为gydF4y2Ba[0,gydF4y2BafgydF4y2Ba年代gydF4y2Ba/ 2)gydF4y2Ba/单位时间和周期gydF4y2Ba[0,gydF4y2BafgydF4y2Ba年代gydF4y2Ba/ 2)gydF4y2Ba周期/单位时间，其中gydF4y2BafgydF4y2Ba年代gydF4y2Ba为抽样率。gydF4y2Ba

数据类型:gydF4y2Ba字符gydF4y2Ba|gydF4y2Ba字符串gydF4y2Ba

重构过程中的不一致性容差，指定为逗号分隔对组成gydF4y2Ba“InconsistencyTolerance”gydF4y2Ba一个正标量。时，重构进程停止gydF4y2Ba规范化的不一致gydF4y2Ba低于公差。gydF4y2Ba

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

阶段初始化，指定为逗号分隔对，由gydF4y2Ba“InitializePhaseMethod”gydF4y2Ba和gydF4y2Ba“零”gydF4y2Ba或gydF4y2Ba“随机”gydF4y2Ba．仅指定其中一个gydF4y2Ba“InitializePhaseMethod”gydF4y2Ba或gydF4y2Ba“InitialPhase”gydF4y2Ba．gydF4y2Ba

  • “零”gydF4y2Ba—函数将阶段初始化为零。gydF4y2Ba

  • “随机”gydF4y2Ba—函数将相位初始化为在间隔内均匀分布的随机数gydF4y2Ba(- - - - - -gydF4y2BaπgydF4y2Ba，gydF4y2BaπgydF4y2Ba］gydF4y2Ba．gydF4y2Ba

数据类型:gydF4y2Ba字符gydF4y2Ba|gydF4y2Ba字符串gydF4y2Ba

初始阶段，指定为逗号分隔对，包括gydF4y2Ba“InitialPhase”gydF4y2Ba和一个值域内的实数矩阵gydF4y2Ba(- - - - - -gydF4y2BaπgydF4y2Ba，gydF4y2BaπgydF4y2Ba］gydF4y2Ba．矩阵的大小必须与gydF4y2Ba年代gydF4y2Ba．仅指定其中一个gydF4y2Ba“InitializePhaseMethod”gydF4y2Ba或gydF4y2Ba“InitialPhase”gydF4y2Ba．gydF4y2Ba

例子:gydF4y2Ba角(gydF4y2BastftgydF4y2Ba（gydF4y2BarandngydF4y2Ba(1000,1)))gydF4y2Ba指定随机信号的短时傅里叶变换的相位。gydF4y2Ba

例子:gydF4y2Ba2*pi*（兰特（标准英尺（兰特（1000,1）））-1/2）gydF4y2Ba指定在间隔内均匀分布的随机相位矩阵gydF4y2Ba(- - - - - -gydF4y2BaπgydF4y2Ba，gydF4y2BaπgydF4y2Ba］gydF4y2Ba．矩阵的大小与随机信号的短时傅里叶变换相同。gydF4y2Ba

．gydF4y2Ba

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

输入时间维度，指定为逗号分隔对组成gydF4y2Ba“InputTimeDimension”gydF4y2Ba和gydF4y2Ba“acrosscolumns”gydF4y2Ba或gydF4y2Ba“downrows”gydF4y2Ba．gydF4y2Ba

  • “acrosscolumns”gydF4y2Ba-该函数假定gydF4y2Ba年代gydF4y2Ba是跨列的频率维度是沿行向下的。gydF4y2Ba

  • “downrows”gydF4y2Ba-该函数假定gydF4y2Ba年代gydF4y2Ba是沿行向下的，频率维度是跨列的。gydF4y2Ba

数据类型:gydF4y2Ba字符gydF4y2Ba|gydF4y2Ba字符串gydF4y2Ba

最大优化迭代次数，指定为逗号分隔对组成gydF4y2Ba“MaxIterations”gydF4y2Ba一个正整数标量。当迭代次数大于时，重构过程停止gydF4y2Ba“MaxIterations”gydF4y2Ba．gydF4y2Ba

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

信号重建算法，指定为逗号分隔对组成gydF4y2Ba“方法”gydF4y2Ba其中之一是:gydF4y2Ba

  • “杯子”gydF4y2Ba-原始的重建算法，由Griffin和Lim提出，在gydF4y2Ba[1]gydF4y2Ba．gydF4y2Ba

  • “fgla”gydF4y2Ba- Perraudin、Balazs和Søndergaard提出的快速Griffin-Lim算法，描述于gydF4y2Ba［2］gydF4y2Ba．gydF4y2Ba

  • “legla”gydF4y2Ba-一种由Le Roux、Kameoka、Ono和Sagayama提出的快速算法，描述于gydF4y2Ba［3］gydF4y2Ba．gydF4y2Ba

数据类型:gydF4y2Ba字符gydF4y2Ba|gydF4y2Ba字符串gydF4y2Ba

相邻段之间重叠样本的数目，指定为逗号分隔对，由gydF4y2Ba“OverlapLength”gydF4y2Ba一个正整数小于gydF4y2Ba“窗口”gydF4y2Ba．成功的信号重建需要gydF4y2Ba“OverlapLength”gydF4y2Ba匹配用于生成STFT幅值的重叠段数。如果你忽略了gydF4y2Ba“OverlapLength”gydF4y2Ba或者将其指定为空，则将其设置为小于或等于窗口长度的75%的最大整数，即默认Hann窗口的96个样本。gydF4y2Ba

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

截断订单gydF4y2Ba“legla”gydF4y2Ba更新规则，指定为逗号分隔的对，由gydF4y2Ba“截断顺序”gydF4y2Ba一个正整数。这个论点只适用于gydF4y2Ba“方法”gydF4y2Ba设置为gydF4y2Ba“legla”gydF4y2Ba并控制该方法每次迭代中更新的相位值的数量。如果不指定,gydF4y2Ba“截断顺序”gydF4y2Ba是用自适应算法确定的。gydF4y2Ba

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

快速Griffin-Lim算法的更新参数，指定为逗号分隔对组成gydF4y2Ba“UpdateParameter”gydF4y2Ba一个正标量。这个论点只适用于gydF4y2Ba“方法”gydF4y2Ba设置为gydF4y2Ba“fgla”gydF4y2Ba并指定该方法的更新规则的参数。gydF4y2Ba

数据类型:gydF4y2Ba单gydF4y2Ba|gydF4y2Ba双gydF4y2Ba
复数的支持:金宝appgydF4y2Ba是的gydF4y2Ba

光谱窗口，指定为逗号分隔对，由gydF4y2Ba“窗口”gydF4y2Ba和一个向量。成功的信号重建需要gydF4y2Ba“窗口”gydF4y2Ba来匹配用来产生STFT幅值的窗口。如果不指定窗口或将其指定为空，函数将使用长度为128的周期性Hann窗口。的长度gydF4y2Ba“窗口”gydF4y2Ba必须大于等于2。gydF4y2Ba

有关可用窗口的列表，请参见gydF4y2Ba窗户gydF4y2Ba．gydF4y2Ba

例子:gydF4y2Ba损害(128年,“周期性”)gydF4y2Ba和gydF4y2Ba(1-cos(2 *π* (128:1:1)/ 128))/ 2gydF4y2Ba两者都指定使用的默认窗口gydF4y2Bastftmag2siggydF4y2Ba．gydF4y2Ba

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

输出参数gydF4y2Ba

全部折叠gydF4y2Ba

重构的时域信号，以矢量形式返回。gydF4y2Ba

信号被重建的时间瞬间，以矢量的形式返回。gydF4y2Ba

重构过程信息，作为包含以下字段的结构返回:gydF4y2Ba

  • ExitFlaggydF4y2Ba- - - - - -终止标志。gydF4y2Ba

    • 的值gydF4y2Ba0gydF4y2Ba表示算法在达到最大迭代次数时停止。gydF4y2Ba

    • 的值gydF4y2Ba1gydF4y2Ba当算法满足相对公差时停止。gydF4y2Ba

  • NumIterationsgydF4y2Ba-迭代的总次数。gydF4y2Ba

  • 不一致gydF4y2Ba-在最后两次迭代之间朝向收敛的平均相对改进。gydF4y2Ba

  • ReconstructedPhasegydF4y2Ba-最后一次迭代的重构阶段。gydF4y2Ba

  • ReconstructedSTFTgydF4y2Ba-上次迭代时的短时傅里叶变换重构。gydF4y2Ba

更多关于gydF4y2Ba

全部折叠gydF4y2Ba

短时傅里叶变换gydF4y2Ba

短时傅里叶变换(STFT)用于分析非平稳信号的频率内容如何随时间变化。gydF4y2Ba

信号的短时傅立叶变换是通过滑动an来计算的gydF4y2Ba分析窗口gydF4y2Ba的长度gydF4y2Ba 米gydF4y2Ba 并计算加窗数据的离散傅里叶变换。窗口每隔一段时间跳过原始信号gydF4y2Ba RgydF4y2Ba 大多数窗口函数在边缘处逐渐变细以避免光谱振铃。如果重叠长度不为零gydF4y2Ba lgydF4y2Ba 时，加窗段重叠补偿窗口边缘的信号衰减。每个加窗段的DFT被添加到一个矩阵中，该矩阵包含每个时间点和频率点的幅值和相位。STFT矩阵的列数为gydF4y2Ba

kgydF4y2Ba ＝gydF4y2Ba ⌊gydF4y2Ba NgydF4y2Ba xgydF4y2Ba −gydF4y2Ba lgydF4y2Ba 米gydF4y2Ba −gydF4y2Ba lgydF4y2Ba ⌋gydF4y2Ba ，gydF4y2Ba

在哪里gydF4y2Ba NgydF4y2Ba xgydF4y2Ba 是原始信号的长度吗gydF4y2Ba xgydF4y2Ba （gydF4y2Ba ngydF4y2Ba ）gydF4y2Ba 和gydF4y2Ba⌊⌋gydF4y2Ba符号表示楼层功能。矩阵的行数等于gydF4y2BaNgydF4y2BaDFTgydF4y2Ba的DFT点的个数，对于中心变换和双边变换和gydF4y2Ba⌊gydF4y2BaNgydF4y2BaDFTgydF4y2Ba/ 2⌋+ 1gydF4y2Ba片面的转换。gydF4y2Ba

STFT矩阵为gydF4y2Ba XgydF4y2Ba （gydF4y2Ba fgydF4y2Ba ）gydF4y2Ba ＝gydF4y2Ba ［gydF4y2Ba XgydF4y2Ba 1gydF4y2Ba （gydF4y2Ba fgydF4y2Ba ）gydF4y2Ba XgydF4y2Ba 2gydF4y2Ba （gydF4y2Ba fgydF4y2Ba ）gydF4y2Ba XgydF4y2Ba 3.gydF4y2Ba （gydF4y2Ba fgydF4y2Ba ）gydF4y2Ba ⋯gydF4y2Ba XgydF4y2Ba kgydF4y2Ba （gydF4y2Ba fgydF4y2Ba ）gydF4y2Ba ］gydF4y2Ba 这样gydF4y2Ba 米gydF4y2Ba 这个矩阵的第Th元素是gydF4y2Ba

XgydF4y2Ba 米gydF4y2Ba （gydF4y2Ba fgydF4y2Ba ）gydF4y2Ba ＝gydF4y2Ba ∑gydF4y2Ba ngydF4y2Ba ＝gydF4y2Ba −gydF4y2Ba ∞gydF4y2Ba ∞gydF4y2Ba xgydF4y2Ba （gydF4y2Ba ngydF4y2Ba ）gydF4y2Ba ggydF4y2Ba （gydF4y2Ba ngydF4y2Ba −gydF4y2Ba 米gydF4y2Ba RgydF4y2Ba ）gydF4y2Ba egydF4y2Ba −gydF4y2Ba jgydF4y2Ba 2gydF4y2Ba πgydF4y2Ba fgydF4y2Ba ngydF4y2Ba ，gydF4y2Ba

在哪里gydF4y2Ba

  • ggydF4y2Ba （gydF4y2Ba ngydF4y2Ba ）gydF4y2Ba -窗函数的长度gydF4y2Ba 米gydF4y2Ba ．gydF4y2Ba

  • XgydF4y2Ba 米gydF4y2Ba （gydF4y2Ba fgydF4y2Ba ）gydF4y2Ba -以时间为中心的加窗数据的DFTgydF4y2Ba 米gydF4y2Ba RgydF4y2Ba ．gydF4y2Ba

  • RgydF4y2Ba -连续dft之间的跳数大小。跃点大小是窗口长度之间的差值gydF4y2Ba 米gydF4y2Ba 重叠长度gydF4y2Ba lgydF4y2Ba ．gydF4y2Ba

短时傅里叶变换的幅值的平方得到gydF4y2Ba光谱图gydF4y2Ba表示函数的功率谱密度。gydF4y2Ba

规范化的不一致gydF4y2Ba

的gydF4y2Ba规范化的不一致gydF4y2Ba度量在连续的优化迭代中对重构过程收敛性的改进。gydF4y2Ba

规范化的不一致性定义为gydF4y2Ba

不一致gydF4y2Ba ＝gydF4y2Ba 为gydF4y2Ba STFTgydF4y2Ba （gydF4y2Ba ISTFTgydF4y2Ba （gydF4y2Ba 年代gydF4y2Ba 美国东部时间gydF4y2Ba ）gydF4y2Ba ）gydF4y2Ba −gydF4y2Ba 年代gydF4y2Ba 美国东部时间gydF4y2Ba 为gydF4y2Ba 为gydF4y2Ba 年代gydF4y2Ba 美国东部时间gydF4y2Ba 为gydF4y2Ba ，gydF4y2Ba

在哪里gydF4y2Ba年代gydF4y2Ba美国东部时间gydF4y2Ba为每次迭代估计的复短时间傅里叶变换，括号中表示矩阵范数，gydF4y2BaSTFTgydF4y2Ba为短时傅里叶变换，gydF4y2BaISTFTgydF4y2Ba表示它的逆矩阵。gydF4y2Bastftmag2siggydF4y2Ba使用MATLABgydF4y2Ba®gydF4y2Ba函数gydF4y2Ba规范gydF4y2Ba来计算矩阵规范。有关STFT及其逆函数的更多信息，请参见gydF4y2Ba短时傅里叶变换gydF4y2Ba和gydF4y2Ba短时傅里叶反变换gydF4y2Ba．gydF4y2Ba

工具书类gydF4y2Ba

[1] Griffin, Daniel W.和Jae S. Lim。“修正短时傅里叶变换的信号估计”。gydF4y2BaIEEE声学、语音和信号处理汇刊gydF4y2Ba．第32卷第2期，1984年4月，236-243页。https://doi.org/10.1109/TASSP.1984.1164317。gydF4y2Ba

[2] 佩罗丁、纳塔娜·埃尔、彼得·巴拉兹和彼得·桑德加德：《一种快速格里芬-林算法》，载于gydF4y2Ba2013 IEEE信号处理在音频和声学中的应用研讨会gydF4y2Ba，纽约州新帕尔茨，2013年10月20日至23日。https://doi.org/10.1109/WASPAA.2013.6701851.gydF4y2Ba

[3] Le Roux, Jonathan, Kameoka Hirokazu, Nobutaka Ono, Shigeki Sagayama。基于谱图一致性的幅度STFT谱图快速信号重建在gydF4y2Ba第十三届数字音频效果国际会议论文集(DAFx-10)gydF4y2Ba，奥地利格拉茨，2010年9月6日至10日。gydF4y2Ba

扩展功能gydF4y2Ba

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

另请参阅gydF4y2Ba

功能gydF4y2Ba

介绍了R2020bgydF4y2Ba

信号处理工具箱文档gydF4y2Ba

金宝app

基于MATLAB的信号处理深度学习gydF4y2Ba

基于MATLAB的信号处理深度学习gydF4y2Ba

下载白皮书gydF4y2Ba