微分方程求解

IBB-EUT · 发表于 2020-8-13 11:09

EDA365欢迎您登录！

您需要登录才可以下载或查看，没有帐号？注册

x

我最近在做四阶龙格库塔算法的微分方程求解，求大神分享代码，感谢！

Uifhjvv · 发表于 2020-8-13 13:10

function yout = runge_kutta(fun,t,yin,varargin)
t_interval = t(2)-t(1);
k1 = t_interval*fun(t(1),yin,varargin);
k2 = t_interval*fun(t(1)+t_interval/2,yin+k1/2varargin);
k3 = t_interval*fun(t(1)+t_interval/2,yin+k2/2,varargin);
k4 = t_interval*fun(t(1)+t_interval,yin+k3,varargin);
yout = yin+(k1+2*k2+2*k3+k4)/6;
end

IBB-EUT · 发表于 2020-8-14 13:50

Uifhjvv 发表于 2020-8-13 13:10
6 ^. R6 k- ~" h% {/ ^2 k- pfunction yout = runge_kutta(fun,t,yin,varargin)
& n5 ]7 t3 y% gt_interval = t(2)-t(1);( T7 _5 p5 I9 |' h
k1 = t_interval*fun(t(1),y ...

谢谢大佬，我先去试试

shuddkk · 发表于 2020-8-17 10:40

楼上正解，，

cichishia · 发表于 2020-8-18 10:15

帐号		自动登录	找回密码
密码			注册

微分方程求解

EDA365欢迎您登录！

点评

浏览过的版块

推荐内容 /1