- 1、本文档共79页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
[理学]数值分析习题解答
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插值法:
xk f(xk) f[x0,x1] f(x0,x1,x2) 0 1 1 2 1 2 4 3/2 1/2
若用分段线性插值法,则插值函数有两个分枝。
第二题:
由插值余项可得:
第三题:
所用的插值数据如下表:
x 0.4000 0.5000 0.6000 0.7000 0.8000 y 0.3894 0.4794 0.5646 0.6442 0.7174 dy 0.9211 0.8776 0.8253 0.7648 0.6967 Newton插值法可得差商表:
xk f(xl) 一阶差商 二阶差商 三阶差商 四阶差商 0.4000 0.3894 NaN NaN NaN NaN 0.5000 0.4794 0.9000 NaN NaN NaN 0.6000 0.5646 0.8760 -0.2400 NaN NaN 0.7000 0.6442 0.8493 -0.2533 -0.1333 NaN 0.8000 0.7174 0.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.7001 0.4016 0.1081 -0.1744 -0.4375 x 0.1000 0.2000 0.3000 0.4000 0.5000 差商表为:
yk f-1(yk) 一阶差商 二阶差商 三阶差商 四阶差商 0.7001 0.1000 NaN NaN NaN NaN 0.4016 0.2000 -0.3350 NaN NaN NaN 0.1081 0.3000 -0.3378 0.0096 NaN NaN -0.1744 0.4000 -0.3431 0.0140 -0.0153 NaN -0.4375 0.5000 -0.3516 0.0198 -0.0186 0.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
文档评论(0)