Discussion :

[S4sha]

J'ai cette commande Matlab à optimiser :

Code :

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

1
2
 
A(:,find(A(1,:)<=0))=[];

où A est une matrice (3,n), et n est assez grand. Le but étant de retirer les colonnes de A dont le premier élément est nul. Cette instruction est répétée un très grand nombre de fois. Matlab me propose d'utiliser "logical" à la place, mais je ne vois pas trop comment. Que puis-je faire ?

Merci !

[rostomus]

Bonjour,

Code :

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

1
2
 
A(:,A(1,:)<=0)=[];

ou bien

Code :

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

1
2
 
A=A(:,A(1,:)>0);

[rostomus]

La deuxieme est plus rapide que la premiere:

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
 
clc
A=randn(3,10);
b=A;
tic
for i=1:10000
    A(:,A(1,:)<=0)=[];
    A=b;
end
toc
tic
for i=1:10000
    A=A(:,A(1,:)>0);
    A=b;
end
toc
clear

ca donne:

Code :

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

1
2
3
4
 
Elapsed time is 0.451000 seconds.
Elapsed time is 0.150000 seconds.
>>

[S4sha]

et pourquoi le A=b ?

[rostomus]

et pourquoi le A=b ?

Bonjour, le programme que je vous ai donné est seulement pour comparer entre les deux methodes en temps d'execution.

[tug83]

J'ai cette commande Matlab à optimiser :

Code :

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

1
2
 
A(:,find(A(1,:)<=0))=[];

où A est une matrice (3,n), et n est assez grand. Le but étant de retirer les colonnes de A dont le premier élément est nul. Cette instruction est répétée un très grand nombre de fois. Matlab me propose d'utiliser "logical" à la place, mais je ne vois pas trop comment. Que puis-je faire ?

Merci !

Tu trouveras en plus du code précédent une solution relative aux boucle for et la notion de logical:
http://www.mathworks.com/support/sol...ution=1-35LRTD

[creativo]

Bonjour,

Une syntaxe différente pour optimiser X2 le temps d'exécution :

Actuel :

Code :

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

1
2
3
4
A=2.458521*randn(3,6665535);
tic
A(:,find(A(1,:)<=0))=[];
toc

Proposition :

Code :

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

1
2
3
4
A=2.458521*randn(3,6665535);
tic
A=A(A(1,:)>0);
toc

Donne sur ma machine les temps d'éxécution suivants :

elapsed_time =

1.7030

elapsed_time =

0.8120

Ce gain suffit-il à résoudre votre problème ?

Bon courage
Créativo