常微分方程的RK4解法

#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#define H 0.1
int main()
{float x,y,r=3;float k1,k2,k3,k4;float fun1(float,float);int i,n;printf("请输入x,y的初值：\n");scanf("%f%f",&x,&y);printf("请输入运算次数：\n");scanf("%d",&n);printf("运算结果为：\n"); for(i=0;i<n;i++){printf("x%-2d=%.2f,\ty%-2d=%f,\ty(%-2d)=%f,\tr(%-2d)=%f",i,x,i,y,i,r,i,r-y);printf("\n"); k1=fun1(x,y);k2=fun1(x+H/2,y+H*k1/2);k3=fun1(x+H/2,y+H*k2/2);k4=fun1(x+H,y+H*k3);x=x+H;y=y+H*(k1+2*k2+2*k3+k4)/6;r=3/(1+x*x*x);}system("pause");return 0;
}float fun1(float x,float y){float f;f=-x*x*y*y;return f;}

常微分方程的RK4解法相关推荐

[常微分方程的数值解法系列五] 龙格-库塔(RK4)法
龙格-库塔法简介基本思想具体方法一阶二阶求解参数特殊二阶三阶高阶步长选择例子在惯性导航以及VIO等实际问题中利用IMU求解位姿需要对IMU测量值进行积分得到需要的位置和姿态,其 ...
二阶龙格库塔公式推导_[常微分方程的数值解法系列五] 龙格-库塔(RK4)法
在惯性导航以及VIO等实际问题中利用IMU求解位姿需要对IMU测量值进行积分得到需要的位置和姿态,其中主要就是求解微分方程.但之前求解微分方程的解析方法主要是应用于一些简单和特殊的微分方程求解中,对于 ...
[常微分方程的数值解法系列三] 改进欧拉法(预估校正法)
改进欧拉法简介预估-校正截断误差例子在惯性导航以及VIO等实际问题中利用IMU求解位姿需要对IMU测量值进行积分得到需要的位置和姿态,其中主要就是求解微分方程.但之前求解微分方程的解析方法主 ...
[常微分方程的数值解法系列二] 欧拉法
欧拉法简介几何意义证明泰勒展开近似求导近似积分近似几种欧拉方式向前欧拉公式向后欧拉公式梯形公式中点公式截断误差求解过程向前欧拉公式例子向前欧拉公式在惯性导航以及VIO ...
[常微分方程的数值解法系列四] 中值法
中值法简介具体步骤截断误差例子在惯性导航以及VIO等实际问题中利用IMU求解位姿需要对IMU测量值进行积分得到需要的位置和姿态,其中主要就是求解微分方程.但之前求解微分方程的解析方法主要是应 ...
【数理知识】《数值分析》李庆扬老师-第9章-常微分方程初值问题数值解法
第8章回到目录无第9章-常微分方程初值问题数值解法 9.1 引言利普希茨 (Lipschitz) 条件 / 利普希茨常数定理1 解的存在唯一性定理定理2 解对初值依赖的敏感性 9.2 简单 ...
二阶常微分方程的数值解法（中心差分法和有限体积法）
二阶常微分方程的数值解法(中心差分法和有限体积法) 这里我们介绍中心差分法和有限体积法求解方程. 题目: 用差分法的中心差分格式和有限体积法求解两点边值问题 u′′−α(2x−1)u′−2αu=0,0 ...
一阶常微分方程的数值解法（二阶显式、隐式 Adams 公式及 Milne 方法）
一阶常微分方程的数值解法这里我们介绍二阶显式.隐式 Adams 公式及 Milne 方法求解方程. 题目: 对初值问题 u ′ = u − t 2 , 0 ≤ t ≤ 1 , u ( 0 ) = 0 ...
常微分方程数值解的c语言程序,常微分方程的数值解法一阶常微分方程数值解的C语言编程实现.doc...
常微分方程的数值解法一阶常微分方程数值解的C语言编程实现导读:就爱阅读网友为您分享以下"一阶常微分方程数值解的C语言编程实现"资讯,希望对您有所帮助,感谢您对92的支持! 一阶 ...

常微分方程的RK4解法

常微分方程的RK4解法相关推荐

最新文章

热门文章