Pendolo semplice con matlab

di il
1 risposte

Pendolo semplice con matlab

Salve a tutti ragazzi, sono nuovo del forum e sono anche un novello programmatore di matlab.
Per dirla tutta sono diciamo obbligato a imparare per motivi universitari.
Comunque volevo chiedere il vostro aiuto per un problema che ho con lo sviluppo del programma di un esercizio nel quale mi sono impuntanto.
-------------------------------------------------------------
Allora l'esercizio proposto è questo:

INTEGRAZIONE EQUAZIONI DIFFERENZIALI
-Studiare l'evoluzione nel tempo di un sitema non lineare (il pendolo semplice) a partire da diverse condizioni iniziali.

Allo scopo si ricorda che il pendolo semplice è retto dalle seguenti equazioni

x1° = x2
x2° = -(g/l)*sin(x1)

se non c'è smorzamento, altrimenti le equazioni sono

x1° = x2
x2° = -(g/l)*sin(x1)-(ß/M)*x2

In ogni caso x1= ? e x2= ?° dove ? è l'angolo che il braccio del pendolo fa rispetto alla posizione di riposo. Con l e M si sono indicate rispettivamente la lunghezza e la massa del pendolo, mentre ß rappresenta il coefficiente di attrito viscoso (assumere valori numerici per il calcolo).
Si può anche linearizzare il sistema attorno al punto di equilibrio e valutare la differenza di comportamento dei due sistemi quando le condizioni iniziali permettono effettivamente la linearizzazione e quando ciò non è più valido.

COMANDI DI BASE: ode23, function, figure.
---------------------------------------------------------
Ho provato ad abbozzare qualcosa che non so innanzitutto se è giusto o meno, in più ci sono pure degli errori che nn riesco a trovare che matlab mi segnala continuamente:
m = 10;
b = 15;
g = 9.81;
L = 1;

xbar = [pi/4 0];

ubar = m*g;

function xdot=equaz(t,x);
f=20;
B=(b/m);
k=(g/L);
xdot(1)=x(2);
xdot(2)=-B*x(2)-k*x(1)+f;
xdot=xdot';

ode23('equaz',[tinizio,tfine],[valore iniz 1,val iniz 2])
Aiutatemi per favore ci sto uscendo di testaaaaaa....

1 Risposte

Devi accedere o registrarti per scrivere nel forum
1 risposte