matlab大神,求解微分方程组代码
来源:学生作业帮 编辑:搜搜考试网作业帮 分类:综合作业 时间:2024/07/28 08:44:09
matlab大神,求解微分方程组代码
求代码
![matlab大神,求解微分方程组代码](/uploads/image/z/4353308-44-8.jpg?t=matlab%E5%A4%A7%E7%A5%9E%2C%E6%B1%82%E8%A7%A3%E5%BE%AE%E5%88%86%E6%96%B9%E7%A8%8B%E7%BB%84%E4%BB%A3%E7%A0%81)
lumda=1;
mu=0.3;
fun=@(t,is)[lumda*is(2)*is(1)-mu*is(1);-lumda*is(2)*is(1)];
[T,Y]=ode45(fun,[0 100],[0.02 0.98]); %时间区间个可以选择
subplot(211),plot(T,Y);
xlabel('t');
legend('i','s');
subplot(212),plot(Y(:,2),Y(:,1));
xlabel('s');
ylabel('i');
![](http://img.wesiedu.com/upload/c/77/c7798b317b13365de787e5eef8a153f6.jpg)
mu=0.3;
fun=@(t,is)[lumda*is(2)*is(1)-mu*is(1);-lumda*is(2)*is(1)];
[T,Y]=ode45(fun,[0 100],[0.02 0.98]); %时间区间个可以选择
subplot(211),plot(T,Y);
xlabel('t');
legend('i','s');
subplot(212),plot(Y(:,2),Y(:,1));
xlabel('s');
ylabel('i');
![](http://img.wesiedu.com/upload/c/77/c7798b317b13365de787e5eef8a153f6.jpg)