主要内容gydF4y2Ba

stftgydF4y2Ba

短时傅里叶变换gydF4y2Ba

所有的页面崩溃gydF4y2Ba

语法gydF4y2Ba

s = stft (x)gydF4y2Ba
s = stft (x, fs)gydF4y2Ba
s = stft (x, ts)gydF4y2Ba
s = stft (gydF4y2Ba___gydF4y2Ba、名称、值)gydF4y2Ba
(s、f) = stft (gydF4y2Ba___gydF4y2Ba)gydF4y2Ba
f (s, t) = stft (gydF4y2Ba___gydF4y2Ba)gydF4y2Ba
stft (gydF4y2Ba___gydF4y2Ba)gydF4y2Ba

描述gydF4y2Ba

年代gydF4y2Ba= stft (gydF4y2BaxgydF4y2Ba)gydF4y2Ba返回gydF4y2Ba短时傅里叶变换gydF4y2Ba(STFT)gydF4y2BaxgydF4y2Ba。gydF4y2Ba

例子gydF4y2Ba

年代gydF4y2Ba= stft (gydF4y2BaxgydF4y2Ba,gydF4y2BafsgydF4y2Ba)gydF4y2Ba返回的STFTgydF4y2BaxgydF4y2Ba使用采样率gydF4y2BafsgydF4y2Ba。gydF4y2Ba

年代gydF4y2Ba= stft (gydF4y2BaxgydF4y2Ba,gydF4y2BatsgydF4y2Ba)gydF4y2Ba返回的STFTgydF4y2BaxgydF4y2Ba使用样品的时间gydF4y2BatsgydF4y2Ba。gydF4y2Ba

例子gydF4y2Ba

年代gydF4y2Ba= stft (gydF4y2Ba___gydF4y2Ba,gydF4y2Ba名称,值gydF4y2Ba)gydF4y2Ba使用名称-值对参数指定附加选项。选项包括FFT窗口和长度。这些参数可以被添加到任何以前的输入语法。gydF4y2Ba

例子gydF4y2Ba

(gydF4y2Ba年代gydF4y2Ba,gydF4y2BafgydF4y2Ba)= stft (gydF4y2Ba___gydF4y2Ba)gydF4y2Ba返回频率gydF4y2BafgydF4y2BaSTFT的评估。gydF4y2Ba

例子gydF4y2Ba

(gydF4y2Ba年代gydF4y2Ba,gydF4y2BafgydF4y2Ba,gydF4y2BatgydF4y2Ba)= stft (gydF4y2Ba___gydF4y2Ba)gydF4y2Ba返回的时候STFT评估。gydF4y2Ba

stft (gydF4y2Ba___gydF4y2Ba)gydF4y2Ba没有输出参数的块的大小STFT在当前图窗口。gydF4y2Ba

例子gydF4y2Ba

全部折叠gydF4y2Ba

打开生活的脚本gydF4y2Ba

生成两秒的电压控制振荡器输出,由正弦信号采样控制在10 kHz。gydF4y2Ba

fs = 10 e3;t = 0:1 / fs: 2;x = vco(罪(2 *π* t) (0.1 - 0.4) * fs, fs);gydF4y2Ba

计算和情节的STFT的信号。使用256 Kaiser窗长度和形状参数gydF4y2Ba βgydF4y2Ba =gydF4y2Ba 5gydF4y2Ba 。指定的长度重叠220样品和DFT长度为512点。情节的STFT默认colormap和视图。gydF4y2Ba

stft (x, fs,gydF4y2Ba“窗口”gydF4y2Ba凯瑟(256 5),gydF4y2Ba“OverlapLength”gydF4y2Ba,220,gydF4y2Ba“FFTLength”gydF4y2Ba,512)gydF4y2Ba

图包含一个坐标轴对象。坐标轴对象与标题短时傅里叶变换包含一个类型的对象的形象。gydF4y2Ba

改变视图显示STFT瀑布阴谋。设置colormapgydF4y2Ba飞机gydF4y2Ba。gydF4y2Ba

colormap视图(-45、65)gydF4y2Ba飞机gydF4y2Ba

