在下面的代码加几行程序，实现生成频谱图load ('ecg'); Nstd = 0.2; NR = 500; MaxIter = 5000; [modes its]=ceemdan(ecg,0.2,500,5000); t=1:length(ecg); [a b]=size(modes); figure; subplot(a+1,1,1); plot(t,ecg);% the ECG signal is in the first row of the subplot ylabel('ECG') set(gca,'xtick',[]) axis tight; for i=2:a subplot(a+1,1,i); plot(t,modes(i-1,:)); ylabel (['IMF ' num2str(i-1)]); set(gca,'xtick',[]) xlim([1 length(ecg)]) end; subplot(a+1,1,a+1) plot(t,modes(a,:)) ylabel(['IMF ' num2str(a)]) xlim([1 length(ecg)])

可以加入以下代码来实现生成频谱图：

for i=1:a
    figure;
    [f, P] = periodogram(modes(i,:), [], [], 1000); %计算信号的功率谱密度
    plot(f, P);
    title(['IMF ' num2str(i) ' Spectrum']);
    xlabel('Frequency (Hz)');
    ylabel('Power/Frequency (dB/Hz)');
end

这段代码会计算每个IMF信号的功率谱密度，并绘制出频谱图。可以通过观察频谱图来分析信号的频域特征。