ggpdf.m

来自「基于小波变换的特征检索算法」· M 代码 · 共 15 行

function y = ggpdf(x, a, b)
% GGPDF	Generalized Gaussian probability density function
%	Y = GGPDF(X, A, B) returns the generalized Gaussian density function
%	with parameters A and B, at the values in X.
%
%	Y = (B/(2*A*gamma(1/B))) * exp(-(X/A)^B)

%   Return NaN if the arguments are outside their respective limits.
if (a <= 0 | b <= 0)
    tmp = NaN;
    y = tmp(ones(size(x)));    
else
    y = log(b) - log(2) - log(a) - gammaln(1/b) - (abs(x).^b) ./ (a^b);
    y = exp(y);
end