参考学习资料 建模第一讲微分方程简介.docVIP

第一讲微分方程简介

一、微分方程的概念

我们先看下面两个例子

例1（求曲线方程问题）已知曲线上任意一点处的切线斜率为，且此曲线过点，求此曲线方程。

解设所求曲线方程为，根据导数的几何意义，有

或（1）

对上式两边积分，得（2）

由于所求曲线过点，因此（2）式应满足条件，代入（2）式得

即所求曲线方程为

例2（运动问题）设质量为的物体，在时刻时自由下落，在空气中受到的阻力与物体的下落速度成正比，求物体下落时的位移与时间的关系式。

解设为物体下落时的位移，

于是，物体下落的速度为，加速度为

空气阻力为

由牛顿第二定律，可以列出方程

即

或（3）

于是，问题归结为求满足（3）式的未知函数。

为了简便起见，现在我们只考虑物体在真空中下落，即没有空气阻力，

这时，方程（3）变为

（4）

对上式两边做两次积分，得

考虑物体下落时的初始状态，即，初速度，将这两个条件代入上述两式，得

于是，自由落体运动的位移公式为

定义（1）

含有未知函数的导数或微分的\o等式等式称为微分方程。未知函数是一元函数的微分方程称为常微分方程；未知函数为多元函数的微分方程称为偏微分方程。

定义（2）

微分方程中所出现的未知函数的导数或微分的最高阶数称为微分方程的阶。满足微分方程的函数称为微分方程的解。求一个微分方程的解的过程称为解微分方程。

定义（3）

如果常微分方程的解中含有独立的任意常数（“独立”指不能将不同的常数合并），且独立的任意常数的个数与方程的阶相同，则这样的解称为方程的通解；不含任意常数的解称为特解。

由通解确定特解时，通常需要一些函数值和导数值，这些已知的函数值和导数值称为方程的初值或初始条件，把微分方程和其初值合到一起，称为微分方程的初始问题。

二、几种特殊类型的常微分方程的解法：

在一阶常微分方程中，有一类比较简单且常见的方程——可分离变量的方程，它的一般形式为

（分别是的连续函数）

求解方法如下：

（1）分离变量

（2）两边积分

（3）求出积分，得通解

其中，分别是的原函数，为任意常数。

例3求微分方程满足初始条件的特解。

解分离变量

两边积分，得

即

得通解为

将代入通解中，得

于是，所求方程满足初始条件的特解为。

三、用MATLAB求解常微分方程:

1.微分方程的解析解

求微分方程（组）的解析解命令：

dsolve(‘方程1,方程2，...方程n’,‘初始条件’，‘自变量’)

记号：在表达微分方程时，用字母D表示求微分，D2，D3等表示求高阶微分。任何D后所跟的字母为因变量，自变量可以指定或由系统规则选定为缺省。

例如，微分方程应表达为：

例4求常微分方程的通解。

解输入命令：dsolve(‘Dy=2’，’x’)

ans=

2*x+C1

该微分方程的通解为

其默认的独立变量为t。

例5求常微分方程的通解。

解输入命令：dsolve(‘D2y-2*Dy+5*y=exp(x)*sin(2*x)’,‘x’)

例6求常微分方程满足初始条件的特解。

解输入命令：dsolve(‘D2y-y=4*x*exp(x)’,‘y(0)=0,Dy(0)=1’,‘x’)

ans=

exp(x)-exp(-x)+(-1+x)*x*exp(x)

该微分方程的特解为

例7求常微分方程组的通解。

解输入命令：[x,y]=dsolve(‘Dx=y+1,Dy=x+1’)

例如，求上例满足初始条件的特解。

输入命令：[x,y]=dsolve(‘Dx=y+1,Dy=x+1’,’x(0)=-2,y(0)=0’,’t’)

作业：

解微分方程

2.求微分方程满足初始条件的特解。

3.解微分方程

4.求微分方程的通解。

5.求方程满足初始条件的特解