第四章 微分方程与差分方程方法.docVIP

第四章 微分方程与差分方程方法 第一节 微分方程模型 我们在数学分析中所研究的函数,是反映客观现实世界运动过程中量与量之间的一种关系,但我们在构造数学模型时,遇到的大量实际问题往往不能直接写出量与量之间的关系,却能比较容易地建立这些变量和它们的导数(或微分)间的关系式,这种联系着自变量、未知函数及其导数(或微分)的关系式称为微分方程. §4.1.1 微分方程简介 这一节,我们将介绍关于微分方程的一些基本概念. 一、微分方程的阶数 首先我们具体的来看一个微分方程的例子. 例4-1 物体冷却过程的数学模型 将某物体放置于空气中,在时刻,测量得它的温度为,10分钟后测量得温度为.我们要求决定此物体的温度和时间的关系,并计算20分钟后物体的温度.这里我们假定空气的温度保持为. 解:根据物理学中的牛顿冷却定律可知,热量总是从温度高的物体向温度低的物体传导;一个物体的温度变化速度与这一物体的温度与其所在介质温度的差值成正比.设物体在时刻的温度为,则温度的变化速度可以用来表示.我们得到描述物体温度变化的微分方程 (4.1.1) 其中是比例常数. 方程(4.1.1)中含有未知函数及它的一阶导数,这样的方程,我们称为一阶微分方程.微分方程中出现的未知函数最高阶导数的阶数称为微分方程的阶数. 方程 (4.1.2) 中未知函数最高阶导数的阶数是三阶,则方程(4.1.2)称为三阶微分方程. 二、常微分方程与偏微分方程 如果在微分方程中,自变量的个数只有一个,我们称这种微分方程为常微分方程;自变量的个数为两个或两个以上的微分方程称为偏微分方程. 方程 (4.1.3) 就是偏微分方程的例子,其中是未知函数,、、都是自变量.而方程(4.1.1)(4.1.2)都是常微分方程的例子. 三、线性与非线性微分方程 如果阶常微分方程 (4.1.4) 的左端为关于未知函数及其各阶导数的线性组合,则称该方程为线性微分方程,否则称为非线性方程.一般的阶线性微分方程具有形式 (4.1.5) 其中是关于的已知函数. 当时,称(4.1.5)为阶齐次线性微分方程;当时,称(4.1.5)为阶非齐次线性微分方程. 例如,方程(4.1.2)是三阶线性微分方程.而方程 (4.1.6) 是一阶齐次非线性方程. 四、微分方程初边值问题 我们把含有个独立的任意常数的解 称为阶微分方程(4.1.4)的通解.为了确定微分方程一个特定的解,我们通常给出这个解所必须满足的条件,这就是定解条件.常见的定解条件是初始条件.阶微分方程(4.1.4)的初始条件是指如下的个条件: 求解微分方程满足初始条件的解,称为初值问题.求解初值问题的过程,就是通过初始条件确定通解中的常数,从而求得满足初始条件的特解的过程. 若方程所给出的定解条件,既有自变量初始时刻的值,也有自变量取终值时的值,则称该问题为边值问题. §4.1.2微分方程的求解及Matlab实现 线性微分方程和低阶特殊微分方程往往可以通过解析解的方法求解,但一般的非线性微分方程是没有解析解的,即使可以求得解析解,参数也很难确定,需要用数值解的方式求解.具体的求解方法很多,本节我们主要介绍如何利用matlab来求解微分方程的解析解和数值解. §4.1.2.1微分方程的解析解 一、基本理论 有些简单的常微分方程可用一些技巧,如分离变量法、积分因子法、常数变异法、降阶法等,化为可直接积分的方程从而求得显示解.例如,一阶常系数线性常微分方程 可化为 两边积分可得的通解为 . 一般的常系数线性微分方程 (4.1.7) 的解满足叠加原理,即方程(4.1.7)的通解是对应齐次方程 (4.1.8) 的通解与非齐次方程(4.1.7)的一个特解的和. 一阶常系数线性常微分方程总可用这一思路求得显示解. 高阶线性常系数微分方程可用特征根法求对应齐次微分方程(4.1.8)的通解.(4.1.8)对应的代数特征方程 若方程的特征根均可以求出,且两两相异,则方程(4.1.7)的解可以表示为 其中,为待定系数,是满足方程(4.1.8)的一个特解,这个特解可以通过常数变异法,比较系数法得到.有重