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

用Matlab……求谢谢各位大神了... 用Matlab……求 谢谢各位大神了 展开
 我来答
lhmhz
高粉答主

2018-12-29 · 专注matlab等在各领域中的应用。
lhmhz
采纳数:7264 获赞数:17014

向TA提问 私信TA
展开全部

用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的曲线图

推荐律师服务: 若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询

为你推荐:

下载百度知道APP,抢鲜体验
使用百度知道APP,立即抢鲜体验。你的手机镜头里或许有别人想知道的答案。
扫描二维码下载
×

类别

我们会通过消息、邮箱等方式尽快将举报结果通知您。

说明

0/200

提交
取消

辅 助

模 式