Bonsoir à tous, je souhaite faire varier le t afin d'obtenir plusieurs valeurs estimées de fiabilité R_chap. Je tente de le faire avec la boucle sur t mais je n'y arrive pas .
Quelqu'un a une idée pour calculer R_chap pour plusieurs valeurs de temps ex t_f=0, 100, 500 ? Merci.
Mon code :
Code : Sélectionner tout - Visualiser dans une fenêtre à part
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 clear all; close all; clc; lambda = 1e-4; % t_f=0:1000-1; R_chap=[]; N=10^3; % nb simulations % Systeme hydraulique for i=1:N % boucle de simulation for t=0:100:1000 u=rand(11,1); % On gégère T1 & T2 va exponentielle %NIVEAU1 N1 T1(i)=-(1/lambda)*log(u(1)); % instant de défaillance composant 1; T2(i)=-(1/lambda)*log(u(2)); V1(i)=-(1/lambda)*log(u(3)); P0(i)=-(1/lambda)*log(u(4)); T3(i)=-(1/lambda)*log(u(5)); %NIVEAU2 N2 T4(i)=-(1/lambda)*log(u(6)); V3(i)=-(1/lambda)*log(u(7)); T7(i)=-(1/lambda)*log(u(8)); %NIVEAU3 T5(i)=-(1/lambda)*log(u(9)); V2(i)=-(1/lambda)*log(u(10)); T6(i)=-(1/lambda)*log(u(11)); N1(i)=min([T1(i) T2(i) V1(i) P0(i) T3(i)]); % instant de défaillance du sous-syst série N2(i)=min([T4(i) V3(i) T7(i)]); % instant de défaillance du sous-syst redondant N3(i)=min([T5(i) V2(i) T6(i)]); N23(i)=max([N2(i) N3(i)]); X(i)=min([N2(i) N23(i)]); %construction de la fct indicatrice if X(i)>= t else Ind(i)=0; end end Ind(i)=1; % Xj MTTF=sum(X)/N; % Fiabilité estimée R_chap(t)=sum(Ind)/N %histogram(X) end
Partager