Question

求解一个MATLAB程序及画图(如何在MATLAB中写出求解程序) 该题初始条件N(0)=3000,E(0)=5,A(0)=100，方程：

该题的微分方程组为：
dN/dt=0.885N(1-N/(1.25*10^9))-4*10(-3)*E*N
dE/dt=5*A-2*E
dA/dt=4.8*10^(-7)*N-0.274*A
该图像应该是个增函数然后再一定时间后稳定，中间似乎有一段是阻尼震荡，可是我的图像一直不对，希望大神可以把源程序贴一下，这个是模拟生物系统的函数

Answer 1

参考代码：

ode = inline('[0.885*y(1)*(1-y(1)/(1.25*10^9))-4*10^(-3)*y(2)*y(1); 5*y(3)-2*y(2); 4.8*10^(-7)*y(1)-0.274*y(3)]','t','y');
ode45(ode,[0 100],[3000 5 100])

追问

能把所有的程序发下吗，加上plot

追答

这个代码已包含画图的部分，ode* 系列函数不带参数输出的时候会自动绘图。