光谱图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___gydF4y2Baspectrumtype)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

光谱图(gydF4y2Ba___gydF4y2Bafreqloc)gydF4y2Ba

描述gydF4y2Ba

例子gydF4y2Ba

年代gydF4y2Ba=光谱图(gydF4y2BaxgydF4y2Ba)gydF4y2Ba返回gydF4y2Ba短时傅里叶变换gydF4y2Ba(STFT)的输入信号gydF4y2BaxgydF4y2Ba。每一列的gydF4y2Ba年代gydF4y2Ba包含一个短期的估计,time-localized频率的内容gydF4y2BaxgydF4y2Ba。大小的平方gydF4y2Ba年代gydF4y2Ba被称为gydF4y2Ba光谱图gydF4y2Ba的时频表示gydF4y2BaxgydF4y2Ba[1]gydF4y2Ba。gydF4y2Ba

例子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年代gydF4y2Ba,gydF4y2BawgydF4y2Ba,gydF4y2BatgydF4y2Ba)=光谱图(gydF4y2Ba___gydF4y2Ba)gydF4y2Ba返回一个向量的归一化频率,gydF4y2BawgydF4y2Ba一个向量的时间瞬间,gydF4y2BatgydF4y2Ba,STFT的计算。这个语法可以包括从以前的语法输入参数的任意组合。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返回STFT在指定的归一化频率gydF4y2BawgydF4y2Ba。gydF4y2BawgydF4y2Ba必须至少有两个元素,否则函数解释它吗gydF4y2BanfftgydF4y2Ba。gydF4y2Ba

例子gydF4y2Ba

(gydF4y2Ba年代gydF4y2Ba,gydF4y2BafgydF4y2Ba,gydF4y2BatgydF4y2Ba)=光谱图(gydF4y2BaxgydF4y2Ba,gydF4y2Ba窗口gydF4y2Ba,gydF4y2BanoverlapgydF4y2Ba,gydF4y2BafgydF4y2Ba,gydF4y2BafsgydF4y2Ba)gydF4y2Ba返回STFT在指定的频率周期性gydF4y2BafgydF4y2Ba。gydF4y2BafgydF4y2Ba必须至少有两个元素,否则函数解释它吗gydF4y2BanfftgydF4y2Ba。gydF4y2Ba

例子gydF4y2Ba

(gydF4y2Ba___gydF4y2Ba,gydF4y2BapsgydF4y2Ba)=光谱图(gydF4y2Ba___gydF4y2Ba,gydF4y2BaspectrumtypegydF4y2Ba)gydF4y2Ba还返回一个矩阵,gydF4y2BapsgydF4y2Ba,比例的光谱图gydF4y2BaxgydF4y2Ba。gydF4y2Ba

如果您指定gydF4y2BaspectrumtypegydF4y2Ba作为gydF4y2Ba“psd”gydF4y2Ba的每一列gydF4y2BapsgydF4y2Ba包含一个估计功率谱密度(PSD)的窗口的部分。gydF4y2Ba
如果您指定gydF4y2BaspectrumtypegydF4y2Ba作为gydF4y2Ba“权力”gydF4y2Ba的每一列gydF4y2BapsgydF4y2Ba包含一个视窗化段功率谱的估计。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___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没有输出参数的阴谋gydF4y2BapsgydF4y2Ba在当前图窗口分贝。gydF4y2Ba

例子gydF4y2Ba

光谱图(gydF4y2Ba___gydF4y2Ba,gydF4y2BafreqlocgydF4y2Ba)gydF4y2Ba指定的轴绘制频率。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

图包含一个坐标轴对象。坐标轴对象包含样本,ylabel归一化频率(空白乘以πr d i n s / s m p l e)包含一个类型的对象的形象。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`光谱图gydF4y2Ba`功能和STFT的定义gydF4y2Ba

打开生活的脚本gydF4y2Ba

生成一个信号,由复值凸二次啁啾取样在600 Hz 2秒。的唧唧声有一个初始频率250赫兹和最后一个50赫兹的频率。gydF4y2Ba

fs = 6 e2;ts = 0:1 / fs: 2;x =唧唧声(ts 250, ts(结束),50岁,gydF4y2Ba“二次”gydF4y2Ba0,gydF4y2Ba“凸”gydF4y2Ba,gydF4y2Ba“复杂”gydF4y2Ba);gydF4y2Ba

光谱图gydF4y2Ba函数gydF4y2Ba

使用gydF4y2Ba光谱图gydF4y2Ba函数计算STFT的信号。gydF4y2Ba

把信号分成部分,gydF4y2Ba $米gydF4y2Ba =gydF4y2Ba 49gydF4y2Ba$ 样品长。gydF4y2Ba
指定gydF4y2Ba $lgydF4y2Ba =gydF4y2Ba 11gydF4y2Ba$ 毗邻的重叠部分的样本。gydF4y2Ba
丢弃最后,短段。gydF4y2Ba
每一段Bartlett窗口的窗口。gydF4y2Ba
评估每一段的离散傅里叶变换gydF4y2Ba ${NgydF4y2Ba}_{DFTgydF4y2Ba} =gydF4y2Ba 1024年gydF4y2Ba$ 点。默认情况下,gydF4y2Ba光谱图gydF4y2Ba计算双边变换的复值信号。gydF4y2Ba

M = 49;L = 11;g = bartlett (M);Ndft = 1024;f (s, t) =光谱图(x, g, L, Ndft, fs);gydF4y2Ba

使用gydF4y2BawaterplotgydF4y2Ba函数来计算并显示谱图,定义为STFT的大小的平方。gydF4y2Ba

waterplot (s, t)gydF4y2Ba

图包含一个坐标轴对象。坐标轴对象包含频率(赫兹),ylabel时间(s)包含一个补丁类型的对象。gydF4y2Ba

STFT的定义gydF4y2Ba

计算的STFTgydF4y2Ba ${NgydF4y2Ba}_{xgydF4y2Ba}$ 使用定义样本信号。将信号划分为gydF4y2Ba $⌊gydF4y2Ba \frac{{NgydF4y2Ba}_{xgydF4y2Ba} - - - - - -gydF4y2Ba lgydF4y2Ba}{米gydF4y2Ba - - - - - -gydF4y2Ba lgydF4y2Ba} ⌋gydF4y2Ba$ 重叠的部分。窗口每一部分并评估其离散傅里叶变换gydF4y2Ba ${NgydF4y2Ba}_{DFTgydF4y2Ba}$ 点。gydF4y2Ba

(之后,~)=缓冲区(1:长度(x), M, L,gydF4y2Ba“nodelay”gydF4y2Ba);X = fft (X(之后)。* g, Ndft);gydF4y2Ba

计算STFT的时间和频率范围。gydF4y2Ba

找到时间值,除以时间向量重叠的部分。时间值的中点段,每段视为一个区间打开处于这个区间的低端。gydF4y2Ba
找到频率值,指定一个尼奎斯特间隔关闭在零频率和开放在低端。gydF4y2Ba

