Discussion :

[nidjzaf]

bonsoir

ça fait une semaine que je travaille sur ce programme mais j n'arrive pas à comprendre où est le problème donc s'il vous plait aidez moi à trouver la faute et aussi comment est ce que je dois écrire sous matlab que l'entrée u doit être inférieure a 180
NB:je travaille sur matlab 6.5
et merci d'avance pour vos réponse

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
function xdot = pendul(t, x)
m =0.1;                        
mc =1;                          
M =m+mc;
g=9.8;                    
l=0.5;                              
 
d=4/3*l*M-m*l*(cos(x(1)))^2;
F=M*g*sin(x(1))-m*l*(x(2))^2*cos(x(1))*sin(x(1))/d;
G=cos(x(1))/d;
 
u=180;
xdot(1)= x(2);
xdot(2)=F+G*u
x0=[-pi/60 0];
tspan=[0 20];
[t x]=ode45('pendule',tspan,x0);
plot(t,u(:,1));

[FR119492]

Salut!
Tout d'abord, on ne résout pas une équation ou un système différentiel, on l'intègre.

En outre, il nous serait plus facile de te dépanner si tu nous disais quel est ton système différentiel et comment tu vois que la réponse n'est pas correcte.
Jean-Marc Blanc

[salseropom]

Salut et bienvenue sur ce forum.

Ton programme me paraît très bizarre : il n'a qu'une seule fonction où tout est mélangé. Voici ce que je te propose (non testé) :

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
 
function xdot = pendul(t, x, u)
 
xdot = zeros(2,1); % xdot doit etre un vecteur colonne
 
m =0.1;                        
mc =1;                          
M =m+mc;
g=9.8;                    
l=0.5;                              
 
d=4/3*l*M-m*l*(cos(x(1)))^2;
F=M*g*sin(x(1))-m*l*(x(2))^2*cos(x(1))*sin(x(1))/d;
G=cos(x(1))/d;
 
u=180;
xdot(1)= x(2);
xdot(2)=F+G*u

et voici ton programme principal : main.m

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
 
 
u = 180;
 
if(u > 180)
  error('u doit etre inferieur ou egal a 180')
end
 
x0=[-pi/60 ; 0]; % x0 doit etre un vecteur colonne
tspan=[0 20];
[t x]=ode45(@pendule,tspan,x0,[],u);
plot(t,u(:,1));

[nidjzaf]

Salut et bienvenue sur ce forum.

Ton programme me paraît très bizarre : il n'a qu'une seule fonction où tout est mélangé. Voici ce que je te propose (non testé) :

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
 
function xdot = pendul(t, x, u)
 
xdot = zeros(2,1); % xdot doit etre un vecteur colonne
 
m =0.1;                        
mc =1;                          
M =m+mc;
g=9.8;                    
l=0.5;       
 
 
 
 
 
d=4/3*l*M-m*l*(cos(x(1)))^2;
F=M*g*sin(x(1))-m*l*(x(2))^2*cos(x(1))*sin(x(1))/d;
G=cos(x(1))/d;
 
u=180;
xdot(1)= x(2);
xdot(2)=F+G*u

et voici ton programme principal : main.m

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
 
 
u = 180;
 
if(u > 180)
  error('u doit etre inferieur ou egal a 180')
end
 
x0=[-pi/60 ; 0]; % x0 doit etre un vecteur colonne
tspan=[0 20];
[t x]=ode45(@pendule,tspan,x0,[],u);
plot(t,u(:,1));

merci salseropom
le programme n'a pas marché
j'ai toujours le méme probleme (input argument 'x' is undefined.)

[phryte]

Bonjour.
Le programme marche bien avec :

Code :

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

1
2
[t x]=ode23(@pendul,tspan,x0);
plot(t,x);

attention la function est pendul (sans le e)

[salseropom]

merci salseropom
le programme n'a pas marché
j'ai toujours le méme probleme (input argument 'x' is undefined.)

Salut,
effectivement, phryte a raison, c'est pendul (sans le e). Ensuite plot(t,u(:,1)) ne veut rien dire car ta sortie c'est x et non u.

Je n'ai pas matlab sur moi. A quelle ligne ça plante. Si, comme dit phryte, ca marche avec ode23, ça doit marcher avec odeXX (ode45, ode15s, ode23tb etc...) car la syntaxe est la même. Après, qui dit algo différents dit résultats différents (ou légèrement différents)

Si ça plante toujours, peux-tu faire un copier - coller de la sortie Matlab ?

[FR119492]

Salut!

ce que je dois faire c'est de tracer le trajectoire (théta) d'un pendule invérse je voulé vous ecrire son equation différentiel mais je ne sais comment le faire dans ce forum ,elle est écrite dans le programme que j'ai écrit (xdot(2) )

Je reformule le problème:
"Ecrire les équations qui régissent le comportement d'un pendule inverse, puis les intégrer en n'utilisant que du papier et un crayon."

Est-tu capable, en faisant abstraction du temps nécessaire, de résoudre ce problème?
Jean-Marc Blanc

[nidjzaf]

Salut!
Tout d'abord, on ne résout pas une équation ou un système différentiel, on l'intègre.

En outre, il nous serait plus facile de te dépanner si tu nous disais quel est ton système différentiel et comment tu vois que la réponse n'est pas correcte.
Jean-Marc Blanc

désolé pour le retard
ce que je dois faire c'est de tracer le trajectoire (théta) d'un pendule invérse je voulé vous ecrire son equation différentiel mais je ne sais comment le faire dans ce forum ,elle est écrite dans le programme que j'ai écrit (xdot(2) )