Discussion :

[abel413]

Bonjour à tous,

J'ai le problème suivant : on se donne une matrice de la forme suivante :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
X=
1	3
10	12
17	24
29	35
41	46
51	58
60	65
73	79
82	84
91	99
102	104
112	116
118	123

Le but est de réécrire le résultat sous la forme d'une seule colonne dans un vecteur Y et de la façon suivante : sur le vecteur X on regarde chaque ligne, par exemple la première : 1 3 ; puis dans X, je dois écrire en colonne, par pas de 1, toutes les valeurs entre le début et la fin d'une ligne bornes comprises, ce qui donnerait ici : 1 2 3 écrit en colonne. De même pour 10 12, ce qui donnerait avec cet exemple le début suivant :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
Y=
 
1
2
3
10
11
12
17
18
19
20
...

J'ai réussit à le faire à l'aide de la fonction suivante :

Code :

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

1
2
3
4
5
6
7
8
function Y=creer(X1_Imp1_liquide_melange)
Y=[];
 
%donne les événements
 
    for i=1:size(X1_Imp1_liquide_melange,1)
        Y=[Y;(X1_Imp1_liquide_melange(i,1):X1_Imp1_liquide_melange(i,2))'];
    end

Mais je cherche un moyen de le réaliser directement sans boucle for.

EDIT:

Pour résumé, mon problème commence déjà par l’interrogation suivante : si on se donne l'intervalle [1 3], comment écrire tous les nombres de cet intervalle par pas de 1 (afin d'obtenir [1 2 3]) sans faire de boucle ?

[Jerome Briot]

Tu devrais pouvoir accélérer le code (dans le cas de matrice de grande dimension) car tu peux déterminer la taille finale du vecteur Y et donc faire une préallocation préalable :

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
X = [1	3
    10	12
    17	24
    29	35
    41	46
    51	58
    60	65
    73	79
    82	84
    91	99
    102	104
    112	116
    118	123];
 
d = diff(X,[],2);
nY = sum(d)+size(X,1);
Y = zeros(nY,1);
 
k = 1;
for i = 1:size(X,1)
    Y(k:k+d(i),1) = X(i,1):X(i,2);    
    k = k+d(i)+1;
end

A tester et à comparer à ta solution.

[abel413]

Merci beaucoup Dut pour cet effort. C'est effectivement plus rapide de la sorte malgré la boucle for !

[Jean Dumoncel]

Bonjour,

Une solution sans boucle :

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
clear
X = [1	3
10	12
17	24
29	35
41	46
51	58
60	65
73	79
82	84
91	99
102	104
112	116
118	123];
 
Y = cellfun(@(Y) Y(1):Y(2), num2cell(X,2),'UniformOutput', false).';
Y = cell2mat(Y).';

[abel413]

Merci à vous et désolé de répondre seulement maintenant.

magelan, j'aimerais bien tester ta variante mais comme je ne comprends pas tout à fait ce que tu as écrit, j'ai voulu le tester tel quel :

Code :

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

1
2
3
4
function Y=creer_sans_boucle(X1_Imp1_liquide_melange_evenements)
 
Y = cellfun(@(Y) Y(1):Y(2), num2cell(X1_Imp1_liquide_melange_evenements,2),'UniformOutput', false).';
Y = cell2mat(Y).';

Mais je reçois l'erreur suivante :

Code :

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

1
2
3
4
Too many inputs.
 
Error in ==> creer_sans_boucle at 3
Y = cellfun(@(Y) Y(1):Y(2), num2cell(X1_Imp1_liquide_melange_evenements,2),'UniformOutput', false).';

Apparement le cellfun de ma version matlab (7.0.1.24704 (R14) Service Pack 1) n'accepte pas autant de paramètres. Le problème est le même avec le code que tu m'as proposé tel quel ci-dessus.

EDIT:

En parcourant le forum, il est possible qu'il y ait un conflit de versions matlab et qu'avec la mienne, 7.0.1.24704 (R14) Service Pack 1 cela ne fonctionne pas. Y a-t-il malgré tout un moyen de faire marcher les deux lignes avec ma version ?

[abel413]

Quoi que je fasse, il y a toujours un problème : too many ou too few arguments et ce genre de chose. Au moins, magelan, le code que tu m'as proposé plus haut fonctionne-t-il tel quel sur ta version de matlab ?

[Jerome Briot]

