Problème de simplification

**Alecks** · 28/04/2010, 20h18

Bonsoir à tous.

J'ai un petit problème de simplification avec une des opérations de mon m-file et j'espère que vous pouvez m'aider.

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
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
 
%%vitesse
i=1; %%counter
wob=2; %%vitesse angulaire constante
 
for Th = 35:35  
    Th=Th/180*pi;
    syms B vax wab;
 
    Wab=[0 0 wab];
    Wob=[0 0 wob];
    l=0.45;
    r=0.1;
    rab=[l*cos(B) l*sin(B) 0];
    rob=[-r*cos(Th) r*sin(Th) 0];
 
    X1=cross(Wab,rab);
    X2=cross(Wob,rob);
 
    eq1=-vax+X1(1)+X2(1);
    eq2=X1(2)+X2(2);
    eq3=asin((r*sin(Th)+r)/l)-B;
 
    S=solve(eq1,eq2,eq3,'wab','vax','B');
 
    Beta=(180*double(S.B))/pi;
    B=double(S.B);
    vax=double(S.vax);
    wab=double(S.wab);
 
    %% store vax
     vax_table(i) = vax;
 
    % %% accélération
    syms Aax aphAB;
    AphAB=[0 0 aphAB];
 
    Wab=[0 0 wab];
    Wob=[0 0 wob];
 
    Y1=cross(AphAB, rab);
    Y2=cross(Wab,cross(Wab,rab));
    Y3=cross(Wob,cross(Wob,rob));
 
    eq4=-Aax+Y1(1)+Y2(1)+Y3(1);
    eq5=Y1(2)+Y2(2)+Y3(2);
 
 
    S=solve(eq4,eq5,'aphAB','Aax');
 
    aphAB=(S.aphAB);
    Aax=(S.Aax)
    %Aax_table(i)=ans
 
    i=i+1;
end

En exécutant ce code, on obtient un nombre très laid. J'ai essayé toutes les fonctions de simplication de je connais et il n'y a rien à faire. En fait, il n'y a pas de simplification à faire parce que tout est connu. D'ailleurs, si je copie/coller l'expression au long dans la console, il en ressort un chiffre réduit.

Mon problème est que j'ai besoin que ça se fasse tout seul dans mon script.

Quelle est la manière de faire?

Merci.

**Jean Dumoncel** · 28/04/2010, 21h44

Bonjour,

Code :

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

1
2
3
4
5
6
7
>> help ans
 ANS Most recent answer.
    ANS is the variable created automatically when expressions
    are not assigned to anything else. ANSwer.
 
    Reference page in Help browser
       doc ans

ans est donc utilisé lorsque tu effectues un calcul sans lui assigner de variable. Dans ton code, tous les calculs sont "rangés" dans des variables, la valeur qui est dans ans provient d'un de tes calculs précédent qui n'apparaît pas dans ton programme, donc je doute que tu veuille utiliser ce résultat.

C'est Aax ton résultat que tu trouves très laid? Tu utilises des variables symboliques, c'est notmal que le résultat contienne des variables symboliques et non un nombre. Si tu ne veux pas obtenir un résultat symbolique, pourquoi utilises-tu des variables symboliques?

Autre chose, dans ta boucle tu fais :

Code :

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

for Th = 35:35

donc Th varie de 35 à ... 35, il ne prendra donc que la valeur 35, c'est ce que tu souhaites?

**Alecks** · 28/04/2010, 22h23

Dans ma boucle, je veux faire 0 à 360 en fait, mais pour l'instant je teste, donc je fait une boucle unique pour la vitesse d'exécution.

Ensuite, si j'utilise des variables symboliques, c'est parce que je résous un système à 2 inconnus dont Aax :

S=solve(eq4,eq5,'aphAB','Aax')

Par la suite, mes deux variables ont leur valeurs mais elle ne sont pas réduites.

Voilà ce qui apparaît dans la console de Matlab si tu exécute mon script :

Code :

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

1
2
3
4
5
6
 
Aax =
 
-1/6338253001141147007483516026880*(430714436108490670451780542401*sin(B)^2+1454189027634344751696581754880*sin(B)+430714436108490670451780542401*cos(B)^2-2076797161242243356021573550080*cos(B))/cos(B)
 
>>

La variable B est donné dans le script donc elle est connu alors je devrais avoir un chiffre et non une expression comme celle-là. Si je la copie/colle dans la console, j'obtiens :

Code :

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

1
2
3
4
5
6
7
8
>> -1/6338253001141147007483516026880*(430714436108490670451780542401*sin(B)^2+1454189027634344751696581754880*sin(B)+430714436108490670451780542401*cos(B)^2-2076797161242243356021573550080*cos(B))/cos(B)

ans =

    0.1695

>>

C'est ce nombre réduit que je veux avoir mais sans le faire manuellement.

Y'a-t-il un moyen de détruire les sym maintenant que Aax à une valeur?
Est-ce que ça peut être ça qui cause le problème?

**Jean Dumoncel** · 28/04/2010, 22h43

Ah, oui, je n'avais pas vu que tu attribuais une valeur à B, dsl.

Donc pour obtenir la valeur numérique de Aax, une méthode (cela force juste B à prendre une valeur numérique) :

Code :

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