1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83
| % expert 1
p_d1=0.001+(0.004-0.001)*rand(1000,1); %génération de 1000 valeurs de p_d1
p_d2=0.001+(0.004-0.001)*rand(1000,1); %génération de 1000 valeurs de p_d2
p_d3=0.0009+(0.002-0.0009)*rand(1000,1); %génération de 1000 valeurs de p_d3
p_d4=0.0009+(0.002-0.0009)*rand(1000,1); %génération de 1000 valeurs de p_d4
p_d5=0.009+(0.03-0.009)*rand(1000,1); %génération de 1000 valeurs de p_d5
p_d6=0.004+(0.007-0.004)*rand(1000,1); %génération de 1000 valeurs de p_d6
p_d7=0.004+(0.007-0.004)*rand(1000,1); %génération de 1000 valeurs de p_d7
p_d8=0.009+(0.03-0.009)*rand(1000,1); %génération de 1000 valeurs de p_d8
p_g1=p_d1.*p_d2; %probabilité echec gate 1
p_g5=p_d3.*p_d4; %probabilité echec gate 5
p_g3=1-(1-p_g5).*(1-p_d5); %probabilité echec gate 3
p_g6=p_d6.*p_d7; %probabilité echec gate 6
p_g4=1-(1-p_g6).*(1-p_d8); %probabilité echec gate 4
p_g2=p_g3.*p_g4; %probabilité echec gate 2
p_TE=p_g1.*p_g2; %probabilité echec TOP EVENT
figure(1);
hist(p_TE);
a=mean(p_TE)
% expert 2
p1_d1=0.001+(0.003-0.001)*rand(1000,1); %génération de 1000 valeurs de p1_d1
p2_d2=0.001+(0.003-0.001)*rand(1000,1); %génération de 1000 valeurs de p2_d2
p3_d3=0.003+(0.006-0.003)*rand(1000,1); %génération de 1000 valeurs de p3_d3
p4_d4=0.003+(0.006-0.003)*rand(1000,1); %génération de 1000 valeurs de p4_d4
p5_d5=0.008+(0.02-0.008)*rand(1000,1); %génération de 1000 valeurs de p5_d5
p6_d6=0.002+(0.006-0.002)*rand(1000,1); %génération de 1000 valeurs de p6_d6
p7_d7=0.002+(0.006-0.002)*rand(1000,1); %génération de 1000 valeurs de p7_d7
p8_d8=0.008+(0.02-0.008)*rand(1000,1); %génération de 1000 valeurs de p8_d8
p1_g1=p1_d1.*p2_d2; %probabilité echec gate 1
p5_g5=p3_d3.*p4_d4; %probabilité echec gate 5
p3_g3=1-(1-p5_g5).*(1-p5_d5); %probabilité echec gate 3
p6_g6=p6_d6.*p7_d7; %probabilité echec gate 6
p4_g4=1-(1-p6_g6).*(1-p8_d8); %probabilité echec gate 4
p2_g2=p3_g3.*p4_g4; %probabilité echec gate 2
p7_TE=p1_g1.*p2_g2; %probabilité echec TOP EVENT
y = dsadf('expinv',1000,p_TE);
figure(4)
dscdf(p_TE);
figure(2);
hist(p7_TE);
a=mean(p7_TE)
%moyenne des deux experts
p_d1=0.001+(0.004-0.001)*rand(1000,1); %génération de 1000 valeurs de p_d1
p_d2=0.001+(0.004-0.001)*rand(1000,1); %génération de 1000 valeurs de p_d2
p_d3=0.0009+(0.002-0.0009)*rand(1000,1); %génération de 1000 valeurs de p_d3
p_d4=0.0009+(0.002-0.0009)*rand(1000,1); %génération de 1000 valeurs de p_d4
p_d5=0.009+(0.03-0.009)*rand(1000,1); %génération de 1000 valeurs de p_d5
p_d6=0.004+(0.007-0.004)*rand(1000,1); %génération de 1000 valeurs de p_d6
p_d7=0.004+(0.007-0.004)*rand(1000,1); %génération de 1000 valeurs de p_d7
p_d8=0.009+(0.03-0.009)*rand(1000,1); %génération de 1000 valeurs de p_d8
p1_d1=0.001+(0.003-0.001)*rand(1000,1); %génération de 1000 valeurs de p1_d1
p2_d2=0.001+(0.003-0.001)*rand(1000,1); %génération de 1000 valeurs de p2_d2
p3_d3=0.003+(0.006-0.003)*rand(1000,1); %génération de 1000 valeurs de p3_d3
p4_d4=0.003+(0.006-0.003)*rand(1000,1); %génération de 1000 valeurs de p4_d4
p5_d5=0.008+(0.02-0.008)*rand(1000,1); %génération de 1000 valeurs de p5_d5
p6_d6=0.002+(0.006-0.002)*rand(1000,1); %génération de 1000 valeurs de p6_d6
p7_d7=0.002+(0.006-0.002)*rand(1000,1); %génération de 1000 valeurs de p7_d7
p8_d8=0.008+(0.02-0.008)*rand(1000,1); %génération de 1000 valeurs de p8_d8
p_g1=p_d1.*p_d2; %probabilité echec gate 1
p_g5=p_d3.*p_d4; %probabilité echec gate 5
p_g3=1-(1-p_g5).*(1-p_d5); %probabilité echec gate 3
p_g6=p_d6.*p_d7; %probabilité echec gate 6
p_g4=1-(1-p_g6).*(1-p_d8); %probabilité echec gate 4
p_g2=p_g3.*p_g4; %probabilité echec gate 2
p_TE=p_g1.*p_g2; %probabilité echec TOP EVENT
p1_g1=p1_d1.*p2_d2; %probabilité echec gate 1
p5_g5=p3_d3.*p4_d4; %probabilité echec gate 5
p3_g3=1-(1-p5_g5).*(1-p_d5); %probabilité echec gate 3
p6_g6=p6_d6.*p7_d7; %probabilité echec gate 6
p4_g4=1-(1-p6_g6).*(1-p_d8); %probabilité echec gate 4
p2_g2=p3_g3.*p4_g4; %probabilité echec gate 2
p7_TE=p1_g1.*p2_g2; %probabilité echec TOP EVENT
G1=(p_g1+p1_g1)/2%=p_d1.*p_d2; %probabilité echec gate 1 pour les deux experts
G2=(p_g5+p5_g5)/2
G3=(p_g3+p3_g3)/2
G4 = (p_g6+p6_g6)/2
G5 = (p_g4+p4_g4)/2
G6 = (p_g2+p2_g2)/2
G7 = (p_TE+p7_TE)/2
figure(3);
dscdf(p_TE)
hist(G7);
a=mean(G7) |
Partager