hhmspectrum1.m

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

function psd = hhmspectrum1(A,f,DEFSPL)

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

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

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

psd=fliplr(psd);