litir1.m

此文件包括蒙特卡罗方法的基本思想及应用

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