数学建模中的差分法.pptVIP

西北大学数学系西北大学数学系差分方程模型则差称为函数的差分，也称一阶差分，记为，即设函数，记为，当t取遍非负整数时，函数值可以排成一个数列，1差分：一差分方程同理，可定义三阶差分等。二阶及二阶以上的差分称为高阶差分。差分的性质：0102二阶差分差分方程：含有未知函数及表示未知函数的几个时期值的符号的方程。如01差分方程的阶方程中含未知函数附标的最大值与最小值的差数。不同形式的方程可互相转化。0212543如果一个函数代入差分方程后，方程两端恒等，则称此函数为该方程的解。例一阶线性方程有解类似于微分方程可定义初值条件，特解等。123454差分方程的解满足方程贰对于一阶线性常系数差分方程壹即平衡点为肆的解，称为上方程的平衡点。叁二差分方程平衡点的稳定性03方程（2）的解可表示为02方程（1）平衡点的稳定性问题可转化为下面方程零点的稳定性。01由于04可得到下面的稳定性结论。当且仅当时，方程（2）的平衡点（零点）是稳定的，从而方程（1）的平衡点稳定。对于n维向量常数矩阵A构成的方程组称其为一阶常系数线性齐次差分方程组。结论1若r(A)1，则其平衡点是全局渐近稳定；若r(A)1，则其平衡点是不稳定的；若r(A)=1，稳定性不确定。称为矩阵A的谱，称为矩阵A的谱半径。A的特征根的集合2对于二阶线性常系数差分方程平衡点为为了得到（4）零点的稳定性我们求解方程（4）。写出特征方程解出特征值通解为其中常数由初值条件确定。当且仅当时，方程（4）的平衡点是稳定的。结论2非齐次线性方程（5）的稳定性可转化为齐次方程（4）来研究。对于n阶线性方程平衡点稳定的条件是特征根3一阶非线性差分方程平衡点通过求解方程而得到。研究稳定性的方法之一是研究其对应的线性部分的稳定性。将方程（6）的右端在点作泰勒展开只取一次项，（6）近似为也是（7）平衡点。当01结论302时，方程（6）与（7）03平衡点的稳定性相同。04结论405当06时，方程（7）平衡点是稳定的；07当08时，方程（7）平衡点是不稳定的。09常微分方程向差分方程转化（数值解）Euler方法求初值问题的近似解。只要给定就可求得先把自变量所在的区间n等分；例1从出发并取，求下列初值问题的近似解。01解02继续下去，自变量使用等间隔值，并生成其n个值，令01步数n可任意大，但n太大，会有误差积累。02优点：容易编程计算。03例2从出发并取，求下列初值问题的近似解。01解02解01Malthus模型的离散形式02例3对于方程组的情形，Euler方法同样可用。先把自变量所在的区间n等分；步数n可任意大，但n太大，会有误差积累。12用Euler法求出前三次逼近，初始条件为解第一组点：对捕食模型第一组点：第二组点：继续下去，就可生成数值解表。01如果方程组为自治系统，在相平面上就可得到近似的轨线图。02第三组点：01.对捕食模型02.用Euler法，在相平面上画出轨线的近似图，观察其变化情况。上机练习1：213Runge-kutta型方法也是用来求初值问题的近似解，但比Euler方法收敛更快。先把自变量所在的区间n等分；Euler方法01常取02单步方法。03差分形式Logistic模型离散化变形令b=r+101.02.03.04.四差分方程模型举例1）平衡点及其稳定性求平衡点：平衡点为根据稳定的条件当b1，只有一个稳定的平衡点，当1b3，变成不稳定的，而成为稳定的，当b3时，是不稳定的平衡点。西北大学数学系西北大学数学系