Discussion :

[mansgueg]

Bonjour,

Je trouve que matlab est très mal adapté pour travailler avec de très longues boucles, rendant la compilation de certain programme très longue (quelques heures pour certain aglo), auriez vous une astuce pour éviter de perdre se temps précieux ?

Merci

[duf42]

Bonsoir,

Il existe de nombreuses fonctions dans MATLAB qui permettent de se passer des boucles et d'optimiser le code, par contre il est difficile de définir un cas général, tu peux soit montrer un code pour qu'on puisse te dire ce qu'on en pense ou alors chercher sur le forum, de nombreux sujets traitent de ce genre de cas de figure.

Duf

[tubaas]

salut

la compilation de certain programme

Matlab ne compile pas mais est un interpréteur.
sinon pour répondre à ta question, en effet Matlab est un langage de programmation matricielle et n'est pas du tout optimisé pour les boucles.
il te faut donc bien construire tes matrices, utiliser les méthodes d'indexages (voir la FAQà ce sujet), ainsi que les fonctions FIND, REPMAT, RESHAPE qui sont souvent utiles.

edit: tu as d'ailleurs ce tutoriel (en lien depuis la FAQ) qui explique la gestion des matrices.

[kmaniche]

Bonsoir,

D'ordre général, pour optimiser son programme, passe par la vectorisation des boucles quand c'est possible.

Mais avant tout, il faut connaitre la répartition du temps d'exécution du programme par instruction. Pour cela, on a recours le plus souvent au profiler.

Code :

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

>> doc profiler

[mansgueg]

Mon code est tout simple pourtant

Code :

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

1
2
3
4
5
 
     for i=1:1:nombre_echantillons_5cm
          profil(i,1)=interp1(distance_clio,profil_clio,0.05*i,'nearest');
 
     end;

et cela prend 1 heure environs en sachant que nombre_echantillons_5cm peut valoir 500000 ...

[kmaniche]

Code :

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

1
2
3
4
 
for i=1:1:nombre_echantillons_5cm
               profil(i,1)=interp1(distance_clio,profil_clio,0.05*i,'nearest');
          end

qui peut être écrit comme :

Code :

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

1
2
3
 
i = 1:nombre_echantillons_5cm ;
profil(:,1)=interp1(distance_clio, profil_clio, 0.05*i,'nearest');

Et si ça prend énormément de temps c'est normal vu le nombre d'échantillons que tu as.

Question, Pourquoi ne pas utiliser un pas de 10 par exemple d'un échantillons à un autre. la différence ne se voit pas à mon avis.

Code :

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

1
2
3
i = 1:10:nombre_echantillons_5cm ;
profil(:,1)=interp1(distance_clio, profil_clio, 0.05*i,'nearest');

=> temps/10

[mansgueg]

Ce paramètre est très important et me permet d'être justement précis, je vais prendre mon mal en patience et laisser tourner le PC 24h/24, merci de vos réponses