主要内容gydF4y2Ba

光谱图gydF4y2Ba

使用短时傅里叶变换光谱图gydF4y2Ba

所有的页面崩溃gydF4y2Ba

语法gydF4y2Ba

s =光谱图(x)gydF4y2Ba
s =光谱图(x,窗口)gydF4y2Ba
s =光谱图(x,窗口,noverlap)gydF4y2Ba
s =光谱图(x,窗口、noverlap nfft)gydF4y2Ba
[s w t] =光谱图(gydF4y2Ba___gydF4y2Ba)gydF4y2Ba
f (s, t) =光谱图(gydF4y2Ba___gydF4y2Bafs)gydF4y2Ba
[s w t] =光谱图(x,窗口、noverlap w)gydF4y2Ba
f (s, t) =光谱图(x,窗口、noverlap f, f)gydF4y2Ba
(gydF4y2Ba___gydF4y2Ba、ps) =光谱图(gydF4y2Ba___gydF4y2Ba)gydF4y2Ba
(gydF4y2Ba___gydF4y2Ba)=光谱图(gydF4y2Ba___gydF4y2Ba“重新分配”)gydF4y2Ba
(gydF4y2Ba___gydF4y2Ba、ps、fc tc) =光谱图(gydF4y2Ba___gydF4y2Ba)gydF4y2Ba
(gydF4y2Ba___gydF4y2Ba)=光谱图(gydF4y2Ba___gydF4y2Bafreqrange)gydF4y2Ba
(gydF4y2Ba___gydF4y2Ba)=光谱图(gydF4y2Ba___gydF4y2Ba、名称、值)gydF4y2Ba
(gydF4y2Ba___gydF4y2Ba)=光谱图(gydF4y2Ba___gydF4y2Baspectrumtype)gydF4y2Ba
光谱图(gydF4y2Ba___gydF4y2Ba)gydF4y2Ba
光谱图(gydF4y2Ba___gydF4y2Bafreqloc)gydF4y2Ba

描述gydF4y2Ba

例子gydF4y2Ba

年代gydF4y2Ba=光谱图(gydF4y2BaxgydF4y2Ba)gydF4y2Ba返回输入信号的短时傅里叶变换,gydF4y2BaxgydF4y2Ba。每一列的gydF4y2Ba年代gydF4y2Ba包含一个短期的估计,time-localized频率的内容gydF4y2BaxgydF4y2Ba。gydF4y2Ba

年代gydF4y2Ba=光谱图(gydF4y2BaxgydF4y2Ba,gydF4y2Ba窗口gydF4y2Ba)gydF4y2Ba使用gydF4y2Ba窗口gydF4y2Ba把信号分成段和执行窗口。gydF4y2Ba

例子gydF4y2Ba

年代gydF4y2Ba=光谱图(gydF4y2BaxgydF4y2Ba,gydF4y2Ba窗口gydF4y2Ba,gydF4y2BanoverlapgydF4y2Ba)gydF4y2Ba使用gydF4y2BanoverlapgydF4y2Ba毗邻的重叠部分的样本。gydF4y2Ba

例子gydF4y2Ba

年代gydF4y2Ba=光谱图(gydF4y2BaxgydF4y2Ba,gydF4y2Ba窗口gydF4y2Ba,gydF4y2BanoverlapgydF4y2Ba,gydF4y2BanfftgydF4y2Ba)gydF4y2Ba使用gydF4y2BanfftgydF4y2Ba采样点计算离散傅里叶变换。gydF4y2Ba

(gydF4y2Ba年代gydF4y2Ba,gydF4y2BawgydF4y2Ba,gydF4y2BatgydF4y2Ba)=光谱图(gydF4y2Ba___gydF4y2Ba)gydF4y2Ba返回一个向量的归一化频率,gydF4y2BawgydF4y2Ba一个向量的时间瞬间,gydF4y2BatgydF4y2Ba计算的谱图。这个语法可以包括从以前的语法输入参数的任意组合。gydF4y2Ba

例子gydF4y2Ba

(gydF4y2Ba年代gydF4y2Ba,gydF4y2BafgydF4y2Ba,gydF4y2BatgydF4y2Ba)=光谱图(gydF4y2Ba___gydF4y2Ba,gydF4y2BafsgydF4y2Ba)gydF4y2Ba返回一个向量的周期性的频率,gydF4y2BafgydF4y2Ba表示的采样率,gydF4y2BafsgydF4y2Ba。gydF4y2BafsgydF4y2Ba必须输入第五gydF4y2Ba光谱图gydF4y2Ba。输入采样率和仍在使用前可选参数的默认值,指定这些参数为空,gydF4y2Ba[]gydF4y2Ba。gydF4y2Ba

例子gydF4y2Ba

(gydF4y2Ba年代gydF4y2Ba,gydF4y2BawgydF4y2Ba,gydF4y2BatgydF4y2Ba)=光谱图(gydF4y2BaxgydF4y2Ba,gydF4y2Ba窗口gydF4y2Ba,gydF4y2BanoverlapgydF4y2Ba,gydF4y2BawgydF4y2Ba)gydF4y2Ba返回指定的归一化频率的谱图gydF4y2BawgydF4y2Ba。gydF4y2Ba

(gydF4y2Ba年代gydF4y2Ba,gydF4y2BafgydF4y2Ba,gydF4y2BatgydF4y2Ba)=光谱图(gydF4y2BaxgydF4y2Ba,gydF4y2Ba窗口gydF4y2Ba,gydF4y2BanoverlapgydF4y2Ba,gydF4y2BafgydF4y2Ba,gydF4y2BafsgydF4y2Ba)gydF4y2Ba返回在指定的周期频率谱图gydF4y2BafgydF4y2Ba。gydF4y2Ba

(gydF4y2Ba___gydF4y2Ba,gydF4y2BapsgydF4y2Ba)=光谱图(gydF4y2Ba___gydF4y2Ba)gydF4y2Ba还返回一个矩阵,gydF4y2BapsgydF4y2Ba,包含估计功率谱密度(PSD)或每段的功率谱。gydF4y2Ba

例子gydF4y2Ba