图包含一个坐标轴对象。坐标轴对象与标题短时傅里叶变换包含一个类型的对象的表面。gydF4y2Ba

打开生活的脚本gydF4y2Ba

生成一个二次啁啾取样1 kHz 2秒。瞬时频率为100赫兹gydF4y2Ba tgydF4y2Ba =gydF4y2Ba 0gydF4y2Ba 和十字架200 HzgydF4y2Ba tgydF4y2Ba =gydF4y2Ba 1gydF4y2Ba 第二。gydF4y2Ba

ts = 0:1/1e3:2;f0 = 100;f1 = 200;x =唧唧声(ts, f0 1 f1,gydF4y2Ba“二次”gydF4y2Ba[],gydF4y2Ba“凹”gydF4y2Ba);gydF4y2Ba

计算和显示二次唧唧喳喳的STFT 1毫秒的时间。gydF4y2Ba

d =秒(1 e - 3);赢得=汉明(100gydF4y2Ba“周期”gydF4y2Ba);stft (x, d,gydF4y2Ba“窗口”gydF4y2Ba,赢了,gydF4y2Ba“OverlapLength”gydF4y2Ba,98,gydF4y2Ba“FFTLength”gydF4y2Ba,128);gydF4y2Ba

图包含一个坐标轴对象。坐标轴对象与标题短时傅里叶变换包含一个类型的对象的形象。gydF4y2Ba

打开生活的脚本gydF4y2Ba

生成一个信号采样5 kHz 4秒。信号由一组脉冲减少持续时间由区域的振动振幅和频率波动增加的趋势。画出信号。gydF4y2Ba

fs = 5000;t = 0:1 / fs: 4 - 1 / f;x = besselj(0600 *(罪(2 *π* (t + 1)。^ 3/30)。^ 5));情节(t, x)gydF4y2Ba

图包含一个坐标轴对象。坐标轴对象包含一个类型的对象。gydF4y2Ba

计算单边、双边和集中短时傅里叶变换的信号。在所有情况下,使用202 -样本Kaiser窗形状系数gydF4y2Ba βgydF4y2Ba =gydF4y2Ba 10gydF4y2Ba 窗口信号段。显示用于计算每个变换的频率范围。gydF4y2Ba

