• 妲己爱主人
    了解作者
  • matlab
    开发工具
  • 29KB
    文件大小
  • rar
    文件格式
  • 0
    收藏次数
  • 1 积分
    下载积分
  • 3
    下载次数
  • 2021-01-24 11:28
    上传日期
龙格-库塔(Runge-Kutta)方法是一种在工程上应用广泛的高精度单步算法。见常微分方程初值问题数值解法。
四阶龙格-库塔算法.rar
  • 四阶龙格-库塔算法
  • Adams.jpg
    33.6KB
  • RF4.jpg
    16.1KB
  • rkftry.m
    766B
  • rk4.m
    585B
内容介绍
function rkf4 clc; clear all; close all; t0=0; t1=1; n=10000; ht=(t1-t0)/n; t=t0:ht:t1; w=zeros(1,n+1); for i=2:n+1 k1=ht*myf(t(i-1),w(i-1)); k2=ht*myf(t(i-1)+ht/4,w(i-1)+k1/4); k3=ht*myf(t(i-1)+3*ht/8,w(i-1)+3*k1/32+9*k2/32); k4=ht*myf(t(i-1)+12*ht/13,w(i-1)+1932*k1/2197-7200*k2/2197+7296*k3/2197); k5=ht*myf(t(i-1)+ht,w(i-1)+439*k1/216-8*k2+3680*k3/513-845*k4/4104); w(i)=w(i-1)+25*k1/216+1408*k3/2565+2197*k4/4104-k5/5; end plot(t,w, 'b','LineWidth', 2);hold on title('Adams Four-Order Pridictor-Corrector'); function f=myf(t,y) kk=6.22*10^(-19); n1=2*10^(3); n2=n1; n3=3*10^(3); f=kk*(n1-y/2)^2*(n2-y/2)^2*(n3-3*y/4)^3;
评论
    相关推荐
    • RKG.zip
      求解常微分方程组初值问题的龙格库塔法分析及其C代码
    • MATLAB解矩阵微分方程.zip
      本文先简要介绍了矩阵的积分和...接着通过弹簧-质量块儿-阻尼器系统的例子,阐述了求解常微分方程数值解的常用方法;最后,在前两部分的基础上,总结了用MATLAB求解形如_X(t) = f(X(t); t) 矩阵微分方程数值解的方法。
    • sanwanjufa.zip
      数值分析龙格库塔法和三弯矩法求值,是用来求解非线性常微分方程的解的重要的一类隐式或显式迭代法。
    • 龙格库塔四阶:经典四阶龙格-库塔法-matlab开发
      数值分析中,Runge-Kutta 方法是一类用于求解常微分方程近似解的迭代方法。 在这里,对于 e = 0.1 的偏心率,实现了从 t0 = 0 到 t = 86400(s) 的归一化二体问题的积分。
    • 数值分析欧拉法、改进的欧拉法、4阶龙格库塔法 MATLAB算法
      数值解与理论解对比可知,四阶龙格-库塔法的精度已经很高,用它来解一般常微分方程已经足够了。 有程序运行说明 数值解与理论解对比可知,四阶龙格-库塔法的精度已经很高,用它来解一般常微分方程已经足够了。 有...
    • 数值计算方法与C语言工程函数
      10.3变步长龙格-库塔法 10.4改进中点法 10.5外推法 10.6二阶守恒方程组 10.7刚性方程组 第十一章 常微分方程(组)的边值问题 11.1引言 11.2单向打靶法 11.3双向打靶法 11.4松驰法 11.5应用实例 第十二...
    • 四阶龙格库塔法(Runge-Kutta)求解常微分方程的Matlab程序及案例.rar
      原创开发,四阶龙格库塔法(Runge-Kutta)求解常微分方程的Matlab程序及案例,内含自定义Matlab函数、丰富的演示实例和详细的说明文档,简单易用。
    • 未命名文件夹 2.zip
      这是一篇讲解运用mathematica进行编程的四阶龙格库塔法的实例,感觉还不错
    • 龙格塔库算法.rar
      龙格库塔法就是用来解决此类常微分方程数值解的一种算法,其求解的精度较高,可以在已知方程导数和初值信息,利用计算机仿真时应用,省去求解微分方程的复杂过程。同时在工程实际问题中也有着广泛的应用。
    • matlabcnhelp.rar
      matlab中文帮助很难找的,快速下载