Discussion :

[antoineleb]

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 :

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

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 :

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

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

[Invité]

Bonjour,

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

Voir 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.

[antoineleb]

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

[duf42]

Bonjour,

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

Duf

[ol9245]

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 ), 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 :

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

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 ;

[antoineleb]

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 :

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

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 :

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

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 :

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

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

[Invité]

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

Code :

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

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 :

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

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 :

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

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

[antoineleb]

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

Bonne journée