作业帮 > 综合 > 作业

急!求下面matlab代码每行都是什么意思?最好每行给个注释?

来源:学生作业帮 编辑:作业帮 分类:综合作业 时间:2024/05/24 13:57:02
急!求下面matlab代码每行都是什么意思?最好每行给个注释?
在没有人类介入的生物圈中,只考虑野兔和狐狸这两种动物的生存和发展.假设这一环境中有足够多的食物供野兔生存,而狐狸只是以野兔为食物.方程中x,y项反映了两生物种群的相互作用对于某一种数量变化的影响.因为狐狸的存在对野兔不利,所以-a0 .
取参数a=0.01,b=0.02 ,自变量t的变化范围(0, 15),初值条件为
x(0)=20,y(0)=20
function Lotka_Volterra
t_end=15;
x0=[20;20];
[t,x]=ode45(@LV_model,[0,t_end],x0)
plot(t,x);
legend('野兔','狐狸')
function y=LV_model(t,x)
a=0.01;b=0.02;
y=[x(1)-a*x(1)*x(2);-x(2)+b*x(1)*x(2)];
这个代码简单,先说模型:
--------------------------------------------------------------
function y=LV_model(t,x)
a=0.01;b=0.02;
y=[x(1)-a*x(1)*x(2);-x(2)+b*x(1)*x(2)];
--------------------------------------------------------
x(1)表示兔子的当前个数,x(2)表示狐狸的当前个数.
y(1)=x(1)-a*x(1)*x(2) 表示兔子增长的变化趋势,y(2)=-x(2)+b*x(1)*x(2)表示狐狸增长的变化趋势.a,b分别表示变化速率系数:a表示随着t的增加,兔子的变化速率是与当前兔子个数和兔子狐狸个数乘积的a系数相关;b系数同理.
再说主函数
-------------------------------------------------------
function Lotka_Volterra
t_end=15; %[0,15]t的变化区间
x0=[20;20]; %初始兔子和狐狸的个数
[t,x]=ode45(@LV_model,[0,t_end],x0) %通过ode45(常微分方程数值求解)求解模型,得到t和x的相应变化情况,即在t=0时,兔子和狐狸的个数是多少,t=0.1时,兔子和狐狸的个数是多少...
plot(t,x); %绘制t兔子和狐狸变化对应关系图
legend('野兔','狐狸')%给不同曲线标上所代表的是什么变量的变化,x(1),x(2)分别代表兔子和狐狸
---------------------------------------------------------