希尔伯特变换和瞬时频率
希尔伯特变换为单组分信号估计信号的瞬时频率。单组分信号在时频平面上由一个描述“岭”。的set of monocomponent signals includes single sinusoids and signals like chirps.
生成一个啁啾取样1 kHz两秒钟。最初指定唧唧声频率是100 Hz,一秒钟后增加到200赫兹。
fs = 1000;t = 0:1 / fs: 2 - 1 / f;y =唧唧声(t、100、1200);
估计的谱图使用短时傅里叶变换实现啁啾pspectrum
函数。描述的信号是一个峰值频率在每个时间点。
pspectrum (y, fs,的谱图)
计算分析信号并区分其相位测量瞬时频率。按比例缩小的导数收益率有意义的估计。
z =希尔伯特(y);instfrq = fs /(2 *π)* diff(打开(角(z)));clf情节(t(2:结束),instfrq) ylim ([0 f / 2])
的instfreq
函数计算并显示瞬时频率在一个步骤。
instfreq (y, fs,“方法”,希尔伯特的)
方法失败时,信号不是单组分。
产生一笔两个60赫兹和90赫兹的频率正弦信号,采样两秒的1023赫兹。计算和绘制光谱图。每个时间点显示了两个组件的存在。
fs = 1023;t = 0:1 / fs: 2 - 1 / f;x =罪(2 *π* 60 * t) +罪(2 *π* 90 * t);pspectrum (x, fs,的谱图)yticks ([90])
计算分析信号及其相区分。放大该地区封闭正弦信号的频率。预测分析信号的瞬时频率正弦信号频率的平均值。
z =希尔伯特(x);instfrq = fs /(2 *π)* diff(打开(角(z)));情节(t(2:结束),instfrq) ylim([90])包含(“时间(s)”)ylabel (的频率(赫兹))
的instfreq
函数还估计平均水平。
instfreq (x, fs,“方法”,希尔伯特的)
估计这两个频率作为时间的函数,使用光谱图
功率谱密度和tfridge
跟踪两个山脊。在tfridge
,指定惩罚改变频率为0.1。
(s、f tt) = pspectrum (x, fs,的谱图);numcomp = 2;(冰箱,~,lr) = tfridge (s、f, 0.1,“NumRidges”,numcomp);pspectrum (x, fs,的谱图)举行在plot3 (tt、冰箱、abs (s (lr)),“线宽”,4)从yticks ([90])