全球旧事资料 分类
16897x3086696x266375x118699(2)根据原散点构建拉格朗日插值多项式,对给出的xi的数值进行插值得到yi,在下图可以看到利用拉格朗日插值多项式插值的图像与曲线拟合度很好。yiColum
s1through134238404149473507422436011127921768318162617911811021020333198565403584Colum
s14through176107947957099377881023713
122977
f3第二题
31题目用雅格比法与高斯-赛德尔迭代法解下列方程组Axb,研究其收敛性,上机验证理论分析是否正确,比较它们的收敛速度,观察右端项对迭代收敛有无影响。1A行分别为A1621A2142A3314;b1324Tb2100200345T2A行分别为A110808A208108A308081;b1321Tb25010T3A行分别为A113A271;b46T32结果分析在以下的程序中,默认初始值为0,最大迭代次数N为50,精度为000001(1)根据题意可得方程:
对此方程分别采用两种方法,通过运行jacobim与gausssidelm程序进行收敛性分析:Jacobi:满足误差限的迭代近似解为:07273需要迭代18次;Gausssidel:满足误差限的迭代近似解为:07273需要迭代11次。08081025250808102525
Jacobi:满足误差限的迭代近似解为:363636需要迭代24次;Gausssidel:满足误差限的迭代近似解为:需要迭代16次。
207071140404207071140404
363636
通过以上计算分析,此方程的解收敛。根据迭代次数的比较,可以得到Gausssidel迭代法要比Jacobi收敛性更快,更好。同时,上述两个方程系数阵都相同,但右端不同,用同样方法不同右端比较,得到结论,右端对迭代收敛有影响。
f(2)
用jacobi到最大迭代次数50次时仍不能收敛。用Gausssidel满足误差限的迭代近似解为:57692迭代16次。0769242308需要
用jacobi到最大迭代次数50次时仍不能收敛。用Gausssidel满足误差限的迭代近似解为:326923要迭代43次。通过上述计算分析,用jacobi迭代时,两个方程都不能收敛,用Gaussssidel迭代时,可以得到收敛的近似解,既,Gausssidel方法具有良好的收敛性。(3)76923423077需
对此方程组,无论jacobi还是Gausssidel迭代法都不能得到收敛的近似解。可以推断是给出的行列式无法满足这两种方法的收敛条件。系数阵既不是行对角占优矩阵,也不是列占优阵,因此用这两种方法,得不到收敛的近似解。
f4第三题
41题目用Ru
geKutta4阶算法对初值问题y20y,y01按不同步长求解,用于观察稳定区间的作用,推荐两种步长h0102。注:此方程的精确解为:ye20x42计算结果分析取不同步长得到求解的微分方程的r
好听全球资料 返回顶部