用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
知道答主
回答量:0
采纳率:0%
帮助的人:0
展开全部

 图片:


程序:

 

    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-04-02 广告
正弦振动多用于找出产品设计或包装设计的脆弱点。看在哪一个具体频率点响应最大(共振点);正弦振动在任一瞬间只包含一种频率的振动,而随机振动在任一瞬间包含频谱范围内的各种频率的振动。由于随机振动包含频谱内所有的频率,所以样品上的共振点会同时激发... 点击进入详情页
本回答由富港检测技术(东莞)有限公司_提供
推荐律师服务: 若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询
?>

为你推荐:

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

类别

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

说明

0/200

提交
取消

辅 助

模 式