Matlab四阶龙格库塔法求解常微分方程

Question

用Matlab……求 谢谢各位大神了

Answer 1

用Matlab四阶龙格库塔法求常微分方程可以按照以下方法去实现。

1、首先建立自定义微分方程函数

function f = ode_fun(x,y)

f=y+2*x/y^2;

end

2、然后用四阶龙格库塔法求其数值解

figure(2)

y0=[1]; %初值y（0）=1

h=0.1;

a=0;

b=5;

[x,y] = runge_kutta(@(x,y)ode_fun(x,y),y0,h,a,b);

disp('        x        y')

A=[x',y']

plot(x,y,'LineWidth',1.5),grid on

xlabel('x'),ylabel('y(x)');

3、绘制y-x的曲线图