Discussion :

[cedji]

Bonjour,
Je dispose des équations et données suivantes :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
omega= 2*pi*f;
lambdaI=j*omega*(phi -(2/3)*eta);
muI=j*omega*eta;
 
K1= sqrt((lambdaR+2*muR)^2 + (lambdaI+2*muI)^2);
K2= sqrt(lambdaR+2*muR+K1);
K3= sqrt(muR^2+muI^2);
K4= sqrt(muR+K3);
 
vL=sqrt(2/density)*K1/K2
vT=sqrt(2/density)*K3/K4

A l'aide de mesures, je peux connaitre vL, vT,lambdaR et muR.
Je souhaite déterminer eta et phi.

Je pensais utiliser la fonction solve mais je ne sais pas de quelle manière car il y a beaucoup d'expressions!
Quelqu'un aurait il une piste à me donner svp ?
J'utilise Matlab 6.1

Merci de votre aide !

[cedji]

Bon j'ai trouvé une idée qui est un peu longue à réaliser... Je compte resoudre les systèmes 2 à 2 pour remonter à la fin aux paramètres inconnus.
la première étape serait de resoudre le sytème suivant :

Code :

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

1
2
K4= sqrt(muR+K3);
vT=sqrt(2/density)*K3/K4

A partir de ce système je veux déterminer l'expression de K3 en fonction de density et vT. Je l'ai fait à la main mais bon j'aimerai retrouver le calcul par matlab car j'aurai besoin de le calculer pour différentes valeurs de density et vT.

J'ai réussi à obtenir l'expression de K3 avec la fonction solve mais je n'arrive pas à remplacer les constantes 'density' et 'vT' par les valeurs que je souhaite...
Qqn pourrait-il m'aider ??
Merci d'avance

[cedji]

oups en fait j'ai parlé trop vite. Je n'ai pas réussi à obtenir l'expression de K3 en fonction des autres paramètres...
Voila le code que j'ai écrit :

Code :

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

S=solve('K4= sqrt(muR+K3)','vT=sqrt(2/density)*K3/K4','K3','K4')

et voila le message d'erreur que je recois :

Code :

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

1
2
3
4
5
6
Warning: Explicit solution could not be found.
> In C:\MATLAB6p1\toolbox\symbolic\solve.m at line 136
 
S =
 
[ empty sym ]

Quel est le problème ?
Car meme lorsque je donne des valeurs aux variables density, muR et vT avant la résolution de l'equation ca ne marche pas non plus...

[Myrne]

Ce n'est pas un message d'erreur mais une annonce pour te dire qu'aucune solution n'a été trouvée au système.

Ton couple d'équations n'est pas particulièrement compliqué, si j'étais toi je poserai les équations tranquillement sur le papier pour voir s'il est possible de trouver une solution. Tu vas probablement arriver à un système du second degré, je suppose que tu connais les intervalles de valeur envisageables pour tes paramètres, tu dois donc pouvoir dire s'il y a des solutions à ton équation ou non.

EDIT : J'ai posé ça vite fait, j'en arrive à

Code :

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

K4² - K4 (vT.sqrt(density/2)) - muR = 0

à toi de dire si ce truc est faisable avec tes paramètres.

[FR119492]

Salut!
Pourquoi ne pas utiliser la méthode de Newton-Raphson?
Jean-Marc Blanc

[membreComplexe12]

Salut!
Pourquoi ne pas utiliser la méthode de Newton-Raphson?
Jean-Marc Blanc

moi j'aurais fais ceci aussi. De plus il existe peut etre une fonction toute prete pour cela?

sinon je pense que ce n'est pas très difficile à implementer dans ton cas