Discussion :

[DoubleZ]

Bonjour

Je débute en programmation MATLAB et ai besoin d'un petit de pouce sur mon problème.

Je cherche à générer un arbre. A partir d'un critère (condition if) à chaque pas de temps, un nœud peut en générer un certain nombre d'autres (dans mon exemple, ce nombre est nb_quantiles que j'ai fixé à 4) si le critère est rempli et on poursuit le calcul sur ces nouveaux nœuds. Dans le cas contraire, le nombre de nœuds reste identique et le calcul se poursuit.

J'ai un énorme souci de temps de calcul. J'utilise "beaucoup" de boucles for et if imbriquées et je pense que cela vient de là. Notamment, le nombre de noeuds générés peut vite "exploser" (4 puissance le nombre de fois où le critère est rempli). J'ai essayé de réduire les tailles de vecteurs (avec reshape) sans succès. Je n'ai pas trouvé d'autres moyens de réaliser cet arbre sur le forum. Auriez-vous une idée à me proposer?



Voici mon bout de code (une belle usine à gaz)

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
 
%****  M et V grandeurs aléatoires  ********
 
M=(.001+rand(1,2501,1));            % Masses aléatoires (borne inf;centre;borne sup)
V=(1+2*rand(1,2501,1));             % Vitesses aléatoires (borne inf;centre;borne sup)
Ecinet=transpose(0.5*M.*(V.*V));    % Energies cinétiques aléatoires
nb_ab=1;                            % Initialisation compteur d'abrasion
 
for i=2:size(Ecinet,1)
    if Ecinet(i)>4.1                % Condition d'abrasion
        nb_ab=nb_ab+1;              % Mise à jour numéro abrasion
    end
end
 
Cab_proba_quantiles=zeros((nb_quantiles-1)^(nb_ab),1); % On affecte au vecteur Cab_proba_quantiles la taille maximum due à l'abrasion
Clast_ab=zeros((nb_quantiles-1)^(nb_ab),1);            % On affecte au vecteur Clast_ab qui stocke les valeurs à la dernière abrasion connue
limite=(nb_quantiles-1)^nb_ab;
 
while limite>10000
    nb_ab=nb_ab-1;
    limite=(nb_quantiles-1)^(nb_ab);
end
 
tab=1;      % Réinitialisation instant abrasion
nb_ab=1;    % Réinitialisation numéro abrasion
 
figure(17)
 
for t=2:size(Ecinet,1)
    if Ecinet(t)>4.1  % Condition d'abrasion
       Clast_ab=Cab_proba_quantiles;
       for i=1:(nb_quantiles-1)^(nb_ab)
           Cab_proba_quantiles((1:(nb_quantiles-1))+(i-1)*(nb_quantiles-1),1)=Clast_ab(i,1)+Cproba_quantiles(1:nb_quantiles-1,t-tab+2);
           % Changement de variable et multiplication de cas dus à l'abrasion
           if size(Cab_proba_quantiles,1)>10000
               Cab_proba_quantiles=reshape(Cab_proba_quantiles,limite,[]);
           end
       end
       Clast_ab=Cab_proba_quantiles; % Stockage des valeurs au moment de dernière abrasion
       nb_ab=nb_ab+1; % Actualisation nombre d'abrasion
       tab=t;         % Actualisation instant d'abrasion
    else
       for i=1:(nb_quantiles-1)^(nb_ab-1)
           Cab_proba_quantiles((1:(nb_quantiles-1))+(i-1)*(nb_quantiles-1),1)=Clast_ab((1:(nb_quantiles-1))+(i-1)*(nb_quantiles-1),1)+Cproba_quantiles(1:nb_quantiles-1,t-tab+1);
           % Changement de variable 
       end 
    end
 
    rep=find(Cab_proba_quantiles==0); % Recherche des éléments nuls du vecteur
    Cab_proba_quantiles(rep)=[];      % Remplacement des éléments nuls du vecteur par rien (=diminution dimension vecteur)
 
    hold on            
   % subplot(3,3,e)  % Affichage des différents scénarios
    plot(t,Cab_proba_quantiles)        
end

Merci d'avance pour vos contributions

[Jerome Briot]

Tu peux déjà éviter d'utiliser find dans le bloc de code suivant :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
rep = Cab_proba_quantiles==0; % Recherche des éléments nuls du vecteur
Cab_proba_quantiles(rep) = [];      % Remplacement des éléments nuls du vecteur par rien (=diminution dimension vecteur)

Voir la FAQ : Qu'est-ce que l'indexage logique ?

Pour le reste, commence par analyser ton code à l'aide du profiler : Profile to Improve Performance