Download probabrs.procedur
Back to the list
1 : * PROBABRS PROCEDUR AM 99/07/13 21:15:06 3631
2 : 'DEBPROC' PROBABRS MuR*'FLOTTANT' SigR*'FLOTTANT'
3 : rR*'FLOTTANT' bR*'FLOTTANT'
4 : MuS*'FLOTTANT' SigS*'FLOTTANT'
5 : rS*'FLOTTANT' bS*'FLOTTANT';
6 : *|=====================================================================
7 : *|
8 : *| << OBJET >> :
9 : *|
10 : *| Procedure permettant le calcul de la probabilite de defaillance
11 : *| idealisee :
12 : *| Prob(R> :
21 : *| -- ARGUMENTS OBLIGATOIRES DANS TOUS LES CAS --
22 : *|
23 : *| MuR, MuS : objets de type REEL designant la valeur moyenne
24 : *| de R et S.
25 : *|
26 : *| SigR, SigS : objets de type REEL designant l'ecart-type
27 : *| de R et S.
28 : *|
29 : *| rR, rS : objets de type REEL designant la valeur du coefficient
30 : *| de symétrie de R et S.
31 : *|
32 : *| bR, bS : objets de type REEL designant la valeur du coefficient
33 : *| d'aplatissement de R et S.
34 : *|
35 : *| << SORTIE >> :
36 : *|
37 : *| Pf : objet de type REEL designant la probabilite de défaillance
38 : *|
39 : *|=====================================================================
40 : *
41 : * ______________________
42 : *
43 : * DECLARATIONS GENERALES
44 : * ______________________
45 : *
46 : SI = 'MOT' 'SI' ; SINON = 'MOT' 'SINON' ; FINSI = 'MOT' 'FINSI' ;
47 : REPETER = 'MOT' 'REPETER'; QUITTER = 'MOT' 'QUITTER'; FIN = 'MOT' 'FIN';
48 : EXISTE = 'MOT' 'EXISTE' ; EGA = 'MOT' 'EGA' ;
49 : ERREUR = 'MOT' 'ERREUR' ;
50 : ENTIER = 'MOT' 'ENTIER ' ;
51 : FLOTTANT = 'MOT' 'FLOTTANT' ;
52 : MOT = 'MOT' 'MOT ' ;
53 : * _____
54 : *
55 : * TESTS
56 : * _____
57 : *
58 : * _________________
59 : *
60 : * DEBUT DES CALCULS
61 : * _________________
62 : *
63 : *
64 : Pf =0.;
65 : oldS = 0.;
66 : FR = 0.;
67 : NPAS = 600;
68 : tabmin = PROG (MuS-(6*SigS)) (MuR-(6*SigR));
69 : tabmax = PROG (MuS+(6*SigS)) (MuR+(6*SigR));
70 : binf = mini tabmin;
71 : bsup = maxi tabmax;
72 : h = (bsup-binf)/NPAS;
73 : REPE BOUC NPAS;
74 : absx = binf + (&BOUC*h);
75 : valR = PROB MuR SigR rR bR absx;
76 : valS = PROB MuS SigS rS bS absx;
77 : fS = (valS-oldS)/h;
78 : FR = valR;
79 : Pf = Pf + (h*fS*FR);
80 : oldS = valS;
81 : FIN BOUC;
82 : 'FINPROC' Pf;
83 :