跪求matlab大神指导!>> function f=sirfun(t,x)lmda=1;mu=0.3;f(1)=lmda*x(1)*x(2)-mu*x(1);f(2)=-lmda*x(1)*x(2);f=f(:);%作解曲线图ts=[0,40];x0=[0.02 0.98];[t,x]=ode45(@sirfun,ts,x0);s=x(:,2);i=x(:,1);r=1-x(:,1)-x(:,2);plot(t,s,t,i,t,r);legend('
来源:学生作业帮助网 编辑:作业帮 时间:2024/06/21 15:02:45
![跪求matlab大神指导!>> function f=sirfun(t,x)lmda=1;mu=0.3;f(1)=lmda*x(1)*x(2)-mu*x(1);f(2)=-lmda*x(1)*x(2);f=f(:);%作解曲线图ts=[0,40];x0=[0.02 0.98];[t,x]=ode45(@sirfun,ts,x0);s=x(:,2);i=x(:,1);r=1-x(:,1)-x(:,2);plot(t,s,t,i,t,r);legend('](/uploads/image/z/9905613-69-3.jpg?t=%E8%B7%AA%E6%B1%82matlab%E5%A4%A7%E7%A5%9E%E6%8C%87%E5%AF%BC%21%3E%3E+function+f%3Dsirfun%28t%2Cx%29lmda%3D1%3Bmu%3D0.3%3Bf%281%29%3Dlmda%2Ax%281%29%2Ax%282%29-mu%2Ax%281%29%3Bf%282%29%3D-lmda%2Ax%281%29%2Ax%282%29%3Bf%3Df%28%3A%29%3B%25%E4%BD%9C%E8%A7%A3%E6%9B%B2%E7%BA%BF%E5%9B%BEts%3D%5B0%2C40%5D%3Bx0%3D%5B0.02+0.98%5D%3B%5Bt%2Cx%5D%3Dode45%28%40sirfun%2Cts%2Cx0%29%3Bs%3Dx%28%3A%2C2%29%3Bi%3Dx%28%3A%2C1%29%3Br%3D1-x%28%3A%2C1%29-x%28%3A%2C2%29%3Bplot%28t%2Cs%2Ct%2Ci%2Ct%2Cr%29%3Blegend%28%27)
跪求matlab大神指导!>> function f=sirfun(t,x)lmda=1;mu=0.3;f(1)=lmda*x(1)*x(2)-mu*x(1);f(2)=-lmda*x(1)*x(2);f=f(:);%作解曲线图ts=[0,40];x0=[0.02 0.98];[t,x]=ode45(@sirfun,ts,x0);s=x(:,2);i=x(:,1);r=1-x(:,1)-x(:,2);plot(t,s,t,i,t,r);legend('
跪求matlab大神指导!
>> function f=sirfun(t,x)
lmda=1;mu=0.3;
f(1)=lmda*x(1)*x(2)-mu*x(1);
f(2)=-lmda*x(1)*x(2);
f=f(:);
%作解曲线图
ts=[0,40];x0=[0.02 0.98];
[t,x]=ode45(@sirfun,ts,x0);
s=x(:,2);i=x(:,1);
r=1-x(:,1)-x(:,2);
plot(t,s,t,i,t,r);
legend('易感染者S','已感染者I','移出者R')
程序报错
function f=sirfun(t,x)
|
Error:Function definitions are not permitted at the prompt or in scripts.
请问是什么错误啊?应该怎么修改?
跪求matlab大神指导!>> function f=sirfun(t,x)lmda=1;mu=0.3;f(1)=lmda*x(1)*x(2)-mu*x(1);f(2)=-lmda*x(1)*x(2);f=f(:);%作解曲线图ts=[0,40];x0=[0.02 0.98];[t,x]=ode45(@sirfun,ts,x0);s=x(:,2);i=x(:,1);r=1-x(:,1)-x(:,2);plot(t,s,t,i,t,r);legend('
function f=sirfun(t,x)
lmda=1;mu=0.3;
f(1)=lmda*x(1)*x(2)-mu*x(1);
f(2)=-lmda*x(1)*x(2);
f=f(:);
要保存到 sirfun.m 文件中.
再在命令窗口中输入
>>%作解曲线图
ts=[0,40];x0=[0.02 0.98];
[t,x]=ode45(@sirfun,ts,x0);
s=x(:,2);i=x(:,1);
r=1-x(:,1)-x(:,2);
plot(t,s,t,i,t,r);
legend('易感染者S','已感染者I','移出者R')
就好了.
看点入门书即可.