su_shi1 2023-02-18 11:14 采纳率: 100%
浏览 57
已结题

关于#matlab#的问题:matlab微分方程模型,说我索引超出数组元素的数目(2),程序无法执行

matlab微分方程模型,说我索引超出数组元素的数目(2)。求大家帮我看看我哪里写错了。

img

img

  • 写回答

5条回答 默认 最新

  • A1358722560 2023-02-18 14:40
    关注

    输入变量必须是3维的,但是你的是二维的,维数不对
    t与y都必须是三维的才行

    clc;
    
    R=[0.1,0.2,0.3];
    A=[0,0.05,0.08; 0.1,0,0.18;0.06,0.12,0];
    B=[2,1.8,1.4];
    P0=0.8;
    D=[0.02,0.01,0.05];
    [t, y]=ode45(@(t,y)eq3(t,y,R,A,B,P0,D),[0,12,15],[100,100,100]);
    
    function dy=eq3(t,y,R,A,B,P0,D)%说明微分变里是3维的
        dy=zeros(3,1);
        %微分方程组
        dy(1)=R(1,1)*y(1)*(1-A(2,1)*y(2)-A(3,1)*y(3))-B(1,1)*P0*(1+sin(pi/6*t))-D(1,1)*y(1);
        dy(2)=R(1,2)*y(2)*(1-A(1,2)*y(1)-A(3,2)*y(3))-B(1,2)*P0*(1+sin(pi/6*t))-D(1,2)*y(2);
        dy(3)=R(1,3)*y(3)*(1-A(1,3)*y(1)-A(2,3)*y(2))-B(1,3)*P0*(1+sin(pi/6*t))-D(1,3)*y(3);
    end
    
    
    

    如果问题得到解决的话请点 采纳~~

    本回答被题主选为最佳回答 , 对您是否有帮助呢?
    评论 编辑记录
查看更多回答(4条)

报告相同问题?

问题事件

  • 系统已结题 2月26日
  • 已采纳回答 2月18日
  • 赞助了问题酬金15元 2月18日
  • 创建了问题 2月18日

悬赏问题

  • ¥50 matlab,文献复现
  • ¥20 寻找Flink CDC指导,有偿咨询
  • ¥15 cruise如何仿真不同温度下车辆的经济性?
  • ¥20 扑克牌游戏的完整代码
  • ¥88 求帮安装Python架构的vnpy
  • ¥20 逻辑回归法分析品牌知名度和购买意向
  • ¥15 华为云主机-异常连接
  • ¥50 vue 成绩列表 输入成绩,并统计各分数段人数
  • ¥20 请教怎么绕过这个卡密系统 大家交流一下
  • ¥15 设计数字电路得到要求结果