找回密码
 -注 册-
搜索
热搜: MTBF GJB MIL FMEA
查看: 2731|回复: 1

菜鸟求教,请高手指点,帮我看看程序

[复制链接]
发表于 2008-10-17 08:56:28 | 显示全部楼层 |阅读模式
这是我的用于可靠性威布尔参数估计的计算程序:是用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的,这种情况下又要将数据做何处理比较好呢?

小小菜鸟请求大侠们予以指点!不圣感激!!!
发表于 2009-4-24 21:24:30 | 显示全部楼层
MATLAB自带Weibull参数估计的函数:weibfit()
还有,你的似然函数写得不对
回复

使用道具 举报

您需要登录后才可以回帖 登录 | -注 册-

本版积分规则

QQ|Archiver|手机版|小黑屋|可靠性网 ( 粤ICP备14066057号 )

GMT+8, 2025-4-20 21:45

Powered by Discuz! X3.5

© 2001-2023 Discuz! Team.

快速回复 返回顶部 返回列表