西安交通大学2015年计算方法上机作业.doc
《西安交通大学2015年计算方法上机作业.doc》由会员分享,可在线阅读,更多相关《西安交通大学2015年计算方法上机作业.doc(30页珍藏版)》请在咨信网上搜索。
1、计算方法上机作业学院:能动学院班级:硕4011班P81页第四题求20阶三对角方程组的解,其中:,算法分析:追赶法的算法组织如下:(1) 输入三对角矩阵和右端向量;(2) 将压缩为四个一维数组,将分解矩阵压缩为三个一维数组(3) 对做分解(也可以用分解)导出追赶法的计算步骤如下:(4) 回代求解 (5) 停止,输出结果 结果:将结果带入方程可以验证结果完全正确。第五题算法分析:将原始线性代数方程组改写为的形式,其中为的矩阵函数。于是可以得到迭代格式:,此即为迭代法的迭代格式。如果在计算时,将已经算出的分量立即代换对应分量,则得到迭代法的迭代格式。1. 迭代法的算法组织如下:(1) 给出迭代格式(
2、2) 给出迭代初始向量、允许误差和最大迭代次数(3) 按照迭代格式进行迭代,直至达满足迭代停止条件(4) 停止,输出结果2. 迭代法的算法组织如下:(1) 给出迭代格式(2) 给出迭代初始向量、允许误差和最大迭代次数(3) 按照迭代格式,并且将已经算出的分量立即代换对应分量进行迭代,直至达满足迭代停止条件(4) 停止,输出结果结果:将结果带入方程验证的确满足误差1e-3P131页第一题已知,对a 计算函数在点处的值b 求插值数据点的插值多项式和三次样条插值多项式c 对,计算和相应的函数值d 计算,解释所得到结果。算法分析1. 插值插值多项式是实质是插值多项式的零次式和一次式的推广,一般形式为:
3、2. 三次样条插值在节点分成的每个小区间上利用其节点处的二次导数值(即弯矩值)进行线性插值,再在此区间上积分两次,利用节点处的函数值求得两个积分常数,从而得到三次样条函数的表达式。要确定三次样条插值函数,用到的光滑性,利用三次样条函数在节点处的一阶导数连续得到满足个方程:其中三弯矩方程组只有个方程,不能确定未知量。为此,可以给出区间的端点处的一阶导数值。由和可以导出其中联合先前得到的三弯矩方程组和两个补充方程得到恰定的三弯矩方程组。它是严格对角占优的三对角方程组,可以用追赶法求解。结果:a 直接计算即可得出对于n=5有对于n=10有对于n=20有b对于n=5有牛顿插值多项式样条插值函数(分别对
4、于各个小区间)为:对于n=10有牛顿插值函数p = (9094987740384*X9 + 5272349642029869237892763424*X8 - 16938813146112*X7 - 10950865478705002051580730624*X6 + 10150693091136*X5 + 7491941821973715378406714008*X4 + 1915628554944*X3 - 2014100801013926045821422321*X2 + 322192441744*X + 210052147079480949741593257)/243810615467
5、456022706126848 pp = 3.73*10(-14)*X9 + 21.6*X8 - 6.95*10(-14)*X7 - 44.9*X6 + 4.16*10(-14)*X5 + 30.7*X4 + 7.86*10(-15)*X3 - 8.26*X2 + 1.32*10(-15)*X + 0.862对于n=10有样条插值函数:对于n=20有牛顿插值函数- 6.82*10(-10)*X19 - 25770.0*X18 + 2.52*10(-9)*X17 + 95660.0*X16 - 3.79*10(-9)*X15 - 1.47*105*X14 + 3.02*10(-9)*X13 +
6、1.21*105*X12 - 1.46*10(-9)*X11 - 58660.0*X10 + 4.3*10(-10)*X9 + 17220.0*X8 - 6.6*10(-11)*X7 - 3066.0*X6 + 4.06*10(-12)*X5 + 328.0*X4 - 3.19*10(-13)*X3 - 21.6*X2 + 3.08*10(-15)*X + 0.993对于n=20有样条插值函数: - 5.25*X3 - 14.37*X2 - 12.98*X - 3.827 1.723*X3 + 4.461*X2 + 3.963*X + 1.256 0.07504*X3 + 0.5066*X2
7、+ 0.7996*X + 0.4127 0.9738*X3 + 2.394*X2 + 2.121*X + 0.721 1.552*X3 + 3.436*X2 + 2.746*X + 0.846 3.552*X3 + 6.434*X2 + 4.245*X + 1.096 5.726*X3 + 9.044*X2 + 5.289*X + 1.235 12.54*X3 + 15.17*X2 + 7.127*X + 1.419 - 32.79*X3 - 12.02*X2 + 1.688*X + 1.056 - 89.07*X3 - 28.91*X2 - 2.345*10(-15)*X + 1.0 89.
8、07*X3 - 28.91*X2 - 6.395*10(-15)*X + 1.0 32.79*X3 - 12.02*X2 - 1.688*X + 1.056 - 12.54*X3 + 15.17*X2 - 7.127*X + 1.419 - 5.726*X3 + 9.044*X2 - 5.289*X + 1.235 - 3.552*X3 + 6.434*X2 - 4.245*X + 1.096 - 1.552*X3 + 3.436*X2 - 2.746*X + 0.846 - 0.9738*X3 + 2.394*X2 - 2.121*X + 0.721 - 0.07504*X3 + 0.506
9、6*X2 - 0.7996*X + 0.4127 - 1.723*X3 + 4.461*X2 - 3.963*X + 1.256 5.25*X3 - 14.37*X2 + 12.98*X - 3.827c 容易计算各点函数值 YK = 0.0400 0.0416 0.0433 0.0451 0.0471 0.0491 0.0513 0.0536 0.0561 0.0588 0.0617 0.0648 0.0681 0.0716 0.0755 0.0796 0.0841 0.0890 0.0943 0.1000 0.1063 0.1131 0.1206 0.1289 0.1379 0.1479
10、0.1590 0.1712 0.1848 0.2000 0.2169 0.2358 0.2571 0.2809 0.3077 0.3378 0.3717 0.4098 0.4525 0.5000 0.5525 0.6098 0.6711 0.7353 0.8000 0.8621 0.9174 0.9615 0.9901 1.0000 0.9901 0.9615 0.9174 0.8621 0.8000 0.7353 0.6711 0.6098 0.5525 0.5000 0.4525 0.4098 0.3717 0.3378 0.3077 0.2809 0.2571 0.2358 0.2169
11、 0.2000 0.1848 0.1712 0.1590 0.1479 0.1379 0.1289 0.1206 0.1131 0.1063 0.1000 0.0943 0.0890 0.0841 0.0796 0.0755 0.0716 0.0681 0.0648 0.0617 0.0588 0.0561 0.0536 0.0513 0.0491 0.0471 0.0451 0.0433 0.0416 0.0400 对于n=5 可以计算相应插值函数的值Nk= -0.0496 -0.1257 -0.1906 -0.2449 -0.2891 -0.3239 -0.3497 -0.3672 -0.
12、3769 -0.3793 -0.3750 -0.3643 -0.3479 -0.3263 -0.2997 -0.2688 -0.2340 -0.1957 -0.1542 -0.1101 -0.0636 -0.0152 0.0347 0.0859 0.1379 0.1905 0.2434 0.2962 0.3487 0.4005 0.4515 0.5014 0.5499 0.5968 0.6419 0.6850 0.7260 0.7646 0.8008 0.8342 0.8649 0.8927 0.9174 0.9391 0.9576 0.9728 0.9846 0.9932 0.9983 1.
13、0000 0.9983 0.9932 0.9846 0.9728 0.9576 0.9391 0.9174 0.8927 0.8649 0.8342 0.8008 0.7646 0.7260 0.6850 0.6419 0.5968 0.5499 0.5014 0.4515 0.4005 0.3487 0.2962 0.2434 0.1905 0.1379 0.0859 0.0347 -0.0152 -0.0636 -0.1101 -0.1542 -0.1957 -0.2340 -0.2688 -0.2997 -0.3263 -0.3479 -0.3643 -0.3750 -0.3793 -0
14、.3769 -0.3672 -0.3497 -0.3239 -0.2891 -0.2449 -0.1906 -0.1257 -0.0496 Sk=0.0383 0.0378 0.0371 0.0363 0.0355 0.0348 0.0344 0.0343 0.0347 0.0357 0.0373 0.0397 0.0430 0.0472 0.0526 0.0592 0.0672 0.0765 0.0874 0.1000 0.1143 0.1302 0.1475 0.1660 0.1857 0.2063 0.2277 0.2497 0.2721 0.2948 0.3175 0.3403 0.3
15、628 0.3849 0.4064 0.4273 0.4472 0.4661 0.4837 0.5000 0.5147 0.5279 0.5396 0.5497 0.5582 0.5652 0.5706 0.5745 0.5768 0.5776 0.5768 0.5745 0.5706 0.5652 0.5582 0.5497 0.5396 0.5279 0.5147 0.5000 0.4837 0.4661 0.4472 0.4273 0.4064 0.3849 0.3628 0.3403 0.3175 0.2948 0.2721 0.2497 0.2277 0.2063 0.1857 0.
16、1660 0.1475 0.1302 0.1143 0.1000 0.0874 0.0765 0.0672 0.0592 0.0526 0.0472 0.0430 0.0397 0.0373 0.0357 0.0347 0.0343 0.0344 0.0348 0.0355 0.0363 0.0371 0.0378 0.0383 对于n=20 可以计算相应插值函数的值Nk= 8.2065 7.7749 4.8826 2.1401 0.3516 -0.4925 -0.6818 -0.5282 -0.2658 -0.0308 0.1204 0.1843 0.1839 0.1492 0.1065 0
17、.0729 0.0560 0.0558 0.0680 0.0865 0.1057 0.1218 0.1334 0.1407 0.1454 0.1497 0.1557 0.1648 0.1779 0.1948 0.2149 0.2373 0.2611 0.2861 0.3123 0.3403 0.3713 0.4066 0.4476 0.4950 0.5492 0.6094 0.6741 0.7408 0.8063 0.8669 0.9189 0.9589 0.9841 0.9927 0.9841 0.9589 0.9189 0.8669 0.8063 0.7408 0.6741 0.6094
18、0.5492 0.4950 0.4476 0.4066 0.3713 0.3403 0.3123 0.2861 0.2611 0.2373 0.2149 0.1948 0.1779 0.1648 0.1557 0.1497 0.1454 0.1407 0.1334 0.1218 0.1057 0.0865 0.0680 0.0558 0.0560 0.0729 0.1065 0.1492 0.1839 0.1843 0.1204 -0.0308 -0.2658 -0.5282 -0.6818 -0.4925 0.3516 2.1401 4.8826 7.7749 8.2065 Sk=0.039
19、0 0.0403 0.0423 0.0446 0.0471 0.0494 0.0516 0.0539 0.0563 0.0588 0.0616 0.0647 0.0680 0.0716 0.0755 0.0796 0.0841 0.0890 0.0943 0.1000 0.1063 0.1131 0.1206 0.1289 0.1379 0.1479 0.1589 0.1712 0.1848 0.2000 0.2170 0.2359 0.2572 0.2810 0.3077 0.3375 0.3711 0.4089 0.4517 0.5000 0.5541 0.6127 0.6742 0.73
20、72 0.8000 0.8606 0.9152 0.9594 0.9891 1.0000 0.9891 0.9594 0.9152 0.8606 0.8000 0.7372 0.6742 0.6127 0.5541 0.5000 0.4517 0.4089 0.3711 0.3375 0.3077 0.2810 0.2572 0.2359 0.2170 0.2000 0.1848 0.1712 0.1589 0.1479 0.1379 0.1289 0.1206 0.1131 0.1063 0.1000 0.0943 0.0890 0.0841 0.0796 0.0755 0.0716 0.0
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 西安交通大学 2015 计算方法 上机 作业
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【天****】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【天****】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。