[计算机数值分析]四阶龙格-库塔经典格式解常微分方程的初值问题

龙格-库塔方法的设计思想：

四阶龙格-库塔方法的经典格式:

程序设计框图：

例：设取步长h=0.2，从x=0到x=1用四阶经典格式解决以下常微分方程的初值问题。

运行示例：

源码：

#include<iostream>
#include<iomanip>
double f(double x, double y);
double f1(double x);using namespace std;int main(void)
{double x0, y0, h;   //x0,y0为初值;h为步长double x1, y1;   //步进后的新值double K1, K2, K3, K4;   //四阶经典格式中的斜率int N;   //N为步数cout << "请输入初值:";cin >> x0 >> y0;cout << "请输入步长:";cin >> h;cout << "请输入步数:";cin >> N;//输出提示信息int i = 1;cout << "\t" << setw(10) << "xn" << "\t" << setw(10) << "yn" << "\t" << setw(10) << "y(xn)" << endl;for (i; i <= N; i++){x1 = x0 + h;   //离散点//求斜率K1 = f(x0, y0);K2 = f(x0 + h / 2, y0 + h / 2 * K1);K3 = f(x0 + h / 2, y0 + h / 2 * K2);K4 = f(x1, y0 + h * K3);//求新值y1 = y0 + h / 6 * (K1 + 2 * K2 + 2 * K3 + K4);//输出本次步进后的离散点数据cout << i << "\t" << setw(10) << x1 << "\t" << setw(10) << y1 << "\t" << setw(10) << f1(x1) << endl;//把新的值赋给旧值，开启下一趟循环x0 = x1;y0 = y1;}return 0;
}double f(double x, double y)   //欧拉格式中的x，y的函数关系式，即f(xn,yn)
{double result;result = y - 2 * x / y;return result;
}double f1(double x)   //实际函数解析式
{double result;result = sqrt(1 + 2 * x);return result;
}

[计算机数值分析]四阶龙格-库塔经典格式解常微分方程的初值问题相关推荐

数值分析—四阶龙格库塔python复现
数值分析-四阶龙格库塔python复现 python入门,复现<数值分析>中的欧拉.龙格库塔,这里以经典四阶龙格库塔来示例. 四阶龙格库塔公式 y_(n+1)=y_n+h/6 (K_1+2 ...
数值分析C++实现用四阶龙格-库塔(Runge-Kutta)方法求解常微分方程初值问题
问题:用四阶龙格-库塔(Runge-Kutta)方法求解常微分方程初值问题. 算法描述算法的程序框图: 具体算法: (1)读取a,b,n,f (2)计算步长h = (b-a)/n,x0=a,y0=f ...
龙格库塔公式法解微分方程组初值问题实例
用四阶龙格库塔公式(取h=0.1)解下列微分方程组初值问题: { y 1 ′ = 120 − 2 y 1 + 2 y 2 ( 0 ≤ x ≤ 1.0 ) , y 2 ′ = 2 y 1 − 5 y 2 ...
[计算机数值分析]改进欧拉格式解常微分方程的初值问题
欧拉方法yn+1 = yn +hf(xn,yn) 是一种显式算法,计算量小,但精度很低.梯形方法yn+1 = yn + h / 2 *[f(xn,yn) + f(xn+1,yn+1)] 虽然提高了精度 ...
二阶水箱 matlab 四阶龙格库塔,请问这个二阶常微分方程组用龙格四阶库塔法怎么编写...
要是不用MATLAB自带的ode45函数也可以网上下载一个4阶龙格库塔算法来代替附上一个仅供参考 function y=DELGKT4_rungekuta(f,h,a,b,y0,varvec) ...
四阶龙格库塔方程（Rungekutta）解二阶常微分方程组并计算船舶在迎浪下的纵摇埀荡耦合运动方程-附Matlab代码
今年年初的时候给师姐做了DDPG算法的船舶减横摇控制算法,师姐还有想法要让我把纵摇-埀荡两个自由度的减摇也做出来,这个任务归我了.实际上不管是多少个自由度的减摇,其实都需要解运动方程,当初做单自由度横 ...
四阶龙格库塔方程解二阶常微分方程组并计算船舶在迎浪下的纵摇埀荡耦合运动方程-附Python代码
0 写在前面这篇博客是在将我上一篇博客的matlab代码移植到python中,应为后续要开展深度强化学习下的船舶减摇研究,总的来说还是在python上做这项工作比较合适.具体方程的解法和背景不在赘述 ...
一阶欧拉近似matlab,MATLAB改进欧拉法与四阶龙格-库塔求解一阶常微分方程.doc
MATLAB改进欧拉法与四阶龙格-库塔求解一阶常微分方程姓名:樊元君学号:2012200902 日期:2012.11.06 一.实验目的掌握改进欧拉法与四阶龙格-库塔求解一阶常微分方程的初值问题 ...
四阶龙格库塔法的基本思想_四阶龙格库塔实验报告.docx
四阶龙格库塔实验报告三.四阶Runge-Kutta法求解常微分方程一.龙格库塔法的思想根据第九章的知识可知道,Euler方法的局部截断误差是,而当用Euler方法估计出再用梯形公式进行校正,即采用改 ...

[计算机数值分析]四阶龙格-库塔经典格式解常微分方程的初值问题

[计算机数值分析]四阶龙格-库塔经典格式解常微分方程的初值问题相关推荐

最新文章

热门文章