(gydF4y2Ba___gydF4y2Ba)=光谱图(gydF4y2Ba___gydF4y2Ba“重新分配”)gydF4y2Ba抽调每个PSD或功率谱估计的能量中心的位置。如果你的信号包含距离时间或光谱组件,那么该选项生成一个尖锐的谱图。gydF4y2Ba

例子gydF4y2Ba

(gydF4y2Ba___gydF4y2Ba,gydF4y2BapsgydF4y2Ba,gydF4y2Ba足球俱乐部gydF4y2Ba,gydF4y2BatcgydF4y2Ba)=光谱图(gydF4y2Ba___gydF4y2Ba)gydF4y2Ba返回两个矩阵,gydF4y2Ba足球俱乐部gydF4y2Ba和gydF4y2BatcgydF4y2Ba包含中心的频率和时间,每个PSD的能量或功率谱估计。gydF4y2Ba

(gydF4y2Ba___gydF4y2Ba)=光谱图(gydF4y2Ba___gydF4y2Ba,gydF4y2BafreqrangegydF4y2Ba)gydF4y2Ba返回PSD或功率谱估计在指定的频率范围gydF4y2BafreqrangegydF4y2Ba。有效的选择gydF4y2BafreqrangegydF4y2Ba是gydF4y2Ba“单向的”gydF4y2Ba,gydF4y2Ba双侧的gydF4y2Ba,gydF4y2Ba“中心”gydF4y2Ba。gydF4y2Ba

例子gydF4y2Ba

(gydF4y2Ba___gydF4y2Ba)=光谱图(gydF4y2Ba___gydF4y2Ba,gydF4y2Ba名称,值gydF4y2Ba)gydF4y2Ba使用名称-值对参数指定附加选项。选项包括最低阈值和输出时间维度。gydF4y2Ba

(gydF4y2Ba___gydF4y2Ba)=光谱图(gydF4y2Ba___gydF4y2Ba,gydF4y2BaspectrumtypegydF4y2Ba)gydF4y2Ba返回PSD估计如果gydF4y2BaspectrumtypegydF4y2Ba被指定为gydF4y2Bapsd的gydF4y2Ba并返回功率谱估计gydF4y2BaspectrumtypegydF4y2Ba被指定为gydF4y2Ba“权力”gydF4y2Ba。gydF4y2Ba

例子gydF4y2Ba

光谱图(gydF4y2Ba___gydF4y2Ba)gydF4y2Ba没有输出参数的光谱图块在当前图窗口。gydF4y2Ba

例子gydF4y2Ba

光谱图(gydF4y2Ba___gydF4y2Ba,gydF4y2BafreqlocgydF4y2Ba)gydF4y2Ba指定的轴绘制频率。gydF4y2Ba

例子gydF4y2Ba

全部折叠gydF4y2Ba

打开生活的脚本gydF4y2Ba

生成gydF4y2Ba NgydF4y2Ba xgydF4y2Ba =gydF4y2Ba 1gydF4y2Ba 0gydF4y2Ba 2gydF4y2Ba 4gydF4y2Ba 样品的一个信号,由一个正弦曲线。归一化频率的正弦曲线gydF4y2Ba 2gydF4y2Ba πgydF4y2Ba /gydF4y2Ba 5gydF4y2Ba rad /样本和gydF4y2Ba 4gydF4y2Ba πgydF4y2Ba /gydF4y2Ba 5gydF4y2Ba rad /样品。更高的频率正弦信号的另一个正弦信号的振幅的10倍。gydF4y2Ba

N = 1024;n = 0: n - 1;w0 =π/ 2 * 5;x =罪(w0 * n) + 10 *罪(2 * w0 * n);gydF4y2Ba

使用函数计算短时傅里叶变换的违约。绘制光谱图。gydF4y2Ba

s =光谱图(x);光谱图(x,gydF4y2Ba“桠溪”gydF4y2Ba)gydF4y2Ba

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

重复计算。gydF4y2Ba

  • 将信号分成几个部分的长度gydF4y2Ba ngydF4y2Ba 年代gydF4y2Ba cgydF4y2Ba =gydF4y2Ba ⌊gydF4y2Ba NgydF4y2Ba xgydF4y2Ba /gydF4y2Ba 4gydF4y2Ba 。gydF4y2Ba 5gydF4y2Ba ⌋gydF4y2Ba 。gydF4y2Ba

  • 窗口的部分使用汉明窗。gydF4y2Ba

  • 指定50%重叠的部分。gydF4y2Ba

  • 计算FFT,使用gydF4y2Ba 马克斯gydF4y2Ba (gydF4y2Ba 2gydF4y2Ba 5gydF4y2Ba 6gydF4y2Ba ,gydF4y2Ba 2gydF4y2Ba pgydF4y2Ba )gydF4y2Ba 点,gydF4y2Ba pgydF4y2Ba =gydF4y2Ba ⌈gydF4y2Ba 日志gydF4y2Ba 2gydF4y2Ba ngydF4y2Ba 年代gydF4y2Ba cgydF4y2Ba ⌉gydF4y2Ba 。gydF4y2Ba

验证这两种方法给出相同的结果。gydF4y2Ba

Nx =长度(x);nsc =地板(Nx / 4.5);11月=地板(nsc / 2);nff = max (256 2 ^ nextpow2 (nsc));t =光谱图(x,汉明(nsc), 11月nff);maxerr = max (abs (abs (t (:)) abs (s (:))))gydF4y2Ba
maxerr = 0gydF4y2Ba

把信号分成8组相等的长度,有50%的重叠部分。指定相同的FFT长度在前面的步骤。计算短时傅里叶变换,并验证它给相同的结果与前面的两个过程。gydF4y2Ba

ns = 8;ov = 0.5;lsc =地板(Nx / (ns - (ns-1) * ov));t =光谱图(x, lsc地板(ov * lsc), nff);maxerr = max (abs (abs (t (:)) abs (s (:))))gydF4y2Ba
maxerr = 0gydF4y2Ba
打开生活的脚本gydF4y2Ba

