四阶龙格库塔算法计算一阶常微分方程组MATLAB

四阶龙格库塔算法是一种常用的数值解常微分方程组的方法。下面是使用MATLAB实现四阶龙格库塔算法计算一阶常微分方程组的代码：

function \[x, Y\] = runge_kutta(fun1, x0, y0, h, n)
% fun1是一个函数句柄，表示常微分方程组的右侧函数
% x0是初始点的x坐标
% y0是初始点的y坐标
% h是步长
% n是迭代次数

x = zeros(1, n+1);
Y = zeros(3, n+1);
x(1) = x0;
Y(:, 1) = y0;

for i = 1:n
    K1 = fun1(x(i), Y(:, i));
    K2 = fun1(x(i) + 1/2*h, Y(:, i) + 1/2*K1*h);
    K3 = fun1(x(i) + 1/2*h, Y(:, i) + 1/2*K2*h);
    K4 = fun1(x(i) + h, Y(:, i) + K3*h);
    
    x(i+1) = x(i) + h;
    Y(:, i+1) = Y(:, i) + 1/6*(K1 + 2*K2 + 2*K3 + K4)*h;
end

end

在这段代码中，`fun1`是一个函数句柄，表示常微分方程组的右侧函数。`x0`和`y0`是初始点的坐标，`h`是步长，`n`是迭代次数。函数返回计算得到的`x`和`Y`。

你可以根据自己的常微分方程组的具体形式，修改`fun1`函数的定义，然后调用`runge_kutta`函数进行计算。
#### 引用[.reference_title]
- *1* *2* [4 详解matlab实现龙格库塔算法求解复杂常微分方程组](https://blog.csdn.net/Linhua009900/article/details/127376412)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^insert_down28v1,239^v3^insert_chatgpt"}} ] [.reference_item]
[ .reference_list ]