数值计算与MATLAB方法课后答案.doc

第一章习题 1. 序列 满足递推关系 ，取 及 试分别计算 ，从而说明递推公式对于计算是不稳定的。 ? n 1 1 0.01 0.0001 2 0.01 0.0001 0.000001 3 0.0001 0.000001 0 4 0.000001 0 10-10 5 0 10-10 ? ? n 1 1.000001 0.01 0.000099 2 0.01 0.000099 -0 3 0.000099 -0 -0 4 -0 -0 -1.0001 5 -0 -1.0001 ? 初始相差不大，而 却相差那么远，计算是不稳定的。 2. 取y0=28，按递推公式 ，去计算y100，若取 （五位有效数字），试问计算y100将有多大误差？y100中尚留有几位有效数字？ 解：每递推一次有误差 因此，尚留有二位有效数字。 ? 3．函数 ，求f(30)的值。若开方用六位函数表，问求对数时误差有多大？若改用另一等价公式 计算，求对数时误差有多大？ 设z=ln(30-y)， ，y*=29.9833， |E(y)| 0.5(10-4 z*=ln(30-y*)=ln(0.0167)=-4.09235 若改用等价公式 设z=-ln(30+y)， ，y*=29.9833， |E(y)| 0.5(10-4 z*=-ln(30+y*)=-ln(59.9833)=-4.09407 ? 4．下列各数都按有效数字给出，试估计f的绝对误差限和相对误差限。 1)????? f=sin[(3.14)(2.685)] 设f=sin xy x*=3.14, E(x) 0.5(10-2, y*=2.685, E(y) 0.5(10-3, sin(x*y*)=0.838147484, cos(x*y*)=-0.545443667 |2.685((-0.5454) (0.5(10-2+3.14(-0.5454) (0.5(10-3| =0.81783(10-20.82(10-2 |Er(f)| 0.82(10-2/0.83815 0.9894(10-210-2 ? 2)????? f=(1.56)3.414 设f = xy , x*=1.56, E(x) 0.5(10-2, y*=3.414, E(y) 0.5(10-3, |3.414(1.562.414(0.5(10-2+1.563.414(ln1.56(0.5(10-3| =|3.414(2.925518039(0.5(10-2 +4.563808142(0.444685821(0.5(10-3| =0.051 |Er(f)| 0.051/1.563.414=0.0112 ? 5．计算 ，利用下列等式计算，哪一个得到的结果最好，为什么？ ? 6．下列各式怎样计算才能减少误差？ ? 7. 求方程x2-56x+1=0的二个根，问要使它们具有四位有效数字， 至少要取几位有效数字？如果利用伟达定理，(又该取几位有效数字呢？ 解一: 若要取到四位有效数字, 如果利用伟达定理, 解二: 由定理二,欲使x1,x2有四位有效数字,必须使 由定理一知，(至少要取 7位有效数字。 如果利用伟达定理, 由定理一知，(至少要取 4位有效数字。 ? 8．求近似数x* 的幂(x*)n 的相对误差公式。 解1：y*=(x*)n, dy*=d(x*)n=n(x*)n-1dx*, E(y)( n(x*)n-1E(x) 解2: Er(y) (d lny=d lnxn = nd lnx ( nEr(x) ? 9．设y=lgx，证明x* 的常用对数的绝对误差不超过x* 的相对误差的一半，反之，y*的反对数的相对误差不超过y* 的绝对误差的2.5倍。 ? 10．简化求多项式 的值的运算式。 第二章?????? 非线性方程求根 习题2-1 1. 试寻找f(x)= x 3+6.6 x 2-29.05 x +22.64=0的实根上下界,及正根 所在的区间，区间长度取