生成一个二次唧唧声,gydF4y2BaxgydF4y2Ba,1千赫采样2秒。短促声波的频率为100赫兹最初和跨越200 HzgydF4y2BatgydF4y2Ba= 1。gydF4y2Ba

t = 0:0.001:2;x =唧唧声(t, 100, 1200,gydF4y2Ba“二次”gydF4y2Ba);gydF4y2Ba

计算和显示的光谱图gydF4y2BaxgydF4y2Ba。gydF4y2Ba

  • 128年长度将信号分成几个部分,有窗的汉明窗。gydF4y2Ba

  • 指定120个样本之间的重叠的部分。gydF4y2Ba

  • 计算谱在gydF4y2Ba ⌊gydF4y2Ba 1gydF4y2Ba 2gydF4y2Ba 8gydF4y2Ba /gydF4y2Ba 2gydF4y2Ba +gydF4y2Ba 1gydF4y2Ba ⌋gydF4y2Ba =gydF4y2Ba 6gydF4y2Ba 5gydF4y2Ba 频率和gydF4y2Ba ⌊gydF4y2Ba (gydF4y2Ba lgydF4y2Ba egydF4y2Ba ngydF4y2Ba ggydF4y2Ba tgydF4y2Ba hgydF4y2Ba (gydF4y2Ba xgydF4y2Ba )gydF4y2Ba - - - - - -gydF4y2Ba 1gydF4y2Ba 2gydF4y2Ba 0gydF4y2Ba )gydF4y2Ba /gydF4y2Ba (gydF4y2Ba 1gydF4y2Ba 2gydF4y2Ba 8gydF4y2Ba - - - - - -gydF4y2Ba 1gydF4y2Ba 2gydF4y2Ba 0gydF4y2Ba )gydF4y2Ba ⌋gydF4y2Ba =gydF4y2Ba 2gydF4y2Ba 3gydF4y2Ba 5gydF4y2Ba 时间箱。gydF4y2Ba

光谱图(x, 128120128年,1 e3)gydF4y2Ba

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

布莱克曼窗取代汉明窗。减少重叠到60样本。情节的时间轴,使其值增加从上到下。gydF4y2Ba

光谱图(x,布莱克曼(128),60128年,1 e3) ax = gca;斧子。YDir =gydF4y2Ba“反向”gydF4y2Ba;gydF4y2Ba

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

打开生活的脚本gydF4y2Ba

计算和显示每个段的二次唧唧喳喳的PSD起价100 Hz,跨越200 HzgydF4y2BatgydF4y2Ba= 1秒。指定1 kHz的采样率,一段长度为128样本,和120个样本的重叠。使用128 DFT分和默认的汉明窗。gydF4y2Ba

fs = 1000;t = 0:1 / fs: 2;x =唧唧声(t, 100, 1200,gydF4y2Ba“二次”gydF4y2Ba);fs光谱图(x, 128120128年,gydF4y2Ba“桠溪”gydF4y2Ba)标题(gydF4y2Ba二次唧唧喳喳的gydF4y2Ba)gydF4y2Ba

图包含一个坐标轴对象。坐标轴对象与标题二次唧唧声包含一个类型的对象的形象。gydF4y2Ba

计算和显示每一部分的PSD线性啁啾取样1 kHz,开始在直流和十字架150 HzgydF4y2BatgydF4y2Ba= 1秒。指定一个段的长度256样品和250个样本的重叠。使用默认的汉明窗DFT和256点。gydF4y2Ba

x =唧唧声(t, 0, 1150);fs光谱图(x, 256250256年,gydF4y2Ba“桠溪”gydF4y2Ba)标题(gydF4y2Ba“线性啁啾”gydF4y2Ba)gydF4y2Ba

图包含一个坐标轴对象。坐标轴对象与标题线性啁啾包含一个类型的对象的形象。gydF4y2Ba

计算和显示每个段的PSD的对数啁啾取样1 kHz,始于20赫兹和跨越60赫兹gydF4y2BatgydF4y2Ba= 1秒。指定一个段的长度256样品和250个样本的重叠。使用默认的汉明窗DFT和256点。gydF4y2Ba

1 x =唧唧声(t, 20日,60岁,gydF4y2Ba“对数”gydF4y2Ba);光谱图(x, 256250年,[],fs,gydF4y2Ba“桠溪”gydF4y2Ba)标题(gydF4y2Ba对数唧唧喳喳的gydF4y2Ba)gydF4y2Ba

图包含一个坐标轴对象。与标题对数坐标轴对象啁啾包含一个类型的对象的形象。gydF4y2Ba

使用频率轴的对数刻度。光谱图变成一条线。gydF4y2Ba

甘氨胆酸ax =;斧子。YScale =gydF4y2Ba“日志”gydF4y2Ba;gydF4y2Ba

图包含一个坐标轴对象。坐标轴对象与标题对数唧唧声表面包含一个类型的对象。gydF4y2Ba

打开生活的脚本gydF4y2Ba

使用gydF4y2Ba光谱图gydF4y2Ba函数来衡量和跟踪信号的瞬时频率。gydF4y2Ba

生成一个二次啁啾取样1 kHz两秒钟。最初指定唧唧声以便其频率100赫兹,一秒钟后增加到200赫兹。gydF4y2Ba

fs = 1000;t = 0:1 / fs: 2 - 1 / f;y =唧唧声(t, 100, 1200,gydF4y2Ba“二次”gydF4y2Ba);gydF4y2Ba

估计的频谱使用短时傅里叶变换实现啁啾gydF4y2Ba光谱图gydF4y2Ba函数。100年长度将信号分成几个部分,有窗的汉明窗。指定80个样本之间的重叠的部分和评估频谱gydF4y2Ba ⌊gydF4y2Ba 1gydF4y2Ba 0gydF4y2Ba 0gydF4y2Ba /gydF4y2Ba 2gydF4y2Ba +gydF4y2Ba 1gydF4y2Ba ⌋gydF4y2Ba =gydF4y2Ba 5gydF4y2Ba 1gydF4y2Ba 频率。gydF4y2Ba

fs光谱图(y, 100年,80100年,gydF4y2Ba“桠溪”gydF4y2Ba)gydF4y2Ba

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

