跪求matlab大神指导!错在哪里 如何修改啊?>> figure;for s0=[0.9 0.7 0.5 0.3 0.1] x0=[1-s0,s0]; [t,x]=ode45(@sirfun,ts,x0); s=x(:,2);i=x(:,1); plot(s,i);hold on;endxlabel('S');ylabel('R')axis([0 1 0 1])plot([0 1],[1,0],'r')tit
来源:学生作业帮助网 编辑:作业帮 时间:2024/07/12 02:15:51
![跪求matlab大神指导!错在哪里 如何修改啊?>> figure;for s0=[0.9 0.7 0.5 0.3 0.1] x0=[1-s0,s0]; [t,x]=ode45(@sirfun,ts,x0); s=x(:,2);i=x(:,1); plot(s,i);hold on;endxlabel('S');ylabel('R')axis([0 1 0 1])plot([0 1],[1,0],'r')tit](/uploads/image/z/9905614-70-4.jpg?t=%E8%B7%AA%E6%B1%82matlab%E5%A4%A7%E7%A5%9E%E6%8C%87%E5%AF%BC%21%E9%94%99%E5%9C%A8%E5%93%AA%E9%87%8C+%E5%A6%82%E4%BD%95%E4%BF%AE%E6%94%B9%E5%95%8A%3F%3E%3E+figure%3Bfor+s0%3D%5B0.9+0.7+0.5+0.3+0.1%5D++++x0%3D%5B1-s0%2Cs0%5D%3B++++%5Bt%2Cx%5D%3Dode45%28%40sirfun%2Cts%2Cx0%29%3B+++++s%3Dx%28%3A%2C2%29%3Bi%3Dx%28%3A%2C1%29%3B++++plot%28s%2Ci%29%3Bhold+on%3Bendxlabel%28%27S%27%29%3Bylabel%28%27R%27%29axis%28%5B0+1+0+1%5D%29plot%28%5B0+1%5D%2C%5B1%2C0%5D%2C%27r%27%29tit)
跪求matlab大神指导!错在哪里 如何修改啊?>> figure;for s0=[0.9 0.7 0.5 0.3 0.1] x0=[1-s0,s0]; [t,x]=ode45(@sirfun,ts,x0); s=x(:,2);i=x(:,1); plot(s,i);hold on;endxlabel('S');ylabel('R')axis([0 1 0 1])plot([0 1],[1,0],'r')tit
跪求matlab大神指导!错在哪里 如何修改啊?
>> figure;
for s0=[0.9 0.7 0.5 0.3 0.1]
x0=[1-s0,s0];
[t,x]=ode45(@sirfun,ts,x0);
s=x(:,2);i=x(:,1);
plot(s,i);hold on;
end
xlabel('S');ylabel('R')
axis([0 1 0 1])
plot([0 1],[1,0],'r')
title('相平面图')
hold off;
程序报错
? Error using ==> feval
Undefined command/function 'sirfun'.
Error in ==> funfun\private\odearguments at 110
f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0.
Error in ==> ode45 at 173
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn,
请问错在哪里如何修改?
跪求matlab大神指导!错在哪里 如何修改啊?>> figure;for s0=[0.9 0.7 0.5 0.3 0.1] x0=[1-s0,s0]; [t,x]=ode45(@sirfun,ts,x0); s=x(:,2);i=x(:,1); plot(s,i);hold on;endxlabel('S');ylabel('R')axis([0 1 0 1])plot([0 1],[1,0],'r')tit
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 文件中.
再在命令窗口中输入
>>figure;
for s0=[0.9 0.7 0.5 0.3 0.1]
x0=[1-s0,s0];
[t,x]=ode45(@sirfun,ts,x0);
s=x(:,2);i=x(:,1);
plot(s,i);hold on;
end
xlabel('S');ylabel('R')
axis([0 1 0 1])
plot([0 1],[1,0],'r')
title('相平面图')
hold off;
就好了.
看点入门书即可.