m0_62959064 2021-10-14 10:44 采纳率: 62.5%
浏览 500
已结题

MATLAB欧拉法求解微分方程

img

  • 写回答

1条回答 默认 最新

  • joel_1993 2021-10-14 11:07
    关注

    你好!代码供参考

    function main %主函数
    odefun = @(t,y) y^2 + t;
    y0 = 1;
    tspan = [0, 0.5];
    [t,y] = euler(odefun, tspan, y0);
    plot(t,y)
    end
    function [t, y] = euler(odefun, tspan, y0)
    % odefun ode函数
    % tspan: 求解时间范围
    % x0:初值
    t = tspan(:);
    if(numel(tspan)==2)
        t = linspace(tspan(1), tspan(2), 1001)';
    end
    dy0 = odefun(t(1), y0);
    y = zeros(numel(t), numel(dy0));
    y(1,:) = y0(:);
    
    for i = 2:1:numel(t)
        y(i,:) = y(i-1,:) + odefun(t(i-1),y(i-1,:))*(t(i)-t(i-1));
    end
    end
    
    

    img

    本回答被题主选为最佳回答 , 对您是否有帮助呢?
    评论

报告相同问题?

问题事件

  • 系统已结题 10月23日
  • 已采纳回答 10月15日
  • 创建了问题 10月14日

悬赏问题

  • ¥15 python使用pulp线性优化时报错
  • ¥15 为什么我的uibot导入py模块出错呀。py文件放在了uibot里对应的python文件夹了,卸了重安也不行
  • ¥15 开源或低价数据中台哪个最好
  • ¥15 arduino编程出现字符串疑似覆盖现象
  • ¥15 我的b站在没有碰到屏幕的情况下偶尔会自动跳出进度条,就像在屏幕上点了一下一样,但我并没有点。而且视频进度并没有变。这可能是什么原因造成的?
  • ¥30 STK matlab python仿真
  • ¥15 关于IMageEnView 图标定位问题
  • ¥20 求解答(matlab)
  • ¥30 ffmpeg库使用过程中遇到的问题
  • ¥15 pyqt5 中python如何通过Qtwebchannel主动发消息给web前端