matlab笔记：一元微分方程求解

为什么要在这里写这东西。。还不是怕自己忘了。。

主要是ode45的用法。。这个函数用来解决微分方程组

先上代码

tspan=[1:0.01:10]；
X0=[0;15]；
delta=@(x,y) [y(2);-29*y(1)-4*y(2)]；
[T,X]=ode45(delta,tspan,X0)；
plot(T,X(:,1))；
plot(T,X(:,2))；

先声明一点，y(1)表示原函数，y(2)才表示导数，y(3)表示二阶导。。以此类推

然后解释一下代码。。

第一行就是定义一个函数规则（叫function handle？有兴趣以后了解一下），感觉就是类似于c++里面的虚函数吧，仅仅是为了传参而生。。

他这个函数就是返回一个列向量，内容可以自己定义

在被ode45调用时，这个列向量delta代表的含义是

y(2)=delta(1)

y(3)=delta(2)

y(4)=delta(3)

.....

有时候我们对一些导数是没什么需求的，例如上例只有一个微分方程，所以只表示成了y(3)=xxx的形式。。

而对y(2)的方程已经是形式上的存在，令y(2)=y(2)该方程就废了。。

说说再ode45的返回值和参数。。

第一个参数代表函数名，也可以是定义了函数的文件名，不过为了定义一个函数专门写个文件感觉没什么必要就是。

第二个参数是一个行向量，可以算是定义域了吧，就是在求解微分方程自变量的取值范围

第三个参数是一个列向量，表示初值，就是x=tspan(1)时各函数值取值为多少什么的。。

然后就是返回值，第一个返回值就是tspan'

第二个返回值是一个矩阵

X(i,j)表示x取T(i)时的j阶导数

所以可以看出ode45并没有真正求解微分方程，他只是把特定点（给定的tspan）的导数值求出来而已。。

不过只要tspan各数间距够小，就能用plot画出图像。。对求解问题也是非常有帮助。。

matlab笔记：一元微分方程求解相关推荐

MATLAB实现一元三次方程求解/盛金公式
MATLAB实现一元三次方程求解/盛金公式一元三次方程求解中,1945年卡尔丹诺把冯塔纳的三次方程求根公式发表出来,但该公式形式比较复杂,直观性也较差.1989年范盛金对一元三次方程求解进行了深入的 ...
Maple笔记2--常微分方程求解
转需看原文地址:Maple笔记2--常微分方程求解作者:Lionel 来源:网络论坛转载(VB资料库) 常微分方程求解微分方程求解是数学研究与应用的一个重点和难点. Maple能够显式或隐式地解析 ...
maple 解代数方程组得多项式_Maple笔记2--常微分方程求解
来源:网络论坛转载(VB资料库) 常微分方程求解微分方程求解是数学研究与应用的一个重点和难点. Maple能够显式或隐式地解析地求解许多微分方程求解. 在常微分方程求解器dsolve中使用了一些传统 ...
matlab二阶非线性微分方程求解,二阶非线性常微分方程的打靶法.ppt
二阶非线性常微分方程的打靶法二阶非线性常微分方程的打靶法计算思路主要分为以下五步: 给定容许误差ε,迭代初始值γ1,对k=1,2,...做: (1)用四阶Runge-Kutta 方法求解初值问, ...
matlab差分方程和微分方程,实验二微分方程与差分方程模型matlab求解
实验二微分方程与差分方程模型matlab求解实验二: 微分方程与差分方程模型 Matlab 求解一.实验目的[1] 掌握解析.数值解法,并学会用图形观察解的形态和进行解的定性分析:[2] 熟悉 MA ...
【MATLAB】符号数学计算（七）：符号微积分、符号微分方程求解、符号代数方程求解
一.符号表达式的极限 limit(F,x,a):求当时,符号表达式F的极限. limit(F,a):符号表达式F采用默认自变量(可由函数findsym求得),该函数求F的自变量趋于a时的极限值. li ...
matlab输入数据作方程,用MATLAB函数编写并求解微分方程
实验八:用MATLAB 函数编写并求解微分方程一.实验原理为了对连续系统进行方针,首先需要建立其数学模型,然后利用计算机求这些数学模型,从而得出数学模型的数值解.由于连续系统是通过微分方程老建模的 ...
matlab方程求解的实验报告,用matlab对微分方程求解实验报告.doc
PAGE PAGE 1 o <高等数学>上机作业(三) 课程 <高等数学> 上机内容微分方程求解成绩姓名专业班级学号教学班指导教师上机 ...
matlab行星运动轨迹仿真动画,利用Matlab可视化功能实现微分方程求解行星运动轨迹...
利用Matlab可视化功能实现微分方程求解行星运动轨迹 1.背景在物理学璀璨的发展史上,物理学家花了很长时间研究我们头顶浩瀚的星空,试图探究星星的运行模式,以及地球自身的运动模式.其中不乏像亚里士多 ...

matlab笔记：一元微分方程求解

matlab笔记：一元微分方程求解相关推荐

最新文章

热门文章