这是我的用于可靠性威布尔参数估计的计算程序:是用matlab编的。
>>x=[134.6,139.7,142,139.1,137.8,133.7,150.1,159.4,138.3,145.2,163.6,138.1,138.5,146.9,145.4]
>>symsmn;
e=(m/n)*(x/n).^(m-1)
>>k=exp(-(x/n).^m)
>>F=log(prod(e.*k))
>>F1=diff(F,’m’)
>>F2=diff(F,’n’)
>>[m,n]=solve(‘F1=0′,’F2=0’)
**********
结果是:Warning:Explicitsolutioncouldnotbefound.
>Insolveat140
m=
[emptysym]
n=
[]
请问:1、这是不是说明F1=0和F2=0这个方程组无解,如果遇到这种情况怎么将数据处理才能得到结果?
2、我还遇到过结果是负值的情况,按威布尔分布的参数的定义是应大于0的,这种情况下又要将数据做何处理比较好呢?
小小菜鸟请求大侠们予以指点!不圣感激!!!:P
MATLAB自带Weibull参数估计的函数:weibfit()
还有,你的似然函数写得不对