查看单个帖子
旧 2008-10-17, 08:44   #1
whyhzj
初级会员
 
注册日期: 2008-09-03
帖子: 9
声望力: 0
whyhzj 正向着好的方向发展
问题 菜鸟求救,请大侠指点

这是我的用于可靠性威布尔参数估计的计算程序:
>> 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、这是不是说明F1=0和F2=0这个方程组无解,如果遇到这种情况怎么将数据处理才能得到结果?
2、我还遇到过结果是负值的情况,按威布尔分布的参数的定义是应大于0的,这种情况下又要将数据做何处理比较好呢?

小小菜鸟请求大侠们予以指点!不圣感激!!!
whyhzj 当前离线   回复时引用此帖