随机数生成器= [gydF4y2Ba“单向的”gydF4y2Ba“双侧”gydF4y2Ba“中心”gydF4y2Ba];gydF4y2Ba为gydF4y2Ba随机数生成器kj = 1:长度()选择= {gydF4y2Ba“窗口”gydF4y2Ba凯瑟(202、10)gydF4y2Ba“FrequencyRange”gydF4y2Ba随机数生成器(kj)};[~ f] = stft (x, fs,选择{:});次要情节(长度(随机数生成器),1,kj) stft (x, fs,选择{}):标题(sprintf (gydF4y2Ba“% s”: [% 5.3 f, % 5.3 f] kHz的gydF4y2Ba随机数生成器(kj), (f (1) f(结束)]/ 1000)gydF4y2Ba结束gydF4y2Ba

图包含3轴对象。轴与标题“单向的”对象1:[0.000,2.500]kHz包含一个类型的对象的形象。轴与标题“双侧”对象2:[0.000,4.975]kHz包含一个类型的对象的形象。坐标轴对象3标题“中心”:[-2.475,2.500]kHz包含一个类型的对象的形象。gydF4y2Ba

重复计算,但现在Kaiser窗的长度更改为203,奇数。的gydF4y2Ba双侧的gydF4y2Ba不改变频率间隔。另外两个在较高频率间隔变得开放。gydF4y2Ba

为gydF4y2Ba随机数生成器kj = 1:长度()选择= {gydF4y2Ba“窗口”gydF4y2Ba凯瑟(203、10)gydF4y2Ba“FrequencyRange”gydF4y2Ba随机数生成器(kj)};[~ f] = stft (x, fs,选择{:});次要情节(长度(随机数生成器),1,kj) stft (x, fs,选择{}):标题(sprintf (gydF4y2Ba“% s”: [% 5.3 f, % 5.3 f] kHz的gydF4y2Ba随机数生成器(kj), (f (1) f(结束)]/ 1000)gydF4y2Ba结束gydF4y2Ba

图包含3轴对象。轴与标题“单向的”对象1:[0.000,2.488]kHz包含一个类型的对象的形象。轴与标题“双侧”对象2:[0.000,4.975]kHz包含一个类型的对象的形象。坐标轴对象3标题“中心”:[-2.488,2.488]kHz包含一个类型的对象的形象。gydF4y2Ba

打开生活的脚本gydF4y2Ba

生成一个三通道信号组成的三个不同啁啾取样1千赫一秒钟。gydF4y2Ba

  1. 第一个通道由一个凹二次线性调频与瞬时频率100赫兹gydF4y2Bat =gydF4y2Ba0和跨越300 HzgydF4y2Bat =gydF4y2Ba1秒。它有一个初始阶段等于45度。gydF4y2Ba

  2. 第二通道由一个凸二次啁啾与瞬时频率100赫兹gydF4y2Bat =gydF4y2Ba0和跨越500 HzgydF4y2Bat =gydF4y2Ba1秒。gydF4y2Ba

  3. 第三通道由一个对数线性调频与瞬时频率300赫兹gydF4y2Bat =gydF4y2Ba0和跨越500 HzgydF4y2Bat =gydF4y2Ba1秒。gydF4y2Ba

计算多道信号使用一个周期的STFT汉明窗的长度128和50个样本的重叠长度。gydF4y2Ba

fs = 1 e3;t = 0:1 / fs: 1 - 1 / f;x =[唧唧声(t, 100, 1300,gydF4y2Ba“二次”gydF4y2Ba45岁的gydF4y2Ba“凹”gydF4y2Ba);唧唧声(t, 100, 1500,gydF4y2Ba“二次”gydF4y2Ba[],gydF4y2Ba“凸”gydF4y2Ba);唧唧声(t, 300, 1500,gydF4y2Ba“对数”gydF4y2Ba)];F (S, T) = stft (x, fs,gydF4y2Ba“窗口”gydF4y2Ba128年,汉明(gydF4y2Ba“周期”gydF4y2Ba),gydF4y2Ba“OverlapLength”gydF4y2Ba,50);gydF4y2Ba

想象每个通道的STFT瀑布阴谋。控制轴的行为使用helper函数gydF4y2BahelperGraphicsOptgydF4y2Ba。gydF4y2Ba

瀑布(F T abs (S (:,: 1))) helperGraphicsOpt (1)gydF4y2Ba

图包含一个坐标轴对象。坐标轴对象标题输入通道:1包含一个补丁类型的对象。gydF4y2Ba

瀑布(F T abs (S (:,: 2))) helperGraphicsOpt (2)gydF4y2Ba

图包含一个坐标轴对象。坐标轴对象标题输入通道:2包含一个对象类型的补丁。gydF4y2Ba

瀑布(F T abs (S (:,:, 3))) helperGraphicsOpt (3)gydF4y2Ba

图包含一个坐标轴对象。坐标轴对象标题输入通道:3包含一个对象类型的补丁。gydF4y2Ba

这个helper函数设置当前轴的外观和行为。gydF4y2Ba

函数gydF4y2Ba甘氨胆酸helperGraphicsOpt (ChannelId) ax =;斧子。XDir =“反向”gydF4y2Ba;斧子。ZLim = 30 [0];ax.Title。字符串=(gydF4y2Ba的输入通道:gydF4y2Banum2str (ChannelId)];ax.XLabel。字符串=gydF4y2Ba的频率(赫兹)gydF4y2Ba;ax.YLabel。字符串=gydF4y2Ba的时间(秒)gydF4y2Ba;斧子。视图= 45 [30];gydF4y2Ba结束gydF4y2Ba

输入参数gydF4y2Ba

全部折叠gydF4y2Ba

输入信号,指定为一个向量,矩阵或MATLABgydF4y2Ba®gydF4y2Ba时间表gydF4y2Ba。gydF4y2Ba

请注意gydF4y2Ba

如果你反gydF4y2Ba年代gydF4y2Ba使用gydF4y2BaistftgydF4y2Ba和想要的结果是相同的长度gydF4y2BaxgydF4y2Ba的价值,gydF4y2Ba长度(长度(x) -noverlap) /((窗口)-noverlap)gydF4y2Ba必须是一个整数。gydF4y2Ba

  • 如果输入有多个渠道,指定gydF4y2BaxgydF4y2Ba一个矩阵,每一列对应一个频道。gydF4y2Ba

  • 对时间表的输入,gydF4y2BaxgydF4y2Ba必须包含均匀增加有限行乘以。如果时间表有缺失或重复的时间点,你可以修理它使用的技巧gydF4y2Ba清洁时间表与失踪、重复或不均匀gydF4y2Ba。gydF4y2Ba

  • 对于多通道输入时间表,指定gydF4y2BaxgydF4y2Ba作为一个时间表与单个变量包含一个矩阵或时间表与多个变量每个都包含一个列向量。所有的变量必须具有相同的精度。gydF4y2Ba

每个通道的gydF4y2BaxgydF4y2Ba必须有一个长度大于或等于窗口长度。gydF4y2Ba

例子:gydF4y2Ba唧唧声(0:1/4e3:2,250, 1500,“二次”)gydF4y2Ba指定一个单通道的唧唧声。gydF4y2Ba

例子:gydF4y2Ba时间表(兰德(5,2)' SampleRate ', 1)gydF4y2Ba随机变量指定了一个双通道采样1赫兹为4秒。gydF4y2Ba

例子:gydF4y2Ba时间表(兰德(5,1),兰德(5、1),“SampleRate”, 1)gydF4y2Ba随机变量指定了一个双通道采样1赫兹为4秒。gydF4y2Ba

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

采样率,指定为一个积极的标量。这种说法只适用于当gydF4y2BaxgydF4y2Ba是一个向量或矩阵。gydF4y2Ba

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

样品时间,指定为一个gydF4y2Ba持续时间gydF4y2Ba标量。这种说法只适用于当gydF4y2BaxgydF4y2Ba是一个向量或矩阵gydF4y2Ba

例子:gydF4y2Ba秒(1)gydF4y2Ba是一个gydF4y2Ba持续时间gydF4y2Ba标量代表1秒时间连续信号样本的区别。gydF4y2Ba

数据类型:gydF4y2Ba持续时间gydF4y2Ba

名称-值参数gydF4y2Ba

指定可选的双参数作为gydF4y2BaName1 = Value1,…,以=家gydF4y2Ba,在那里gydF4y2Ba的名字gydF4y2Ba参数名称和吗gydF4y2Ba价值gydF4y2Ba相应的价值。名称-值参数必须出现在其他参数,但对的顺序无关紧要。gydF4y2Ba

R2021a之前,用逗号来分隔每一个名称和值,并附上gydF4y2Ba的名字gydF4y2Ba在报价。gydF4y2Ba

例子:gydF4y2Ba“窗口”,汉明(100),“OverlapLength”, 50岁' FFTLength ', 128gydF4y2Bawindows 100 -样本数据使用汉明窗,和50个样本之间的重叠相邻段和128点FFT。gydF4y2Ba

光谱窗口中,指定为一个向量。如果你不指定窗口或指定它为空,长度为128的函数使用一个损害窗口。的长度gydF4y2Ba“窗口”gydF4y2Ba必须大于或等于2。gydF4y2Ba

可用窗口的列表,请参阅gydF4y2Ba窗户gydF4y2Ba。gydF4y2Ba

例子:gydF4y2Ba损害(N + 1)gydF4y2Ba和gydF4y2Ba(1-cos(2 *π* (0:N) / N)) / 2gydF4y2Ba都指定一个损害窗口长度gydF4y2BaNgydF4y2Ba+ 1。gydF4y2Ba

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

重叠的样本数量,指定为一个正整数的长度小于gydF4y2Ba“窗口”gydF4y2Ba。如果您省略gydF4y2Ba“OverlapLength”gydF4y2Ba或指定它是空的,它被设置为最大的整数不到75%的窗口长度,这是默认的96个样本损害窗口。gydF4y2Ba

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

数量的DFT点,指定为一个正整数。的值必须大于或等于窗口长度。如果输入信号的长度小于DFT长度,用0填充数据。gydF4y2Ba

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

STFT频率范围,指定为gydF4y2Ba“中心”gydF4y2Ba,gydF4y2Ba双侧的gydF4y2Ba,或gydF4y2Ba“单向的”gydF4y2Ba。gydF4y2Ba

  • “中心”gydF4y2Ba——集中STFT计算一个双边。如果gydF4y2Ba“FFTLength”gydF4y2Ba是偶数,那么gydF4y2Ba年代gydF4y2Ba计算的时间间隔gydF4y2Ba(-gydF4y2BaπgydF4y2Ba,gydF4y2BaπgydF4y2Ba]gydF4y2Barad /样品。如果gydF4y2Ba“FFTLength”gydF4y2Ba是奇数,那么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——计算一个双边STFT间隔gydF4y2Ba[0,2gydF4y2BaπgydF4y2Ba)gydF4y2Barad /样品。如果你指定时间信息,然后间隔gydF4y2Ba[0,gydF4y2BafgydF4y2Ba年代gydF4y2Ba)gydF4y2Ba周期/单位时间。gydF4y2Ba

  • “单向的”gydF4y2Ba——计算一个片面的STFT。如果gydF4y2Ba“FFTLength”gydF4y2Ba是偶数,那么gydF4y2Ba年代gydF4y2Ba计算的时间间隔gydF4y2Ba[0,gydF4y2BaπgydF4y2Ba]gydF4y2Barad /样品。如果gydF4y2Ba“FFTLength”gydF4y2Ba是奇数,那么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,gydF4y2BastftgydF4y2Ba输出值的积极尼奎斯特范围和不节约总功率。gydF4y2Ba

例如,看到的gydF4y2BaSTFT频率范围gydF4y2Ba。gydF4y2Ba

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

输出时间维度,指定为gydF4y2Ba“acrosscolumns”gydF4y2Ba或gydF4y2Ba“downrows”gydF4y2Ba。将这个值设置为gydF4y2Ba“downrows”gydF4y2Ba如果你想要的时间维度gydF4y2Ba年代gydF4y2Ba的行和列频率维度。将这个值设置为gydF4y2Ba“acrosscolumns”gydF4y2Ba如果你想要的时间维度gydF4y2Ba年代gydF4y2Ba在列和频率维度的行。这个输入是忽略了如果没有输出参数的函数被调用。gydF4y2Ba

输出参数gydF4y2Ba

全部折叠gydF4y2Ba

短时傅里叶变换,作为一个矩阵或三维数组返回。时间增加的列gydF4y2Ba年代gydF4y2Ba和频率增加的行。第三个维度,如果存在,对应于输入通道。gydF4y2Ba

  • 如果信号gydF4y2BaxgydF4y2Ba有gydF4y2BaNgydF4y2BaxgydF4y2Ba样品,然后gydF4y2Ba年代gydF4y2Ba有gydF4y2BakgydF4y2Ba列,gydF4y2BakgydF4y2Ba=⌊(gydF4y2BaNgydF4y2BaxgydF4y2Ba- - - - - -gydF4y2BalgydF4y2Ba)/ (gydF4y2Ba米gydF4y2Ba- - - - - -gydF4y2BalgydF4y2Ba)⌋gydF4y2Ba,gydF4y2Ba米gydF4y2Ba的长度是gydF4y2Ba“窗口”gydF4y2Ba,gydF4y2BalgydF4y2Ba是gydF4y2Ba“OverlapLength”gydF4y2Ba地板,⌊⌋符号表示函数。gydF4y2Ba

  • 的行数gydF4y2Ba年代gydF4y2Ba等于指定的值gydF4y2Ba“FFTLength”gydF4y2Ba。gydF4y2Ba

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

频率STFT的评估,作为一个向量返回。gydF4y2Ba

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

瞬间的时间,作为一个向量返回。gydF4y2BatgydF4y2Ba包含对应的时间值数据的中心部分用于计算短时功率谱估计。gydF4y2Ba

  • 如果一个采样率gydF4y2BafsgydF4y2Ba提供,那么向量包含时间值在几秒钟内。gydF4y2Ba

  • 如果样品的时间gydF4y2BatsgydF4y2Ba提供,那么向量是一系列时间用相同的时间格式作为输入。gydF4y2Ba

  • 如果没有提供时间信息,那么向量包含示例数据。gydF4y2Ba

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

更多关于gydF4y2Ba

全部折叠gydF4y2Ba

短时傅里叶变换gydF4y2Ba

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

STFT的信号是通过滑动一个计算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符号表示函数的地板上。矩阵的行数等于gydF4y2BaNgydF4y2BaDFTgydF4y2BaDFT点的数量,集中和双边变换gydF4y2Ba⌊gydF4y2BaNgydF4y2BaDFTgydF4y2Ba/ 2⌋+ 1gydF4y2Ba片面的转换。gydF4y2Ba

STFT矩阵是由gydF4y2Ba XgydF4y2Ba (gydF4y2Ba fgydF4y2Ba )gydF4y2Ba =gydF4y2Ba (gydF4y2Ba XgydF4y2Ba 1gydF4y2Ba (gydF4y2Ba fgydF4y2Ba )gydF4y2Ba XgydF4y2Ba 2gydF4y2Ba (gydF4y2Ba fgydF4y2Ba )gydF4y2Ba XgydF4y2Ba 3gydF4y2Ba (gydF4y2Ba fgydF4y2Ba )gydF4y2Ba ⋯gydF4y2Ba XgydF4y2Ba kgydF4y2Ba (gydF4y2Ba fgydF4y2Ba )gydF4y2Ba ]gydF4y2Ba 这样gydF4y2Ba 米gydF4y2Ba 这个矩阵的元素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 - DFT的窗口的数据集中的时候了gydF4y2Ba 米gydF4y2Ba RgydF4y2Ba 。gydF4y2Ba

  • RgydF4y2Ba ——跳连续阶之间的大小。跳大小窗口长度的区别gydF4y2Ba 米gydF4y2Ba 和重叠长度gydF4y2Ba lgydF4y2Ba 。gydF4y2Ba

大小的平方STFT的收益率gydF4y2Ba光谱图gydF4y2Ba的功率谱密度函数的表示形式。gydF4y2Ba

完美的重建gydF4y2Ba

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

  • 输入大小——如果你反转的输出gydF4y2BastftgydF4y2Ba使用gydF4y2BaistftgydF4y2Ba和想要的结果是相同的长度作为输入信号gydF4y2BaxgydF4y2Ba的价值,gydF4y2Ba kgydF4y2Ba =gydF4y2Ba (gydF4y2Ba lgydF4y2Ba egydF4y2Ba ngydF4y2Ba ggydF4y2Ba tgydF4y2Ba hgydF4y2Ba (gydF4y2Ba xgydF4y2Ba )gydF4y2Ba - - - - - -gydF4y2Ba ngydF4y2Ba ogydF4y2Ba vgydF4y2Ba egydF4y2Ba rgydF4y2Ba lgydF4y2Ba 一个gydF4y2Ba pgydF4y2Ba )gydF4y2Ba (gydF4y2Ba lgydF4y2Ba egydF4y2Ba ngydF4y2Ba ggydF4y2Ba tgydF4y2Ba hgydF4y2Ba (gydF4y2Ba wgydF4y2Ba 我gydF4y2Ba ngydF4y2Ba dgydF4y2Ba ogydF4y2Ba wgydF4y2Ba )gydF4y2Ba - - - - - -gydF4y2Ba ngydF4y2Ba ogydF4y2Ba vgydF4y2Ba egydF4y2Ba rgydF4y2Ba lgydF4y2Ba 一个gydF4y2Ba pgydF4y2Ba )gydF4y2Ba 必须是一个整数。gydF4y2Ba

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

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

引用gydF4y2Ba

[1]Mitra Sanjit K。gydF4y2Ba数字信号处理:一个基于计算机的方法gydF4y2Ba。第二版。纽约:麦格劳-希尔,2001年。gydF4y2Ba

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

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

扩展功能gydF4y2Ba

版本历史gydF4y2Ba

介绍了R2019agydF4y2Ba

另请参阅gydF4y2Ba

功能gydF4y2Ba

对象gydF4y2Ba

主题gydF4y2Ba