tbuf = ts(之后);色=意味着(tbuf(2:最终,:));弗林特= 0:fs / Ndft: fs-fs / Ndft;gydF4y2Ba

比较的输出gydF4y2Ba光谱图gydF4y2Ba的定义。显示谱图。gydF4y2Ba

maxdiff = max (max (abs (sx)))gydF4y2Ba

maxdiff = 0gydF4y2Ba

waterplot (X,弗林特,色调)gydF4y2Ba

图包含一个坐标轴对象。坐标轴对象包含频率(赫兹),ylabel时间(s)包含一个补丁类型的对象。gydF4y2Ba

函数gydF4y2Bawaterplot (s, t)gydF4y2Ba%瀑布图的谱图gydF4y2Ba瀑布(f t abs (s)”。^ 2)组(gca, XDir =gydF4y2Ba“反向”gydF4y2Ba、查看= 50[30])包含(gydF4y2Ba“频率(赫兹)”gydF4y2Ba)ylabel (gydF4y2Ba“时间(s)”gydF4y2Ba)gydF4y2Ba结束gydF4y2Ba

比较gydF4y2Ba`光谱图gydF4y2Ba`和gydF4y2Ba`stftgydF4y2Ba`功能gydF4y2Ba

打开生活的脚本gydF4y2Ba

生成一个信号组成的啁啾取样1.4 kHz 2秒。短促声波的频率减少线性测量期间从600赫兹到100赫兹。gydF4y2Ba

fs = 1400;x =唧唧声(0:1 / fs: 2600、2100);gydF4y2Ba

stftgydF4y2Ba违约gydF4y2Ba

计算信号的STFT使用gydF4y2Ba光谱图gydF4y2Ba和gydF4y2BastftgydF4y2Ba功能。使用默认值gydF4y2BastftgydF4y2Ba功能:gydF4y2Ba

把信号分成128 -样本段和窗口每一段周期损害窗口。gydF4y2Ba
指定96个样本之间的重叠的部分。这长度相当于75%的窗口长度。gydF4y2Ba
指定128 DFT点和中心STFT在零频率,频率表示的赫兹。gydF4y2Ba

验证两个结果是相等的。gydF4y2Ba

M = 128;g =损害(M,gydF4y2Ba“周期”gydF4y2Ba);L = 0.75 * M;Ndft = 128;[sp, fp, tp) =光谱图(x, g, L, Ndft, fs,gydF4y2Ba“中心”gydF4y2Ba);[s、f t] = stft (x, fs);dff = max (max (abs (sp-s)))gydF4y2Ba

dff = 0gydF4y2Ba

使用gydF4y2Ba网gydF4y2Ba函数绘制的两个输出。gydF4y2Ba

nexttile网(tp、fp、abs (sp)。^ 2)标题(gydF4y2Ba“声谱图”gydF4y2Ba)(2),轴gydF4y2Ba紧gydF4y2Banexttile网(t、f、abs (s) ^ 2)标题(gydF4y2Ba“stft”gydF4y2Ba)(2),轴gydF4y2Ba紧gydF4y2Ba

图包含2轴对象。坐标轴对象1标题声谱图包含一个类型的对象的表面。坐标轴对象2标题stft包含一个类型的对象的表面。gydF4y2Ba

光谱图gydF4y2Ba违约gydF4y2Ba

重复的计算使用默认值gydF4y2Ba光谱图gydF4y2Ba功能:gydF4y2Ba

把信号划分成段的长度gydF4y2Ba $米gydF4y2Ba =gydF4y2Ba ⌊gydF4y2Ba {NgydF4y2Ba}_{xgydF4y2Ba} /gydF4y2Ba 4gydF4y2Ba 。gydF4y2Ba 5gydF4y2Ba ⌋gydF4y2Ba$ ,在那里gydF4y2Ba ${NgydF4y2Ba}_{xgydF4y2Ba}$ 是信号的长度。汉明窗窗口每一部分。gydF4y2Ba
指定50%重叠部分。gydF4y2Ba
计算FFT,使用gydF4y2Ba $马克斯gydF4y2Ba (gydF4y2Ba 256年gydF4y2Ba,gydF4y2Ba {2gydF4y2Ba}^{⌈gydF4y2Ba {日志gydF4y2Ba}_{2gydF4y2Ba} 米gydF4y2Ba ⌉gydF4y2Ba})gydF4y2Ba$ 点。只有积极的归一化频率计算光谱图。gydF4y2Ba

M =地板(长度(x) / 4.5);g =汉明(M);L =地板(M / 2);Ndft = max (256 2 ^ nextpow2 (M));(sx,外汇,tx) =光谱图(x);(圣,英国《金融时报》,tt) = stft (x,窗口= g, OverlapLength = L,gydF4y2Ba…gydF4y2BaFFTLength = Ndft FrequencyRange =gydF4y2Ba“单向的”gydF4y2Ba);dff = max (max (sx-st))gydF4y2Ba

dff = 0gydF4y2Ba

使用gydF4y2BawaterplotgydF4y2Ba函数绘制的两个输出。把频率轴gydF4y2Ba $πgydF4y2Ba$ 在这两种情况下。为gydF4y2BastftgydF4y2Ba输出,样本数量除以有效采样率,gydF4y2Ba $2gydF4y2Ba πgydF4y2Ba$ 。gydF4y2Ba

图nexttile waterplot (sx, fx /π,tx)标题(gydF4y2Ba“声谱图”gydF4y2Ba)nexttile waterplot(圣,《金融时报》/π,tt /(2 *π))标题(gydF4y2Ba“stft”gydF4y2Ba)gydF4y2Ba

$图包含2轴对象。坐标轴对象1标题声谱图,包含频率/ \π,ylabel样本包含一个补丁类型的对象。坐标轴对象2标题stft,包含频率/ \π,ylabel样本包含一个补丁类型的对象。gydF4y2Ba$

