dist_value.m

Hilbert-Huang Transform的完整源码

% This is a utility program being called by "significance.m".
%
%   function PDF = dist_value(yPos, yBar, nDof)
%

function PDF = dist_value(yPos, yBar, nDof)

%   function PDF = dist_value2(yPos, yBar, nDof)
%
%   PDF: a normalized output array 
%   yPos: An input array at which PDF values are calculated
%	yBar: The expected value of yPos 
%	nDof: The number of degree of freedom
%
%
% References can be found in the "Reference" section.
%
% The code is prepared by Zhaohua Wu. For questions, please read the "Q&A" section or
% contact
%   [e-mail]
%
ylen = length(yPos);
eBar = exp(yBar);

evalue=exp(yPos);

for i=1:ylen,
    tmp1 = evalue(i)/eBar-yPos(i);
    tmp2 = -tmp1*nDof*eBar/2;
    tmp3(i) = 0.5*nDof*eBar*log(nDof) + tmp2;
end

rscale = max(tmp3);

tmp4 = tmp3 - rscale;

PDF= exp(tmp4);