Hilbert转换估计仅用于单一组分信号的信号的瞬时频率。单一组件信号通过单个“脊的时频平面描述。该组单一组织信号包括单个正弦波和像啁啾等信号。
生成在1 kHz上采样的啁啾,两秒钟。指定Chirp,因此其频率最初是100 Hz,一秒后增加到200 Hz。
FS = 1000;T = 0:1 / FS:2-1 / FS;y = chirp(t,100,1,200);
使用在中实施的短时傅里叶变换来估计啁啾的谱图pspectrum.
功能。信号通过每个时间点的单个峰值频率很好地描述。
pspectrum(y,fs,'谱图')
计算分析信号并区分其阶段以测量瞬时频率。缩放的衍生物产生了有意义的估计。
z =希尔伯特(y);instfrq = fs /(2 * pi)* diff(未包装(角度(z)));CLF绘图(T(2:end),instfrq)ylim([0 fs / 2])
这instfreq.
功能计算并在一步中显示瞬时频率。
instfreq(y,fs,'方法'那'希尔伯特')
当信号不是单一组分时,该方法失败。
产生频率60 Hz和90 Hz的两个正弦曲线的总和,在1023Hz上采样两秒钟。计算和绘制频谱图。每个时间点都显示了两个组件的存在。
FS = 1023;T = 0:1 / FS:2-1 / FS;X = SIN(2 * PI * 60 * T)+ SIN(2 * PI * 90 * T);pspectrum(x,fs,'谱图')yticks([60 90])
计算分析信号并区分其相位。放大围绕正弦曲线频率的区域。分析信号预测是正弦频率的平均值的瞬时频率。
z = hilbert(x);instfrq = fs /(2 * pi)* diff(未包装(角度(z)));绘图(T(2:结束),instfrq)ylim([60 90])xlabel('时间'')ylabel('频率(Hz)')
这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([60 90])