首页 文章

在MATLAB中绘制二维希尔伯特谱

提问于
浏览
2

我想创建一个希尔伯特谱的二维图 . 我想要的是时间与频率的关系图,其中信号的幅度由图中的颜色变化表示 .

我所做的是this,但我需要它看起来像this .

先感谢您

编辑:

我有使用希尔伯特黄变换生成的振幅和瞬时频率随时间变化的值 . 假设我们只使用第一个IMF,我们有100个样本 . 那么我们拥有的是:

instantaneous_frequency = [f1 f2 f3 ... f100]
    instantaneous_amplitude = [a1 a2 a3 ... a100]
    time = [t1 t2 t3 ... t100]

我需要的是一种方法来绘制它们像第二个图像,而不使用频谱图功能,因为它使用STFT,我已经应用了HHT .

制作第一个图的代码是:

Time_Window = 1:101;
signal= rand(1,101);
hilb = hilbert(signal);
inst_amp = abs(hilb);
inst_th = angle(hilb);
inst_freq = diff(a_inst_th)/(1/256)/(2*pi); %instantaneous frequency 
%inst_freq = remove_outliers(inst_freq,Time_Window(1:end-1)); 
inst_freq(end+1) = inst_freq(end); % This is done due to diff()
plot(Time_Window*4,inst_freq,'k.','MarkerSize',5)

1 回答

  • 0

    已经有一个MATLAB函数来生成spectrogram,它看起来很像你想要的 . 看看是否能解决您的问题 .

相关问题