三阶及四阶Runge-Kutta法

1. 实验目的

能运用Matlab编程实现Runge-Kutta法与求解微分方程初值问题的数值解；
能用图像来比较数值解与精确解；
熟悉Matlab编程环境。

2. 实验内容

实验题目

用三阶、四阶Runge-Kutta方法求常微分方程初值问题

的数值解，步长h=0.2,0.5,1，并与真解进行比较。

2.实验原理

三阶Runge-Kutta方法：

3. 实验步骤

步1. 用dsolve函数求初值问题的解析解；

步2. 用ezplot函数做出解析解的图像；

步3. 用Matlab编写三阶、四阶Runge-Kutta法求解初值问题的通用程序；

步4. 用编写的程序求初值问题的数值解；

步5. 用plot函数将数值解显示在图像里，并与解析解进行比较。

4. 结果与分析

分别取步长h为1、0.5、0.2可以得到下列图像

5. 实验总结

通过上面的实验我们能够得到四阶Runge-Kutta与三阶Runge-Kutta法能够较快且较好的得到比较精确的解。

function [u,t]=Runge_Kutta3(t0,u0,T,h)
% t表示区间[t0,T]上步长为h节点构成的向量
% u表示由3阶Runge-Kutta法得到的数值解构成的向量
t=[t0:h:T];
n=length(t);
u=zeros(1,n);
u(1)=u0;
for i=2:nK1=t(i-1)*u(i-1);K2=(t(i-1)+0.5*h)*(u(i-1)+0.5*h*K1);K3=t(i)*(u(i-1)-h*K1+2*h*K2);u(i)=u(i-1)+h*(K1+4*K2+K3)/6;
endfunction [u,t]=Runge_Kutta4(t0,u0,T,h)
% t表示区间[t0,T]上步长为h节点构成的向量
% u表示由4阶Runge-Kutta法得到的数值解构成的向量
t=[t0:h:T];
n=length(t);
u=zeros(1,n);
u(1)=u0;
for i=2:nK1=t(i-1)*u(i-1);K2=(t(i-1)+0.5*h)*(u(i-1)+0.5*h*K1);K3=(t(i-1)+0.5*h)*(u(i-1)+0.5*h*K2);K4=t(i)*(u(i-1)+h*K3);u(i)=u(i-1)+h*(K1+2*K2+2*K3+K4)/6;
end
clc
clear all
u_true=dsolve('Du=t*u','u(0)=1','t');
t0=0; T=2;h=0.5;u0=1;kmax=1;
[u2,t2]=Runge_Kutta3(t0,u0,T,h);
[u3,t3]=Runge_Kutta4(t0,u0,T,h);
ezplot(u_true,[t0,T])
hold on
plot(t2,u2,'bs')
hold on
plot(t3,u3,'bo')
grid on
hold off
legend('u true'，'RK3法','RK4法')

三阶及四阶Runge-Kutta法相关推荐

四阶魔方玩法总结V1.0
四阶魔方玩法总结V1.0 1.引言今写此文,我主要是为了方便自己再次玩其魔方的时候,可以快速的想起,避免又从头学起.毕竟自己学会的,理解的,写出来的东西,再次玩魔方的时候,仅仅是回顾和追忆的过程,不 ...
12面魔方公式图解法_【高级篇】（三）三阶魔方CFOP高级玩法之——F2L
一.F2L这一步要干什么 1.先了解一下"棱角对"和"槽位"的概念棱角对:即由一个棱块和一个角块构成,是F2L的基本单元(共四组) 槽位:给"棱角对 ...
三阶魔方大中小魔公式_三阶魔方的入门玩法教程|魔方玩法|魔方视频教程|魔方公式图解|--想成为魔方高手就来魔方乐园吧...
深圳魔方爱好者可以到实体店跟魔方老师学魔方了!查询店铺地址在正式开始学习魔方还原前,请进入网站的首页,抽上10分钟的时间学习一些魔方的基础知识,这样会很有利于我们下面的学习,在还原的过程中有碰到疑问 ...
魔方cfop公式软件_【高级篇】（三）三阶魔方CFOP高级玩法之——F2L
一.F2L这一步要干什么 1.先了解一下"棱角对"和"槽位"的概念棱角对:即由一个棱块和一个角块构成,是F2L的基本单元(共四组) 槽位:给"棱角对 ...
魔方cfop公式软件_【高级篇】（四）三阶魔方CFOP高级玩法之——OLL
一.OLL这一步需要做什么? OLL这一步是最容易理解的,到了这一步,前两层已经全部复原了那不管顶面出现什么情况,都能用一个公式将顶面黄色朝向全部调整正确(顺序不管) 一步到位的代价就是--这步骤有 ...
4阶显式Runge－Kutta法解常微分方程的通用程序--python实现
对于常微分方程,RK方法速度快,精度高,代码简单,是最为实用的数值方法之一.RK方法很简单,类似梯形法,RK法也是根据前一步点的值推算后一步点.具体算法见以下链接 https://wenku.baid ...
欧拉法，改进欧拉法吗，四阶龙格库拉法 c++代码
参考 https://blog.csdn.net/deepsprings/article/details/109850766 #include<iostream> #include< ...
三阶魔方七步还原法公式备忘录
魔方公式备忘转动符号图解魔方七步公式: 1.底面十字还原 2.底角还原 3.中间层还原上棱到左棱 U'L'U'LUFUF' 上棱到右棱 URUR'U'F'U'F 4.顶面十字循环做FRUR'U ...
matlab龙格库塔法求通解,基于matlab及龙格库塔法求解布拉修斯方程.doc
基于matlab及龙格库塔法求解布拉修斯方程 Runge-Kutta法求解布拉修斯解摘要薄剪切层方程主要有三种解法,即相似解,非相似条件下对偏微分方程组的数值解和近似解.布拉修斯解是布拉修斯于19 ...

三阶及四阶Runge-Kutta法

三阶及四阶Runge-Kutta法相关推荐

最新文章

热门文章