disp_hhms_new.m

hilbert-huang 变换分析一段脉搏波 运行maihht即可

function disp_hhms_new(psd,DEFSPL)

fs=100;

figure
interval=0.5*fs/DEFSPL;
xvalue(1)=interval;
yvalue=psd;

for i=2:200
    xvalue(i)=xvalue(i-1)+interval;
end
%换算成dB%%%%%%%%%%%%%%%%%%%%%%%%%
% for j=1:length(yvalue)
%     if yvalue(i)==0
%         yvalue(i)=0;
%     else
%         yvalue(i)=10*log10(yvalue(i));
%     end
% end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
plot(xvalue,yvalue);
xlabel(['frequency(Hz)'])
ylabel(['Amplitude'])
title('Hilbert-Huang marginal spectrum')