希尔伯特变换和瞬时频率

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])

也可以看看

|

相关话题