Transformer une String en fonction

Version imprimable

07/02/2008, 17h16
Jevil

Transformer une String en fonction

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:

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.
07/02/2008, 17h44
Jerome Briot
Tu peux peut être utiliser INLINE :
Code:

1 2 str='x+1'; f=inline(str,'x');
07/02/2008, 18h17
Jevil

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).
07/02/2008, 20h43
Jerome Briot
Peut être quelque chose comme ceci ?
Code:

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)
07/02/2008, 21h27
Jevil
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:

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.
07/02/2008, 21h33
Caro-Line
Ca dépend à quel niveau tu regardes le temps.
Petit test tout bête :
Code:

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 :

Citation:

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 ?

En essayant la même manipulation avec mon code :

Code:

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 :

Citation:

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

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

07/02/2008, 23h07
Jevil

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

Citation:

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 :

Citation:

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

Je ne comprends pas à quoi c'est dû.
07/02/2008, 23h10
Jerome Briot
Le problème semble bien provenir de EVAL (encore une fois)

En testant les deux codes suivant :
Code:

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:

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 :

Citation:

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

Et pour testB :

Citation:

??? 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" :koi:

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.

Essaie ceci :

Code:

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 :?

07/02/2008, 23h39
Jevil

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 ?
07/02/2008, 23h41
Jerome Briot

Citation:

Envoyé par Jevil

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

Essaie ceci :

Code:

clear subset
07/02/2008, 23h46
Jevil

Ca fonctionne ! :D

Merci beaucoup à toi et à Caro !

Bonne soirée !