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
悬赏问题
- ¥15 装 pytorch 的时候出了好多问题,遇到这种情况怎么处理?
- ¥20 IOS游览器某宝手机网页版自动立即购买JavaScript脚本
- ¥15 手机接入宽带网线,如何释放宽带全部速度
- ¥30 关于#r语言#的问题:如何对R语言中mfgarch包中构建的garch-midas模型进行样本内长期波动率预测和样本外长期波动率预测
- ¥15 ETLCloud 处理json多层级问题
- ¥15 matlab中使用gurobi时报错
- ¥15 这个主板怎么能扩出一两个sata口
- ¥15 不是,这到底错哪儿了😭
- ¥15 2020长安杯与连接网探
- ¥15 关于#matlab#的问题:在模糊控制器中选出线路信息,在simulink中根据线路信息生成速度时间目标曲线(初速度为20m/s,15秒后减为0的速度时间图像)我想问线路信息是什么