Question

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

我所做的是ceil()，但我需要它看起来像this。

提前谢谢

编辑：

我有使用Hilbert Huang变换生成的振幅和瞬时频率随时间变化的值。我们假设我们只使用第一个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)

Answer 1

已经有一个MATLAB函数可以产生spectrogram，它看起来很像你想要的。看看是否能解决您的问题。

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

1 个答案: