Discussion :

[sakona]

Bonjour,

Je voudrais trouver l'angle alpha , solution de l'équation suivante:

10*Fm*(W*sin(alpha))-(10*c*(b/cos(alpha))+((W*cos(alpha)-u*(b/cos(alpha)))*tan(Phi)))=0

W,c,b,u et Phi sont connus.

J'ai essayé avec ce code:

Code :

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

1
2
3
4
5
6
7
8
alpha=0;
 
while 10*Fm*(W*sin(alpha))-(10*c*(b/cos(alpha))+((W*cos(alpha)-u*(b/cos(alpha)))*tan(Phi)))~=0
 
    alpha=alpha+1;
end
 
Sol=alpha;

Je me bloque devant cette erreur:

??? Error using ==> alpha
Too many output arguments.

Merci de m'aider à résoudre ce problème.

[le fab]

salut

alpha est une fonction matlab
prend un autre nom pour ta variable (monAlpha par exemple)

Fabien

[sakona]

Merci pour votre réponse

[doinier]

salut,

pourquoi ne pas remplacer alpha par j, et commencer ton incrémentation par 1 au lieu de 0, i.e:

Code :

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

1
2
3
j=1
(....)
j=j+1

par exemple.

[Dam2227]

Salut Sakona.

Tu vas aussi avoir un autre problème. C'est que ton équation ne va pas forcément s'annuler pour des valeurs de alpha (ou monalpha à présent) entières. il faut donc que tu trouves quand la fonction s'approchera au plus près de 0 et ton alpha sera une valeur approchée à +/-1. Étant donné que tu cherches un angle, alors tu peux définir monalpha sur l'intervalle 0:359, évalue ta fonction avec ce vecteur monalpha puis cherche le minimum des valeurs absolues de ta fonction. Enfin, dans Matlab, les fonctions sin cos etc fonctionnent avec des radiants (0 -> 2pi), pas des degrés. Un exemple dans lequel j'essaye de trouver le minimum de la fonction cos:

Code :

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

1
2
3
4
5
 
monalpha = 0.5:1:359.5;
res = cos(monalpha*2*pi/360);
[val,pos] = min(abs(res));
monalpha(pos)

Après, il y a des histoires de symétrie alors à toi de choisir la valeur qui correspond à ton problème.

Enfin, une méthode peut-être plus adaptée et plus robuste pour ton problème est la méthode de dichotomie: http://fr.wikipedia.org/wiki/M%C3%A9thode_de_dichotomie

[sakona]

Bonjour Dam2227

Effectivement, j'essaye de trouver une valeur proche de 0.
Avec tes propositions importantes, je chercherai la bonne valeur de alpha.

[le fab]

j est également un mot clef matlab !
à éviter

[Jerome Briot]

j est également un mot clef matlab !

« j » n'est pas à proprement parler un mot-clé, pas plus que « i ».

Il faut néanmoins connaitre l’ambiguïté qui existe avec la notation complexe : Peut-on utiliser « i » comme indice de boucle dans MATLAB ?