Y a-t-il malgré tout un moyen de faire marcher les deux lignes avec ma version ?

Non.

[abel413]

Dut,

Dans ce cas existe-t-il quelque chose que je puisse faire du même genre que propose magelan (en prenant par exemple l'exemple du vecteur X donné dans le code de magelan) ? Le but est d'éviter toute boucle quand bien même ta variante est assez rapide.

[Jerome Briot]

Pourquoi vouloir à tout prix supprimer la boucle ?

[abel413]

J'ai tout simplement pour consigne, dans le cadre de mon travail, d'éviter toutes les boucles. Je pense que la variante de magelan est peut-être plus rapide que celle avec boucle. En fait, je pars de l'à priori, que dans le cas général un code sans boucles est plus rapide qu'un avec.

Je viens de tester le code avec un collègue et pour lui la ligne avec cellfun passe mais pas celle avec cell2mat tandis que chez moi celle avec cellfun ne fonctionne pas. C'est donc bien clairement la version de matlab qui est problématique. Lui a une édition 2012.

[Jerome Briot]

J'ai tout simplement pour consigne, dans le cadre de mon travail, d'éviter toutes les boucles. Je pense que la variante de magelan est peut-être plus rapide que celle avec boucle. En fait, je pars de l'à priori, que dans le cas général un code sans boucles est plus rapide qu'un avec.

Un simple test :

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
X = [1	3
    10	12
    17	24
    29	35
    41	46
    51	58
    60	65
    73	79
    82	84
    91	99
    102	104
    112	116
    118	123];
 
X = repmat(X,200,1);
 
tic
 
d = diff(X,[],2);
nY = sum(d)+size(X,1);
Y = zeros(nY,1);
 
k = 1;
for i = 1:size(X,1)
    Y(k:k+d(i),1) = X(i,1):X(i,2);
    k = k+d(i)+1;
end
 
toc
 
tic
Y2 = cellfun(@(Y) Y(1):Y(2), num2cell(X,2),'UniformOutput', false).';
Y2 = cell2mat(Y2).';
toc
 
all(Y(:)==Y2(:))

exécuté plusieurs fois sur ma machine :

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
Elapsed time is 0.005882 seconds.
Elapsed time is 0.029674 seconds.
 
ans =
 
     1
 
>> Untitled3
Elapsed time is 0.010957 seconds.
Elapsed time is 0.020074 seconds.
 
ans =
 
     1
 
>> Untitled3
Elapsed time is 0.005929 seconds.
Elapsed time is 0.021568 seconds.
 
ans =
 
     1
 
>> Untitled3
Elapsed time is 0.015660 seconds.
Elapsed time is 0.021886 seconds.
 
ans =
 
     1
 
>> Untitled3
Elapsed time is 0.007153 seconds.
Elapsed time is 0.023044 seconds.
 
ans =
 
     1
 
>> Untitled3
Elapsed time is 0.010017 seconds.
Elapsed time is 0.018972 seconds.
 
ans =
 
     1

La version avec la boucle est plus rapide sur ma machine.

[abel413]

Je te remercie Dut pour ce test, c'est déjà une bonne nouvelle dans ce cas. Je vais donc conserver la version avec boucle. De toute façon, je n'ai pas le choix pour l'instant et sa rapidité me satisfait. Le problème reste la fonction sort je crois mais je n'ai pas reçu de réponse dans le sujet dédié que tu avais créé pour.

[abel413]

Bonjour,

Je reviens à ce sujet car entre temps j'ai effectué d'autres recherches la solution suivante m'a été suggérée :

Code :

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

mcolon(X(:,1),X(:,2))

En utilisant les fonctions suivantes : http://www.mathworks.com/matlabcentr...exchange/29854

Le problème c'est que je ne peux pas tester ceci en raison du message d'erreur suivant que j'obtiens :

Code :

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

1
2
3
4
5
6
7
??? At compilation "mcolsettol" was determined to be variable. As a variable,
it is uninitialized, but function "mcolsettol" is now in context.
In MATLAB 7, you cannot use the same name in a single function as
both a function and a variable.
 
Error in ==> mcolon at 24
tol = mcolsettol();

A nouveau ma version de Matlab semble problématique ainsi que la fonction mcolsettol de ce packetage. Qu'en pensez-vous ? Au moins si quelqu'un pouvait simplement tester cette fonction si elle marche chez lui et faire le test de rapidité de Dut.