litir1.asv

来自「此文件包括蒙特卡罗方法的基本思想及应用」· ASV 代码 · 共 30 行

function Rguji=litiR1(t,thetaa1,thetaa2,thetab1,thetab2,mm)
%t 是要求统存寿
%thetaa1 是元件A1的数学期要
%thetaa2 是元件A2的数学期要
%thetaa1 是元件A1的数学期要
%thetab2 是元件B2的数学期要
%mm 是随机实验次数
thetaa1
proguji=frq/mm
frq=0;
randnum1 = exprnd(thetaa1,1,mm);
randnum2 = exprnd(thetaa2,1,mm);
randnum3 = exprnd(thetab1,1,mm);
randnum4 = exprnd(theta1,1,mm);
for ii=1:mm
    if (randnum1(1,ii)>t)&(randnum2(1,ii)>t)
        pass1=1;
    else
        pass1=0;
    end 
    if (randnum3(1,ii)>t)&(randnum4(1,ii)>t)
        pass2=1;
    else
        pass2=0;
    end 
    if (pass1+pass2)>=1
        frq=frq+1;
    end    
end    
Rguji=frq/mm