数值计算与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(x0.5′10-2, y*=2.685, E(y0.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.838150.9894′10-210-2 ? 2????? f=(1.563.414 设f = xy , x*=1.56, E(x0.5′10-2, y*=3.414, E(y0.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的二个根，问要使它们具有四位有效数字， 至少要取几位有效数字？如果利用伟达定理，D又该取几位有效数字呢？ 解一: 若要取到四位有效数字, 如果利用伟达定理, 解二: 由定理二,欲使x1,x2有四位有效数字,必须使 由定理一知，D至少要取7位有效数字。 如果利用伟达定理, 由定理一知，D至少要取 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 x2-29.05 x +22.64=0的实根上下界,及正根 所在的区间，区间长度取1。 解:由笛卡儿符号规则知，f(x=0可能有二个正根或无正根 f(-x= -x 3+6.6 x2+29.05 x +22.64=0 即x 3 -6.6 x2-29.05 x-22.64=0 f(-x=0有一个正根，因此，f(x=0有一个负根。 由定理2-3， f(x=0的正根上界 f(x=0的负根下界 x 0 1 2 3 4 5 6 6.39 f(x + + - + + + +