第二十四讲 一阶常微分方程组

一,常微分方程组:

  • 必须满足:自变量只有1个,因变量有多个
  • 一阶形式:\frac{dx}{dt}=f(x,y,t)\frac{dy}{dt}=g(x,y,t),x和y是因变量,t是自变量

二,什么样的常微分方程组是线性方程组?

  • 必须满足:x和y以线性形式出现
  • 例如:\frac{dx}{dt}=a(t)x+b(t)y+f_{1}(t)\frac{dy}{dt}=c(t)x+d(t)y+f_{2}(t),a、b、c和d可以是t的函数

三,什么样的线性方程组是常系数方程组?

  • 必须满足:a、b、c和d都是常数
  • 例如:\frac{dx}{dt}=ax+by+f_{1}(t)\frac{dy}{dt}=cx+dy+f_{2}(t)

四,什么样的常系数方程组是齐次的?

  • 必须满足:f_{1}(t)=0f_{2}(t)=0
  • 例如:\frac{dx}{dt}=ax+by\frac{dy}{dt}=cx+dy

五,初始条件:

  • 需要初始条件是为了:求出通解中的任意常数,从而得到特解
  • 需要初始条件的个数=通解中任意常数的个数
  • 通解中任意常数的个数=方程组总的阶数
  • 例如:\frac{dx}{dt}=f(x,y,t)\frac{dy}{dt}=g(x,y,t){x}'{y}'一共两阶,表示有两个任意常数
  • 初始条件:x(t_{0})=x_{0}(位置),y(t_{0})=y_{0}(位置)

六,例题:

  • 如果微分方程是一个模型,那么微分方程组就是两个以上的模型套在一起
  • 假设在一个盛满冰水混合物(0摄氏度)的烧杯里放进一个刚煮熟的鸡蛋(热的,可以看见鸡蛋里面的蛋白和蛋黄),如图:
  • 水的温度是T_{e}(t),蛋黄的温度是T_{1},蛋白的温度是T_{2}
  • 问题:求蛋黄和蛋白的温度随时间的变化率\frac{dT_{1}}{dt}\frac{dT_{2}}{dt}
  • 初始条件:T_{1}(0)=40T_{2}(0)=45
  • 建立模型:
  1. 根据牛顿的传导定律:\frac{dT_{1}}{dt}=a(T_{2}-T_{1}),a表示蛋黄和蛋白间的传导常数
  2. 同理:\frac{dT_{2}}{dt}=a(T_{1}-T_{2})+b(T_{e}-T_{2}),b表示水和蛋白间的传导常数
  3. 整理成标准方程组:{T_{1}}'=-aT_{1}+aT_{2}{T_{2}}'=aT_{1}-(a+b)T_{2}+bT_{e}
  4. 如果是沸水中煮鸡蛋,水温的变化函数是:T_{e}=100e^{-kt},100表示起始温度,e^{-kt}表示温度以指数规律冷却。但这里是冰水混合物,假设水温恒定为T_{e}=0(这么做其实是为了将方程组简化)
  5. 假设传导函数:a=2b=3
  6. 方程组变为:{T_{1}}'=-2T_{1}+2T_{2}①,{T_{2}}'=2T_{1}-5T_{2}
  • 换元法解方程组:
  1. T_{1}表示T_{2}T_{2}=\frac{​{T_{1}}'+2T_{1}}{2}
  2. 将结果代入②式:{(\frac{​{T_{1}}'+2T_{1}}{2})}'=2T_{1}-5(\frac{​{T_{1}}'+2T_{1}}{2})
  3. 整理:{({T_{1}}'+2T_{1})}'=4T_{1}-5({T_{1}}'+2T_{1})
  4. 整理:{T_{1}}''+7{T_{1}}'+6T_{1}=0
  5. 数学难度守恒定律:消去了T_{2},但是T_{1}变成了二阶,说明消元后剩下单个方程的阶数等于原方程组阶数之和
  6. 顺便一提:只有当系数7和6都为正数时,{T_{1}}''+7{T_{1}}'+6T_{1}=0这个系统才是稳定的(解随时间趋于无穷大而趋于0),例如,可以把7看成是阻尼常数,6是弹性常数,如果阻尼常数为0或为负,则系统会一直振荡下去
  • 解齐次方程{T_{1}}''+7{T_{1}}'+6T_{1}=0
  1. 求特征方程:r^{2}+7r+6=0
  2. 解特征根:r_{1}=-1r_{1}=-6
  3. 得通解:T_{1}=c_{1}e^{-t}+c_{2}e^{-6t}
  4. T_{2}的通解:T_{2}=\frac{​{T_{1}}'+2T_{1}}{2}=\frac{(c_{1}e^{-t}+c_{2}e^{-6t})'+2(c_{1}e^{-t}+c_{2}e^{-6t})}{2}
  5. 整理:T_{2}=\frac{1}{2}c_{1}e^{-t}-2c_{2}e^{-6t}
  6. 代入初始条件:T_{1}(0)=40T_{2}(0)=45
  7. 方程组变为:40=c_{1}+c_{2}①,45=\frac{c_{1}}{2}-2c_{2}
  8. 解出系数:c_{1}=50c_{2}=-10
  9. 代入通解,得特解:T_{1}=50e^{-t}-10e^{-6t}T_{2}=25e^{-t}+20e^{-6t}

七,自治微分方程组的几何意义:

  • 自治微分方程组:\frac{dx}{dt}=f(x,y)\frac{dy}{dt}=g(x,y),等式右边不含自变量t
  • 参数曲线:x=x(t)y=y(t),如图(红色线,有无数条):
  • 解的含义:{x}'(t){y}'(t)表示参数曲线上,某一时刻的变化率(因为无数个变化率加在一起就是曲线,所以解实际上是参数曲线)
  • 自治微分方程组的几何意义:无数变化率组成的斜率向量场,如图(蓝色箭头):
  • 含有两个一阶自治方程的方程组等于一个速度矢量场
  • 6
    点赞
  • 26
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值