非刚性常微分方程的数值解法通常会用四阶龙格库塔算法，其matlab函数对应ode45。 对于dy/dx = f(x,y)，y(0)=y0。 其四阶龙格库塔公式如下： 对于通常计算，四阶已经够用，四阶以 … Continue reading MATLAB 龙格库塔法
非刚性常微分方程的数值解法通常会用四阶龙格库塔算法，其matlab函数对应ode45。 对于dy/dx = f(x,y)，y(0)=y0。 其四阶龙格库塔公式如下： 对于通常计算，四阶已经够用，四阶以 … Continue reading MATLAB 龙格库塔法