lingjian.m

来自「lingjian.m-----蒙特卡罗方法 lingjian.m使用零」· M 代码 · 共 37 行

chengben=[0 25 0;0 50 20;200 50 20;500 100 50;0 0 50;100 25 10;100 25 0];
rongcha=[0.05 0.1 0.1 0.1 0.1 0.1 0.05];
x10=0.1;x20=0.3;x30=0.1;x40=0.1;x50=1.5;x60=16;x70=0.75;
NC=0;NF=0;
for n=1:1000
    x1=normal(x10,(rongcha(1)*x10/3)^2,1);
    x2=normal(x20,(rongcha(2)*x20/3)^2,1);
    x3=normal(x30,(rongcha(3)*x30/3)^2,1);
    x4=normal(x40,(rongcha(4)*x40/3)^2,1);
    x5=normal(x50,(rongcha(5)*x50/3)^2,1);
    x6=normal(x60,(rongcha(6)*x60/3)^2,1);
    x7=normal(x70,(rongcha(7)*x70/3)^2,1);
    y=174.42*(x1/x5)*(x3/(x2-x1))^0.85*sqrt((1-2.62*(1-0.36*(x4/x2)^(-0.56))^1.5*(x4/x2)^1.16)/(x6*x7));
    if (y>=1.2&y<=1.4)|(y>1.6&y<=1.8)
        NC=NC+1;
    end
    if y<1.2|y>1.8
        NF=NF+1;
    end
end
Q=NC+NF*9;
C=0;
for i=1:7
    if rongcha(i)==0.01
        C=C+chengben(i,1);
    end
    if rongcha(i)==0.05
        C=C+chengben(i,2);
    end
    if rongcha(i)==0.1
        C=C+chengben(i,3);
    end
end
fprintf('Q=%f\n',Q);
fprintf('C=%f\n',C);
Z=Q+C