第六讲方程求解.docVIP

第六讲 方程求解 线性方程组解法 消去法 迭代法 非线性方程解法 主要有因式分解法、图形放大法、数值迭代逼近法，这里介绍数值迭代逼近法。 求实根的对分区间法 设是区间上的连续函数，已知在中有的一个根，为具体起见，不妨假设，，则根据连续函数的性质，方程在内必有一实根。求此实根的对分区间法的步骤如下： 计算在区间中点的值， 判断：若，则计算停止。否则若与异号，则根位于中，以代替；若与异号，则根位于中，以代替。 重复步骤（1）和（2），直到区间缩小到容许误差范围之内。此时，区间中点即可作为所要求的根。 迭代法 若给定方程，把它改写成下面等价的形式（要求为连续函数），例如，可令，再做迭代 当给定初值以后，就计算出一个序列，如果序列收敛，，则有 也就是说是方程的根，也叫做的不动点。 什么情况能保证迭代序列的收敛呢？ 定理：把方程改写成，如果满足Lipschitz条件：即对任意的和都有： ， 其中，是与，无关的正常数，称为Lipschitz常数（简称李氏常数）。若，则迭代 收敛到准确解，且有估计 一般地，给定充分小的正数，如果则停止计算。 越接近于0，收敛越快。 迭代收敛的加速 如能获得较小的，可使迭代速度更快。 一种加速迭代过程： 牛顿法（切线法） 解非线性方程的牛顿法是把非线性方程线性化的一种近似方法。把在点附近展开成泰勒级数 取其线性部分作为非线性方程的近似方程，则有 设，则其解为 再把在点附近展开成泰勒级数，类似可得 这样，得到牛顿法的一个迭代序列 牛顿法有十分明显的几何意义，当选取初值以后，过作的切线，其切线方程为 求此切线方程和轴的交点，即得 因此，此法也叫切线法。 割线法 选定曲线上的两个点和，过这两点作一条直线，则直线方程为 当时，设直线与轴的交点是，则 用作为曲线与轴的交点的近似值，可以预期比，更接近于解，重复上述过程，其迭代公式为： 三、非线性方程组解法 牛顿迭代法 解非线性方程的牛顿法是把非线性问题线性化，对非线性问题线性化的方法也可用来解决方程组的问题。 以一个二阶方程组为例来介绍， 设此方程组的一组初始近似值为，把和都在附近用二元泰勒展开，并取其线性部分，得到方程组 只要系数矩阵的行列式 则方程组的解可以写成 四、MATLAB实现及引例求解 命令: 1、解任意函数方程和方程组 solve, 2、解多项式方程roots(p), 3、解线性方程组 sym(A)\sym(b),A\b,inv(A)*b,pinv(A*b) 解非线性方程组 fsolve 五、范例、实验 3