Transformer une String en fonction

**Jevil** · 07/02/2008, 17h16

Bonjour,

Je souhaite transformer une String du genre 'x+1' en fonction.

Pour être plus clair, voici le contexte :
J'ai une interface dans laquelle l'utilisateur saisit une formule ('x+1') que je récupère donc dans une variable formule de type String. Je souhaite ensuite appliquer cette formule à une variable y en fonction de x. J'ai d'abord penser utiliser la fonction eval de Matlab :

Code :

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

 y = eval(formule)

mais étant donné que j'appelle cette fonction dans une boucle et que la formule peut être complexe, le temps de traitement est très long (à cause de eval()).

Pour résumer différemment, est-il possible de créer une variable f de type function_handle et de valeur @(x)x+1 à partir de la String 'x+1' ?

En espérant avoir été assez clair, merci pour vos réponses.

**Jerome Briot** · 07/02/2008, 17h44

Tu peux peut être utiliser INLINE :

Code :

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

1
2
str='x+1';
f=inline(str,'x');

**Jevil** · 07/02/2008, 18h17

Merci je ne connaissais pas,
mais le temps de traitement est aussi long qu'avec la fonction eval (et beaucoup plus long que lorsque l'on rentre la formule en dur).

**Jerome Briot** · 07/02/2008, 20h43

Peut être quelque chose comme ceci ?

Code :

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

1
2
3
4
5
6
7
8
9
10
 
f='x+1';
 
g=eval(['@(x) ' f]);
 
for n=1:3
 
   g(n)
 
end

(merci Caro pour l'idée)

**Jevil** · 07/02/2008, 21h27

Ca fontionne et c'est plus rapide que le inline.
En revanche ça reste toujours un peu long et dans ce cas je ne comprends pas pourquoi le traitement est plus long que si on rentrait directement :

Code :

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

1
2
3
4
5
6
7
g=@(x) x+1;
 
for n=1:3
 
   g(n)
 
end

puisque des deux manières, on a bien @(x) x+1 dans g.

**Caro-Line** · 07/02/2008, 21h33

Ca dépend à quel niveau tu regardes le temps.
Petit test tout bête :

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
clear all
g=@(x) x+1;
tic
for n=1:30
   g(n);
end
toc
 
clear all
f='x+1';
g=eval(['@(x) ' f]);
tic
for n=1:30
   g(n);
end
toc

Moi j'obtiens :

Elapsed time is 0.002189 seconds.
Elapsed time is 0.002059 seconds.

(en le faisant plusieurs fois bien sûr)
Donc à peu près équivalent.

Tu n'as pas d'autres choses dans ta boucle qui pourraient être optimisées ?

**Jevil** · 07/02/2008, 21h58

En essayant la même manipulation avec mon 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
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
clear all
sizeX = 400;
sizeY = 400;
im = zeros(sizeX,sizeY);
im = im2double(im);
 
f='z+1';
formule=eval(['@(z) ' f]);
tic
for k=-sizeX/2+1:sizeX/2
   for l=-sizeY/2+1:sizeY/2  
       a = k/(sizeX/4);
       b = l/(sizeY/4);
       r = sqrt(a*a+b*b);
       z = a + i*b;
       n = 0;  
       while r < 10 && n <= 10
          z = formule(z);
          a = real(z);
          b = imag(z);
          r = sqrt(a*a + b*b);
          n = n+1; 
       end
       im(k+sizeX/2,l+sizeY/2) = min( abs(a) , abs(b));
   end
end
toc
 
clear all
sizeX = 400;
sizeY = 400;
im = zeros(sizeX,sizeY);
im = im2double(im);
 
formule=@(z)z+1;
tic
for k=-sizeX/2+1:sizeX/2
   for l=-sizeY/2+1:sizeY/2  
       a = k/(sizeX/4);
       b = l/(sizeY/4);
       r = sqrt(a*a+b*b);
       z = a + i*b;
       n = 0;  
       while r < 10 && n <= 10
          z = formule(z);
          a = real(z);
          b = imag(z);
          r = sqrt(a*a + b*b);
          n = n+1; 
       end
       im(k+sizeX/2,l+sizeY/2) = min( abs(a) , abs(b));
   end
end
toc

J'obtiens :

Elapsed time is 43.046651 seconds.
Elapsed time is 1.463195 seconds.

Ca donne une grosse différence sur une boucle plus importante.

**Jevil** · 07/02/2008, 23h07

C'est bizarre, quand je place ce morceau de code dans une fonction et que je l'appelle, j'obtiens :

Elapsed time is 42.617504 seconds.
Elapsed time is 1.298563 seconds.

Et quand j'en fait un script (ou que je l'appelle directement dans la Command Window), j'obtiens :

