基于MATLAB实现固定步长的四阶龙格库塔方法求解常微分方程初值问题

最新推荐文章于 2024-08-21 11:02:45 发布

原创

最新推荐文章于 2024-08-21 11:02:45 发布 · 4.3k 阅读

35 ·

CC 4.0 BY-SA版权

文章标签：

#matlab

本文介绍使用四阶龙格-库塔(RK)方法求解特定微分方程的MATLAB代码实现，并通过实例展示如何计算一暂态电路问题的数值解。同时，给出了该问题的精确解析解，并比较了数值解与解析解的差异。

代码

构造函数

function [yn,xn] = RKorder4(f,a,b,h,y0)
%a，b为上下界；h为步长；y0为初值。
%返回yn、xn向量，并绘制曲线
yn=zeros(1,((b-a)/h)+1);
yn(1)=y0;
xn=a:h:b;
for i=1:((b-a)/h)
    k1=h*f(xn(i),yn(i));
    k2=h*f(xn(i)+h/2,yn(i)+k1/2);
    k3=h*f(xn(i)+h/2,yn(i)+k2/2);
    k4=h*f(