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

Question

在没有人类介入的生物圈中，只考虑野兔和狐狸这两种动物的生存和发展。假设这一环境中有足够多的食物供野兔生存，而狐狸只是以野兔为食物。方程中x,y项反映了两生物种群的相互作用对于某一种数量变化的影响。因为狐狸的存在对野兔不利，所以-a<0 ；而野兔的存在对狐狸是有利的，所以b>0 。
取参数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)];

Answer 1

这个代码简单，先说模型：
--------------------------------------------------------------
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)分别代表兔子和狐狸
---------------------------------------------------------