Elapsed time is 25.860350 seconds.
Elapsed time is 25.854857 seconds.

Je ne comprends pas à quoi c'est dû.

**Jerome Briot** · 07/02/2008, 23h10

Le problème semble bien provenir de EVAL (encore une fois)

En testant les deux codes suivant :

Code :

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

1
2
3
4
5
6
7
8
function testA
 
formule=@(x)x+1;
 
y=0;
for n=1:1E15
    y=y+formule(n);
end

Code :

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

1
2
3
4
5
6
7
8
function testB
 
formule=eval('@(x)x+1;');
 
y=0;
for n=1:1E6
    y=y+formule(n);
end

Et en les interrompant pendant leur exécution avec un CTRL+C, on obtient :

Pour testA :

??? Operation terminated by user during ==> testA at 7

Et pour testB :

??? Operation terminated by user during ==> testB>@(x)x+1

In ==> testB at 7
y=y+formule(n);

Il semble qu'il y ait comme une sorte de récursivité autour de EVAL, même dans ce cas particulier. La fonction testB serait donc appelée plusieurs fois afin d'évaluer "formule"

C'est d'ailleurs le fonctionnement normal de EVAL. Son argument pouvant être modifié entre deux appels, MATLAB ne peut le précompilé. D'où un temps d'exécution plus important.

**Jerome Briot** · 07/02/2008, 23h21

Essaie ceci :

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
56
57
58
59
function test

sizeX = 400;
sizeY = 400;
im = zeros(sizeX,sizeY);
im = im2double(im);

formule='z+1';
fid=fopen('subtest.m','wt');
    fprintf(fid,'function z=subtest(z)\nz=%s;',formule);
fclose(fid);

tic
for k=-sizeX/2+1:sizeX/2
   for l=-sizeY/2+1:sizeY/2  
       a = k/(sizeX/4);
       b = l/(sizeY/4);
       r = sqrt(a*a+b*b);
       z = a + i*b;
       n = 0;  
       while r < 10 && n <= 10
          z = subtest(z);
          a = real(z);
          b = imag(z);
          r = sqrt(a*a + b*b);
          n = n+1; 
       end
       im(k+sizeX/2,l+sizeY/2) = min( abs(a) , abs(b));
   end
end
toc

clear all

sizeX = 400;
sizeY = 400;
im = zeros(sizeX,sizeY);
im = im2double(im);

formule=@(z)z+1;
tic
for k=-sizeX/2+1:sizeX/2
   for l=-sizeY/2+1:sizeY/2  
       a = k/(sizeX/4);
       b = l/(sizeY/4);
       r = sqrt(a*a+b*b);
       z = a + i*b;
       n = 0;  
       while r < 10 && n <= 10
          z = formule(z);
          a = real(z);
          b = imag(z);
          r = sqrt(a*a + b*b);
          n = n+1; 
       end
       im(k+sizeX/2,l+sizeY/2) = min( abs(a) , abs(b));
   end
end
toc

ça commence à être tiré par les cheveux

**Jevil** · 07/02/2008, 23h39

Oui j'avais déjà essayer cette astuce, elle fonctionne mais le problème est que j'aimerais pouvoir modifier ma formule pendant l'exécution du programme. Or une fois que le fichier subset.m, et donc la fonction subset, sont créés, le fait de remodifier le fichier ne modifie pas la fonction : la première fonction subset en fonction de la première formule saisie est conservée pendant toute la durée de l'exécution.

Y aurait-il moyen de réinitialiser le programme pour que la fonction subset soit "effacée" lorsque l'on saisit une nouvelle formule ?