西南交通大学研究数值分析上机实习报告.docVIP

数值分析上机实习报告要求 1．应提交一份完整的实习报告。具体要求如下： （1）要有封面，封面上要标明姓名、学号、专业和联系电话； （2）要有序言，说明所用语言及简要优、特点，说明选用的考量； （3）要有目录，指明题目、程序、计算结果，图标和分析等内容所在位置，作到信息简明而完全； （4）要有总结，全方位总结机编程计算的心得体会； （5）尽量使报告清晰明了，一般可将计算结果、图表及对比分析放在前面，程序清单作为附录放在后面，程序中关键部分要有中文说明或标注，指明该部分的功能和作用。 2．程序需完好保存到期末考试后的一个星期，以便老师索取用于验证、询问或质疑部分内容。 3．认真完成实验内容，可以达到既学习计算方法又提高计算能力的目的，还可以切身体会书本内容之精妙所在，期间可以得到很多乐趣。 4．拷贝或抄袭他人结果是不良行为，将视为不合格。 5．报告打印后按要求的时间提交给任课老师。 上机实习必须在规定的时间内完成，可要求在考前或考后一个星期内提交。不合格者和不交者不通过 1. 已知：a=-5,b=5, 以下是某函数 f(x)的一些点（xk,yk）, 其中xk=a+0.1(k-1) ,k=1,..,101 xk=a+0.1k, 请用插值类方法给出函数f(x)的一个解决方案和具体结果。并通过实验考虑下列问题 实验前分析： 所给的节点一共有101个，用Lagrange插值发最高可以做次数100的插值多项式做低次插值是需要进行分组，间插值区域分成互不重叠的区间。分组时须按区间位置高低的顺序依次划分，有可能最后余下的区间节点数不足够做n次的Lagrange插值多项式，需要特殊处理。分完组后在分得的区间套用插值公式得到：,那么最后的插值多项式是分段的多项式： 实验结果： 下面列出其中某些次的插值多项式 的系数 1 0.0000 0.0000 0.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 1.0000 0.0000 0.0000 0.5000 8.2083 53.866 177.62 289.37 189.48 -0.2500 -3.7500 -22.254 -64.297 -94.672 -54.204 6.6667 76.792 351.44 798.62 896.04 397.66 -37.833 -336.08 -1177.9 -2027.1 -1707.0 -558.15 125.50 792.17 1935.5 2261.1 1239.2 248.96 -284.00 -1071.1 -1458.3 -829.80 -166.75 -6.1346 427.33 536.08 25.083 -131.88 0.0694 10.000 -427.33 536.08 -25.083 -131.88 -0.0694 10.000 284.00 -1071.1 1458.3 -829.80 166.74 -6.1346 -125.50 792.17 -1935.5 2261.1 -1239.2 248.96 37.833 -336.08 1177.9 -2027.1 1707.0 -558.15 -6.6667 76.792 -351.44 798.62 -896.04 397.66 0.2500 -3.7500 22.254 -64.298 94.672 -54.204 -0.5000 8.2083 -53.867 177.63 -289.37 189.48 0.0000 0.0000 0.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 1.0000 0.0000 0.0000 Ln(x)的次数n越高，逼近f(x)的程度越好？ 答：不是，比如 当时只在内收敛，在这区间外是发散的。 高次插值收敛性如何？ 答：高阶倒数存在时，高次插值是收敛的，但不一定 如何选择等距插值多项式次数 ？ 答：等距插值采用牛顿插值公式，可以选用牛顿前插公式或牛顿后插公