SXJM第五部分微分与差分方程培训教程文件.ppt

方程(1)的相轨线是一族包含平衡点A( ) 的封闭曲线. x y o k2/c k1/b A 2. 平衡点A的实际意义 记 T=t1-t0 , 称T 周期,将原方程 中的第二个方程改写为 两边从t0 到 t1 积分，得 结论：食用鱼和食肉鱼的平衡量恰为它们的数量在一个周期内的平均值. o x(食用鱼) y A 3. D′A ncona现象的解释 为考察捕鱼业对两种鱼类的影响, 引入捕捞 能力系数ε，将方程(1)改写为 方程(3)的平衡点为A′( ),由于捕 捞能力系数的引进，食用鱼的平均量增大, 而 食肉鱼的平均量则减少了. Volterra原理：为了减少强者,只需捕获弱者. * * t0=0;t1=15; %设置区域 y0=[20 20]; %给定初值 [t,y]=ode23(‘lot1’,t0,t1,y0); %求解 plot(t,y) %绘图 function z=lot1(x,y) z(1,:)=y(1)-0.01*y(1).*y(2); z(2,:)=-y(2)+0.02*y(1).*y(2); 求解常微分方程: (1)定义函数; (2)调用ode23 可化为一阶常微分方程组求解。 例如，二阶常微分方程初值问题 引进新的变量,令z=y’，即可将上述二阶方程化为如下的一阶方程组的初值问题: 高阶常微分方程的数值解法： 例 求下列高阶微分方程的数值解： 解：显然 假设 则 即二阶问题化为微分方程组的初值问题。 考虑方程： 结合下述三种边界条件之一： 边界问题的解法：打靶法、有限差分法 (7.4)式中 它们分别称为第一、第二、第三边界问题。 七.二阶常微分方程边值问题数值方法 1、打靶法 基本思路：将边值问题转化为初值问题考虑。或者说适当选择初始值使初值问题的解满足边值条件。然后用求解初值问题的任一种有效的数值方法求解。 以第一边界条件为例 考虑边值问题： 取y0=a , 考虑初值问题 待定，由数值解法求解得到在 处的解 ， 使 ，这里 为给定允许误差界，就停止迭代改进, 输出作为数值解。 求解非线性方程 若 , 则得所求数值解。当 该方程可用二分法、正割法或Newton法等来求解。若求得 对第二类边界问题，也可转化为考虑初值问题(7.5)，取 对第三类边界问题，仍可转化为考虑初值问题(7.5)，取 以 为待定参数。 ，以 为待定参数。 2、有限差分法 则离散化成差分方程 将区间[a,b]进行等分： 设在 处的数值解为 。 用中心差分近似微分，即 二阶中心差商 对应的边界条件也离散成 第一边界问题： 第二边界问题： 第三边界问题： 其中 为已知函数， 若 是 的线性函数时，f 可写成 由常微分方程的理论知，通过变量替换总可以消去 则近似差分方程成离散差分方程为 其中 方程中的 项，不妨设变换后的方程为 将以上方程合并同类项整理得方程组 其中只要 ，则方程组的系数矩阵为弱对角占优的三对角阵， 其中 。 而且还有误差估计： 方程组为三对角线方程组，可以用追赶法求解。 若 不是 的线性函数时，所得方程组是非线性组，可以用解非线性方程组的方法求解。例如，可用简单迭代法、Newton迭代法求解. 例: 用导数的近似计算公式将下列边值问题化为线性方程组的求解,试写出用Seidel迭代法求解的算法(不具体求解) 在x1=0.1, x2=0.2,…, x9=0.9,处的数值解. 解: h=0.1, xj=jh. 设 yj=y(xj) , 由条件, y0=0, y10=1. 代入方程得. 于是, 可得方程组: 5.4 微分方程的定性分析 随着科学技术的发展，常微分方程定性分析 在各个学科领域已成为必不可少的数学工具， 也是数学建模的必备基础理论. 一. 微分方程定性理论的基本任务和主要研究方法 极少情况下，能够用初等函数或初等函数的积分表示微分方程的解. 求微分方程的数值解 解决方法 对微分方程进行定性分析 一般提法：不去积