可靠性模拟.doc

设两系统都是由 4 个元件组成,每个元件的寿命服从参数为θ的指数分布，每个元件 是否正常工作相互独立.两系统的连接方式如下图所示，求两系统寿命大于T=100的概率. function Rguji=litiR1(t,thetaa1,thetaa2,thetab1,thetab2,mm) %t 是要求统存寿 %thetaa1 是元件A1的数学期要 %thetaa2 是元件A2的数学期要 %thetab1 是元件B1的数学期要 %thetab2 是元件B2的数学期要 %mm 是随机实验次数 frq=0; randnuma1 = exprnd(thetaa1,1,mm); randnuma2 = exprnd(thetaa2,1,mm); randnumb1 = exprnd(thetab1,1,mm); randnumb2 = exprnd(thetab2,1,mm); for ii=1:mm if (randnuma1(1,ii)t)(randnuma2(1,ii)t) pass1=1; else pass1=0; end if (randnumb1(1,ii)t)(randnumb2(1,ii)t) pass2=1; else pass2=0; end if (pass1+pass2)=1 frq=frq+1; end end Rguji=frq/mm function Rguji=litiR2(t,thetaa1,thetaa2,thetab1,thetab2,mm) %t 是要求统存寿 %thetaa1 是元件A1的数学期要 %thetaa2 是元件A2的数学期要 %thetab1 是元件B1的数学期要 %thetab2 是元件B2的数学期要 %mm 是随机实验次数 frq=0; randnuma1 = exprnd(thetaa1,1,mm); randnuma2 = exprnd(thetaa2,1,mm); randnumb1 = exprnd(thetab1,1,mm); randnumb2 = exprnd(thetab2,1,mm); for ii=1:mm if (randnuma1(1,ii)t)|(randnumb1(1,ii)t) pass1=1; else pass1=0; end if (randnuma2(1,ii)t)|(randnumb2(1,ii)t) pass2=1; else pass2=0; end if (pass1*pass2)==1 frq=frq+1; end end Rguji=frq/mm设两系统都是由 4 个元件组成,每个元件的寿命服从期望为100的指数分布，每个元件是否正常工作相互独立.两系统的连接方式如下图所示，求系统寿命大于T=30的概率(用随机投点法) 程序: frq=0;mm=10000; A1=exprnd(100,1,mm); A2=exprnd(100,1,mm); B1=exprnd(100,1,mm); B2=exprnd(100,1,mm); for i=1:mm if((A1(i)30A2(i)30)|(B1(i)30B2(i)30)) frq=frq+1; end end gailv=frq/mm 结果: gailv =7.985000000000000e-001 所以系统寿命大于T=30的概率为79.85%。