__________________
Senast redigerad av Tom.Of.Finland 2017-06-17 kl. 21:04.
Senast redigerad av Tom.Of.Finland 2017-06-17 kl. 21:04.
Citat:
Jag vet inte vad du menar att den är analog/tidskontinuerlig, den är självfallet diskret i tiden precis som alla ode lösare är, det finns ingen annan typ av numerisk lösare. Problemet är att den är adaptiv, eftersom du lägger på ett slumpmässigt fel så kommer den göra tidssteget extremt litet över hela lösningsintervallet vilket kommer göra att det tar väldigt lång tid. Sedan så anpassas inte standardavvikelsen efter tidssteget på ett korrekt sätt heller, vilket gör att det går åt skogen. Du måste helt enkelt använda dig av någon annan lösare än ode45.
Okej. ODE45 må vara diskret, men matematiken är analog.
Jag har hittat en lösare: lsim
Citat:
Hejsan, jag snår tråden lite men misstänker att det finns en del kunskap i den här träden för att kunna besvara en imo enklare fråga.
Vilken ODE ska man använda till system av diff. ekv som är first order linary diff eq?
Jag använder just nu ode45 till det och det verkar funka. ode15s funkar också men den är lite långsammare.
Jag skriver ner systemet lika gärna
ds/dt = −k1se + k2c
de/dt = −k1se + k2c + k3c
dc/dt = k1se − k2c − k3c
dp/dt = k3c
Har googlat som faen men hittar inget.
Citat:
Hejsan, jag snår tråden lite men misstänker att det finns en del kunskap i den här träden för att kunna besvara en imo enklare fråga.
Vilken ODE ska man använda till system av diff. ekv som är first order linary diff eq?
Jag använder just nu ode45 till det och det verkar funka. ode15s funkar också men den är lite långsammare.
Jag skriver ner systemet lika gärna
ds/dt = −k1se + k2c
de/dt = −k1se + k2c + k3c
dc/dt = k1se − k2c − k3c
dp/dt = k3c
Har googlat som faen men hittar inget.
Vilken ODE ska man använda till system av diff. ekv som är first order linary diff eq?
Jag använder just nu ode45 till det och det verkar funka. ode15s funkar också men den är lite långsammare.
Jag skriver ner systemet lika gärna
ds/dt = −k1se + k2c
de/dt = −k1se + k2c + k3c
dc/dt = k1se − k2c − k3c
dp/dt = k3c
Har googlat som faen men hittar inget.
I regel börjar man alltid med ODE45. Fungerar inte det så testa ODE23S. Men skit i ODE-lösare! Bygg en själv! Du har mer frihet över den då. Lyssna här!
Jag skulle rekommendera dig att bygga om detta till en tillståndsmodell och simulera detta med lsim från GNU Octave med kontrolpaketet "Control".
Jag skulle också rekommendera dig att bygga detta system. Detta är skrivet i C-kod.
Citat:
n = 1000; iterationer
dt = 0.1; sekunder
u = 10; Insignal t.ex. 10 volt
for(int i = 0; i < n; i++)
{
dx = A*x + B*u;
y[i] = C*x;
T[i] = t;
x = x + dx*t; // Uppdatera nästa tillstånd
}
dt = 0.1; sekunder
u = 10; Insignal t.ex. 10 volt
for(int i = 0; i < n; i++)
{
dx = A*x + B*u;
y[i] = C*x;
T[i] = t;
x = x + dx*t; // Uppdatera nästa tillstånd
}
Plotta sedan T och y.
Citat:
Hejsan, jag snår tråden lite men misstänker att det finns en del kunskap i den här träden för att kunna besvara en imo enklare fråga.
Vilken ODE ska man använda till system av diff. ekv som är first order linary diff eq?
Jag använder just nu ode45 till det och det verkar funka. ode15s funkar också men den är lite långsammare.
Jag skriver ner systemet lika gärna
ds/dt = −k1se + k2c
de/dt = −k1se + k2c + k3c
dc/dt = k1se − k2c − k3c
dp/dt = k3c
Har googlat som faen men hittar inget.
Vilken ODE ska man använda till system av diff. ekv som är first order linary diff eq?
Jag använder just nu ode45 till det och det verkar funka. ode15s funkar också men den är lite långsammare.
Jag skriver ner systemet lika gärna
ds/dt = −k1se + k2c
de/dt = −k1se + k2c + k3c
dc/dt = k1se − k2c − k3c
dp/dt = k3c
Har googlat som faen men hittar inget.
Skillnaden på ode45 och ode15s är att ode45 bygger på en Runge Kutta metod, medan ode15s används för styva problem. Detta betyder ungefär att ode15s löser ett ekvationssystem i varje tidssteg, detta gör i regel att man får en stabilare lösning men det blir långsammare. Däremot att använda en ode45 på ett styvt problem kommer inte fungera bra, om den klarar av att lösa systemet korrekt så kommer den ha kortat ned tidssteget så mycket så att ode15s ändå kommer att vara snabbare än ode45.
Eftersom ode45 verkar fungera för dig så är det den du bör använda dig av.
Skapa ett konto eller logga in för att kommentera
Du måste vara medlem för att kunna kommentera