用Matlab软件求常微分方程的数值解[t，x].PPTVIP

差分方程模型 在实际中许多变量是离散变化的，如人口.商品件数.生产周期等，而离散的运算具有可操作性，差分方程正是联系连续变量与离散的一座桥梁（如摩尔.库仑）。差分方程主要用来解决离散型问题。 设游泳者的速度大小和方向均不随时间变化，即令 而流速 ，其中 u 和 v 为常数, ? 为游泳者和x 轴正向间的夹角。 于是游泳者的路线 (x(t), y(t)) 满足 如果方程组有解，则有 由引理可知，最佳路径必为连接起点与终点的直线！所以， 若已知L, H, v, T (1).以H = 1160 m, L = 1000 m, v= 1.89 m/s 和 第一名成绩T=848 s 代入得 ? =117.46，u=1.54 m/s。 (2).以H = 1160 m, L = 1000 m, v= 1.89 m/s 和 u=1.5 m/s代入代入得 ? =121.85，T=910s，即15分10秒。 一定有解吗？ （3）能够到达的条件： 几何解释 如果游泳者始终以和岸边垂直的方向（y轴正向）游，则有u= vH/L≈2.19 m/s。 情形1：游泳者速度不可能这么快，因此永远游不到终点, 被冲到终点的下游去了。情形2：游泳者速度有可能这么快，因此可以游到终点。 男子 1500 米自由泳世界记录为 14分41秒66, 其平均速度为1.7 m/s。 查资料！ 能够成功到达终点的选手的速度限制 2002年： H = 1160 m, L = 1000 m, v= 1.89 m/s，需要u 1.43 m/s。 1934年： 竞渡的直线距离为5000 m, 垂直距离仍为H = 1160 m, 则L=4864 m,仍设v= 1.89 m/s， 游泳者的速度只要满足 u 0.44 m/s，就可以选到合适的角度游到终点.。(游 5000米很多人可以做到) (4) 由对称性及偏角引理 在本节我们将介绍差分方程的求解方法以及在实际中的应用。 差分方程是处理离散性问题的一种方法，在实际中有广泛的应用，它和微分方程是相辅相成的，对连续性数据的处理也可以采用离散的方法，这时也可以用差分方程。 １．差分方程的平衡点及其稳定性 设有未知序列 ，称 为k阶差分方程，若有 ，满足： 则称 是差分方程 的解 包含k个任意常数的解称为 的通解 为已知条件时，称其为 的初始条件 通解中的任意常数都有初始条件确定的解称为 的特解 形如： 称为k阶线形差分方程，其中 为已知的系数，且 若差分方程 中 的，则称差分方程 为k阶齐次线性差分方程，否则，称为k阶非齐次线性差分方程 若有常数a是差分方程的 解，即 则称a是差分方程的平衡点 若 已知，则形如 的差分方程的解可以在计算机上实现，下面给出一些理论上容易求解的特殊差分方程的解及其简单应用。 一阶常系数线形差分方程： 的通解为： 二阶常系数线形差分方程： 当 时，它有一特解 ，当 时，且 时，它有一特解 ，不管哪种情形， 是方程 的平衡点，设方程 的特征方程 其解为： （1）当 是两个不同实根时，方程 的通解为： （2）当 是两个相同实根时，方程 的通解为： （3）当 是一对共轭复根时， 方程 的通解为： 易知，当且仅当特征方程的任意特征根 时，平衡点 是稳定的 二阶方程的上述结果可以推广到k阶线形方程，即k阶线性方程平衡点稳定的条件是特征方程的根 均满足： （即均在复平面上的单位圆内） 差分方程的求解方法和微分方程相似，为什么会这样呢？？？ 下面仅简单介绍一阶非线性差分方程： 的平衡点的稳定性 其平衡点 由代数方程 解出，为分析其平衡点的稳定性，将方程 的右端 在 点作