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’
$\frac{dh(1)}{dt}=h(2)$
$\frac{dh(2)}{dt}=h(3)$
$\frac{dh(3)}{dt}=-3h(1)h(3)+2h(2{)}^2-h(4)$
$\frac{dh(4)}{dt}=h(5)$
$\frac{dh(5)}{dt}=-2.1h(1)h(5)$

四、边值问题

五、刚性方程组

求解下列刚性方程组

当我们用ode45求解一个方程发现x数量大，计算速度慢时，说明是刚性方程组，可以改用ode15s求解。

六、建模实验

1.夏天把开有空调的室内一支读数为 20 的温度计放到户外，10min 后读数为25.2， 再过 10min 后读数为 28.52，假设温度计值的增速和温差成正比，请推算户外温度。
设温度计读数为T，户外温度为c，则可以根据题意列出方程
T’=k(c-T),T(0)=20
将方程中的每个变量都设成符号向量，用dsolve求解

求出方程以后，用函数代替符号向量

解得c=34.3830 k=0.0449
2.某公司生产一种耐用消费品，市场占有率为 5%时开始做广告，一段时间的市场跟踪调查后，该公司发现：单位时间内购买人口百分比的相对增长率与当时还没有买的百分比成正比，且估得此比例系数为 0.5.
（1） 建立该问题的数学模型，分别求其解析解和数值解，并作比较；
（2） 厂家要做多少时间广告，可以使市场占有率达到 80%？

设市场占有率 x, 增速 x’, 相对增长率 x’/x. 建立微分方程
%x’/x=0.5*(1-x), 等价于 x’=0.5x(1-x) ，初始值 x(0)=0.05



求市场占有率在80%时的时间