跟踪线性调频频率通过寻找与最高能量在时频脊gydF4y2Ba ⌊gydF4y2Ba (gydF4y2Ba 2gydF4y2Ba 0gydF4y2Ba 0gydF4y2Ba 0gydF4y2Ba - - - - - -gydF4y2Ba 8gydF4y2Ba 0gydF4y2Ba )gydF4y2Ba /gydF4y2Ba (gydF4y2Ba 1gydF4y2Ba 0gydF4y2Ba 0gydF4y2Ba - - - - - -gydF4y2Ba 8gydF4y2Ba 0gydF4y2Ba )gydF4y2Ba ⌋gydF4y2Ba =gydF4y2Ba 9gydF4y2Ba 6gydF4y2Ba 时间点。覆盖的瞬时频率谱图。gydF4y2Ba

(~ f t, p) =光谱图(y, 100、80100 fs);(冰箱,~,lr) = tfridge (p f);持有gydF4y2Ba在gydF4y2Baplot3 (t、冰箱、abs (p (lr)),gydF4y2Ba“线宽”gydF4y2Ba,4)gydF4y2Ba从gydF4y2Ba

图包含一个坐标轴对象。坐标轴对象包含2图像类型的对象,线。gydF4y2Ba

打开生活的脚本gydF4y2Ba

生成512个样本的唧唧声呈现正弦变化频率的内容。gydF4y2Ba

N = 512;n = 0: n - 1;x = exp (1 j *π*罪(8 * n / n) * 32);gydF4y2Ba

计算为中心的双边短时傅里叶变换的唧唧声。把信号分成32-sample段与16-sample重叠。指定64 DFT点。绘制光谱图。gydF4y2Ba

(标量、fs、ts) =光谱图(64年x, 32岁,16日,gydF4y2Ba“中心”gydF4y2Ba);64年光谱图(x, 32岁,16日,gydF4y2Ba“中心”gydF4y2Ba,gydF4y2Ba“桠溪”gydF4y2Ba)gydF4y2Ba

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

获得相同的结果通过计算64平均间隔频率间隔的光谱图gydF4y2Ba (gydF4y2Ba - - - - - -gydF4y2Ba πgydF4y2Ba ,gydF4y2Ba πgydF4y2Ba ]gydF4y2Ba 。的gydF4y2Ba“中心”gydF4y2Ba选择是没有必要。gydF4y2Ba

fintv = -π+π/ 32:π/ 32:π;(向量阵线、电视)=光谱图(fintv x, 32岁,16日);fintv光谱图(x, 32岁,16日,gydF4y2Ba“桠溪”gydF4y2Ba)gydF4y2Ba

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

打开生活的脚本gydF4y2Ba

生成一个线性调频信号在1千赫采样2秒。最初指定唧唧声以便其频率100赫兹,1秒后增加到200赫兹。gydF4y2Ba

Fs = 1000;t = 0:1 / Fs: 2;y =唧唧声(t, 100, 1200,gydF4y2Ba“二次”gydF4y2Ba);gydF4y2Ba

估计信号的重新分配谱图。gydF4y2Ba

  • 128年长度将信号分成几个部分,与Kaiser窗窗口的形状参数gydF4y2Ba βgydF4y2Ba =gydF4y2Ba 1gydF4y2Ba 8gydF4y2Ba 。gydF4y2Ba

  • 指定120个样本之间的重叠的部分。gydF4y2Ba

  • 计算谱在gydF4y2Ba ⌊gydF4y2Ba 1gydF4y2Ba 2gydF4y2Ba 8gydF4y2Ba /gydF4y2Ba 2gydF4y2Ba ⌋gydF4y2Ba =gydF4y2Ba 6gydF4y2Ba 5gydF4y2Ba 频率和gydF4y2Ba ⌊gydF4y2Ba (gydF4y2Ba lgydF4y2Ba egydF4y2Ba ngydF4y2Ba ggydF4y2Ba tgydF4y2Ba hgydF4y2Ba (gydF4y2Ba xgydF4y2Ba )gydF4y2Ba - - - - - -gydF4y2Ba 1gydF4y2Ba 2gydF4y2Ba 0gydF4y2Ba )gydF4y2Ba /gydF4y2Ba (gydF4y2Ba 1gydF4y2Ba 2gydF4y2Ba 8gydF4y2Ba - - - - - -gydF4y2Ba 1gydF4y2Ba 2gydF4y2Ba 0gydF4y2Ba )gydF4y2Ba ⌋gydF4y2Ba =gydF4y2Ba 2gydF4y2Ba 3gydF4y2Ba 5gydF4y2Ba 时间箱。gydF4y2Ba

光谱图(y,凯瑟(128年,18),120128年,Fs,gydF4y2Ba“重新分配”gydF4y2Ba,gydF4y2Ba“桠溪”gydF4y2Ba)gydF4y2Ba

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

打开生活的脚本gydF4y2Ba

生成一个线性调频信号在1千赫采样2秒。最初指定唧唧声以便其频率100赫兹,1秒后增加到200赫兹。gydF4y2Ba

Fs = 1000;t = 0:1 / Fs: 2;y =唧唧声(t, 100, 1200,gydF4y2Ba“二次”gydF4y2Ba);gydF4y2Ba

估计时变功率谱密度(PSD)的信号。gydF4y2Ba

  • 128年长度将信号分成几个部分,与Kaiser窗窗口的形状参数gydF4y2Ba βgydF4y2Ba =gydF4y2Ba 1gydF4y2Ba 8gydF4y2Ba 。gydF4y2Ba

  • 指定120个样本之间的重叠的部分。gydF4y2Ba

  • 计算谱在gydF4y2Ba ⌊gydF4y2Ba 1gydF4y2Ba 2gydF4y2Ba 8gydF4y2Ba /gydF4y2Ba 2gydF4y2Ba ⌋gydF4y2Ba =gydF4y2Ba 6gydF4y2Ba 5gydF4y2Ba 频率和gydF4y2Ba ⌊gydF4y2Ba (gydF4y2Ba lgydF4y2Ba egydF4y2Ba ngydF4y2Ba ggydF4y2Ba tgydF4y2Ba hgydF4y2Ba (gydF4y2Ba xgydF4y2Ba )gydF4y2Ba - - - - - -gydF4y2Ba 1gydF4y2Ba 2gydF4y2Ba 0gydF4y2Ba )gydF4y2Ba /gydF4y2Ba (gydF4y2Ba 1gydF4y2Ba 2gydF4y2Ba 8gydF4y2Ba - - - - - -gydF4y2Ba 1gydF4y2Ba 2gydF4y2Ba 0gydF4y2Ba )gydF4y2Ba ⌋gydF4y2Ba =gydF4y2Ba 2gydF4y2Ba 3gydF4y2Ba 5gydF4y2Ba 时间箱。gydF4y2Ba

