习题一: 1.下列各数都是经过四舍五入得到的近似数,试指出它们有几位有效数字以及它们的绝对误差限,相对误差限。 ⑵ ⑷ ⑹ 解:四舍五入得到的近似数的误差都不超过其末位的半个单位,从其最后一位数字开始到前面第一个非零数字为止的所有数字,均是有效数字。 故有效数字的位数分别为2,5,2 由有效数字与绝对误差限的关系:   分别为-1,2,3 解得绝对误差限分别为:    由有效数字与相对误差限的关系:    2.为使下列各数的近似值的相对误差限不超过,问各近似值分别应取几位有效数字?  解:由定理1知,。 由于,所以。 欲使 解得: 至少应取,即3位有效数字。 4.计算,取,利用下列等价表达式计算,哪一个效果最好?为什么? (1); (2); (3); (4) 解:第(3)个计算的结果最好。 由原则三:“避免两相近的数相减”,可排除(2)、(4), 由原则五:“简化计算步骤,减少运算次数”,可排除(1)。 7.利用等式变换使下列表达式的结果比较精确。 (1), ; (2), ; (3), ; (4),。 解:(1)(方法一)(方法二) (2)  (3)  (4)  习题二: 2.方程在区间中有一实根,若用二分法求此根,使其误差不超过,问应将区间对分几次?并用二分法求此根。 解:由题意,。若要,则  解得,故取,即将区间二分六次即可。 n      0 3(-) 4(+) 3.5 +  1 3 3.5 3.25 -  2 3.25 3.5 3.375 -  3 3.375 3.5 3.4375 +  4 3.375 3.4375 3.40625 +  5 3.375 3.40625 3.390625 -  6 3.390625 3.40625 3.3984375 -   7.下面是求的两个迭代格式: (1) (2) 求它们的收敛阶。 解 (1) (*1)  是的根。 我们将(*1)式进行变形,得 方程两边对进行求导,得  再对上式求导,得  故所给迭代格式为二阶收敛。 (2) (*2) , 是的根。 我们将(*2)式进行变形,得  方程两边对进行求导,得   对上式求导,得   再对上式求导,得  故所给迭代格式为三阶收敛。 8.设充分接近方程的某个根,给定迭代函数,其中   (1) 试证至少有三阶收敛速度。 证明:因为初值充分接近根,所以我们可以用定理3来证明。 是方程的根,。并且  我们将(1)式进行变形,得  对上式求一阶导数  解得: 则有。 对上式再求导数  解得: 则有。 故所给迭代格式至少三阶收敛。 10.设,应用牛顿迭代法分别求与之根,从而导出的两种迭代格式,并求,取,。 解 (1) 其牛顿迭代格式为: 。 (2) 其牛顿迭代格式为:  在上述两种迭代格式中分别取,,,并且要求,由此求得的近似值(4.21716421)。 习题三 2.用列主元素法解下列方程组 (1) (2) 解 (1)考查方程的增广矩阵   回代后求解得  (2)考查方程的增广矩阵   回代后求解得  观察可得:上述两个方程的系数矩阵都是严格对角占优矩阵,所以在用列主元素法求解方程时,选主元与不选主元的结果是一样的。 3.用追赶法解下列方程组。 (1)  解:对方程组的增广矩阵进行顺序消元,且每步都将主元系数化为1,  5.求第3题中系数矩阵的分解,并用此分解法解对应的线性方程组。 解   用紧凑格式法: 故  求解方程组,得  求解方程组,得  7.给定,,求及。 解:       解得:   13.设,计算的条件数。 解:  ,。 故 由于,又因为是对称矩阵,所以 ,  由 解得,故。 同理可得 故。 15.对下列方程组考察用雅可比迭代法与高斯-塞德尔迭代法是否收敛?若收敛,写出其迭代格式;若不收敛,能不能将方程组变形,使之用雅可比迭代法与高斯-塞德尔迭代法时收敛? (2) 解:① 解得: ,故Jacobi法发散。 ② 解得: ,故G-S法发散。 如果将方程组变形为: 这时方程组的系数矩阵为严格对角占优的矩阵,故Jacobi法和G-S法都收敛。 (4) 解:① 当时,带入上式, 得  当时,带入上式, 得  说明在之间有一根,则 故Jacobi法发散。 ②因为系数矩阵是对称矩阵,而  即系数矩阵是对称正定的矩阵,G-S法收敛。 迭代格式:  16.写出方程组 的SOR迭代格式。 解:先写出方程组的G-S迭代格式:  则该方程组的SOR迭代格式为:  习题四 1.已知,求的二次插值多项式。 解:  3.给出函数的数表如下,分别用线性插值与二次插值求的近似值,并估计截断误差。 解:线性插值: 内插的效果比外推的要好, 选取节点,,则   因为,所以   二次插值: 选 则: 得:   9.根据函数值表求四次牛顿插值多项式,并用其计算和的近似值。  1.615 1.634 1.702 1.828 1.921   2.41450 2.46459 2.65271 3.03035 3.34066  解:构造差商表: 1.615 2.41450      1.634 2.46459 2.636316     1.702 2.65271 2.738046 1.496029    1.828 3.03035 2.891315 1.314428 -1.441278   1.921 3.34066 3.026667 1.360108 -0.620644 8.824022    11.已知函数的函数值表,解答下列问题。  0.0 0.1 0.2 0.3 0.4 0.5   1.00 1.32 1.68 2.08 2.52 3.00  (1)试列出相应的差分表, (2)分别写出牛顿向前插值公式和牛顿向后插值公式。 解: (1)构造差分表: 0.0 1.00       0.1 1.32 0.32      0.2 1.68 0.36 0.04     0.3 2.08 0.40 0.04 0    0.4 2.52 0.44 0.04 0 0   0.5 3.00 0.48 0.04 0 0 0  (2)牛顿向前插值公式:  习题六: 3.确定下列求积公式中的代定参数,使其代数精度尽量高,并指明所构造出的求积公式具有的代数精度。 (1) 解:要使代数精度尽量高,则求积公式应对于精确成立  即:  即: 解得: 此时求积公式为:  当时,, , 左边等于右边。故求积公式对于也精确成立。 当时, ,左边不等于右边, 故求积公式对于不精确成立。 所以,当时,求积公式的代数精度最高,代数精度为3。 6.分别用复化梯形法和复化辛甫生法计算积分:,8等分积分区间。 解:函数在各个等分点的函数值分别为:  0      0 0.031128 0.061538 0.090566      1  0.117647 0.142349 0.164384 0.183607 0.2  复化梯形公式:    习题七: 取步长,用欧拉法和改进的欧拉法解下列初值问题。  解:欧拉法: 欧拉公式对于该初值问题的具体格式为:  ① ② ③ ④ 改进的欧拉法对于该初值问题的具体格式:  ①   ②     ③    ④