[Index for rtools]
[Return to Master Index]
[Return to www.radarworks.com]
RT_weibest
(rtools/RT_weibest.m)
Function Synopsis
[a,mu] = RT_weibest(x,fast)
Help text
function [a,mu] = RT_weibest(x, fast)
Estimates the shape parameters A and mean MU of a Weibull distribution
that best fits the data in X.
Note that this is not the `usual' form for the Weibull, which uses a
non-linear scale parameter (b). For simulation, explicit definition
in terms of the mean MU is more meaningful. For theoretical
manipulation of equations, you're better off using b though.
The conversion from b is MU = a^(-b) * gamma(1 + b) which is performed in
RT_weibest.m
A Weibull distribution in Amplitude, is also a Weibull distribution
in Intensity. This `power law invariance' is useful for simulation, but
it causes confusion over whether the `usual' shape parameter (c) refers
to Amplitude _or_ Intensity in radar texts. It most often refers to Amplitude,
but this toolbox uses Intensity as the natural measure, so check your
conventions.
The form used is such that the mean value is MU and larger
A is closer to a Gaussian distribution. A = 1
is equivalent to an exponential distribution.
The method uses an initial least squares fit to Weibull paper.
It is not the Maximum Likelihood Solution, but A is reasonably
accurate (MU is quite poor). To return this estimate, set the
FAST flag.
With FAST unset or 0, a numerical MLE solution is generated (slow)
but this can fail when a is too small (a ~< 0.1).
Source: Direct evaluation of linear approximation, improvement
via a call to RT_weibLL.m
Part of the Matlab Radar Toolbox v0.11 by Glen Davidson at www.radarworks.com
Please acknowledge use of this toolbox in any publication or software.
This header must be retained in any derivative work
Cross-Reference Information
This function calls
This function is called by
Produced by mat2html on Mon Mar 17 16:14:19 2003
Cross-Directory links are: OFF