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

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

我所做的是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)

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值