重节点差商MATLAB函数,第一章习题解答-中国石油大学.doc

第一章习题解答-中国石油大学

《数值分析》习题集

数值分析教学团队 编

中国石油大学(北京)理学院

第一章习题解答

在下列各对数中，x 是精确值 a 的近似值。

试估计x的绝对误差和相对误差。

解：(1)

(2)

(3)

(4)

2、已知四个数：。试估计各近似数的有效位数和误差限，并估计运算和的相对误差限。

解：

由相对误差限公式：

所以有：

3、设精确数ａ>0，x是ａ的近似值，x的相对误差限是0.2，求㏑x的相对误差限。

解：设

4、长方体的长宽高分别为50cm，20cm和10cm，试求测量误差满足什么条件时其表面积的误差不超过1cm2。

解：设

所以，测量误差小于0.00625时其表面积的误差不超过1cm2。

5、设x和y的相对误差为0.001，则xy的相对误差约为多少？

解：由公式：

则有：

xy的相对误差约为0.002.

改变下列表达式，使计算结果更准确。

(1) (2)

(3) (4)

解：(1)

(2)

(3)

(4)

7、计算的近似值，取。利用以下四种计算格式，试问哪一种算法误差最小。

(1) (2)

(3) (4)

解：计算各项的条件数

由计算可知，最小，即第三种算法误差最小。

解：在计算机上计算该级数的是一个收敛的级数。因为随着的增大，会出现大数吃小数的现象。

9、 通过分析浮点数集合F=(10，3，-2，2)在数轴上的分布讨论一般浮点数集的分布情况。

解：浮点数集合F=(10，3，-2，2)在数轴上离原点越近，分布越稠密；离原点越远，分布越稀疏。一般浮点数集的分布也符合此规律。

10、试导出计算积分的递推计算公式，并分析此递推公式的数值稳定性。

解：

分析计算的误差，设初值的误差为，递推过程的舍入误差不计，并记

，则有 。显然，随着计算的递推，误差越来越小，因此递推公式是稳定的。

11、为减少乘除法运算次数，应将下面算式怎样改写？

解：令，则

试推导求函数值的绝对误差和相对误差。

解：

试推导求函数值的条件数。

解：

序列满足递推公式

求计算到的的误差，并讨论计算过程的稳定性。

解：

…

误差逐渐增大，计算不稳定

写出下面Matlab程序所描述的数学表达式。

for j=1: n

for i=1: m

y ( i )= A (i, j)*x(j)+y(i)

end

end

for j=1:n

y=x(j)*A(:,j)+y

end

解：(1)

(2) A为n列的矩阵，x为n行1列的列向量，y为与A有相同行数的列向量

写出下面Matlab程序所描述的数学表达式。

for i=1: m

for j=1: n

A(i,j)= A (i, j)+x(i)*y(j)

end

end

(2) for j=1:m

A(:,j)= A (:, j)+ y(j)* x(:)

end

解：(1)

?(2) ??A为m列的矩阵,y为与A有相同行数的列向量

第二章习题解答

试证明中的子集“上三角阵”对矩阵乘法是封闭的。

证明：设为上三角阵，则

C=AB，则

，即上三角阵对矩阵乘法封闭。

已知矩阵，求A的行空间及零空间N(A)的基。

解：对进行行变换，

，的基为

由Ax=0可得

N(A)的基为

3.已知矩阵，试计算A的谱半径。

解：

4.试证明其中

。

试证明若是内积空间H中不含零向量的正交向量组，则必线性无关。

证明：假设存在使

由于是正交向量组，则

等式两边与作内积得

则必线性无关

计算下列向量的。

(1)，

(2)，

(3)，k为正整数。

解：(1)，

同理(1)，得，

。

在C[a,b]中，试证明：

(1)为f(x)的范数，(2)为f(x)的范数。

证明：(1)

为f(x)的范数。

(2)

是C[a,b]上的一个范数。

在内积空间中给出Cauchy-Schawz不等式，其中内积，A为对称正定矩阵。

解：。

已知向量，求x，y之间的距离。

解：

试计算，其中m，n是正整数，

解：

11.已知，试计算，，，。

12.在上，由构造带权的首1正交多项式，和。

解：

13.给出点集及权，试构造正交函数组，和。

14.已知向量Gauss变换阵将向量x变为。

解：。

15.已知向量x=(1，2，2)T ，y =(0，3，4)T 。试构造Householder阵H使H x为y的倍数，即H x=ky。给出变换阵H和系数k。