菜鸟求救,请大侠指点
 

这是我的用于可靠性威布尔参数估计的计算程序:
>> 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]

>> syms m n;
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: Explicit solution could not be found.
> In solve at 140

m =
[ empty sym ]
n =
[]
请问:1、这是不是说明Ｆ1=0和Ｆ２＝０这个方程组无解，如果遇到这种情况怎么将数据处理才能得到结果？
２、我还遇到过结果是负值的情况，按威布尔分布的参数的定义是应大于０的，这种情况下又要将数据做何处理比较好呢？

小小菜鸟请求大侠们予以指点！不圣感激！！！

回复: 菜鸟求救,请大侠指点
 

不就是求最大似然估计吗，以下代码即可

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];
[mn,mnci]=mle(x,'distribution','weibull')

结果（mn为参数估计结果，mnci的两列分别为相应的置信区间）：

mn =

147.7382 15.3725


mnci =

142.6437 10.7835
153.0146 21.9144

回复: 菜鸟求救,请大侠指点
 

好的,不胜感激!有问题再向您请教!

回复: 菜鸟求救,请大侠指点
 

这个weibull分布拟合我做过！

回复: 菜鸟求救,请大侠指点
 

[QUOTE=xiezhh;18029]不就是求最大似然估计吗，以下代码即可

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];
[mn,mnci]=mle(x,'distribution','we...[/QUOTE]

高手哈，向你学习。