hhmspectrum_new.m

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

function psd = hhmspectrum_new(A,f,DEFSPL)

  t = 1:size(A,2);
  sply = DEFSPL;

lt=length(f(1,:));

psd=[];
psd(sply) = 0;

for i=1:size(f(:,1))
  for j = 1:lt
    ff=floor(f(i,j)*2*(sply-1))+1;
    if ff <= sply % in case f(i,j) > 0.5
      psd(ff)=psd(ff)+A(i,j);
    end
  end  
end

psd=psd';