rdiepo
05-11-2008, 02:16
Ciao a tutti...non capisco una cosa!!!:muro:
Ho un sistema di due equazioni:
dx/dt = y
dy/dt= -x
Si tratta quindi di un sistema lineare della forma;
dz/dt = A*z
dove A=[0 1;-1 0]
dove gli autovalori sono pertanto +j e -j.
Il problema è questo:
Se utilizzo la funzione ode45 e faccio il grafico ottengo un andamento oscillante ma stabile delle due variabili (come dovrebbe essere)
MA se utilizzo il seguente codice(che è quello che devo usare):
ft=30;
dt=0.1;
T=0:dt:ft;
X0=[1 1]';
X=X0;
XArray=X
A=[0 1;-1 0]
for k=dt:dt:ft
X=X+(A*X)*dt;
XArray=[XArray X];
end
plot(T,XArray)
ottengo un andamento delle variabili che oscillano ma instabili (l'ampiezza aumenta all'infinito)
perchè?????? aiutooo
Ho un sistema di due equazioni:
dx/dt = y
dy/dt= -x
Si tratta quindi di un sistema lineare della forma;
dz/dt = A*z
dove A=[0 1;-1 0]
dove gli autovalori sono pertanto +j e -j.
Il problema è questo:
Se utilizzo la funzione ode45 e faccio il grafico ottengo un andamento oscillante ma stabile delle due variabili (come dovrebbe essere)
MA se utilizzo il seguente codice(che è quello che devo usare):
ft=30;
dt=0.1;
T=0:dt:ft;
X0=[1 1]';
X=X0;
XArray=X
A=[0 1;-1 0]
for k=dt:dt:ft
X=X+(A*X)*dt;
XArray=[XArray X];
end
plot(T,XArray)
ottengo un andamento delle variabili che oscillano ma instabili (l'ampiezza aumenta all'infinito)
perchè?????? aiutooo