西安交通大学计算方法上机作业Word格式.docx
- 文档编号:13100297
- 上传时间:2022-10-05
- 格式:DOCX
- 页数:30
- 大小:413.40KB
西安交通大学计算方法上机作业Word格式.docx
《西安交通大学计算方法上机作业Word格式.docx》由会员分享,可在线阅读,更多相关《西安交通大学计算方法上机作业Word格式.docx(30页珍藏版)》请在冰豆网上搜索。
追赶法的算法组织如下:
(1)输入三对角矩阵和右端向量;
(2)将压缩为四个一维数组,将分解矩阵压缩为三个一维数组
(3)对做分解(也可以用分解)导出追赶法的计算步骤如下:
(4)回代求解
(5)停止,输出结果
结果:
将结果带入方程可以验证结果完全正确。
第五题
将原始线性代数方程组改写为的形式,其中为的矩阵函数。
于是可以得到迭代格式:
,此即为迭代法的迭代格式。
如果在计算时,将已经算出的分量立即代换对应分量,则得到迭代法的迭代格式。
1.迭代法的算法组织如下:
(1)给出迭代格式
(2)给出迭代初始向量、允许误差和最大迭代次数
(3)按照迭代格式进行迭代,直至达满足迭代停止条件
(4)停止,输出结果
2.迭代法的算法组织如下:
(3)按照迭代格式,并且将已经算出的分量立即代换对应分量进行迭代,直至达满足迭代停止条件
结果:
将结果带入方程验证的确满足误差<
1e-3
P131页第一题
已知,对
a计算函数在点处的值
b求插值数据点的插值多项式和三次样条插值多项式
c对,计算和相应的函数值
d计算,,解释所得到结果。
算法分析
1.插值
插值多项式是实质是插值多项式的零次式和一次式的推广,一般形式为:
2.三次样条插值
在节点分成的每个小区间上利用其节点处的二次导数值(即弯矩值)进行线性插值,再在此区间上积分两次,利用节点处的函数值求得两个积分常数,从而得到三次样条函数的表达式。
要确定三次样条插值函数,用到的光滑性,利用三次样条函数在节点处的一阶导数连续得到满足个方程:
其中
三弯矩方程组只有个方程,不能确定未知量。
为此,可
以给出区间的端点处的一阶导数值。
由和可以导出
联合先前得到的三弯矩方程组和两个补充方程得到恰定的三弯矩方程组。
它是严格对角占优的三对角方程组,可以用追赶法求解。
a直接计算即可得出
对于n=5有
对于n=10有
对于n=20有
b对于n=5有牛顿插值多项式
样条插值函数(分别对于各个小区间)为:
对于n=10有牛顿插值函数
p=
(9094987740384*X^9+5272349642029869237892763424*X^8-16938813146112*X^7-10950865478705002051580730624*X^6+10150693091136*X^5+7491941821973715378406714008*X^4+1915628554944*X^3-2014100801013926045821422321*X^2+322192441744*X+210052147079480949741593257)/243810615467456022706126848
pp=
3.73*10^(-14)*X^9+21.6*X^8-6.95*10^(-14)*X^7-44.9*X^6+4.16*10^(-14)*X^5+30.7*X^4+7.86*10^(-15)*X^3-8.26*X^2+1.32*10^(-15)*X+0.862
对于n=10有样条插值函数:
对于n=20有牛顿插值函数
-6.82*10^(-10)*X^19-25770.0*X^18+2.52*10^(-9)*X^17+95660.0*X^16-3.79*10^(-9)*X^15-1.47*10^5*X^14+3.02*10^(-9)*X^13+1.21*10^5*X^12-1.46*10^(-9)*X^11-58660.0*X^10+4.3*10^(-10)*X^9+17220.0*X^8-6.6*10^(-11)*X^7-3066.0*X^6+4.06*10^(-12)*X^5+328.0*X^4-3.19*10^(-13)*X^3-21.6*X^2+3.08*10^(-15)*X+0.993
对于n=20有样条插值函数:
-5.25*X^3-14.37*X^2-12.98*X-3.827
1.723*X^3+4.461*X^2+3.963*X+1.256
0.07504*X^3+0.5066*X^2+0.7996*X+0.4127
0.9738*X^3+2.394*X^2+2.121*X+0.721
1.552*X^3+3.436*X^2+2.746*X+0.846
3.552*X^3+6.434*X^2+4.245*X+1.096
5.726*X^3+9.044*X^2+5.289*X+1.235
12.54*X^3+15.17*X^2+7.127*X+1.419
-32.79*X^3-12.02*X^2+1.688*X+1.056
-89.07*X^3-28.91*X^2-2.345*10^(-15)*X+1.0
89.07*X^3-28.91*X^2-6.395*10^(-15)*X+1.0
32.79*X^3-12.02*X^2-1.688*X+1.056
-12.54*X^3+15.17*X^2-7.127*X+1.419
-5.726*X^3+9.044*X^2-5.289*X+1.235
-3.552*X^3+6.434*X^2-4.245*X+1.096
-1.552*X^3+3.436*X^2-2.746*X+0.846
-0.9738*X^3+2.394*X^2-2.121*X+0.721
-0.07504*X^3+0.5066*X^2-0.7996*X+0.4127
-1.723*X^3+4.461*X^2-3.963*X+1.256
5.25*X^3-14.37*X^2+12.98*X-3.827
c容易计算各点函数值
YK=
0.04000.04160.04330.04510.04710.04910.05130.05360.05610.05880.06170.06480.06810.07160.07550.07960.08410.08900.09430.10000.10630.11310.12060.12890.13790.14790.15900.17120.18480.20000.21690.23580.25710.28090.30770.33780.37170.40980.45250.50000.55250.60980.67110.73530.80000.86210.91740.96150.99011.00000.99010.96150.91740.86210.80000.73530.67110.60980.55250.50000.45250.40980.37170.33780.30770.28090.25710.23580.21690.20000.18480.17120.15900.14790.13790.12890.12060.11310.10630.10000.09430.08900.08410.07960.07550.07160.06810.06480.06170.05880.05610.05360.05130.04910.04710.04510.04330.04160.0400
对于n=5可以计算相应插值函数的值
Nk=
-0.0496-0.1257-0.1906-0.2449-0.2891-0.3239-0.3497-0.3672-0.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.01520.03470.08590.13790.19050.24340.29620.34870.40050.45150.50140.54990.59680.64190.68500.72600.76460.80080.83420.86490.89270.91740.93910.95760.97280.98460.99320.99831.00000.99830.99320.98460.97280.95760.93910.91740.89270.86490.83420.80080.76460.72600.68500.64190.59680.54990.50140.45150.40050.34870.29620.24340.19050.13790.08590.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
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 西安交通大学 计算方法 上机 作业