数模中的数值解法.pptVIP

2：常微分方程数值解决 * 1:非线性方程及非线性方程组解决 它等价于:　 非线性方程：设 为非线性函数，方程 称为非线性方程。 重根：若 则称 是方程的 重根或 的 重零点． 方法1 对分区间法 逐次对分区间得区间套： 其中 ，即 均为有根区间，则 是方程（8-1）的根。 若 ，则 内必有方程（8-1）的根，称其为有根区间。 其误差为: 1. 确定有根区间: 2. 确定更小的有根区间: 取 3. 取根的近似值: 用对分区间法求根步骤： ２．逐次对分有根区间求近似根 若 ，则 ；否则 直到满足精度。 1. 确定求根区间 误差公式还可以用于事先估计等分区间的次数 解：∵　 ∴ 为有根区间 故 ， 即为满足精度要求的近似解． 例： 求方程 的一实根。要求误差不超过 （即精确到四位小数）。 方法2　简单迭代法 迭代法（单步）的一般形式： 称为迭代函数。 构造迭代格式: 出发，逐次代入迭代公式产生序列 ，以 近似非线性方程的根. 从初值 （1） 简单迭代法的一般形式及其几何意义　　　　 迭代格式：　　 结论：若 在 处连续，则 是方程的根，即 步骤：取初值 ，逐次代入上式产生序列 ，以 为近似方程 的解，称为简单迭代法。 x y y = x x y y = x x y y = x x y y = x x* x* x* x* y= g(x) y=g(x) y=g(x) y=g(x) x0 p0 x1 p1 ? x0 p0 x1 p1 ? x0 p0 x1 p1 ? x0 p0 x1 p1 ? 　　 　　　　 例2　用简单迭代法求方程 的根，要求精确到六位小数，取 解 （1） 迭代格式　　　 将初值代入得 显然，此迭代序列发散。 迭代格式:　　　 　　将初值代入得 : ∴ （2） 　　迭代格式: 　　代入初值得： 例2的结果表明，对同一方程可构造不同的迭代格式， 产生的迭代序列收敛性也不同。迭代序列的收敛性 取决于迭代函数在方程的根的邻近的性态。 （3） （2）迭代法的收敛条件 （2）存在常数 ，使得对一切 都有 则方程 在 内有唯一的根 且对任意初值 迭代格式： 所产生的序列 均收敛于 定理8.1（压缩映象原理） ，都有 （1）对任意 在区间 上满足条件： 设函数 证明： （2）唯一性： （1）存在性： 定理8.2　 （2）存在常数 ，使得对 ，都有 则方程 在 上有唯一的根 且对 迭代序列 均收敛于 且 （1） （3）　Newton法 当 接近 时，由Taylor公式 若 ，则有 基本思想：将非线性方程线性化，以线性方程的解 逼近非线性方程的解。 连续可导, 是方程 设 的根， 由此导出迭代公式 按此迭代公式求方程 的近似解称为Newton法。 其几何意义：以曲线 在点 的切线与 近似曲线与 轴的交点，故Newton法又称切线法。 轴的交点 x y x* x0 例：用Newton法求方程 的根，取 解： 代入初值得： 迭代公式为 (1) 用差商近似导数 微分方程离散化的方法 若用向后差商近似导数，即 （2）用Taylor多项式近似 §1　Euler方法 Euler方法 的解作为微分方程初值问题的数值解，即 称为Euler方法。 以差分方程初值问题 解：Euler公式为 若取 ，则有 用Euler方法求初值问题 的数值解。 若以 求数值解，称为向后Euler法，这是隐式公式， 一般需用迭代法求解。 §2 改进Euler方法 2.1 梯形公式 求解公式 *