非线性方程的根.doc

非线性方程的根 专业班级：信计111 （1）实验目的： 熟悉Matlab编程； 学习非线性方程求根的方法及程序设计算法 （2）实验题目： 1、迭代函数对收敛性的影响 用迭代法求方程的根 方案1：化方程为等价方程取初值，迭代10次。 方案2;化为等价方程 ,取初值，迭代10次，观察其计算值，并加以分析。 2、初值的选取对迭代法的影响 用牛顿法求方程 在 附近的根。 方案1：使用牛顿法并取 ，由得 迭代10次。 方案2：取 ，使用同样的公式迭代10次，观察比较并分析原因。 3.收敛性与收敛速度的比较 求方程 的全部实根, 方案1：用牛顿法求解； 方案2：用简单迭代法； 取相同迭代初值，比较各方法的收敛速度。 (3)实验内容 1. format long x0=0;y0=0; syms x fx gx fx=((x+1)/2)^(1/3); gx=(x^3)*2-1; j=0; while(j10) fx0=((x0+1)/2)^(1/3); gy0=(y0^3)*2-1; x0=fx0; y0=gy0; j=j+1; x0 y0 end 得到结果x0 = 0.793700525984100 y0 = -1 x0 =0.964361757887056 y0 = -3 x0 =0.994024659401817 y0 = -55 x0 = 0.999003116453725 y0 =-332751 x0 =0.999833825129730 y0 = -7.368652968112150e+016 x0 =0.999972303421195 y0 = -8.001921866539816e+050 x0 =0.999995383882224 y0 =-1.024738174056895e+153 x0 = 0.999999230646445 y0 =-Inf x0 =0.999999871774391 y0 =-Inf x0 =0.999999978629065 y0 =-Inf 2. 初值的选取对迭代法的影响 format long x0=0;y0=1.5; syms x fx fx=x-(x.^3-x-1)/(x.^2*3-1); j=0; while(j10) fx0=x0-(x0.^3-x0-1)/(x0.^2*3-1); x0=fx0 fy0=y0-(y0.^3-y0-1)/(y0.^2*3-1); y0=fy0 j=j+1; end 实验结果 x0 = -1 y0 = 1.347826086956522 x0 =-0.500000000000000 y0 = 1.325200398950907 x0 = -3 y0 = 1.324718173999054 x0 = -2.038461538461538 y0 =1.324717957244790 x0 =-1.390282147216736 y0 = 1.324717957244746 x0 =-0.911611897717927 y0 = 1.324717957244746 x0 =-0.345028496748169 y0 = 1.324717957244746 x0 = -1.427750704027270 y0 =1.324717957244746 x0 =-0.942417912509483 y0 =1.324717957244746 x0 = -0.404949357199379 y0 =1.324717957244746 3.收敛性与收敛速度的比较 1.普通迭代，选用初值0.5 format long c=10.^(-6); x0=0.5; syms x fx fx=(x.^2-exp(x)+2)/3; j=1; fx0=(x0.^2-exp(x0)+2)/3 while(abs((x0-fx0)/fx0)c) x0=fx0; fx0=(x0.^2-exp(x0)+2)/3 j=j+1; end fx0 = 0.200426243099957 fx0 = 0.272749065098375 fx0 = 0.253607156584130 fx0 = 0.258550376264936 x0 = 0.257265636335094 fx0 = 0.257598985162190 fx0 = 0.257512454514832 fx0 = 0.257534913615251 fx0 = 0.257529084167956 fx0 =