2015各种计算方法.pptVIP

各种计算方法 信息学院 熊波 算法重在设计 所谓算法就是计算机上使用的计算方法 众所周知，行列式解法的Gramer(格莱姆)法则原则上可以用来求解线性方程组。用这种方法求解一个n阶方程组，要算n+1个n阶行列式的值,这个计算量是相当惊人的！ 一 方程求根的迭代方法 二 线性方程组的迭代法 * * 比如一个20阶的方程组，用巨型计算机来计算，要连续工作100年，这是完全没有实际意义的，而运用消元技术，一个20阶的线性方程组用普通的计算器也能很快地解出来。因此合理地选择算法是科学计算成败的关键 数学思维的化归策略 算法设计的基础是数学思维 数学家的思维特征：将问题加工变形，直到把它转化归纳为能够解决的问题——化归策略 例1：提供水龙头、水壶、煤气灶与火柴，怎样烧开水？ 烧开水程序分三步 打开水龙头将水壶灌满 用火柴点燃煤气灶 把水壶放在煤气灶上 问题2： 所给条件与1相同，只是水壶已经装满了水，这是该怎样烧开水呢？ 根据生活经验，这个问题比1简单，只用执行2、3步即可 但是数学家的思维方式，解决2更“合理”的方案是：先倒掉壶中的水，将问题2化归为问题1，再按1的方案进行 化归策略是算法设计的基本策略 1 求开方值的迭代公式 相对于加减乘除四则运算来说，开方运算无疑是复杂的。人们自然希望将复杂的开方运算化归为四则运算 给定 ，求开方值 的问题转化为解方程 设给定某个预报值 ，希望借助于某种简单方法确定校正量 ，使校正值 能够比较准确地满足所给方程 假设 是个小量，为简化计算舍去高阶小量 ，令 这是一个关于 的一次方程，据此求出 ，从而得到 反复施行上述过程，即可导出开方公式 从给定的某个初值 出发，利用上式反复迭代，即可获 得满足精度要求的开方值 算法：任给 ，对 ，执行算式（1），直 到偏差 为止（ 为给定精度），最终获得的近 似值 即为所求 始 读入 输出 终 例：用开方算法求 ，设取 解： 的准确值为 ，这里迭代5次得到 的结果 开方公式规定了预报值 与校正值 之间的一种函数关 系 ，这里 称为开方法的迭代函数 开方法的收敛性： 所谓收敛性指：迭代误差 当 时趋近于0 （证明略） 定理：开方法对任意给定初值 均收敛 2 NEWTON法 考虑一般形式的函数方程 设已知它的近似根 ，则自然要求校正值 能 更好的满足所给方程(2) 取其线性主部 ，令 从而关于校正值 有如下计算公式： 这就是著名的NEWTON公式 NEWTON迭代的思想是： 将非线性方程的求根过程线性化 NEWTON公式确定了预报值 与校正值 之间的一种函 数关系 ，这里迭代函数为 Newton法的几何解释： 方程 的根 在几何上解释 为曲线 与 轴交点的横坐 标。设 是根 的某个近似值， 对曲线 上横坐标为 的点 引切线，设该切线与 轴交点 的横坐标记为 ，则这样获得的 即为按Newton法求得的近似 根。由于这种几何背景Newton法 亦称为切线法 Newton法的计算流程 始 输入 输出迭代失败标志 终 输出 输出奇异标志 开方法是Newton法的特殊情况 求开方值 就是求解方程 Newton法的迭代函数为 就是开方法的迭代函数 定义：如果迭代误差 当 时成立 则称迭代过程是 阶收敛的 压缩映象原理 对于一般方程 ，为要使用迭代法，先将它改写成 的形式，其中