matlab实现盖尔圆,[理学]数值分析习题解答.doc

[理学]数值分析习题解答

6.(1)设，试求

(2)设，试求

解

；

4．设，

(1)试对进行PLU分解：；

(2)根据PLU分解求解。

解 (1)

(2)

8．分别用Householder变换法和MGS法对进行QR分解

解　(1) Householder法对A进行QR分解

令，调用算法2.1有，所以

,

故

再令,调用算法2.1得,则

,

故

.

10．

求的最小二乘问题的全部解。

解

因为

所以，齐次方程组的基础解系为

的特解为

所以的LS问题的全部解为

第三章习题解

第一题:

因为有三对数据故可以求次数不超过2的插值多项式。

若用Lagrange插值法：

若用Newton插值法：

xkf(xk)f[x0,x1]f(x0,x1,x2)01121243/21/2

若用分段线性插值法，则插值函数有两个分枝。

第二题：

由插值余项可得：

第三题：

所用的插值数据如下表：

x0.40000.50000.60000.70000.8000y0.38940.47940.56460.64420.7174dy0.92110.87760.82530.76480.6967Newton插值法可得差商表：

xkf(xl)一阶差商二阶差商三阶差商四阶差商0.40000.3894NaNNaNNaNNaN0.50000.47940.9000NaNNaNNaN0.60000.56460.8760-0.2400NaNNaN0.70000.64420.8493-0.2533-0.1333NaN0.80000.71740.8200-0.2667-0.1333-0.0000由此可以写出Newton插值多项式N4(x)=?，并且可以求出

要注意这里求出的应当是次数不超过四次的newton插值多项式。

(2)如果用分段线性插值，则只需要求出0.5789所在的小区间上的线性插值多项式的表达式，再把0.5789代入即求sin0.5789的近似值：

(3)如果用分段三次Hermite插值，则只需要求出0.5789所在的小区间上的三次Hermite插值多项式的表达式，再把0.5789代入即求sin0.5789的近似值：

第四题：

对f(x)求导后易判断f(x)存在反函数，因此就是所要求的根。现在对进行插值，使用如下数据表：

f(x)0.70010.40160.1081-0.1744-0.4375x0.10000.20000.30000.40000.5000差商表为：

ykf-1(yk)一阶差商二阶差商三阶差商四阶差商0.70010.1000NaNNaNNaNNaN0.40160.2000-0.3350NaNNaNNaN0.10810.3000-0.33780.0096NaNNaN-0.17440.4000-0.34310.0140-0.0153NaN-0.43750.5000-0.35160.0198-0.01860.0125利用Newton插值法，可得近似根：

真正的零点为。这种求根方法称为反插值法。

第六题：

解：由于是等距节点，故。再由于函数值都为时尚早，故各阶差商都为0，从而。由(2。16)得三次样条的基本方程组为

(1)由第一边界条件可得另外两个方程

从而由(2。20)建立第一边界的三次样条方程组

求解此方程组得解，再由公式(2。15)就可以建立三次样条的系数矩阵：

因此可以写出第一边界条件下的三次样条函数：

(2)由第二边界条件由第二边界条件及三次样条的基本方程可得

由此可得解，再由公式(2。15)就可以建立三次样条的系数矩阵：

因此可以写出第二边界条件下的三次样条函数：

利用MATLAB的函数csape求出的结果与此相同。

第七题：

(1)第一数据表的插值多项式是二点三次Hermite插值多项式，利用公式(1。23)，插值基函数为：

因此所求插值多项式为：

(2)共有五对插值数据，可以唯确定一个至多四次的插值多项式，由于数据的特殊性采用待定系数法：

由于在x=0处的函数值及导数值均为0，故得从而只有三个系数需要确定：

再由其余三对数据可得如下方程组：

第八题：

插值基函数见课本P60公式(1。7)，要注意插值节点节点是互不相同的。

(1)证：常值函数在节点处的Lagrange插值多项式恰为

因为f(x