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
本回答被题主选为最佳回答 , 对您是否有帮助呢?解决评论 打赏 举报无用 1
悬赏问题
- ¥30 Matlab打开默认名称带有/的光谱数据
- ¥50 easyExcel模板 动态单元格合并列
- ¥15 res.rows如何取值使用
- ¥15 在odoo17开发环境中,怎么实现库存管理系统,或独立模块设计与AGV小车对接?开发方面应如何设计和开发?请详细解释MES或WMS在与AGV小车对接时需完成的设计和开发
- ¥15 CSP算法实现EEG特征提取,哪一步错了?
- ¥15 游戏盾如何溯源服务器真实ip?需要30个字。后面的字是凑数的
- ¥15 vue3前端取消收藏的不会引用collectId
- ¥15 delphi7 HMAC_SHA256方式加密
- ¥15 关于#qt#的问题:我想实现qcustomplot完成坐标轴
- ¥15 下列c语言代码为何输出了多余的空格