Question

用MATLAB对微分方程组进行仿真

微分方程组为dS/dt=2.2831- 0.00088*S*I-0.002228*S
dE/dt=0.00088*S*I-1.33*E
dI/dt=0.7*E-0.600228*I
最好将最后的那个图形搞出来，发到我邮箱940672887@qq.com

Answer 1

clc
f=@(t,x)([2.2831-.00088*x(1)*x(3)-.002228*x(1);.00088*x(1)*x(3)-1.33*x(2);.7*x(2)-.600228*x(3)]);
ode45(f,[0 1],[1 1 1])

Answer 2

[S,E,I]=dsolve('DS=2.2831-0.00088*S*I-0.002228*S','DE=0.00088*S*I-1.33*E','DI=0.7*E-0.600228*I')

追问

我主要是想要这个微分方程组的仿真图形，最好用matlab仿真出来

追答

建立函数：
function xdot=lorenz(t,x)
xdot=[-0.00088*x(2)-0.002228,0,0;0,0.00088*x(1),-1.33;0,-0.600228,0.7]*x+[2.2831;0;0];

命令窗口：
x0=[1 1 1]';
[t,x]=ode45('lorenz',[1,10],x0);

好像解不出来

追问

好的谢谢 我去试试

Answer 3

红人堂还好