2006年《系统仿真学报》增刊征订启示
2006年《系统仿真学报》增刊征订启示 第l8卷第9期
2006年9月系统仿真Sept.,
2006
设,为口,的矩估计,采用以下方法获得参数的矩
估计:
口是方程
上:,?f11),11,
+F(1+2/a,
的解,其中--~和分别为样本均值和样本方差,r(.)为伽
马函数.最后取
=
~-/r(1+1/)(12)
c1,c2的取值依赖于样本大小.由初始值{Ci)及初始区
域,对D=[It,b】及其收缩区域D'=[a',h'】加上如下条
件
日"?0,f=1.2
4仿真算例
某轴承产品寿命服从威布尔分布,对其进行疲劳寿命试
验(单位:h),得40个失效数据为13_31,16.00,17.88, 28.92,32.18,33.0,41.52,42.12,45.60,48-48,51.84,
51.93,51.96,52_35,54.12,55.56,58.08,60.01,66.77,
67.80,68.64,68.64,68.88,79.82,84.12,87.88,91.48,
93.12,94.19,98.03,98.64,105.12,105.84,l13.03,l13.10,
127.92,128.04,159.38,173.40,201.46.试用SNTO算法
进行参数极大似然估计.
基于上述数据,应用SNTO算法求取威布尔分布参数的
极大似然估计.算法初始参数为:=5,C.=10, =0.000001,巩=987,n,:233.
图1描述了极大似然估计后的威布尔分布曲线. 表1给出了矩估计和利用SNTO求得的极大似然估计以 及文献[1]对分法得到的极大似然估计.LLMAX表示对数似 然函数,在相应参数估计的函数值.
表1三种估计比较
环,计算速度优于文献[1]的对分法.SNTO的这种计算效率 优势来源于数论网格取点方法,在保证精度的情形下有效缩 小了搜索范围.
槲
孽
较
时间(h)
图1原始数据和威布尔分布曲线
5结论
仿真算例表明这种方法既保证了计算精度,又具有较高 的算法效率.该算法优点包括:只
要求
对教师党员的评价套管和固井爆破片与爆破装置仓库管理基本要求三甲医院都需要复审吗
函数是连续的而不要 求函数是单峰,可微的,特别适合求解有多个局部极值的似 然函数问题;而且,只需计算似然函数在数论网格点的值, 而无需求其微商值,计算程序对似然函数的依赖较小,具有 较强的通用性.
参考文献-
[1]金星,洪延姬,等.可靠性数据计算及应用[M].北京:国防工业出 版社,2003.
[2]方开泰,王元.数论方法在统计中的应用[M].北京:科学出版社, 1996.
[3]华罗庚,王元.数论在近似
分析
定性数据统计分析pdf销售业绩分析模板建筑结构震害分析销售进度分析表京东商城竞争战略分析
中应用[M].北京:科学出版社, 1978.
『41GWCran.MomentEstimatorsfor山e3一ParameterWeibull 2006年<系统仿真》增刊征订启示
2006年《(系统仿真》出版2期增刊.增刊1出版第六届全国虚拟现实与可视化学术会议
(CCVRV'06)4~138篇论文.定价:260元.增刊2出版中国系统仿真学会第五次全国会员代表
大会暨2006年学术年会收录的294篇优秀论文.定价:400元.欢迎读者定购. 联系电话:88527147汇款地址:北京142信箱13分箱编辑部收邮编:100854 ?2536?
浙公网安备 33021202002483