输出的频率和时间每个PSD的重心估计。设置为0的PSD小于这些元素gydF4y2Ba - - - - - -gydF4y2Ba 3gydF4y2Ba 0gydF4y2Ba dB。gydF4y2Ba

[~,~,~,pxx, fc, tc) =光谱图(y,凯瑟(128年,18),120128年,Fs,gydF4y2Ba…gydF4y2Ba“MinThreshold”gydF4y2Ba,-30);gydF4y2Ba

情节的非零元素功能重心频率和时间。gydF4y2Ba

情节(tc (pxx > 0), fc (pxx > 0),gydF4y2Ba“。”gydF4y2Ba)gydF4y2Ba

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

打开生活的脚本gydF4y2Ba

生成一个信号采样在1024 Hz 2秒。gydF4y2Ba

nSamp = 2048;Fs = 1024;t = (0: nSamp-1) / Fs;gydF4y2Ba

在第一第二,信号由一个400 Hz正弦信号和凹二次唧唧喳喳。指定啁啾对区间中点,以便它是对称的,开始和结束在250赫兹的频率达到至少150 Hz。gydF4y2Ba

t1 = t (1: nSamp / 2);x11 =罪(2 *π* 400 * t1);x12 =唧唧声(t1-t1 (nSamp / 4), 150年,nSamp / Fs, 1750年,gydF4y2Ba“二次”gydF4y2Ba);x1 = x11 + x12;gydF4y2Ba

其余的由两个线性啁啾信号降低频率。一个唧唧喳喳的初始频率250赫兹,减少到100赫兹。其他唧唧喳喳的初始频率400赫兹,减少到250赫兹。gydF4y2Ba

t2 = t (nSamp / 2 + 1: nSamp);x21 =唧唧声(t2,400 nSamp / Fs 100);将=唧唧声(t2,550 nSamp / Fs 250);x2 = x21 +将;gydF4y2Ba

添加高斯白噪声的信号。指定一个20 dB的信噪比。重置的随机数字生成器可重复的结果。gydF4y2Ba

信噪比= 20;rng (gydF4y2Ba“默认”gydF4y2Ba)团体= (x1, x2);团体=团体+ randn(大小(sig)) *性病(sig) / db2mag(信噪比);gydF4y2Ba

计算和信号的谱图。指定一个皇帝窗口形状参数的长度为63gydF4y2Ba βgydF4y2Ba =gydF4y2Ba 1gydF4y2Ba 7gydF4y2Ba 10相邻重叠部分的样本较少,和一个FFT长度为256。gydF4y2Ba

nwin = 63;风= kaiser (nwin, 17);nlap = nwin-10;nfft = 256;光谱图(团体、风、nlap nfft, Fs,gydF4y2Ba“桠溪”gydF4y2Ba)gydF4y2Ba

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

阈值谱图,以便任何元素值小于信噪比设置为零。gydF4y2Ba

光谱图(团体、风、nlap nfft, Fs,gydF4y2Ba“MinThreshold”gydF4y2Ba信噪比,gydF4y2Ba“桠溪”gydF4y2Ba)gydF4y2Ba

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

重新分配每个PSD的位置的估算其能量的中心。gydF4y2Ba

光谱图(团体、风、nlap nfft, Fs,gydF4y2Ba“再分配”gydF4y2Ba,gydF4y2Ba“桠溪”gydF4y2Ba)gydF4y2Ba

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

阈值重新分配谱图,以便任何元素值小于信噪比设置为零。gydF4y2Ba

光谱图(团体、风、nlap nfft, Fs,gydF4y2Ba“再分配”gydF4y2Ba,gydF4y2Ba“MinThreshold”gydF4y2Ba信噪比,gydF4y2Ba“桠溪”gydF4y2Ba)gydF4y2Ba

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

打开生活的脚本gydF4y2Ba

加载一个音频信号,其中包含两个减少啾啾和宽带飞溅的声音。计算短时傅里叶变换。把波形分成400 -样本段与300 -样本重叠。绘制光谱图。gydF4y2Ba

负载gydF4y2Ba长条木板gydF4y2Ba%,类型soundsc (y, Fs)gydF4y2Basg = 400;ov = 300;光谱图(y, sg, ov [], Fs,gydF4y2Ba“桠溪”gydF4y2Ba)colormapgydF4y2Ba骨gydF4y2Ba

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

使用gydF4y2Ba光谱图gydF4y2Ba函数输出信号的功率谱密度(PSD)。gydF4y2Ba

(s、f t, p) =光谱图([],y, sg, ov Fs);gydF4y2Ba

跟踪两声响使用gydF4y2BamedfreqgydF4y2Ba函数。发现越强,低频唧唧喳喳,限制搜索频率高于100赫兹和次宽带声音的开始。gydF4y2Ba

f1 = f > 100;t1 = t < 0.75;m1 = medfreq (p (f1, t1), f (f1));gydF4y2Ba

发现微弱高频唧唧喳喳,限制搜索频率高于2500赫兹,乘以0.3秒和0.65秒之间。gydF4y2Ba

f2 = f > 2500;t2 = t & t < 0.3 > 0.65;m2 = medfreq (p (f2, t2), f (f2));gydF4y2Ba

谱图上叠加的结果。频率值除以1000 kHz来表达它们。gydF4y2Ba

