数值分析试卷习题5.docVIP

第五章 习题解答与问题 一、习题解答 1．求经过A(0，1)，B(1，2)，C(2，3)三个样点的插值多项式 解：令x0 = 0，x1 = 1，x2 = 2，则有f(x0)=1，f(x1)=2，f(x2)=3，由Lagrange二次插值公式 = x+1 2．已知函数的数据如下表 x 0 1 2 3 y 1 3 9 27 试作一个三次插值多项式P3(x),利用P3(x)计算 解：令xk = k (k = 0，1，2，3)，则根据函数表有f(xk)= 3k。构造差商表 x f(x) 0 1 1 3 2 2 9 6 2 3 27 18 6 4/3 根据Newton插值公式 由于被插值函数，故取 x = 1/2，便得 3．已知函数y = f(x)的数据如下表 x – 1 0 1 y -1 0 1 y’ 0 解：由于x=0是二重零点，令H3(x)= x2(a x + b )。又由H3(-1)= - 1，H3(1)=1得方程组 解之：a =1，b = 0 所以，H3(x) = x3。 4．设被插值函数f(x)在区间[x0，x1]上具有2阶连续导数，求证：两点线性插值函数L(x)的误差界满足不等式 证：由拉格朗日插值误差定理，得 令h(x) = (x – x0)(x – x1)，求导数并令其为零，可得极值点x*=0.5×(x0 + x1)。显然这是函数|h(x)|在所讨论区间内的最大值点，故 5．若给出 sin x 在区间[0，(]上的等距结点的数值表，用线性插值计算表中没有的数据时，希望截断误差小于5×10-3，问函数表中自变量数据的步长应取多少为好？ 解：设应取的步长为h，利用第4题的结论，只须 所以，应取h = 0.2。 6．常用对数表中，只能查出 lgN 的值，其中N 是具有形式：a1a2a3.a4 的实数。试研究线性插值计算近似值的误差界。 解：由N 的取值形式可知，对数表中函数自变量取值的步长为 h = 0.1。而 由第4题结论，线性插值计算近似值的误差界 7．证明两点三次Hermite插值余项是 并由此求出分段三次Hermite插值的误差限. 证明：由插值条件知 R(xk)=R’(xk)=0, R(xk+1)=R’(xk+1)=0 取x异于xk和xk+1，(x) = f(x) – H(x)=C(x)(x – xk)2(x – xk+1)2 构造辅助函数， 显然，F(t)有三个零点xk, x, xk+1,由Roll定理知,存在F’(t)的两个零点t0,t1满足xkt0t1xk+1,而xk和xk+1也是F’(x)的零点,故F’(x)至少有四个相异零点. 反复应用Roll定理,得F(4)(t)至少有一个零点设为ξ。由此得 所以，，故 对于等距离节点（步长取为h），分段Hermite插值的误差余项为 8．求一个次数不高于4次的多项式P(x),使它满足:P(0) =0，P’(0) =0，P(1) =1，P’(1) =1，P(2) =1，并写出其余项表达式. 解：由题意 P(x) = x2(ax2 + b x + c )，由插值条件得方程组 求解，得 a =1/4，b= – 3/2 ，c =9/4。所以 插值余项为 9．设，在-5≤x≤5上取 n = 10，按等距结点求分段线性插值函数Ih(x)，计算各结点间中点处的Ih(x)和f(x) 的值，并估计误差。 解：因为 f(xk) = 1/(1+ k2)，(k = -5，…，-1，0，1，…，5)分段线性插值函数为 其中，分段线性插值的基函数为 对于某一节点中点xk0 = (xk + xk+1)/2= k+ 0.5，于是 ， 因为 所以 10．已知二元函数u(x，y)在四个点P1(x1，y1)，P2(x2，y1)，P3(x2，y2)，P4(x1，y2)处的函数值为u1，u2，u3，u4，试推导双线性插值函数 uh(x，y) = a + b x + c y + d xy 的表达式。 解：令，，则利用基函数得 uh(w，v) =(1 – w)(1 – v)u1 + w (1 – v)u2 + w v u3 + (1 – w) v u4 所以，有 二、例题 在代数插值问题中，x0，x1，……，xn是(n+1)个互异的插值结点，由这(n+1)个结点构造拉格朗日插值基函数 l0(x)，l1(x)，……，ln(x) 试证明恒等式：l0(x) + l1(x) + …… + ln(x) = 1。 证明：记P(x)= l0(x) + l1(x) + …… + ln(x) – 1，由于拉格朗日基函数在插值结点处满足 故对所有插值结点有 l0(xk)