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