XIONGBAOBAOYA 2022-03-08 16:12 采纳率: 100%
浏览 66
已结题

用ODE45求解微分方程组出现这样的报错应该怎么修改

img

毕业论文卡在这一步了,有人指点的话,万分感谢。
大概就是一个耦合方程,好像是语句问题,我第一次学习ODE45,不太清楚应该怎么修改。

  • 写回答

2条回答 默认 最新

  • 技术专家团-Joel 2022-03-08 16:37
    关注

    你的L是向量,这是不允许的,你可以把L定义成和时间相关的量,或者你的L就是一个数,这样的ODE才是可以求解的,比方说

    L = @(t) t;
    K0 = @(t) pi/2*L(t);
    K1 = @(t)2*sqrt(2) *K0(t);
    K2 = K1;
    f = @(t,x)[i*K1(t)*x(2); i*K1(t)*x(1)+i*K2(t)*x(3); i*K2(t)*x(2)];
    [t,y] = ode45(f, [0,1], [1;0;0]);
    plot(t,y)
    
    

    或者你就让L为常数1

    
    L = 1;
    K0 = pi/2*L;
    K1 = 2*sqrt(2) *K0;
    K2 = K1;
    f = @(t,x)[i*K1*x(2); i*K1*x(1)+i*K2*x(3); i*K2*x(2)];
    [t,y] = ode45(f, [0,1], [1;0;0]);
    plot(t,y)
    
    
    本回答被题主选为最佳回答 , 对您是否有帮助呢?
    评论
查看更多回答(1条)

报告相同问题?

问题事件

  • 系统已结题 3月19日
  • 已采纳回答 3月11日
  • 创建了问题 3月8日

悬赏问题

  • ¥15 visionmaster启动失败,提示为“机器不满足授权而被禁用”
  • ¥50 用logisim设计16位单时钟周期cpu
  • ¥15 IDEA中圈复杂度如何具体设置
  • ¥50 labview采集不了数据
  • ¥15 请上面代码做什么处理或什么混淆
  • ¥15 英雄联盟自定义房间置顶
  • ¥15 W5500网线插上无反应
  • ¥15 如何用字典的Key,显示在WPF的xaml中
  • ¥15 weautomate读取Excel表格信息然后填写到网页一直报错,如何解决?
  • ¥15 C#如何在Webview2中获取网页验证码