|
|
|
Simulation eines biologischen Systems |
|
| Julia |
Gast
|
 |
Beiträge: ---
|
|
|
|
Anmeldedatum: ---
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
 |
|
Verfasst am: 05.12.2011, 19:04
Titel: Simulation eines biologischen Systems
|
 |
| |
 |
 |
Mein Problem: Simuliere as VZ2 glied min matlab mit hilfe der ode solver.
Wechselwirkung Räuber -Beute:Lotka-Volterra-Gl.
x(punkt)1 = (a1 +b1*x1 +c1*x2)*x1 Beute
x(punkt)2 = a2 +b2x2+c2x1)x2 räuber
x1,2 Anzahl der individuen erster, zeiter art
a1,2 wachstumskoeff. inerhalb der eigenen art
b1,2 wechselwirkung innerhalb der eigenen art
c1,2 wechselwirkung zwischen den arten
folgende zahlenwerte sind zu verwenden: a1=6;b1=-1/2;c1=-3;a2=-3;b2=-3,c2=1,
meine bisherige unvollständige lösung
function []=raeuber_beute()
x1 = 5; %Anfangsbestand Beute
x2 = 1; %Anfangsbestand Räuber
[t,y] = odeEuler(@dgl,[0 30],[x1 x2])
plot(y(:,1),y(:,2),'--g');
function [t,y]=odeEuler(dgl,t_sim,x)
dt = 0.001;
t_akt = t_sim(1);
x_akt = x;
i = 1;
while (t_akt < t_sim(2))
t(i)= t_akt;
y(i, =x_akt;
i =i+1
function [dx]=dgl(t,x)
a1 = 6;
b1 = -0.5;
c1 = -3;
a2 = -3;
b2 = -3;
c2 = 1;
dx1=(a1 + b1*x1 + c1*x2)*x1;
dx2=(a2 + b2*x2 + c2*x1)*x2;
falls mir jemand helfen könnte es zu vervollständigen wäre ich sehr dankbar.
grüße julia
|
|
|
|
|
|
| Andreas Goser |
Forum-Meister
|
|
Beiträge: 3.654
|
|
|
|
Anmeldedatum: 04.12.08
|
|
|
|
Wohnort: Ismaning
|
|
|
|
Version: 1.0
|
|
|
|
|
|
|
Verfasst am: 06.12.2011, 08:25
Titel:
|
|
Ist das wirklich in MATLAB zu lösen und nicht in Simulink? Das ist eigentlich eine perfekte Nutzungs von Simulink.
Andreas
|
|
|
|
| steve |
Ehrenmitglied
|
|
Beiträge: 2.028
|
|
|
|
Anmeldedatum: 03.09.07
|
|
|
|
Wohnort: Wien
|
|
|
|
Version: R2024a
|
|
|
|
|
|
|
Verfasst am: 06.12.2011, 09:10
Titel:
|
|
Moin,
@Andreas: Ja, es wäre sicherlich eine super Anwendung für Simulink - allerdings kenne ich es aus dem Studium auch noch, dass solche Probleme mit MATLAB gelöst werden sollten. Es geht dabei ausschließlich darum die Verwendung der ODE-Solver und die Programmierung von DGL-Systemen zu üben.
@Julia:
1. Verwende bitte zukünftig die Code-Umgebung, dann ist dein Posting viel besser zu lesen.
2. An welcher Stelle hast du denn Probleme/Sorgen? Gibt es Fehlermeldungen? Oder möchtest du gerne eine Komplettlösung haben (machen wir hier äußerst ungerne  )
Gruß
Alex
_________________
>> I told me to.
____________________________________
Matlab Cheat Sheet
goMatlab-Knigge - dran gehalten?!
Schon in den FAQ gesucht?
Ist vielleicht bei den Skripten oder den Tutorials was für dich dabei?
|
|
|
|
|
|
|
Einstellungen und Berechtigungen
|
|
Du kannst Beiträge in dieses Forum schreiben.
Du kannst auf Beiträge in diesem Forum antworten.
Du kannst deine Beiträge in diesem Forum nicht bearbeiten.
Du kannst deine Beiträge in diesem Forum nicht löschen.
Du kannst an Umfragen in diesem Forum nicht mitmachen.
Du kannst Dateien in diesem Forum posten
Du kannst Dateien in diesem Forum herunterladen
|
|
Impressum
| Nutzungsbedingungen
| Datenschutz
| FAQ
|  RSS
Hosted by:
Copyright © 2007 - 2024
goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks
MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.
|
|
