《西安交通大学计算方法上机作业》由会员分享，可在线阅读，更多相关《西安交通大学计算方法上机作业（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.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
