Discussion :

[math003]

Bonjour, je travaille actuellement sur la création d'un programme sous MATLAB.
Cependant je rencontre quelque difficultés, je vous explique mon problème :

En faite je possède 3 équations du type :

a=b*X+c*(-d+e*R)
X=(f-2*Y/g)*h+i
Y=j*Z^3+k*Z^2+m

avec a, b, c, d, e, f, g, h, i, j, k, m des constantes connues
et X, Y, Z, R des variables.

Pour mon problème je cherche R, je dois faire varier Z ce qui me donnera Y puis X et au final R.

Je veux concrètement les valeurs de R pour les différentes valeurs imposées de Z.

Je pense qu'il faut utiliser des boucles imbriquées cependant je n'arrive vraiment pas à rédiger ce programme.

Pouvez-vous m'aider s'il vous plait.
Cordialement
Mathieu

[Merel]

Salut,

Montre nous ce que tu as déjà programmé avec les erreurs ou les étapes où tu bloques. A partir de là, on pourra t'aider ;-).

[math003]

Pour l'instant je n'ai rien fait en gros, j'ai seulement fait la 1ere boucle :

Code :

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

1
2
3
4
5
%%Boucle calcul de To(z+deltaZ/2)
To=Alpha.*Tomax
for z=0:-0.005:-0.2
  To=A.*(z+deltaZ/2)^3+B.*(z+deltaZ/2)^2+D
end

ce qui correspond à la 3eme équation dans mon 1er message.

Seulement mon problème est que je ne sais vraiment pas comment procéder pour résoudre mon problème, c'est à dire qu'elle outil de matlab utiliser.

[Invité]

Bonjour,

Comme tu le dis, il suffit de faire le chemin Z => Y => X => R, donc commence par ré-écrire les équations pour cela (que la première en fait).
Il te suffira par la suite d'écrire quelque chose comme :

Code :

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

1
2
3
4
Z = ...
Y = f(Z);
X = g(Y);
R = h(X);

Préfère écrire ici :

Code :

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

1
2
3
4
z=0:-0.005:-0.2
for i = 1:length(z)
    z(i)
end

[math003]

Désolé mais je suis loin de gérer MATLAB alors à quoi correspond le "i" dans le 2eme encadré ?

Et pour le 1er encadré, je dois définir par la suite f, g, h ??
Et si oui comment ?
Ce sont bien des fonctions ?

[Invité]

i est un indice de tableau qui va varier de 1 au nombre d'éléments de z (=length(z)).
f, g, h désignent bien de façon générique des fonctions, voulant ici simplement dire

A = f(B) <=> A exprimé en fonction B

[math003]

Mais je dois tout de même définir A du type :

A=...... * B ???

et je le fais après avoir dit : A=f(B) ?

Je n'utilises donc plu de boucles ?

[math003]

Pour l'instant j'ai écrit ma 1ere équation :

Rs=[[Pi*Rp^2+Rp^2*((SigmaS-(Es/Ep)*SigmaP)/CoefdepoissonS+Pi)]/[(SigmaS-(Es/Ep)*SigmaP)/CoefdepoissonS+2*Pi-2*Po*Pi]]^(1/2)

Rs=h(SigmaP)
SigmaP=f(To)
To=g(z)

Et pour la suite je suis perdu, il me manque des morceaux