数值分析上机实践.docVIP

1、分别用牛顿法，及牛顿-Steffensen法 求ln(x+sinx)=0的根。初值x0分别取0.1, 1,1.5, 2, 4进行计算。 求sinx=0的根。初值x0分别取1，1.4，1.6, 1.8，3进行计算。 分析其中遇到的现象与问题。 （1）、newton法 (1)、对应M文件为newton1如下： function x=newton1(x0) x=x0-((x0+sin(x0))*log(x0+sin(x0)))/(1+cos(x0)); while (x-x00)|(x-x00) x0=x; x=x0-((x0+sin(x0))*log(x0+sin(x0)))/(1+cos(x0)); end newton1(0.1) ans = 0.5110 newton1(1) ans = 0.5110 1.5、2.4作为初值不适合，计算结果出现复数。 (2)对应M文件为newton2如下： function x=newton2(x0) x=x0-sin(x0)/cos(x0); while (x-x00)|(x-x00) x0=x; x=x0-sin(x0)/cos(x0); end newton2(1) ans = 0 newton2(1.4) ans = 3.1416 newton2(1.6) ans = 31.4159 newton2(1.8) ans = 6.2832 newton2(3) ans = 3.1416 newton2(3.3) ans = 3.1416 newton2(3.9) ans = 3.1416 初始值取值不同，有可能造成迭代公式发散。计算不出正确值，或是计算得出错误的值。 （2）、newton-steffensen加速法 (1)、对应的M文件为NS111如下： function x=NS111(x0) y0=x0-((x0+sin(x0))*log(x0+sin(x0)))/(1+cos(x0)); z0=y0-((y0+sin(y0))*log(y0+sin(y0)))/(1+cos(y0)); x=x0-(y0-x0)*(y0-z0)/(z0-2*y0+x0); while (abs(x-x0)0.00001) x0=x; y0=x0-((x0+sin(x0))*log(x0+sin(x0)))/(1+cos(x0)); z0=y0-((y0+sin(y0))*log(y0+sin(y0)))/(1+cos(y0)); x=x0-(y0-x0)*(y0-z0)/(z0-2*y0+x0); end x=x0 NS111(0.1) x = -10.7295 - 3.1908i ans = -10.7295 - 3.1908i NS111(1) x = 0.5140 ans = 0.5140 NS111(1.5) x = 0.5140 - 0.0000i ans = 0.5140 - 0.0000i NS111(2.4) x = 2.9105e+002 -4.0048e-001i ans = 2.9105e+002 -4.0048e-001i （2）、对应M文件为NS2如下： function x=NS2(x0) y0=x0-sin(x0)/cos(x0); z0=y0-sin(y0)/cos(y0); x=x0-(y0-x0)*(y0-z0)/(z0-2*y0+x0); while (abs(x-x0)0.00001) x0=x; y0=x0-sin(x0)/cos(x0); z0=y0-sin(y0)/cos(y0); x=x0-(y0-x0)*(y0-z0)/(z0-2*y0+x0); end x=x0 NS2(1) x = 0.0311 ans = 0.0311 NS2(1.4) x = -0.0311 ans = -0.0311 NS2(1.6) x = 3.1727 ans = 3.1727 NS2(1.8) x = 9.3937 ans = 9.3937 NS2(3) x =