函数gydF4y2Bawaterplot (s, t)gydF4y2Ba%瀑布图的谱图gydF4y2Ba瀑布(f t abs (s)”。^ 2)组(gca, XDir =gydF4y2Ba“反向”gydF4y2Ba、查看= 50[30])包含(gydF4y2Ba“频率/ \π”gydF4y2Ba)ylabel (gydF4y2Ba“样本”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

图包含一个坐标轴对象。坐标轴对象包含时间(s), ylabel频率(赫兹)包含一个类型的对象的形象。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

图包含一个坐标轴对象。坐标轴对象包含时间(s), ylabel频率(赫兹)包含2图像类型的对象,线。gydF4y2Ba

谱图的复杂信号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

图包含一个坐标轴对象。坐标轴对象包含样本,ylabel归一化频率(空白乘以πr d i n s / s m p l e)包含一个类型的对象的形象。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

图包含一个坐标轴对象。坐标轴对象包含样本,ylabel归一化频率(空白乘以πr d i n s / s m p l e)包含一个类型的对象的形象。gydF4y2Ba

比较gydF4y2Ba`光谱图gydF4y2Ba`和gydF4y2Ba`pspectrumgydF4y2Ba`功能gydF4y2Ba

打开生活的脚本gydF4y2Ba

产生一个信号,由一个压控振荡器和三个高斯原子。信号采样gydF4y2Ba ${fgydF4y2Ba}_{年代gydF4y2Ba} =gydF4y2Ba 2gydF4y2Ba$ kHz 2秒。gydF4y2Ba

fs = 2000;tx = 0:1 / fs: 2;gaussFun = @ (A, x,μf) exp (- (xμ)。^ 2 /(2 * 0.03 ^ 2)。*罪(2 *π* f。* x) *”;s = gaussFun ((1 1 1), tx”[0.1 - 0.65 1], [2 6 2] * 100) * 1.5;x = vco(线性调频(tx + 1。0 tx(结束),3)。* exp (2 * (tx-1) ^ 2), [0.1 - 0.4] * fs, fs);x = s + x”;gydF4y2Ba

短时傅里叶变换gydF4y2Ba

使用gydF4y2BapspectrumgydF4y2Ba函数来计算STFT。gydF4y2Ba

把gydF4y2Ba ${NgydF4y2Ba}_{xgydF4y2Ba}$ 样本信号长度的片段gydF4y2Ba $米gydF4y2Ba =gydF4y2Ba 80年gydF4y2Ba$ 样本,相应的时间分辨率gydF4y2Ba $80年gydF4y2Ba /gydF4y2Ba 2000年gydF4y2Ba =gydF4y2Ba 40gydF4y2Ba$ 毫秒。gydF4y2Ba
指定gydF4y2Ba $lgydF4y2Ba =gydF4y2Ba 16gydF4y2Ba$ 样品或相邻重叠部分的20%。gydF4y2Ba
Kaiser窗窗口每一部分,指定一个泄漏gydF4y2Ba $ℓgydF4y2Ba =gydF4y2Ba 0gydF4y2Ba 。gydF4y2Ba 7gydF4y2Ba$ 。gydF4y2Ba

M = 80;L = 16;路= 0.7;F (S, T) = pspectrum (x, fs,gydF4y2Ba“声谱图”gydF4y2Ba,gydF4y2Ba…gydF4y2BaTimeResolution = M / fs, OverlapPercent = L / M * 100,gydF4y2Ba…gydF4y2Ba泄漏=路);gydF4y2Ba

比较的结果gydF4y2Ba光谱图gydF4y2Ba函数。gydF4y2Ba

直接在样本指定窗口长度和重叠。gydF4y2Ba
pspectrumgydF4y2Ba总是使用Kaiser窗gydF4y2Ba $ggydF4y2Ba (gydF4y2Ba ngydF4y2Ba)gydF4y2Ba$ 。泄漏gydF4y2Ba $ℓgydF4y2Ba$ 和形状系数gydF4y2Ba $βgydF4y2Ba$ 窗外是相关的gydF4y2Ba $βgydF4y2Ba =gydF4y2Ba 40gydF4y2Ba \timesgydF4y2Ba (gydF4y2Ba 1gydF4y2Ba - - - - - -gydF4y2Ba ℓgydF4y2Ba)gydF4y2Ba$ 。gydF4y2Ba
pspectrumgydF4y2Ba总是使用gydF4y2Ba ${NgydF4y2Ba}_{DFTgydF4y2Ba} =gydF4y2Ba 1024年gydF4y2Ba$ 当计算离散傅里叶变换。您可以指定这个号码如果你想计算变换在双边或频率范围。然而,对于片面的变换,默认为真实信号,gydF4y2Ba光谱图gydF4y2Ba使用gydF4y2Ba $1024年gydF4y2Ba /gydF4y2Ba 2gydF4y2Ba +gydF4y2Ba 1gydF4y2Ba =gydF4y2Ba 513年gydF4y2Ba$ 点。此外,您可以指定您想要计算的向量的频率变换,在这个例子中。gydF4y2Ba
如果不能准确划分为一个信号gydF4y2Ba $kgydF4y2Ba =gydF4y2Ba ⌊gydF4y2Ba \frac{{NgydF4y2Ba}_{xgydF4y2Ba} - - - - - -gydF4y2Ba lgydF4y2Ba}{米gydF4y2Ba - - - - - -gydF4y2Ba lgydF4y2Ba} ⌋gydF4y2Ba$ 段,gydF4y2Ba光谱图gydF4y2Ba截断信号而gydF4y2BapspectrumgydF4y2Ba垫与零信号,创建一个额外的部分。输出等效,删除最后一段和最后一个元素的向量。gydF4y2Ba
光谱图gydF4y2Ba返回STFT,其大小的平方是光谱图。gydF4y2BapspectrumgydF4y2Ba返回segment-by-segment功率谱,已方但除以一个因子gydF4y2Ba ${\sumgydF4y2Ba}_{ngydF4y2Ba} ggydF4y2Ba (gydF4y2Ba ngydF4y2Ba)gydF4y2Ba$ 在平方。gydF4y2Ba
对于片面的变换,gydF4y2BapspectrumgydF4y2Ba增加了一个额外的因素2的谱图。gydF4y2Ba

g = kaiser (M, 40 * (1-lk));k =(长度(x) - l) /(马丁);gydF4y2Ba如果gydF4y2Bak ~ =地板(k) =年代(:,1:地板(k));T = T(1:地板(k));gydF4y2Ba结束gydF4y2Baf (s, t) =光谱图(x / (g)和*√(2),g, L, f, f);gydF4y2Ba

使用gydF4y2BawaterplotgydF4y2Ba函数显示谱图由两个函数计算。gydF4y2Ba

次要情节(2,1,1)waterplot (sqrt (S)、F T)标题(gydF4y2Ba“pspectrum”gydF4y2Ba次要情节(2,1,2)waterplot (s、f t)标题(gydF4y2Ba“声谱图”gydF4y2Ba)gydF4y2Ba

图包含2轴对象。与标题pspectrum坐标轴对象1,包含频率(赫兹),ylabel时间(s)包含一个补丁类型的对象。坐标轴对象与标题声谱图2,包含频率(赫兹),ylabel时间(s)包含一个补丁类型的对象。gydF4y2Ba

maxd = max (max (abs (abs (s) ^ 2 s)))gydF4y2Ba

maxd = 2.4419 e-08gydF4y2Ba

功率谱和便利的情节gydF4y2Ba

的gydF4y2Ba光谱图gydF4y2Ba函数有第四个参数对应于segment-by-segment功率谱和功率谱密度。类似的输出gydF4y2BapspectrumgydF4y2Ba,gydF4y2BapsgydF4y2Ba争论已经平方,包括归一化因子gydF4y2Ba ${\sumgydF4y2Ba}_{ngydF4y2Ba} ggydF4y2Ba (gydF4y2Ba ngydF4y2Ba)gydF4y2Ba$ 。一边倒的声音(真正的信号,你仍然必须包括额外的2倍。设置缩放参数的函数gydF4y2Ba“权力”gydF4y2Ba。gydF4y2Ba

(~,~,~,ps) =光谱图(x * sqrt (2), g, L, F, F,gydF4y2Ba“权力”gydF4y2Ba);马克斯(abs (S (:) ps (:)))gydF4y2Ba

ans = 2.4419 e-08gydF4y2Ba

不带输出参数调用时,两个gydF4y2BapspectrumgydF4y2Ba和gydF4y2Ba光谱图gydF4y2Ba图的谱图信号分贝。包括2片面的谱图的因子。设置colormaps情节是一样的。设置gydF4y2BaxgydF4y2Ba限制同样的价值观使可见额外的段的末尾gydF4y2BapspectrumgydF4y2Ba情节。在gydF4y2Ba光谱图gydF4y2Ba情节,显示的频率gydF4y2BaygydF4y2Ba设在。gydF4y2Ba

次要情节(2,1,1)pspectrum (x, fs,gydF4y2Ba“声谱图”gydF4y2Ba,gydF4y2Ba…gydF4y2BaTimeResolution = M / fs, OverlapPercent = L / M * 100,gydF4y2Ba…gydF4y2Ba泄漏=路)标题(gydF4y2Ba“pspectrum”gydF4y2Ba)cc =爬;xl = xlim;次要情节(2,1,2)谱图(x * sqrt (2), g, L, F, F,gydF4y2Ba“权力”gydF4y2Ba,gydF4y2Ba“桠溪”gydF4y2Ba)标题(gydF4y2Ba“声谱图”gydF4y2Ba爬(cc) xlim (xl)gydF4y2Ba

图包含2轴对象。与标题pspectrum坐标轴对象1,包含时间(s), ylabel频率(赫兹)包含一个类型的对象的形象。坐标轴对象与标题声谱图2,包含时间(s), ylabel频率(赫兹)包含一个类型的对象的形象。gydF4y2Ba

函数gydF4y2Bawaterplot (s, t)gydF4y2Ba%瀑布图的谱图gydF4y2Ba瀑布(f t abs (s)”。^ 2)组(gca, XDir =gydF4y2Ba“反向”gydF4y2Ba、查看= 50[30])包含(gydF4y2Ba“频率(赫兹)”gydF4y2Ba)ylabel (gydF4y2Ba“时间(s)”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

使用gydF4y2Ba光谱图gydF4y2Ba函数没有输出参数绘制重新分配谱图。显示的频率gydF4y2BaygydF4y2Ba设在和时间的gydF4y2BaxgydF4y2Ba设在。gydF4y2Ba

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

图包含一个坐标轴对象。坐标轴对象包含时间(s), ylabel频率(赫兹)包含一个类型的对象的形象。gydF4y2Ba

重做情节使用gydF4y2Ba显示亮度图像gydF4y2Ba函数。指定gydF4y2BaygydF4y2Ba设在方向的频率值增加从下到上。添加gydF4y2Ba每股收益gydF4y2Ba重新分配的声谱图转换成分贝时避免潜在的负无穷大。gydF4y2Ba

[~,fr, tr, pxx] =光谱图(y,凯瑟(128年,18),120128年,fs,gydF4y2Ba…gydF4y2Ba“重新分配”gydF4y2Ba);显示亮度图像(tr、fr、pow2db (pxx + eps)轴gydF4y2BaxygydF4y2Ba包含(gydF4y2Ba“时间(s)”gydF4y2Ba)ylabel (gydF4y2Ba“频率(赫兹)”gydF4y2Ba)colorbargydF4y2Ba

图包含一个坐标轴对象。坐标轴对象包含时间(s), ylabel频率(赫兹)包含一个类型的对象的形象。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

打开生活的脚本gydF4y2Ba

产生一个信号,由一个实值啁啾取样2 kHz 2秒。gydF4y2Ba

fs = 2000;tx = 0:1 / fs: 2;x = vco(线性调频(tx 0 tx(结束),2)。* exp (3 * (tx-1) ^ 2),gydF4y2Ba…gydF4y2Ba(0.1 - 0.4)* fs, fs)。*损害(长度(tx));gydF4y2Ba

双面的光谱图gydF4y2Ba

计算和情节的双边STFT的信号。gydF4y2Ba

把信号分成部分,gydF4y2Ba $米gydF4y2Ba =gydF4y2Ba 73年gydF4y2Ba$ 样品长。gydF4y2Ba
指定gydF4y2Ba $lgydF4y2Ba =gydF4y2Ba 24gydF4y2Ba$ 毗邻的重叠部分的样本。gydF4y2Ba
丢弃最后,短段。gydF4y2Ba
每一段平顶建筑物的窗口的窗口。gydF4y2Ba
评估每一段的离散傅里叶变换gydF4y2Ba ${NgydF4y2Ba}_{DFTgydF4y2Ba} =gydF4y2Ba 895年gydF4y2Ba$ 点,并指出这是一个奇数。gydF4y2Ba

M = 73;L = 24;g = flattopwin (M);Ndft = 895;乃文= ~国防部(Ndft 2);[stwo f t] =光谱图(x, g, L, Ndft, fs,gydF4y2Ba“双侧”gydF4y2Ba);gydF4y2Ba

使用gydF4y2Ba光谱图gydF4y2Ba函数没有输出参数绘制双边谱图。gydF4y2Ba

光谱图(x, g, L, Ndft, fs,gydF4y2Ba“双侧”gydF4y2Ba,gydF4y2Ba“权力”gydF4y2Ba,gydF4y2Ba“桠溪”gydF4y2Ba);gydF4y2Ba

图包含一个坐标轴对象。坐标轴对象包含时间(s), ylabel频率(赫兹)包含一个类型的对象的形象。gydF4y2Ba

使用定义计算双边谱图。将信号划分为gydF4y2Ba $米gydF4y2Ba$ 示例部分与gydF4y2Ba $lgydF4y2Ba$ 毗邻的重叠部分的样本。窗口每一部分并计算其离散傅里叶变换gydF4y2Ba ${NgydF4y2Ba}_{DFTgydF4y2Ba}$ 点。gydF4y2Ba

(之后,~)=缓冲区(1:长度(x), M, L,gydF4y2Ba“nodelay”gydF4y2Ba);Xtwo = fft (x(之后)。* g, Ndft);gydF4y2Ba

计算时间和频率范围。gydF4y2Ba

找到时间值,除以时间向量重叠的部分。时间值的中点段,每段视为一个区间打开处于这个区间的低端。gydF4y2Ba
找到频率值,指定一个尼奎斯特间隔关闭在零频率和开放上端。gydF4y2Ba

tbuf = tx(之后);ttwo =意味着(tbuf(2:最终,:));ftwo = 0: fs / Ndft: fs * (1 - 1 / Ndft);gydF4y2Ba

比较的输出gydF4y2Ba光谱图gydF4y2Ba的定义。使用gydF4y2BawaterplotgydF4y2Ba函数显示谱图。gydF4y2Ba

差别=[马克斯(max (abs (stwo-Xtwo))) max (abs (f-ftwo '))马克斯(abs (t-ttwo)))gydF4y2Ba

差别=gydF4y2Ba1×3gydF4y2Ba10gydF4y2Ba^{-12年gydF4y2Ba}0×0.2274 - 0.0002gydF4y2Ba

图nexttile waterplot (Xtwo ftwo ttwo)标题(gydF4y2Ba“两面,定义”gydF4y2Ba)nexttile waterplot (stwo f t)标题(gydF4y2Ba“两面,谱图函数”gydF4y2Ba)gydF4y2Ba

图包含2轴对象。坐标轴对象1标题两面,定义,包含频率(赫兹),ylabel时间(s)包含一个补丁类型的对象。坐标轴对象2标题两面,谱图功能,包含频率(赫兹),ylabel时间(s)包含一个补丁类型的对象。gydF4y2Ba

集中光谱图gydF4y2Ba

计算信号的集中光谱图。gydF4y2Ba

使用相同的用于双边STFT的时间值。gydF4y2Ba
使用gydF4y2BafftshiftgydF4y2Ba转移函数的零频分量STFT频谱的中心。gydF4y2Ba
对奇值gydF4y2Ba ${NgydF4y2Ba}_{DFTgydF4y2Ba}$ ,频率区间两端开放。对于even-valuedgydF4y2Ba ${NgydF4y2Ba}_{DFTgydF4y2Ba}$ ,低端的频率间隔是开放和关闭的上端。gydF4y2Ba

比较输出并显示谱图。gydF4y2Ba

tcen = ttwo;gydF4y2Ba如果gydF4y2Ba~乃文Xcen = fftshift (Xtwo, 1);fcen = - f / 2 * (1 - 1 / Ndft): fs / Ndft: fs / 2;gydF4y2Ba其他的gydF4y2BaXcen = fftshift (circshift (Xtwo, 1), 1);fcen = (- f / 2 * (1 - 1 / Ndft): fs / Ndft: fs / 2) + fs / Ndft / 2;gydF4y2Ba结束gydF4y2Ba[scen f t] =光谱图(x, g, L, Ndft, fs,gydF4y2Ba“中心”gydF4y2Ba);差别=[马克斯(max (abs (scen-Xcen))) max (abs (f-fcen '))马克斯(abs (t-tcen)))gydF4y2Ba

差别=gydF4y2Ba1×3gydF4y2Ba10gydF4y2Ba^{-12年gydF4y2Ba}0×0.2274 - 0.0002gydF4y2Ba

图nexttile waterplot (Xcen fcen tcen)标题(gydF4y2Ba“集中,定义”gydF4y2Ba)nexttile waterplot (scen f t)标题(gydF4y2Ba“集中,谱图函数”gydF4y2Ba)gydF4y2Ba

图包含2轴对象。坐标轴对象1标题居中,定义,包含频率(赫兹),ylabel时间(s)包含一个补丁类型的对象。坐标轴对象2标题居中,谱图功能,包含频率(赫兹),ylabel时间(s)包含一个补丁类型的对象。gydF4y2Ba

片面的光谱图gydF4y2Ba

计算信号的片面的光谱图。gydF4y2Ba

使用相同的用于双边STFT的时间值。gydF4y2Ba
对奇值gydF4y2Ba ${NgydF4y2Ba}_{DFTgydF4y2Ba}$ 首先,片面的STFT包括gydF4y2Ba $(gydF4y2Ba {NgydF4y2Ba}_{DFTgydF4y2Ba} +gydF4y2Ba 1gydF4y2Ba)gydF4y2Ba /gydF4y2Ba 2gydF4y2Ba$ 双边STFT的行。对于even-valuedgydF4y2Ba ${NgydF4y2Ba}_{DFTgydF4y2Ba}$ 首先,片面的STFT包括gydF4y2Ba ${NgydF4y2Ba}_{DFTgydF4y2Ba} /gydF4y2Ba 2gydF4y2Ba +gydF4y2Ba 1gydF4y2Ba$ 双边STFT的行。gydF4y2Ba
对奇值gydF4y2Ba ${NgydF4y2Ba}_{DFTgydF4y2Ba}$ 、频率间隔是收在零频率在奈奎斯特频率和开放。对于even-valuedgydF4y2Ba ${NgydF4y2Ba}_{DFTgydF4y2Ba}$ ,频率区间两端封闭。gydF4y2Ba

比较输出并显示谱图。为实值信号gydF4y2Ba“单向的”gydF4y2Ba参数是可选的。gydF4y2Ba

语气= ttwo;gydF4y2Ba如果gydF4y2Ba~乃文Xone = Xtwo (1: (Ndft + 1) / 2,:);gydF4y2Ba其他的gydF4y2BaXone = Xtwo (1: Ndft / 2 + 1,);gydF4y2Ba结束gydF4y2Ba词语快捷键= 0:fs / Ndft: fs / 2;(宋,f, t) =光谱图(x, g, L, Ndft, fs);差别=[马克斯(max (abs (sone-Xone))) max (abs (f-fone '))马克斯(abs (t-tone)))gydF4y2Ba

差别=gydF4y2Ba1×3gydF4y2Ba10gydF4y2Ba^{-12年gydF4y2Ba}0×0.1137 - 0.0002gydF4y2Ba

图nexttile waterplot (Xone词语快捷键,语气)标题(gydF4y2Ba“片面的,定义”gydF4y2Ba)nexttile waterplot(宋f t)标题(gydF4y2Ba“片面的,光谱图函数”gydF4y2Ba)gydF4y2Ba

图包含2轴对象。坐标轴对象1标题片面的,定义,包含频率(赫兹),ylabel时间(s)包含一个补丁类型的对象。坐标轴对象2标题片面,谱图功能,包含频率(赫兹),ylabel时间(s)包含一个补丁类型的对象。gydF4y2Ba

函数gydF4y2Bawaterplot (s, t)gydF4y2Ba%瀑布图的谱图gydF4y2Ba瀑布(f t abs (s)”。^ 2)组(gca, XDir =gydF4y2Ba“反向”gydF4y2Ba、查看= 50[30])包含(gydF4y2Ba“频率(赫兹)”gydF4y2Ba)ylabel (gydF4y2Ba“时间(s)”gydF4y2Ba)gydF4y2Ba结束gydF4y2Ba

计算段PSDs和功率谱gydF4y2Ba

打开生活的脚本gydF4y2Ba

的gydF4y2Ba光谱图gydF4y2Ba函数矩阵包含功率谱密度(PSD)或每段作为第四输出的功率谱参数。功率谱等于PSD乘以等效噪声带宽(ENBW)的窗口。gydF4y2Ba

产生一个信号,由一个对数啁啾取样为1秒1 kHz。唧唧喳喳的初始频率400赫兹,减少到10赫兹的测量。gydF4y2Ba

fs = 1000;tt = 0:1 / fs: 1 - 1 / f;y =唧唧声(tt 400 tt(结束),10日gydF4y2Ba“对数”gydF4y2Ba);gydF4y2Ba

段PSDs与采样率和功率谱gydF4y2Ba

把信号分成102 -样本段损害窗口和窗口每一部分。指定12样本之间的重叠相邻段和DFT 1024点。gydF4y2Ba

M = 102;g =损害(M);L = 12;Ndft = 1024;gydF4y2Ba

计算信号的谱图用默认PSD光谱类型。输出的STFT和数组段功率谱密度。gydF4y2Ba

(s、f t, p) =光谱图(y, g, L, Ndft, fs);gydF4y2Ba

重复计算的光谱类型指定为gydF4y2Ba“权力”gydF4y2Ba。输出的STFT和数组段功率谱。gydF4y2Ba

(r, ~, ~ q) =光谱图(y, g, L, Ndft, fs,gydF4y2Ba“权力”gydF4y2Ba);gydF4y2Ba

验证谱图在这两种情况下是相同的。图的谱图用对数刻度的频率。gydF4y2Ba

max (max (abs(年代)。^ 2-abs (r)。^ 2))gydF4y2Ba

ans = 0gydF4y2Ba

瀑布(f t abs (s)”。^ 2)组(gca, XScale =gydF4y2Ba“日志”gydF4y2Ba,gydF4y2Ba…gydF4y2BaXDir =gydF4y2Ba“反向”gydF4y2Ba、查看= 50 [30])gydF4y2Ba

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

验证功率谱等于功率谱密度乘以ENBW窗外。gydF4y2Ba

max (max (abs (qp * enbw (g, fs))))gydF4y2Ba

ans = 1.1102 e-16gydF4y2Ba

验证段功率谱矩阵的谱图成正比。比例系数的平方之和窗口元素。gydF4y2Ba

max (max (abs(年代)。^ 2 q *总和(g) ^ 2))gydF4y2Ba

ans = 3.4694 e-18gydF4y2Ba

段PSDs与归一化频率和功率谱gydF4y2Ba

重复计算,但现在在归一化频率工作。结果是相同的,当你指定的采样率gydF4y2Ba $2gydF4y2Ba πgydF4y2Ba$ 。gydF4y2Ba

(~,~,~,pn) =光谱图(y, g, L, Ndft);[~,~,~,qn] =光谱图(y, g, L, Ndft,gydF4y2Ba“权力”gydF4y2Ba);max (max (abs (qn-pn * enbw (g, 2 *π))))gydF4y2Ba

ans = 1.1102 e-16gydF4y2Ba

在音频信号跟踪啾啾gydF4y2Ba

打开生活的脚本gydF4y2Ba

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

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

图包含一个坐标轴对象。坐标轴对象包含时间(s), ylabel频率(赫兹)包含一个类型的对象的形象。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/1000线宽= 4)情节(t (t2), m2/1000,线宽= 4)gydF4y2Ba从gydF4y2Ba

图包含一个坐标轴对象。坐标轴对象包含时间(s), ylabel频率(赫兹)包含3图像类型的对象,线。gydF4y2Ba

三维谱图可视化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

图包含一个坐标轴对象。坐标轴对象包含时间(s), ylabel频率(赫兹)包含一个类型的对象的形象。gydF4y2Ba

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

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

图包含一个坐标轴对象。坐标轴对象包含时间(s), ylabel频率(赫兹)包含一个对象类型的表面。gydF4y2Ba

输入参数gydF4y2Ba

全部折叠gydF4y2Ba

`xgydF4y2Ba`- - - - - -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`[]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

`noverlapgydF4y2Ba`- - - - - -gydF4y2Ba样本数量的重叠gydF4y2Ba
正整数gydF4y2Ba|gydF4y2Ba`[]gydF4y2Ba`

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

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

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

`nfftgydF4y2Ba`- - - - - -gydF4y2Ba数量的DFT点gydF4y2Ba
正整数标量gydF4y2Ba|gydF4y2Ba`[]gydF4y2Ba`

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

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

`wgydF4y2Ba`- - - - - -gydF4y2Ba归一化频率gydF4y2Ba
向量gydF4y2Ba

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

例子:gydF4y2Baπ。/(24]gydF4y2Ba

`fgydF4y2Ba`- - - - - -gydF4y2Ba周期性的频率gydF4y2Ba
向量gydF4y2Ba

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

`fsgydF4y2Ba`- - - - - -gydF4y2Ba采样率gydF4y2Ba
1赫兹gydF4y2Ba(默认)|gydF4y2Ba积极的标量gydF4y2Ba

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

`freqrangegydF4y2Ba`- - - - - -gydF4y2Ba频率范围为PSD估计gydF4y2Ba
`“单向的”gydF4y2Ba`|gydF4y2Ba`“双侧”gydF4y2Ba`|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

`spectrumtypegydF4y2Ba`- - - - - -gydF4y2Ba功率谱扩展gydF4y2Ba
`“psd”gydF4y2Ba`(默认)|gydF4y2Ba`“权力”gydF4y2Ba`

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

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

`freqlocgydF4y2Ba`- - - - - -gydF4y2Ba频率显示轴gydF4y2Ba
`“xaxis”gydF4y2Ba`(默认)|gydF4y2Ba`“桠溪”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

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

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

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

`MinThresholdgydF4y2Ba`- - - - - -gydF4y2Ba阈值gydF4y2Ba
`负gydF4y2Ba`(默认)|gydF4y2Ba真正的标量gydF4y2Ba

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

`OutputTimeDimensiongydF4y2Ba`- - - - - -gydF4y2Ba输出时间维度gydF4y2Ba
`“acrosscolumns”gydF4y2Ba`(默认)|gydF4y2Ba`“downrows”gydF4y2Ba`

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

输出参数gydF4y2Ba

全部折叠gydF4y2Ba

`年代gydF4y2Ba`——短时傅里叶变换gydF4y2Ba
矩阵gydF4y2Ba

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

如果gydF4y2BaxgydF4y2Ba是一个信号的长度gydF4y2BaNgydF4y2Ba_xgydF4y2Ba,然后gydF4y2Ba年代gydF4y2Ba有gydF4y2BakgydF4y2Ba列,gydF4y2Ba
- kgydF4y2Ba=⌊(gydF4y2BaNgydF4y2Ba_xgydF4y2Ba- - - - - -gydF4y2BanoverlapgydF4y2Ba)/ (gydF4y2Ba窗口gydF4y2Ba- - - - - -gydF4y2BanoverlapgydF4y2Ba)⌋如果gydF4y2Ba窗口gydF4y2Ba是一个标量。gydF4y2Ba
- kgydF4y2Ba=⌊(gydF4y2BaNgydF4y2Ba_xgydF4y2Ba- - - - - -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

`wgydF4y2Ba`-归一化频率gydF4y2Ba
向量gydF4y2Ba

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

`tgydF4y2Ba`——时间瞬间gydF4y2Ba
向量gydF4y2Ba

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

`fgydF4y2Ba`——循环的频率gydF4y2Ba
向量gydF4y2Ba

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

`psgydF4y2Ba`——功率谱密度或功率谱gydF4y2Ba
矩阵gydF4y2Ba

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

如果gydF4y2BaxgydF4y2Ba是真实的,gydF4y2BafreqrangegydF4y2Ba离开不明或一组吗gydF4y2Ba“单向的”gydF4y2Ba,然后gydF4y2BapsgydF4y2Ba包含PSD的片面修改周期图估计或功率谱的每一部分。函数乘以电源2在所有频率除了0和奈奎斯特频率保护总功率。gydF4y2Ba
如果gydF4y2BaxgydF4y2Ba是复数的或者gydF4y2BafreqrangegydF4y2Ba被设置为gydF4y2Ba“双侧”gydF4y2Ba或gydF4y2Ba“中心”gydF4y2Ba,然后gydF4y2BapsgydF4y2Ba包含了两方面的改进的周期图估计的每一部分的PSD或功率谱。gydF4y2Ba
如果你指定一个向量的归一化频率gydF4y2BawgydF4y2Ba或一个向量的周期性的频率gydF4y2BafgydF4y2Ba,然后gydF4y2BapsgydF4y2Ba包含了PSD的修改周期图估计或功率谱的每一部分评估输入频率。gydF4y2Ba

`足球俱乐部gydF4y2Ba`,gydF4y2Ba`tcgydF4y2Ba`——能量的中心频率和时间gydF4y2Ba
矩阵gydF4y2Ba

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

提示gydF4y2Ba

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

引用gydF4y2Ba

[1]Boashash Boualem ed。gydF4y2Ba时间频率信号分析和处理:一个全面的参考gydF4y2Ba。第二版。EURASIP和学术出版社系列在信号和图像处理。阿姆斯特丹和波士顿:学术出版社,2016年。gydF4y2Ba

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

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

[4]奥本海姆,艾伦·V。,R在ald W. Schafer, with John R. Buck.离散时间信号处理gydF4y2Ba。第二版。上台北:Prentice Hall出版社,1999年。gydF4y2Ba

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

扩展功能gydF4y2Ba

高大的数组gydF4y2Ba
计算和数组的行比装入内存。gydF4y2Ba

使用笔记和限制:gydF4y2Ba

输入必须是一个高大的列向量。gydF4y2Ba
的gydF4y2Ba窗口gydF4y2Ba参数必须指定。gydF4y2Ba
OutputTimeDimensiongydF4y2Ba必须始终指定和集吗gydF4y2Ba“downrows”gydF4y2Ba。gydF4y2Ba
的gydF4y2Ba重新分配gydF4y2Ba不支持选项。金宝appgydF4y2Ba
不支持语法没有输出参数。金宝appgydF4y2Ba

有关更多信息,请参见gydF4y2Ba高大的数组gydF4y2Ba。gydF4y2Ba

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

使用笔记和限制:gydF4y2Ba

参数使用名称-值指定参数必须编译时常量。gydF4y2Ba
变量的大小gydF4y2Ba窗口gydF4y2Ba必须是双精度。gydF4y2Ba

线程环境gydF4y2Ba
在后台运行代码使用MATLAB®gydF4y2Ba`backgroundPoolgydF4y2Ba`与并行计算工具箱™或加速代码gydF4y2Ba`ThreadPoolgydF4y2Ba`。gydF4y2Ba

使用笔记和限制:gydF4y2Ba

不支持语法没有输出参数。金宝appgydF4y2Ba
必须均匀间隔的频率向量。gydF4y2Ba

有关更多信息,请参见gydF4y2BaMATLAB函数线程环境中运行gydF4y2Ba。gydF4y2Ba

GPU数组gydF4y2Ba
加速代码运行在一个图形处理单元(GPU)使用并行计算工具箱™。gydF4y2Ba

这个函数完全支持GPU数组。金宝app有关更多信息,请参见gydF4y2Ba运行在GPU MATLAB函数gydF4y2Ba(并行计算工具箱)gydF4y2Ba。gydF4y2Ba

版本历史gydF4y2Ba

之前介绍过的R2006agydF4y2Ba

全部展开gydF4y2Ba

R2023a:gydF4y2Ba可视化输出函数使用gydF4y2Ba创建图gydF4y2Ba住编辑任务gydF4y2Ba

您现在可以使用gydF4y2Ba创建图gydF4y2Ba现场可视化编辑任务的输出gydF4y2Ba光谱图gydF4y2Ba交互。你可以选择不同的图表类型和设置可选的参数。任务也自动生成代码,成为你生活的一部分脚本。gydF4y2Ba

另请参阅gydF4y2Ba

应用程序gydF4y2Ba

信号分析仪gydF4y2Ba

功能gydF4y2Ba

goertzel写到gydF4y2Ba|gydF4y2BaistftgydF4y2Ba|gydF4y2Ba周期图gydF4y2Ba|gydF4y2BapspectrumgydF4y2Ba|gydF4y2BapwelchgydF4y2Ba|gydF4y2BastftgydF4y2Ba|gydF4y2BaxspectrogramgydF4y2Ba

主题gydF4y2Ba

外部网站gydF4y2Ba

傅里叶分析(gydF4y2BaMathWorksgydF4y2Ba教学资源)gydF4y2Ba

光谱图gydF4y2Ba

语法gydF4y2Ba

描述gydF4y2Ba

例子gydF4y2Ba

默认值的光谱图gydF4y2Ba

比较gydF4y2Ba光谱图gydF4y2Ba功能和STFT的定义gydF4y2Ba

比较gydF4y2Ba光谱图gydF4y2Ba和gydF4y2BastftgydF4y2Ba功能gydF4y2Ba

光谱图和瞬时频率gydF4y2Ba

谱图的复杂信号gydF4y2Ba

比较gydF4y2Ba光谱图gydF4y2Ba和gydF4y2BapspectrumgydF4y2Ba功能gydF4y2Ba

重新分配声谱图的二次唧唧声gydF4y2Ba

光谱图与阈值gydF4y2Ba

计算中心和片面的声谱图gydF4y2Ba

计算段PSDs和功率谱gydF4y2Ba

在音频信号跟踪啾啾gydF4y2Ba

三维谱图可视化gydF4y2Ba

输入参数gydF4y2Ba

xgydF4y2Ba- - - - - -gydF4y2Ba输入信号gydF4y2Ba向量gydF4y2Ba

窗口gydF4y2Ba- - - - - -gydF4y2Ba窗口gydF4y2Ba整数gydF4y2Ba|gydF4y2Ba向量gydF4y2Ba|gydF4y2Ba[]gydF4y2Ba

noverlapgydF4y2Ba- - - - - -gydF4y2Ba样本数量的重叠gydF4y2Ba正整数gydF4y2Ba|gydF4y2Ba[]gydF4y2Ba

nfftgydF4y2Ba- - - - - -gydF4y2Ba数量的DFT点gydF4y2Ba正整数标量gydF4y2Ba|gydF4y2Ba[]gydF4y2Ba

wgydF4y2Ba- - - - - -gydF4y2Ba归一化频率gydF4y2Ba向量gydF4y2Ba

fgydF4y2Ba- - - - - -gydF4y2Ba周期性的频率gydF4y2Ba向量gydF4y2Ba

fsgydF4y2Ba- - - - - -gydF4y2Ba采样率gydF4y2Ba1赫兹gydF4y2Ba(默认)|gydF4y2Ba积极的标量gydF4y2Ba

freqrangegydF4y2Ba- - - - - -gydF4y2Ba频率范围为PSD估计gydF4y2Ba“单向的”gydF4y2Ba|gydF4y2Ba“双侧”gydF4y2Ba|gydF4y2Ba“中心”gydF4y2Ba

spectrumtypegydF4y2Ba- - - - - -gydF4y2Ba功率谱扩展gydF4y2Ba“psd”gydF4y2Ba(默认)|gydF4y2Ba“权力”gydF4y2Ba

freqlocgydF4y2Ba- - - - - -gydF4y2Ba频率显示轴gydF4y2Ba“xaxis”gydF4y2Ba(默认)|gydF4y2Ba“桠溪”gydF4y2Ba

名称-值参数gydF4y2Ba

MinThresholdgydF4y2Ba- - - - - -gydF4y2Ba阈值gydF4y2Ba负gydF4y2Ba(默认)|gydF4y2Ba真正的标量gydF4y2Ba

OutputTimeDimensiongydF4y2Ba- - - - - -gydF4y2Ba输出时间维度gydF4y2Ba“acrosscolumns”gydF4y2Ba(默认)|gydF4y2Ba“downrows”gydF4y2Ba

输出参数gydF4y2Ba

年代gydF4y2Ba——短时傅里叶变换gydF4y2Ba矩阵gydF4y2Ba

wgydF4y2Ba-归一化频率gydF4y2Ba向量gydF4y2Ba

tgydF4y2Ba——时间瞬间gydF4y2Ba向量gydF4y2Ba

fgydF4y2Ba——循环的频率gydF4y2Ba向量gydF4y2Ba

psgydF4y2Ba——功率谱密度或功率谱gydF4y2Ba矩阵gydF4y2Ba

足球俱乐部gydF4y2Ba,gydF4y2BatcgydF4y2Ba——能量的中心频率和时间gydF4y2Ba矩阵gydF4y2Ba

更多关于gydF4y2Ba

短时傅里叶变换gydF4y2Ba

提示gydF4y2Ba

引用gydF4y2Ba

扩展功能gydF4y2Ba

高大的数组gydF4y2Ba计算和数组的行比装入内存。gydF4y2Ba

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

线程环境gydF4y2Ba在后台运行代码使用MATLAB®gydF4y2BabackgroundPoolgydF4y2Ba与并行计算工具箱™或加速代码gydF4y2BaThreadPoolgydF4y2Ba。gydF4y2Ba

GPU数组gydF4y2Ba加速代码运行在一个图形处理单元(GPU)使用并行计算工具箱™。gydF4y2Ba

版本历史gydF4y2Ba

R2023a:gydF4y2Ba可视化输出函数使用gydF4y2Ba创建图gydF4y2Ba住编辑任务gydF4y2Ba

另请参阅gydF4y2Ba

应用程序gydF4y2Ba

功能gydF4y2Ba

主题gydF4y2Ba

外部网站gydF4y2Ba

比较gydF4y2Ba`光谱图gydF4y2Ba`功能和STFT的定义gydF4y2Ba

比较gydF4y2Ba`光谱图gydF4y2Ba`和gydF4y2Ba`stftgydF4y2Ba`功能gydF4y2Ba

比较gydF4y2Ba`光谱图gydF4y2Ba`和gydF4y2Ba`pspectrumgydF4y2Ba`功能gydF4y2Ba

`xgydF4y2Ba`- - - - - -gydF4y2Ba输入信号gydF4y2Ba
向量gydF4y2Ba

`窗口gydF4y2Ba`- - - - - -gydF4y2Ba窗口gydF4y2Ba
整数gydF4y2Ba|gydF4y2Ba向量gydF4y2Ba|gydF4y2Ba`[]gydF4y2Ba`

`noverlapgydF4y2Ba`- - - - - -gydF4y2Ba样本数量的重叠gydF4y2Ba
正整数gydF4y2Ba|gydF4y2Ba`[]gydF4y2Ba`

`nfftgydF4y2Ba`- - - - - -gydF4y2Ba数量的DFT点gydF4y2Ba
正整数标量gydF4y2Ba|gydF4y2Ba`[]gydF4y2Ba`

`wgydF4y2Ba`- - - - - -gydF4y2Ba归一化频率gydF4y2Ba
向量gydF4y2Ba

`fgydF4y2Ba`- - - - - -gydF4y2Ba周期性的频率gydF4y2Ba
向量gydF4y2Ba

`fsgydF4y2Ba`- - - - - -gydF4y2Ba采样率gydF4y2Ba
1赫兹gydF4y2Ba(默认)|gydF4y2Ba积极的标量gydF4y2Ba

`freqrangegydF4y2Ba`- - - - - -gydF4y2Ba频率范围为PSD估计gydF4y2Ba
`“单向的”gydF4y2Ba`|gydF4y2Ba`“双侧”gydF4y2Ba`|gydF4y2Ba`“中心”gydF4y2Ba`

`spectrumtypegydF4y2Ba`- - - - - -gydF4y2Ba功率谱扩展gydF4y2Ba
`“psd”gydF4y2Ba`(默认)|gydF4y2Ba`“权力”gydF4y2Ba`

`freqlocgydF4y2Ba`- - - - - -gydF4y2Ba频率显示轴gydF4y2Ba
`“xaxis”gydF4y2Ba`(默认)|gydF4y2Ba`“桠溪”gydF4y2Ba`

`MinThresholdgydF4y2Ba`- - - - - -gydF4y2Ba阈值gydF4y2Ba
`负gydF4y2Ba`(默认)|gydF4y2Ba真正的标量gydF4y2Ba

`OutputTimeDimensiongydF4y2Ba`- - - - - -gydF4y2Ba输出时间维度gydF4y2Ba
`“acrosscolumns”gydF4y2Ba`(默认)|gydF4y2Ba`“downrows”gydF4y2Ba`

`年代gydF4y2Ba`——短时傅里叶变换gydF4y2Ba
矩阵gydF4y2Ba

`wgydF4y2Ba`-归一化频率gydF4y2Ba
向量gydF4y2Ba

`tgydF4y2Ba`——时间瞬间gydF4y2Ba
向量gydF4y2Ba

`fgydF4y2Ba`——循环的频率gydF4y2Ba
向量gydF4y2Ba

`psgydF4y2Ba`——功率谱密度或功率谱gydF4y2Ba
矩阵gydF4y2Ba

`足球俱乐部gydF4y2Ba`,gydF4y2Ba`tcgydF4y2Ba`——能量的中心频率和时间gydF4y2Ba
矩阵gydF4y2Ba

高大的数组gydF4y2Ba
计算和数组的行比装入内存。gydF4y2Ba

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

线程环境gydF4y2Ba
在后台运行代码使用MATLAB®gydF4y2Ba`backgroundPoolgydF4y2Ba`与并行计算工具箱™或加速代码gydF4y2Ba`ThreadPoolgydF4y2Ba`。gydF4y2Ba

GPU数组gydF4y2Ba
加速代码运行在一个图形处理单元(GPU)使用并行计算工具箱™。gydF4y2Ba