weibmen.m

来自「rang doppler imaging and motion compensa」· M 代码 · 共 28 行

function [bh,ah] = weibmen(x);

% Estimation of Weibull parameters using the
% method of Menon.

[n1 n2] = size(x);

n = n1*n2;

x = reshape(x,n,1);

m1 = (1/n)*sum(log(x));    % First moment of ln
m2 = (1/n)*sum(log(x).^2); % Second moment of ln

% Compute b parameter

bh = ( (6/pi^2)*(n/(n-1))*(m2 - m1^2) )^(-0.5);

% Compute a parameter

ah = (exp(bh*m1 + 0.577216))^(-1);


% Note, the MLE method is incorporated in the
% the Matlab statistic toolbox.
% The function name is  weibfit(x)