matlab实验报告求代数方程的近似根.doc

数学实验报告 实验序号： 第二次 日期：2012 年 5月10日 班级 0920861 小组成员姓名 徐易斌；王勇 王康 学号 30 12 33 实验名称：求代数方程的近似根 问题背景描述： 求代数方程的根是最常见的数学问题之一，当是一次多项式时，称为线性方程，否则称之为非线性方程． 当是非线性方程时，由于的多样性，尚无一般的解析解法可使用，但如果对任意的精度要求，能求出方程的近似根，则可以认为求根的计算问题已经解决，至少能满足实际要求． 本实验介绍一些求方程实根的近似值的有效方法，要求在使用这些方法前先确定求根区间，或给出某根的近似值． 实验目的： 了解代数方程求根求解的四种方法：对分法、迭代法、牛顿切线法 掌握对分法、迭代法、牛顿切线法求方程近似根的基本过程。 实验原理与数学模型： 对分法 对分法思想：将区域不断对分，判断根在某个分段内，再对该段对分，依此类推，直到满足精度为止．对分法适用于求有根区间内的单实根或奇重实根． 设在上连续，，即 ，或，．则根据连续函数的介值定理，在内至少存在一点 ，使． 下面的方法可以求出该根： 令，计算； 若，则是的根，停止计算，输出结果． 若 ，则令，，若，则令，；． ……，有、以及相应的． (3) 若 (为预先给定的精度要求)，退出计算，输出结果； 反之，返回(1)，重复(1)，(2)，(3)． 以上方法可得到每次缩小一半的区间序列，在中含有方程的根． 当区间长很小时，取其中点为根的近似值，显然有 以上公式可用于估计对分次数． 2. 迭代法 迭代法的基本思想： 由方程构造一个等价方程 从某个近似根出发，令 ，，这种方法称为迭代法． 若 收敛，即 ， 只要连续，有 即 可知，的极限是的根，也就是的根． 当然，若发散，迭代法就失败． 迭代过程收敛的常用判别标准： 当根区间较小，且对某一，明显小于1时，则迭代收敛 2) 迭代法的加速： a) 松弛法： 若与同是的近似值，则是两个近似值的加权平均，其中称为权重，现通过确定看能否得到加速． 迭代方程是： 其中，令，试确定： 当时，有，即当，时， 可望获得较好的加速效果，于是有松弛法：， b) Altken方法： ，是它的根，是其近似根． 设，，因为 ， 用差商近似代替，有 ， 解出，得 由此得出公式 ； ； ， 3. 牛顿(Newton)法（牛顿切线法） 1) 牛顿法的基本思想： 是非线性方程，一般较难解决，多采用线性化方法． 记： 是一次多项式，用作为的近似方程． 的解为 记为，一般地，记 即为牛顿法公式。 实验所用软件及版本： Matlab 7.0 主要内容（要点）： 分别用对分法、普通迭代法、松弛迭代法、Altken 迭代法、牛顿切法线等5种方法，求方程 的正的近似根，．(建议取 ．) 实验过程记录（含基本步骤、主要程序清单及异常情况记录等）： 1.对分法syms x fx; a=0.001;b=3; fx=0.5*x-sin(x); x=(a+b)/2;k=0; ffx=subs(fx,x,x); if ffx==0; disp([the root is:,num2str(x)]) else disp(k ak bk f(xk)) while abs(ffx)0.0001ab; disp([num2str(k), ,num2str(a), ,num2str(b), ,num2str(ffx)]) fa=subs(fx,x,a);ffx=subs(fx,x,x); if fa*ffx0 b=x; else a=x; end k=k+1;x=(a+b)/2; end disp([num2str(k), ,num2str(a), ,num2str(b), ,num2str(ffx)]) end fprintf(所求的解是：x=%f,迭代步数是：%d/n,x,k) 【调试结果】 0 0.001 3 -0.24728 1 1.5005 3 -0.24728 2 1.5005 2.2502 0.34721 3 1.8754 2.2502 -0.016286 4 1.8754 2.0628 0.15002 5 1.8754 1.9691 0.062824 6 1.8754 1.9222 0.022239 7 1.8754 1.8988 0.0027165 8 1.8871 1.8988 -0.0068499 9 1.8929 1.8988 -0.0