图像缩略图

龙格库塔四阶颂歌

版本1.4.0.0(1.32 KB)通过 犹大小号
解决了使用第四阶Runge Kutta法ODE

421个下载

更新2013年1月16日

查看许可证

该代码定义了现有功能和步长,其可以按要求改变。

P.S:此代码具有比现有可用的代码网上新功能。背后张贴这个非常简单的代码,目的是帮助学生理解概念和解决分配。

引用作为

犹太S(2020)。龙格库塔四阶颂歌(//www.tianjin-qmedu.com/matlabcentral/fileexchange/29851-runge-kutta-4th-order-ode),MATLAB中央文件交换。检索

评论和评分(21

Soumyadeep

它不为Y(0)= 0的初始条件工作

尼迪梅农

亚伯拉罕

做得好!什么有关龙格库塔方法二阶常微分方程代码。这种性质的东西:
d ^ 2Y / DX ^ 2 + 0.6 * DY / DX 0.8y = 0

谢谢

math16

Sumith YD

Dogba Djaze

我怎么可以用RK4方法在MATLAB解决SIR模型?您可以编写代码,请

先生,你能帮助我,那我们如何能运用四阶龙格库塔方法4耦合方程?
DX / DT = -ax - EΩ+ YZ
DY / DT =由+ XZ
DZ / DT = CZ +Fω - XY
Dω/ DT =Dω - GZ
一个= 50,B = -16,C = 10,d = 0.2,E = 10,F = 16,G = 0.5
步长0.001。
问候
法伊兹

@Shahzaib阿西夫非常有用program.JazakAllah

我怎么跑到这个代码/电话?

此函数:f '' '' - F * F ''” + 4 * G = 0
在这里我需要将其插入到该代码?
谢谢

函数RK4(F,A,X0,Y0,h)的

龙格%库塔法4阶
%函数@(X,Y),例如F = @(X,Y)(X + Y);
%A =点多达你得到的结果
%X0 = x的初始条件
%Y0 = Y的初始条件
% 一步的大小

X = X0:H:一个;
Y(1)= Y0;

对于i = 1:(长度(X)-1)

K1 = F(X(i)中,Y(I));
K 2 = F(X(I)+ 0.5 * H,Y(I)+ 0.5 * H * K1);
K3 = F((X(I)+ 0.5 * h)中,(Y(I)+ 0.5 * H * K2));
K4 = F((X(I)+ h)时,(Y(I)+ K3 * H));

Y(I + 1)= Y(I)+(1/6)*(K1 + 2 * K2 + 2 * K3 + K4)*小时;

结束

Y(:)

%Shahzaib阿西(zaibi7402)
%shahzaib.7402@gmail.com

克里斯FUNG

清晰的编码

很好的学习。谢谢。

阿伦

裨踢嗯

优秀作品

我做

优秀的程序,
非常有帮助。

更新

1.4.0.0

刚刚更新

1.1.0.0

没什么

MATLAB版本兼容性
创建了R2010a版本
与任何版本兼容
平台兼容性
视窗 苹果系统 Linux的