持有gydF4y2Ba在gydF4y2Ba情节(t (t1)、m1/1000gydF4y2Ba“线宽”gydF4y2Ba4)图(t (t2)、m2/1000,gydF4y2Ba“线宽”gydF4y2Ba,4)gydF4y2Ba从gydF4y2Ba

图包含一个坐标轴对象。坐标轴对象包含3图像类型的对象,线。gydF4y2Ba

打开生活的脚本gydF4y2Ba

生成两秒信号采样10 kHz。指定信号的瞬时频率作为时间的三角函数。gydF4y2Ba

fs = 10 e3;t = 0:1 / fs: 2;x1 = vco(锯齿(2 *π* t, 0.5), (0.1 - 0.4) * fs, fs);gydF4y2Ba

计算和信号的谱图。使用256 Kaiser窗长度和形状参数gydF4y2Ba βgydF4y2Ba =gydF4y2Ba 5gydF4y2Ba 。指定DFT重叠区到另一个区和512年的220个样本点。情节上的频率gydF4y2BaygydF4y2Ba设在。使用默认colormap和视图。gydF4y2Ba

光谱图(x1,凯瑟(256 5),220512年,fs,gydF4y2Ba“桠溪”gydF4y2Ba)gydF4y2Ba

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

改变视图显示谱图作为一个瀑布。设置colormapgydF4y2Ba骨gydF4y2Ba。gydF4y2Ba

colormap视图(-45、65)gydF4y2Ba骨gydF4y2Ba

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

输入参数gydF4y2Ba

全部折叠gydF4y2Ba

输入信号,指定行或列向量。gydF4y2Ba

例子:gydF4y2Ba因为(π/ 4 * (0:159))+ randn (1160)gydF4y2Ba指定一个正弦信号嵌入在高斯白噪声。gydF4y2Ba

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

窗口中,指定为一个整数或一个行或列向量。使用gydF4y2Ba窗口gydF4y2Ba把信号分成部分:gydF4y2Ba

  • 如果gydF4y2Ba窗口gydF4y2Ba是一个整数,然后呢gydF4y2Ba光谱图gydF4y2Ba分gydF4y2BaxgydF4y2Ba成段的长度gydF4y2Ba窗口gydF4y2Ba和窗户每一部分的汉明窗的长度。gydF4y2Ba

  • 如果gydF4y2Ba窗口gydF4y2Ba是一个矢量,然后呢gydF4y2Ba光谱图gydF4y2Ba分gydF4y2BaxgydF4y2Ba成段长度相同的向量和windows每一部分使用gydF4y2Ba窗口gydF4y2Ba。gydF4y2Ba

如果的长度gydF4y2BaxgydF4y2Ba不能准确划分为一个整数的片段gydF4y2BanoverlapgydF4y2Ba重叠的样本,然后gydF4y2BaxgydF4y2Ba被截断。gydF4y2Ba

如果您指定gydF4y2Ba窗口gydF4y2Ba是空的,那么gydF4y2Ba光谱图gydF4y2Ba使用一个汉明窗,这样gydF4y2BaxgydF4y2Ba分为八段gydF4y2BanoverlapgydF4y2Ba重叠的样本。gydF4y2Ba

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

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

重叠的样本数量,指定为一个正整数。gydF4y2Ba

  • 如果gydF4y2Ba窗口gydF4y2Ba是标量,那么gydF4y2BanoverlapgydF4y2Ba必须小于gydF4y2Ba窗口gydF4y2Ba。gydF4y2Ba

  • 如果gydF4y2Ba窗口gydF4y2Ba是一个矢量,然后呢gydF4y2BanoverlapgydF4y2Ba的长度必须小于gydF4y2Ba窗口gydF4y2Ba。gydF4y2Ba

如果您指定gydF4y2BanoverlapgydF4y2Ba是空的,那么gydF4y2Ba光谱图gydF4y2Ba使用许多生产50%重叠部分。如果未指定的段长度,函数集gydF4y2BanoverlapgydF4y2Ba来gydF4y2Ba⌊gydF4y2BaNgydF4y2BaxgydF4y2Ba/ 4.5⌋gydF4y2Ba,在那里gydF4y2BaNgydF4y2BaxgydF4y2Ba是输入信号的长度吗gydF4y2Ba⌊⌋gydF4y2Ba符号表示函数的地板上。gydF4y2Ba

数量的DFT点,指定为一个正整数标量。如果您指定gydF4y2BanfftgydF4y2Ba是空的,那么gydF4y2Ba光谱图gydF4y2Ba将参数设置为gydF4y2Ba马克斯(256 2gydF4y2BapgydF4y2Ba)gydF4y2Ba,在那里gydF4y2BapgydF4y2Ba=⌈日志gydF4y2Ba2gydF4y2BaNgydF4y2BawgydF4y2Ba⌉gydF4y2Ba,gydF4y2Ba⌈⌉gydF4y2Ba符号表示天花板功能,gydF4y2Ba

  • NgydF4y2BawgydF4y2Ba=gydF4y2Ba窗口gydF4y2Ba如果gydF4y2Ba窗口gydF4y2Ba是一个标量。gydF4y2Ba

  • NgydF4y2BawgydF4y2Ba=gydF4y2Ba长度(gydF4y2Ba窗口gydF4y2Ba)gydF4y2Ba如果gydF4y2Ba窗口gydF4y2Ba是一个向量。gydF4y2Ba

归一化频率,指定为一个向量。gydF4y2BawgydF4y2Ba必须至少有两个元素,否则函数解释它吗gydF4y2BanfftgydF4y2Ba。归一化频率在rad /样品。gydF4y2Ba

例子:gydF4y2Baπ。/(24]gydF4y2Ba

循环频率,指定为一个向量。gydF4y2BafgydF4y2Ba必须至少有两个元素,否则函数解释它吗gydF4y2BanfftgydF4y2Ba。的单位gydF4y2BafgydF4y2Ba指定采样率,gydF4y2BafsgydF4y2Ba。gydF4y2Ba

采样率,指定为一个积极的标量。采样率是单位时间内样品的数量。如果时间的单位是秒,然后采样率在赫兹。gydF4y2Ba

