LMD瞬时幅值频谱图求取？

白兔糖

求取PF1瞬时幅值（Am_d）的频谱图，到底该怎么做啊？有这么一段程序，该怎么改？要做成右图那样的。
[Am_d Fm_d]=lmd_test(x);


fs=120;N=2048;
y=fft(Am_d,N);
mag=abs(y(1:N/2))*2/N;
f=(0:N/2-1)'*fs/N;
plot(f,mag);

左图瞬时幅值，右图是它的频谱图。

左图是瞬时幅值，右图就是它的频谱图了。