2021-01-07 matlab数值分析 常微分方程初边值问题数值解 标准龙格库塔四阶四段公式 欧拉法
matlab数值分析 常微分方程初边值问题数值解 标准龙格库塔四阶四段公式 欧拉法
1、标准龙格库塔四阶四段公式
function y=rk4(fun,a,b,y0,n)
h=(b-a)/n;
y(1)=y0;
for k=1:nx=a+(k-1)*h;k1=h*feval(fun,x,y(k));k2=h*feval(fun,x+h/2,y(k)+k1/2);k3=h*feval(fun,x+h/2,y(k)+k2/2);k4=h*feval(fun,x+h,y(k)+k3);y(k+1)=y(k)+(k1+2*k2+2*k3+k4)/6;
end
函数文件
function u=frk4(x,y)
u=y-2*x/y;
调用程序
y=rk4(@frk4,0,1,1,10)
2 欧拉法
function y=euler(f,a,b,y0,h)
n=(b-a)/h;
y(1)=y0;
for i=1:nx(i)=a+(i-1)*h;y(i+1)=h*feval(f,x(i),y(i));
end
函数文件
function u=feuler(x,y)
u=x^3-y/x;
调用程序
y=euler(@feuler,1,2,0.4,0.2)
2021-01-07 matlab数值分析 常微分方程初边值问题数值解 标准龙格库塔四阶四段公式 欧拉法相关推荐
- Matlab求常微分方程组的数值解
上篇博客介绍了Matlab求解常微分方程组解析解的方法:博客地址 微分方程组复杂时,无法求出解析解时,就需要求其数值解,这里来介绍. 以下内容按照Matlab官方文档提供的方程来展开(提议多看官方文档 ...
- 龙格库塔法matlab求解微分方程组,微分方程组的龙格库塔公式求解matlab版.pdf
微分方程组的龙格库塔公式求解matlab版 微分方程组的龙格-库塔公式求解matlab版 南京大学 王寻 1. 一阶常微分方程组 考虑方程组 y'f x,y,z , y x y ...
- STM32学习笔记2021/01/07
关于stm32的寄存器 引言 STM32F103ZET6共有144个引脚,共有7组每组16个共计112个I/O端口:每个GPIO端口有两个32bit的配置寄存器(GPIOx_CRL,GPIOx_CRH ...
- 0基础快速入门CSS技术栈(5)—图解详细阐述说透CSS的盒子模型(超级重要)、圆角边框、盒子阴影及相关重要的笔试题——css的核心中的核心(附详细案例源码解析过程)2021.01.07更新
文章目录 1. 盒子模型(CSS重点) 1.1 看透网页布局的本质 1.2 盒子模型(Box Model) 1.3 盒子边框(border) 1.3.1 边框综合设置 1.3.2 example01 ...
- MATLAB中ode23函数,龙格库塔函数
今天说一说MATLAB中ode23函数的原理,在网上看了好多,但是不知道是怎么计算的,就知道是那么用的,但是最后结果咋回事不知道,今天来讲一讲是怎么计算的. 首先来个程序: function f=eg ...
- [常微分方程的数值解法系列五] 龙格-库塔(RK4)法
龙格-库塔法 简介 基本思想 具体方法 一阶 二阶 求解参数 特殊二阶 三阶 高阶 步长选择 例子 在惯性导航以及VIO等实际问题中利用IMU求解位姿需要对IMU测量值进行积分得到需要的位置和姿态,其 ...
- 二阶龙格库塔公式推导_[常微分方程的数值解法系列五] 龙格-库塔(RK4)法
在惯性导航以及VIO等实际问题中利用IMU求解位姿需要对IMU测量值进行积分得到需要的位置和姿态,其中主要就是求解微分方程.但之前求解微分方程的解析方法主要是应用于一些简单和特殊的微分方程求解中,对于 ...
- matlab微分方程组边值,matlab求解常微分方程边值问题的方法
matlab求解常微分方程边值问题的方法 Matlab 求解常微分方程边值问题的方法:bvp4c 函数 常微分方程的边值问题,即 boundary value problems ,简称 BVP 问题, ...
- Caputo 分数阶微分方程-慢扩散方程初边值问题基于L1 逼近的空间二阶方法及其Matlab程序实现
2.3.3 Caputo 分数阶一维问题基于 L1 逼近的空间二阶方法 考虑如下时间分数阶慢扩散方程初边值问题 { 0 C D t α u ( x , t ) =
最新文章
- 自回归与非自回归模型如何兼得?预训练模型BANG或许可解
- zabbix系列~ 监控模式
- Chrome浏览器如何不让它缓存?
- python super 理解(四)
- 6-1 求链式表的表长
- 计算机设备管理程序在哪,电脑设备管理在哪,教您如何打开电脑设备管理
- The Model Driven Software Network
- C语言strcat()库函数的实现
- .Net程序员学用Oracle系列(16):访问数据库(ODP.NET)
- simple_html_dom.php 使用 乱码处理作者:gaoming13
- 程序员踩坑之旅:将 75000 行 iOS 原生代码迁移到 Flutter!
- GeoServer发布矢量地图服务
- Office Word 之画肘形箭头
- java正则表达标点符号_用正则表达式去除标点符号
- 人类有两大学习能力,即记忆力和理解力
- Consider injecting the bean as one of its interfaces or forcing the use of CGLib-based proxies 问题解决
- GOT-10k: A Large High-Diversity Benchmark forGeneric Object Tracking in the Wild(论文翻译)
- 如何诊断SQL数据?
- c#+wps插件开发(1)
- ogg高版本到低版本同步