MATLAB 第六章 常微分方程 边值问题 刚性方程组

一个人的价值，不是得到了什么，而是拒绝过什么。

一、一阶常微分方程求解

[t,y]=ode45(odefun,tspan,y0)
依次输入微分方程的函数f（t,y），自变量的范围，y的初值，可以得到微分方程的自变量和因变量的向量
求解下面的微分方程


如果没有输出向量，就直接作出图形。

二、高阶常微分方程求解

对于高阶常微分方程，我们首先应该观察阶数确定向量个数，将每阶的求导依次写出来
设h(1)=y
h(2)=y’
则h(1)‘=h(2)
h(2)’=0.01y’^2-2y+sint=0.01h(2) ^ 2-2h(1)+sint


作出图像

寻找t=5时y的值

三、解常微分方程组

分别对每个方程建立辅助变量，t表示自变量，h(1)=x,h(2)=y

下面再看一个比较复杂的方程组


引入辅助变量，t=$\eta$
h(1)=f
h(2)=f’
h(3)=f’’
h(4)=T
h(5)=T’
d h ( 1 ) d t = h ( 2 ) \frac{dh(1)}{dt}=h(2)