Discussion :

[jbollagnier]

Bonjour !

Je viens d'écrire un programme qui doit me résoudre une équation du type (c'est un example - dans mon programme l'équation est beaucoup plus compliqué et la solution aussi!):

Y^2+A=0

Y est l'inconnue et A une constante

l'equation dans Matlab est donc comme ceci:

------
syms Y;

for i:1..5


a=Y^2-A(i);

solve(a);

end

---------
La réponse est donc la suivante et n'est pas un double!:

sqrt(A(i));-sqrt(A(i))

J'aimerais cependant convertir ces deux solutions en "double" et les mettres dans un vecteur du type

Racine1(i)=sqrt(A(i))
Racine2(i)=-sqrt(A(i))

Quelqu'un saurait-il comment faire?

Merci d'avance !

Bonne journée!

JB

[tug83]

Pour résoudre y²=4, je ferai juste:

Code :

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

solve('y^2=4')

qui retourne

ans =

2
-2

[jbollagnier]

Oui mais ceci est juste un example!

J'ai une multitude d'équations à résoudre et c'est pour cela que je veux convertire les solutions de "sym2 en "double"!

Si j'utilise

Racine1(i)=double(solve(a));

il ne peut pas le faire car solve(a) contient deux racines !

La question est juste de convertir les deux solutions en "double".

J'espere que j'ai été plus clair !

Merci d'avance pour votre aide.

JB

[tug83]

J'espere que j'ai été plus clair !

non!

donne un code MATLAB en exemple ce sera plus simple

[jbollagnier]

Ok voici mon programme:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
clear all;
load stressdata3.mat
 
% données : v, B0, Frate, E, h, vp, Ep, Strain
 
z=length(Ep)
 
A=zeros(z,1);
B=zeros(z,1);
C=zeros(z,1);
D=zeros(z,1);
e=zeros(z,1);
Interp=zeros(z,1);
 
 
Sigma=zeros(z,1);
 
 
syms Z;
 
for i=1:z-1
 
    Interp(i)=Ep(i)*(Time(i+1)-Time(i));
 
    A(i)=((v-1)*B0*Frate(i))/(E*2*h)+(vp(i)-1)/Ep(i); 
    B(i)= 2*vp(i)*(1-v)*B0/(Ep(i)*E*2*h);
    C(i)=(1-v)*B0*Strain(i)/2/h/E;
    D(i)=-(1-v)*B0*2*vp(i)/(2*h*E*Interp(i));
    e(i)=-Frate(i);
 
      M(i)=A(i)*(Z+Sigma(i))/2+B(i)*((Z+Sigma(i))/2)^2+(Z-Sigma(i))/(Time(i+1)-Time(i))*C(i)+(Z-Sigma(i))/(Time(i+1)-Time(i))*D(i)*(Time(i+1)-Time(i))*sum(Sigma)+e(i);
 
 
    Sigma(i)=double(Solve(M(i)));
 
 
 
 
end

M(i) est l'équation à resoudre et Z est l'inconnue.
Je veux maintenant avoir les deux solutions de Z dans un vecteur.

Le programme que je viens de taper me renvoi l'erreur suivante:

??? In an assignment A(I) = B, the number of elements in B and
I must be the same.

Error in ==> Stressnumcalc at 35
Sigma(i)=double(Solve(M(i)));

le probleme provient lors de cette conversion "sym" -> "double"
et j'aimerais donc savoir comment procéder.

J'espere que cette dois cela a été plus clair.

Merci d'avance pour votre réponse!

JB

[tug83]

en gros c'est une erreur de programmation qui est reproductible en faisant:

Code :

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

a(1)=double(solve('y^2=4'))

et dont la solution est dans notre FAQ

[jbollagnier]

j'y vois un peu plus clair:

si je mets racine(i:i+1)=double(solve('y^2=4'))

j'aurais
racine(i)=2
racine(i+1)=-2

mais ce que je voudrais c'est avoir deux vecteurs différents !
du style:

racine1(i),racine2(i)=double(solve('y^2=4'))

si quelqu'un sait comment faire !

Merci !

JB

[Caro-Line]

La FAQ, toujours la FAQ....
Comment créer des variables nommées A1, A2, A3, ...,AN ?

[tug83]

en + c'est la note technique préférée de caro

[jbollagnier]

En effet !

En écrivant simplement:

K=double(Solve(M(i)));
Sigma(i)=K(2);
Sigma2(i)=K(1);

J'obtiens ce que je voulais.

Merci de m'avoir éclairer !!!!!!!


JB

[Caro-Line]

Oui... MAIS ce serait plus propre de ne pas créer 2 variables.

D'utiliser soit une matrice (si tu es sûr d'avoir toujours 2 racines), soit un cell array.
Sur un petit exemple inspiré du 1er.
Cas matrice :

Code :

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

1
2
3
4
5
6
syms Y
for i=1:5
    a=Y^2-2*i;
    b(i,:)=solve(a);
end
resultat=double(b)

resultat(i) te donne pour chaque i tes 2 racines

Cas cell array :

Code :

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

1
2
3
4
5
6
syms Y
for i=1:5
    a=Y^2-2*i;
    b{i}=solve(a);
end
resultat=cellfun(@double,b','UniformOutput',false)

resultat{i} te donne pour chaque i tes racines