频率范围的PSD估计,指定为gydF4y2Ba“单向的”gydF4y2Ba,gydF4y2Ba双侧的gydF4y2Ba,或gydF4y2Ba“中心”gydF4y2Ba。对于实值信号,默认gydF4y2Ba“单向的”gydF4y2Ba。对于复数信号,默认gydF4y2Ba双侧的gydF4y2Ba,并指定gydF4y2Ba“单向的”gydF4y2Ba会导致错误。gydF4y2Ba

  • “单向的”gydF4y2Ba——返回片面的实际输入信号的谱图。如果gydF4y2BanfftgydF4y2Ba是偶数,那么gydF4y2BapsgydF4y2Ba有gydF4y2BanfftgydF4y2Ba/ 2 + 1行和计算时间间隔gydF4y2Ba[0,gydF4y2BaπgydF4y2Ba]gydF4y2Barad /样品。如果gydF4y2BanfftgydF4y2Ba是奇数,那么gydF4y2BapsgydF4y2Ba有(gydF4y2BanfftgydF4y2Ba+ 1)/ 2行和间隔gydF4y2Ba[0,gydF4y2BaπgydF4y2Ba)gydF4y2Barad /样品。如果您指定gydF4y2BafsgydF4y2Ba,然后间隔分别为[0,gydF4y2BafsgydF4y2Ba/ 2)周期/单位时间(0,gydF4y2BafsgydF4y2Ba/ 2)周期/单位时间。gydF4y2Ba

  • 双侧的gydF4y2Ba——返回真实的或复杂的双边谱图信号。gydF4y2BapsgydF4y2Ba有gydF4y2BanfftgydF4y2Ba行和计算时间间隔gydF4y2Ba[0,2gydF4y2BaπgydF4y2Ba)gydF4y2Barad /样品。如果您指定gydF4y2BafsgydF4y2Ba,区间[0,gydF4y2BafsgydF4y2Ba)周期/单位时间。gydF4y2Ba

  • “中心”gydF4y2Ba——返回真实的或复杂的集中双边谱图信号。gydF4y2BapsgydF4y2Ba有gydF4y2BanfftgydF4y2Ba行。如果gydF4y2BanfftgydF4y2Ba是偶数,那么gydF4y2BapsgydF4y2Ba计算的时间间隔gydF4y2Ba(-gydF4y2BaπgydF4y2Ba,gydF4y2BaπgydF4y2Ba]gydF4y2Barad /样品。如果gydF4y2BanfftgydF4y2Ba是奇数,那么gydF4y2BapsgydF4y2Ba计算在gydF4y2Ba(-gydF4y2BaπgydF4y2Ba,gydF4y2BaπgydF4y2Ba)gydF4y2Barad /样品。如果您指定gydF4y2BafsgydF4y2Ba,然后间隔分别为(-gydF4y2BafsgydF4y2Ba/ 2,gydF4y2BafsgydF4y2Ba/ /单位时间和(- 2)周期gydF4y2BafsgydF4y2Ba/ 2,gydF4y2BafsgydF4y2Ba/ 2)周期/单位时间。gydF4y2Ba

功率谱扩展,指定为gydF4y2Bapsd的gydF4y2Ba或gydF4y2Ba“权力”gydF4y2Ba。gydF4y2Ba

  • 省略gydF4y2BaspectrumtypegydF4y2Ba,或指定gydF4y2Bapsd的gydF4y2Ba,返回功率谱密度。gydF4y2Ba

  • 指定gydF4y2Ba“权力”gydF4y2Ba尺度的每个估计PSD等效噪声带宽的窗口。结果是一个在每个频率估计的力量。如果gydF4y2Ba“重新分配”gydF4y2Ba选择是,功能集成了PSD在每个频率的宽度本之前重新分配。gydF4y2Ba

频率显示轴,指定为gydF4y2Ba“xaxis”gydF4y2Ba或gydF4y2Ba“桠溪”gydF4y2Ba。gydF4y2Ba

  • “xaxis”gydF4y2Ba——显示频率gydF4y2BaxgydF4y2Ba设在和时间的gydF4y2BaygydF4y2Ba设在。gydF4y2Ba

  • “桠溪”gydF4y2Ba——显示频率gydF4y2BaygydF4y2Ba设在和时间的gydF4y2BaxgydF4y2Ba设在。gydF4y2Ba

如果你调用这个论点被忽略gydF4y2Ba光谱图gydF4y2Ba输出参数。gydF4y2Ba

名称-值参数gydF4y2Ba

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

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

例子:gydF4y2Ba光谱图(x, 100,“OutputTimeDimension”,“downrows”)gydF4y2Ba分gydF4y2BaxgydF4y2Ba成段的长度100和windows每一部分汉明窗的长度的输出光谱图有时间维度的行。gydF4y2Ba

阈值,指定为逗号分隔组成的gydF4y2BaMinThresholdgydF4y2Ba和一个真正的标量用分贝表示。gydF4y2Ba光谱图gydF4y2Ba这些元素设置为零gydF4y2Ba年代gydF4y2Ba这样10日志gydF4y2Ba10gydF4y2Ba(gydF4y2Ba年代gydF4y2Ba)≤gydF4y2Ba打gydF4y2Ba。gydF4y2Ba

输出时间维度,指定为逗号分隔组成的gydF4y2BaOutputTimeDimensiongydF4y2Ba和gydF4y2BaacrosscolumnsgydF4y2Ba或gydF4y2BadownrowsgydF4y2Ba。将这个值设置为gydF4y2BadownrowsgydF4y2Ba,如果你想要的时间维度gydF4y2Ba年代gydF4y2Ba,gydF4y2BapsgydF4y2Ba,gydF4y2Ba足球俱乐部gydF4y2Ba,gydF4y2BatcgydF4y2Ba行和频率沿列维度。将这个值设置为gydF4y2BaacrosscolumnsgydF4y2Ba,如果你想要的时间维度gydF4y2Ba年代gydF4y2Ba,gydF4y2BapsgydF4y2Ba,gydF4y2Ba足球俱乐部gydF4y2Ba,gydF4y2BatcgydF4y2Ba在列和频率沿行维度。这个输入是忽略了如果没有输出参数的函数被调用。gydF4y2Ba

