[理学]数值分析习题解答
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