Discussion :

[Mr_Quignon]

Bonjour à tous,

Je n'ai pas trouvé de solutions à mon problème (d'où mon premier message sur le forum):

Je voudrais savoir comment passer en argument d'une fonction 3 vecteurs de longueur respective m; n et p, pour que celle ci me renvoie une matrice de dimensions m x n x p?

C'est facile en dimension 1. Ca ne marche pas tout seul en dimension 2, il faut utiliser feval par ex:

Code :

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

1
2
3
4
5
6
function z=f(x,y)
  z=2*x^2+y^2;
endfunction
x=linspace(-1,1,100);
y=linspace(-2,2,200);
z=feval(x,y,f);

J'aimerai la même chose mais avec 3 arguments, du genre :

Code :

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

1
2
3
4
5
6
7
function z=f(x,y,w)
  z=2*x^2+y^2+w^2;
endfunction
x=linspace(-1,1,100);
y=linspace(-2,2,200);
W=linspace(-2,2,300);
z=f(x,y,w);   // Avec du coup z de dimension 100x200x300

J'arrive à quelque chose en parcourant les matrices avec boucles for, mais c'est crado et j'aimerais garder l'aspect calcul vectorisé du logiciel puisqu'il est optimisé pour ça.

Pouvez vous m'aider svp?
Merci pour votre temps!

[Jerome Briot]

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
function z=f1(x,y)
    z=2*x^2+y^2
endfunction
 
function z=f2(x,y)    
    [x,y] = ndgrid(x,y)
    z = 2*x.^2+y.^2
endfunction
 
x = linspace(-1,1,4)
y = linspace(-2,2,5)
 
z1 = feval(x,y,f1)
z2 = f2(x,y)
 
disp(and(z1==z2))

[Mr_Quignon]

Impeccable!

Ca marche aussi à 3 arguments, vu que ndgrid accepte autant d'argument qu'on veut je suppose que ça doit marcher pour n arguments.

Merci beaucoup j'aurai jamais trouvé tout seul...

[Jerome Briot]

ndgrid [...] doit marcher pour n arguments.

D'où le nom de la fonction : n-dimensional grid