求数值近似开平方算法.docVIP

求数值近似开平方算法 200721101008 陈林 摘要：本文分别用牛顿迭代法、蒙特卡洛实验法和样条插值方法求数值的近似开方值，并用数值3开平方为算例，简单比较了三种方法的优劣。 牛顿迭代法 1.原理分析：牛顿迭代法实际上是由不动点迭代法的原理结合泰勒展开式构造迭代公式，所以迭代初值必须满足迭代收敛条件，设，则，由初值条件解得，不妨设，，则牛顿迭代公式为： 。 2.实验程序及注释： x0=input(input x0=); %设定迭代初值，本例为1 er=1; n=0; %定义一个计数器记录循环次数 while er0.001 x=0.5*(x0+3/x0); %牛顿迭代公式 er=abs(x-x0); x0=x; %数值解 n=n+1; end x sqrt3=sqrt(3) %真实数值 n er clear %清除内存变量 3.实验结果及分析： x0 x sqrt3 n er 1 1.7321 1.7321 4 9.2047e-005 因为计算机的截断误差和matlab数值显示格式的原因，显示出来的数值解和真实数值是一样的，但从误差变量er的值我们可以看到其计算的误差。误差在小数点后的第5位，保证了5位的有效数字，即1.7321的迭代结果是准确的。迭代4次就达到了要求的精度，说明了牛顿迭代法的收敛速度很快，可以证明它的平方收敛