Discussion :

[MBAYE BABACAR GUEYE]

Bonjour,
J'ai fait un programme matlab qui roule bien, mais à partir d'un certain nombre d'interation, l'exécution devient de plus en plus lente

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
for i=1:nb
    A=(find(mat==i));
    fprintf('A calcule taille=%i\n',length(A))
    if mod(i,nl)~=0
        numeroCase=nc*(mod(i,nl)-1)+floor(i/nl)+1; 
    else
        numeroCase=(nl-1)*nc + floor(i/nl);
    end
    if (~isempty(A)) 
        subplot(nl,nc,numeroCase);
        fprintf('subplot fait ')
        hist(matDat(A,33),minMois:maxMois)
        fprintf('hist effectue')
    end   
end

je constate la lenteur vient est due au subplot car je la constate avant l'écriture du message 'subplot fait '
NB: nc peut aller jusqu'à 40, nl: 70,nb=nc*nl, et size(matDat )=600 000 x 34, et size(mat)=600 000 x 1

[Jerome Briot]

Tu veux afficher 2800 (70*40) graphiques sur la même fenêtre...

Je doute que cela soit très efficace mais essaie toujours en créant les objets Axes en dehors de la boucle de calcul.

J'ai également optimisé légèrement le 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
k = 1;
for i = 1:nl
    for j = 1:nc
        ax(k) = subplot(nl,nc,k);
        k = k+1;
    end
 
end
 
mois  = minMois:maxMois;
 
for i=1:nb
    %     A=(find(mat==i));
    A = mat==i;
    %     fprintf('A calcule taille=%i\n',sum(A))
    if any(A)
 
        if mod(i,nl)~=0
            numeroCase=nc*(mod(i,nl)-1)+floor(i/nl)+1;
        else
            numeroCase=(nl-1)*nc + floor(i/nl);
        end
 
        %         subplot(nl,nc,numeroCase);
        %         fprintf('subplot fait ')
        hist(ax(numeroCase),matDat(A,33),mois);
        %         fprintf('hist effectue')
    end
end

Ce code est-il plus rapide ?

Sinon, utilise le Profiler de MATLAB pour identifier les lignes gourmandes en temps de calcul :

Code :

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

doc profile

[MBAYE BABACAR GUEYE]

merci DUT,
mon prog a fini de tourner hier (aprés 2 jours,je l'ai moi même arreté), j'ai appliqué les modif que t'as proposé, ça va moins lentement mais j'ai toujours une vitesse d'croissante, i.e dans les premiéres itérations (jusqu'à 200 itérations) ca roule rapidement et ensuite trés lentement (vers les 500iéme itération ).

[Jerome Briot]

Comme je l'ai indiqué dans mon précédent message, tu comptes afficher 2800 objets Axes sur une même fenêtre.

As-tu au moins réfléchi à la taille de chacun de ces objets à l'affichage final ?

[MBAYE BABACAR GUEYE]

La lenteur est due au subplot. J'ai opté pour un affichage sur 4 fenetre mais en gardant la mme topologie.Quelqu'un aurait-il une idée sur cela

[MBAYE BABACAR GUEYE]

Salut,
Finalement j'ai trouvé la solution à mon probleme en utilisant 4 fénétres et en respectant la topologie d'origine