用matlab画微分方程组图像

Ds=-0.5*s*i+0.5*i-0.1*s+0.1;De=1.5*s*i-0.1*e+0.9*e*i;Di=i^2+0.1*e*i-1.6*i。都是对t求导这个微分方... Ds=-0.5*s*i+0.5*i-0.1*s+0.1;De=1.5*s*i-0.1*e+0.9*e*i;Di=i^2+0.1*e*i-1.6*i。
都是对t求导这个微分方程组在matlab里面怎么绘图?初值的话取(0.5,0,0.5)和(0,0,1)就行。求大神!!!
展开
lccever
2014-05-11 · 超过39用户采纳过TA的回答
知道小有建树答主
回答量:67
采纳率:0%
帮助的人:53.3万
展开全部

 图片:


程序:

 

    tspan=[0 20];

    y0=[0.5 0 0.5 ];

    y1=[0 0 1];

    [t0,Y0]=ode45('dif_dsdedi',tspan,y0)

    [t1,Y1]=ode45('dif_dsdedi',tspan,y1)

    figure

    subplot(2,1,1)

    plot(t0,Y0(:,1),'-b',t0,Y0(:,2),'-g',t0,Y0(:,3),'-k')

    grid on

    legend('s','e','i')

    xlabel('t-time')

    title('y0=[0.5 0 0.5 ]')

    subplot(2,1,2)

    plot(t1,Y1(:,1),'-b',t1,Y1(:,2),'-g',t1,Y1(:,3),'-k')

    legend('s','e','i')

    grid on

    xlabel('t-time')

    title('y0=[0 0 1]')


%=====子函数

function dy=dif_dsdedi(t,y)

    %y(1)=s y(2)=e y(3)=i 

    dy(1)=-0.5*y(1)*y(3)+0.5*y(3)-0.1*y(1)+0.1;

    dy(2)=1.5*y(1)*y(3)-0.1*y(2)+0.9*y(2)*y(3);

    dy(3)=y(3)^2+0.1*y(2)*y(3)-1.6*y(3);

    dy=dy(:)

东莞大凡
2024-11-14 广告
标定板认准大凡光学科技,专业生产研发厂家,专业从事光学影像测量仪,光学投影测量仪.光学三维测量仪,光学二维测量仪,光学二维测量仪,光学三维测量仪,光学二维测量仪.的研发生产销售。东莞市大凡光学科技有限公司创立于 2018 年,公司总部坐落于... 点击进入详情页
本回答由东莞大凡提供
推荐律师服务: 若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询

为你推荐:

下载百度知道APP,抢鲜体验
使用百度知道APP,立即抢鲜体验。你的手机镜头里或许有别人想知道的答案。
扫描二维码下载
×

类别

我们会通过消息、邮箱等方式尽快将举报结果通知您。

说明

0/200

提交
取消

辅 助

模 式