Newsgroup: it.scienza.chimica
Titolo: preda predatore superpredatore

Ciao a tutti vi dico il mio problema. ho il seguente sistema di
equazioni differenziali non lineari del terzo ordine:
dx1/dt= x1* (r*(1-(x1/50))-(2/(8+x1)*x2));
dx2/dt=x2*(-1.5+(2*(2/8+x1)*x1)-((2/5+x1)*x3)
dx3/dt=x3*(-0.5+(2*(2/5+x2)*x2)
Dato questo sistema devo trovare in matlab i massimi locali della
prima variabile di stato. una volta trovati li dovei plottare in un
grafico in cui in ordinata metto questi massimi locali in ascissa il
parametro di controllo r che dovrà variare tra 0:3 con un passo pari a
0.1
Come posso risolvere questo?

 Come posso risolvere questo?
Non lo so, ma ti posso dire che sei leggermente OT :-)
-- 
news:it-alt.comp.software.openoffice : Il newsgroup dedicato a
OpenOffice.org, la suite open source di applicazioni per ufficio.
Saluti, VITRIOL


OT...Sarebbe???

 OT...Sarebbe???
Fuori tema. Qui si parla di chimica.
Prova a postare su it.scienza.matematica.



In effetti è OT, comunque quelle equazioni hanno forma analoga alle
equazioni cinetiche di un sistema di reazioni oscillanti.

"raffa" ha scritto:
Quello che ormai ha gia' scritto su diversi ng, sembra
tuttavia senza dare ascolto alle risposte ottenute.
Ti e' gia' stato detto, da me e da altri, che le equazioni
che hai scritto sono sbagliate, sia da un punto di vista
sintattico che semantico. Inoltre per risolvere quel
sistema di equazioni devi assegnare le condizioni
iniziali per le variabili x1, x2, x3. Cosa tu voglia poi fare,
una volta risolto il sistema, non riesco a capirlo da
quello che hai scritto, mi permetto comunque di consigliarti
di cercare _prima_ di chiarirti le idee su cio' che devi fare
e su come lo devi fare, e _poi_, se necessario, chiedere
aiuto su un ng (non qui).
Ciao
-- 
Giorgio Bibbiani

On 16 Giu, 12:21, "Giorgio Bibbiani"

A me sembra che le equazioni non siano sbagliate.
Per questo genere di equazioni le condizioni iniziali, almeno entro
"ragionevoli" intervalli, non sono importantissime. Dopo poche
oscillazioni si raggiunge uno "stato stabile", il "ciclo limite**.
Quindi quello che chiede non è scorretto, la funzione x1(t) avrà una
forma a picchi (i massimi) che avranno valori diversi a seconda del
parametro "r".
Comunque non conosco questo caso specifico, per le considerazioni che
ho fatto sopra avevo in mente il "bruxellatore" che gli assomiglia ma
è un po' più semplice.
**Se però si "sbaglia" a scegliere le condizioni iniziali può capitare
che una specie si estingua o che si entri in un regime caotico.


Le equazioni non sono sbagliate in quanto mi sono state assegnate e le
ho opportunamente controllate con il materiale che mi è stato dato.
Seconda cosa forse non mi sono spiegata bene.Ho Assegnato le
condizioni iniziali al sistema ho settato opportunamente i parametri
ed è stata fatta la simulazione del sistema nel tempo delle tre
variabili di stato e l'andamento di queste nel piano delle fasi.
Credo che l'unico che abbia compreso il problema riguarda l'autore
dell'ultimo post:in sistemi lineari del terzo ordine per determinate
condizoni iniziale si possono innescare comportamenti caotici
il mio obiettivo è fare vedere questo mediante simulazioni in matlab e
mi chiedevo se qualcuno potesse darmi qualche dritta per la simulazione

parametro "r".
...
Quello che chiede io non l'ho capito.
Ciao
-- 
Giorgio Bibbiani

On 16 Giu, 16:14, "Giorgio Bibbiani"

Le condizioni iniziali sono
X01=10;
X02=15;
X03=10
hai ragione per ciò che concerne la scrittura delle formula mancano
parentesi...
dx1/dt= (x1* r*(1-(x1/50)))-((2*x1*x2)/(8+x1));
dx2/dt=(x2*(-1.5))+((2*x1*x2)/(8+x1))-((2*x2*x3)/(5+x1))
dx3/dt=(x3*(-0.5))+((2*x2*x3)/(5+x2))
ecco la forma corretta cosa nn capisci?

ecco la forma corretta cosa nn capisci?
Ho letto il tuo ultimo messaggio, in cui specifichi meglio,
anche se ancora in termini generici, la tua richiesta:
"mi chiedevo se qualcuno potesse darmi qualche
dritta per la simulazione".
Io, a parte quanto gia' detto, non saprei cosa consigliarti.
Ciao
-- 
Giorgio Bibbiani


Va bene ma son dettagli. Quando avrà pronto il programma per
l'integrazione potrà fare tutte le prove che vorrà con i più svariati
set di condizioni iniziali.
La domanda che secondo me andrebbe posta a raffa sarebbe invece se ha
problemi con matlab o se non sa che algoritmo usare per integrare le
equazioni.


Ho problemi cn matlab!! mai usato purtroppo...per la risoluzione del
sistema ho usato il metodo runke kutta..la ode 45


Non conosco Matlab, quindi non ti saprei aiutare.
Prova su it.istruzione.universita.ingegneria e/o
it.discussioni.ingegneria
Oppure
http://www.google.com/search?client=opera&rls=it&q=matlab+runge+kutta&sourceid=opera&ie=utf-8&oe=utf-8