EDA365电子论坛网

标题: 微分方程求解 [打印本页]
作者: IBB-EUT    时间: 2020-8-13 11:09
标题: 微分方程求解
我最近在做四阶龙格库塔算法的微分方程求解，求大神分享代码，感谢！
3 y, X: g. p) f% Q4 I% l' W2 y/ P
作者: Uifhjvv    时间: 2020-8-13 13:10
function yout = runge_kutta(fun,t,yin,varargin), ?- t% v% g4 Y9 I. ~
t_interval = t(2)-t(1);) b" Q) ~# U& b5 |2 Z+ E
k1 = t_interval*fun(t(1),yin,varargin);
; i, j2 [, @( ~- O7 x: B) e: L0 o6 Fk2 = t_interval*fun(t(1)+t_interval/2,yin+k1/2varargin);0 |6 Z  n6 Y/ v9 l# z- R
k3 = t_interval*fun(t(1)+t_interval/2,yin+k2/2,varargin);. g) \3 |* n" d  U! N
k4 = t_interval*fun(t(1)+t_interval,yin+k3,varargin);) i# w: @9 z2 {$ m# l9 w, j+ s
yout = yin+(k1+2*k2+2*k3+k4)/6;( o+ z. b6 v" C+ Y7 ^+ f, W
end
作者: IBB-EUT    时间: 2020-8-14 13:50
Uifhjvv 发表于 2020-8-13 13:10
. r) o7 x/ s- p- G' F; ^function yout = runge_kutta(fun,t,yin,varargin)
4 v0 S: m& i/ e4 H9 v1 xt_interval = t(2)-t(1);3 s+ F* J- h) U0 a! G9 `& J) d
k1 = t_interval*fun(t(1),y ...
9 t  W+ G! ?+ B- w! G) X; \# [
谢谢大佬，我先去试试
$ m. g/ w7 i- S/ ~* R( @+ W! A( `
作者: shuddkk    时间: 2020-8-17 10:40
楼上正解，，
作者: cichishia    时间: 2020-8-18 10:15




欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/) Powered by Discuz! X3.2