m0_62959064 2021-10-22 09:13 采纳率: 62.5%
浏览 194
已结题

编写MATLAB函数,用欧拉法求解微分方程

img

  • 写回答

1条回答 默认 最新

  • joel_1993 2021-10-22 10:30
    关注

    你好,代码供参考:

    function eulerMain %主函数
    
    odefun = @(t,y) [y^2+t];% 构建微分方程
    y0 = [1]; % 设置初值y(0)=1
    tspan = [0, 0.8];% 微分方程求解时间跨度
    [t,y] = euler(odefun, tspan, y0);%利用欧拉法求odefun
    plot(t,y)%画图
    xlabel('t')
    ylabel('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月30日
  • 已采纳回答 10月22日
  • 创建了问题 10月22日

悬赏问题

  • ¥100 连续两帧图像高速减法
  • ¥15 组策略中的计算机配置策略无法下发
  • ¥15 如何绘制动力学系统的相图
  • ¥15 对接wps接口实现获取元数据
  • ¥20 给自己本科IT专业毕业的妹m找个实习工作
  • ¥15 用友U8:向一个无法连接的网络尝试了一个套接字操作,如何解决?
  • ¥30 我的代码按理说完成了模型的搭建、训练、验证测试等工作(标签-网络|关键词-变化检测)
  • ¥50 mac mini外接显示器 画质字体模糊
  • ¥15 TLS1.2协议通信解密
  • ¥40 图书信息管理系统程序编写