输出参数gydF4y2Ba

全部折叠gydF4y2Ba

短时傅里叶变换,作为一个矩阵返回。时间增加的列gydF4y2Ba年代gydF4y2Ba和频率增加行,从0开始。gydF4y2Ba

  • 如果gydF4y2BaxgydF4y2Ba是一个信号的长度gydF4y2BaNgydF4y2BaxgydF4y2Ba,然后gydF4y2Ba年代gydF4y2Ba有gydF4y2BakgydF4y2Ba列,gydF4y2Ba

    • kgydF4y2Ba=⌊(gydF4y2BaNgydF4y2BaxgydF4y2Ba- - - - - -gydF4y2BanoverlapgydF4y2Ba)/ (gydF4y2Ba窗口gydF4y2Ba- - - - - -gydF4y2BanoverlapgydF4y2Ba)⌋如果gydF4y2Ba窗口gydF4y2Ba是一个标量。gydF4y2Ba

    • kgydF4y2Ba=⌊(gydF4y2BaNgydF4y2BaxgydF4y2Ba- - - - - -gydF4y2BanoverlapgydF4y2Ba)/ (gydF4y2Ba长度(gydF4y2Ba窗口gydF4y2Ba)gydF4y2Ba- - - - - -gydF4y2BanoverlapgydF4y2Ba)⌋如果gydF4y2Ba窗口gydF4y2Ba是一个向量。gydF4y2Ba

  • 如果gydF4y2BaxgydF4y2Ba是真实的,gydF4y2BanfftgydF4y2Ba是偶数,那么gydF4y2Ba年代gydF4y2Ba有(gydF4y2BanfftgydF4y2Ba/ 2 + 1)行。gydF4y2Ba

  • 如果gydF4y2BaxgydF4y2Ba是真实的,gydF4y2BanfftgydF4y2Ba是奇数,那么gydF4y2Ba年代gydF4y2Ba有(gydF4y2BanfftgydF4y2Ba+ 1)/ 2行。gydF4y2Ba

  • 如果gydF4y2BaxgydF4y2Ba很复杂,gydF4y2Ba年代gydF4y2Ba有gydF4y2BanfftgydF4y2Ba行。gydF4y2Ba

请注意gydF4y2Ba

当gydF4y2BafreqrangegydF4y2Ba被设置为gydF4y2Ba“单向的”gydF4y2Ba,gydF4y2Ba光谱图gydF4y2Ba输出gydF4y2Ba年代gydF4y2Ba价值观的积极的奈奎斯特范围和不节约总功率。gydF4y2Ba

年代gydF4y2Ba不影响gydF4y2Ba“重新分配”gydF4y2Ba选择。gydF4y2Ba

归一化频率,作为一个向量返回。gydF4y2BawgydF4y2Ba长度相等的行数的gydF4y2Ba年代gydF4y2Ba。gydF4y2Ba

瞬间的时间,作为一个向量返回。的时间值gydF4y2BatgydF4y2Ba对应于每一段的中点。gydF4y2Ba

循环频率,作为一个向量返回。gydF4y2BafgydF4y2Ba长度相等的行数的gydF4y2Ba年代gydF4y2Ba。gydF4y2Ba

功率谱密度(PSD)或功率谱,作为一个矩阵返回。gydF4y2Ba

  • 如果gydF4y2BaxgydF4y2Ba是真实的,gydF4y2BafreqrangegydF4y2Ba被设置为gydF4y2Ba“单向的”gydF4y2Ba,然后gydF4y2BapsgydF4y2Ba包含PSD的片面修改周期图估计或功率谱的每一部分。函数乘以电源2在所有频率除了0和奈奎斯特频率保护总功率。gydF4y2Ba

  • 如果gydF4y2BaxgydF4y2Ba是复杂的,如果gydF4y2BafreqrangegydF4y2Ba被设置为gydF4y2Ba双侧的gydF4y2Ba或gydF4y2Ba“中心”gydF4y2Ba,或者如果你指定一个向量的频率,然后gydF4y2BapsgydF4y2Ba包含了两方面的改进的周期图估计的每一部分的PSD或功率谱。gydF4y2Ba

能量的中心频率和时间,返回相同大小的矩阵如短时傅里叶变换。如果你不指定采样率,然后的元素gydF4y2Ba足球俱乐部gydF4y2Ba返回作为归一化频率。gydF4y2Ba

提示gydF4y2Ba

如果一个短时傅里叶变换有0,转换为分贝导致不能绘制的负无穷大。为了避免这种潜在的困难,gydF4y2Ba光谱图gydF4y2Ba增加了gydF4y2Ba每股收益gydF4y2Ba短时傅里叶变换的时不带输出参数调用它。gydF4y2Ba

引用gydF4y2Ba

[1]奥本海姆,艾伦·V。,Ronald W. Schafer, and John R. Buck.离散时间信号处理gydF4y2Ba。第二版上。台北:普伦蒂斯霍尔,1999年。gydF4y2Ba

[2]拉宾,劳伦斯·R。,罗纳德·w·谢弗。gydF4y2Ba数字语音信号的处理gydF4y2Ba。恩格尔伍德悬崖,台北:普伦蒂斯·霍尔出版社,1978年。gydF4y2Ba

[3]Chassande-Motin,埃里克,弗朗索瓦钻,帕特里克Flandrin。“重新分配。“在gydF4y2Ba时频分析:概念和方法gydF4y2Ba。由弗朗茨Hlawatsch编辑和弗朗索瓦钻。伦敦:ISTE /约翰·威利和儿子,2008年。gydF4y2Ba

[4]费洛浦用,肖恩。凯莉弗茨。“time-corrected算法计算瞬时频率(重新分配)谱图,与应用程序”。gydF4y2Ba美国声学学会杂志》上gydF4y2Ba。119卷,2006年1月,页360 - 371。gydF4y2Ba

扩展功能gydF4y2Ba

版本历史gydF4y2Ba

之前介绍过的R2006agydF4y2Ba

另请参阅gydF4y2Ba

应用程序gydF4y2Ba

功能gydF4y2Ba

主题gydF4y2Ba