Fonction à deux variables

Bonjour à tous,

J'essaie de tracer une fonction à deux variables f(x,y) sur MATLAB et je n'y arrive pas... J'ai une fonction dont l'expression change en fonction du signe de x, la voici (l et ang sont des paramètres que je ne fais pas varier) :

Code:

---

1 2 3 4 5 6 7 8 9 10 11

function [d] = beam_length(x,y,l,ang) if x>0 d=sqrt(l^2+(x^2+y^2)+2*l*sqrt(x^2+y^2)*sin(acos(x./sqrt(x^2+y^2))*(1-2*ang/pi))); end if x==0 d=sqrt(l^2+(x^2+y^2)+2*l*sqrt(x^2+y^2)*sin(pi/2*(1-2*ang/pi))); end if x<0 d=sqrt(l^2+(x^2+y^2)+2*l*sqrt(x^2+y^2)*sin((pi-acos(-x./sqrt(x^2+y^2)))*(1-2*ang/pi))); end end

---

Je souhaite donc la représentation graphique de d=beam_length(x,y) en 3D. Pour cela je crée ma grille avec meshgrid et ensuite j'applique ma fonction:

Code:

---

1 2	[x,y]=meshgrid(-0.5:0.01:0.5); r=beam_length(x,y,1.414,pi/4)

---

Error in ==> beam_length at 2
if x>0

??? Output argument "d" (and maybe others) not assigned during call to "beam_length.m (beam_length)".

Voila donc le problème, je pense qu'il fait le test de signe sur le vecteur x au lieu de le faire sur l'élément mais je n'arrive pas à changer cela.

Merci beaucoup pour votre aide

Bonjour,

Merci à l'avenir de mettre le message d'erreur en entier.

Voir :faq: Output argument "***" (and maybe others) not assigned during call to "***"
Ici x est une matrice, et x>0, x==0 et x<0 sont vrais si et seulement si TOUS les éléments sont vrais, ainsi aucune des conditions if n'est vraie. Il en résulte que d n'est jamais fixé, d'où l'erreur.

Merci Jérôme pour cette explication.

Dans ce cas comment je fais pour contourner cette spécificité ? As-tu une piste à me proposer ?

Merci

Bonjour,

Les fonctions ALL et/ou ANY pourraient peut-être t'être utile ici.

Duf

1/ si tu veux que ta fonction length_beam travaille correctement avec des arguments tableaux (ce qui est une très bonne pratique de programmation Matlab :ccool:), il faut que tu utilises les opérations membre à membre (.* ./ .^) en remplacement des opérations matricielles (* / ^).

2/ pour tes test, quand tu teste x, tu teste une matrice et tu veux récupérer de quoi 'arroser' les bons x pour le calcul. On fai ça comme ça :

Code:

---

1 2 3 4

xplus = (x>0) ; d = zeros(taille finale) ; d(xplus) = mon_calcul_pour_les_x>0 ; d(~xplus) = mon_calcul_pour_les_x<=0 ;

---

Bonjour Olivier et merci de ta réponse.

Je ne comprends pas bien ton code du coup ça ne fonctionne pas tout le temps qu'en j'essaie sur Matlab. Voila la 1ère fonction que j'ai écris pour le tester:

Code:

---

1 2 3 4 5 6 7 8 9

function [ d ] = test(x,y) xplus = (x>0) ; d = zeros(101,101) ; d(xplus) = 7 ; d(~xplus) = 4; end   >> [x,y]=meshgrid(-0.5:0.01:0.5); >> r=test(x,y)

---

La ça fonctionne sans problème, par contre si je fais intervenir x ou y voici la réponse:

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12

function [ d ] = test(x,y) xplus = (x>0) ; d = zeros(101,101) ; d(xplus) = (x)+7 ; d(~xplus) = 4; end   ??? In an assignment A(I) = B, the number of elements in B and  I must be the same.   Error in ==> test at 4 d(xplus) = x+7 ;

---

J'arrive donc au même problème comment lui faire comprendre que je veux qu'il prenne les éléments de x et pas la matrice ! De même si j'utilise .*

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13

function [ d ] = test(x,y) xplus = (x>0) ; d = zeros(101,101) ; d(xplus) = x.*y ; d(~xplus) = 4; end   >> r=test(x,y) ??? In an assignment A(I) = B, the number of elements in B and  I must be the same.   Error in ==> test at 4 d(xplus) = x.*y ;

---

Merci encore pour votre aide

Sauf dans le cas d'un scalaire à assigner, le nombre d'éléments à gauche et à droite se doit d'être le même.

Code:

---

1 2	d(xplus) = 7 ; d(~xplus) = 4;

---

Ici 7 et 4 étant des scalaires, pas de problème : tous les éléments sont assignés à la valeur du scalaire.

Code:

---

1 2	d(xplus) = (x)+7 ; d(xplus) = x.*y ;

---

Là, d(plus) correspond aux éléments positifs stricts de x tandis que (x)+7 et x.*y contiennent autant d'éléments que x.
D'où l'erreur In an assignment A(I) = B, the number of elements in B and I must be the same..
Ainsi donc :

Code:

---

1 2	d(xplus) = x(xplus)+7 ; d(xplus) = x(xplus).*y(xplus) ;

---

Merci beaucoup à tous ma fonction fonctionne désormais parfaitement.

Bonne journée