数值分析课程设计-Read.docVIP

数值分析课程设计 理学院2004级信科专业 2004376330张卓鹏 课程设计的目的和意义： 《课程设计》是数值分析的同步课程，是《数值分析》的上机实习课。 《数值分析》课程中构造了各种有效的算法和有效公式，同学们通过上机作课程设计，学习揣摩这些算法的思想和构造，评价算法的优劣性。 通过上机，可以提高我们运用数学软件编程解决问题的能力，为今后从事科学计算和软件开发打下良好的基础。 课程设计的题目： 多项式插值的振荡现象 设计目的： 通过对多项式插值现象的观察,了解多项式的次数与逼近效 的关系，提高同学们分析实验结果的能力。 问题提出： 考虑在一个固定区间上用插值逼近一个函数。显然，Lagrange插值中使用的节点越多，插值多项式的次数就越高。我们自然关心插值多项式增加时，Ln(x)是否也更加靠近被逼近的函数。龙格(Runge)给出的一个例子是极著名并富有启发性的。设区间[-1,1]上的函数 区间[-考虑设计1,1]的一个等距划分，节点为 则拉个朗日插值多项式为 其中的li(x),i=0,1,2,…,n是n次Lagrange插值基函数。 设计要求： 1．选择不断增大的分点数n=2,3,… *画出原函数f(x)及插值多项式函数Ln(x)在[-1,1]上的图像； *给出每一次逼近的最大误差； *比较并分析实验结果。 2．选择其它函数，例如定义在区间[-5,5]上的函数。 重复上述实验看其结果如何。 区间[a,b]上切比雪夫(Chebychev)点的定义为 以x1,x2,…,xn+1为插值节点构造上述各函数的Lagrange插值多项式，比较其结果。 设计过程： 已知函数f(x)在n+1个点x0,x1,…,xn处的函数值为 y0,y1,…,yn 。求一n次多项式函数Pn(x)，使其满足： Pn(xi)=yi,i=0,1,…,n. 解决此问题的拉格朗日插值多项式公式如下 其中Li(x) 为n次多项式： （1）.在MATLAB6.5中输入函数 当取不同的分点数n时，所得图象与原函数图象对比如下： n=2时；最大误差为： Max[L(x)－f(x)]=0.6462 n=3时，max[L(X)-f(X)]=0.7070； n=6时，max[L(X)-f(X)]=0.6169； n=8时，max[L(X)-f(X)]=1.0452； n=10时；最大误差为：max[L(X)-f(X)]=1.9156； n=20时，最大误差为：max[L(X)-f(X)]=58.5855。 图象分析： 从图中可以看出当插值节点很少时，插值的误差很大，插值图象与原图象没有很好重叠在一起，而当随着插值的节点增加，中间能很好的重叠，但是两边出现很大误差，随着n值的增多，总体上分散的越厉害，最大误差也逐渐增加，在n＝3时，最大误差为0.6462，但到了n＝20时，已经变成了58.5855，这种随着节点数增多依然不能很好的接近被插值函数的现象称为“龙格现象”，亦称为多项式插值的振荡现象。因此通过增加节点数从而提高插值多项式的次数来逼近被插函数是不可取的。 （2）.对于定义在定义在区间[-5,5]上的函数h(x)=x/(1+x^4); 输入函数后取节点数n=4，8，10；得到图象如下 n=4时，max[L(X)-h(X)]=0.4020； n=8时，max[L(X)-h(X)]=0.1708； n=10时，max[L(X)-h(X)]=0.1092。 图象分析： 从图象可以看出随着插值节点数的增加出现异常的摆动，中间能较好的接近原函数，但两边却出现很大的误差。 (3).对定义在（-5，5）上的函数 g(x)=arctanx 在matlab里输入相应的lagrange插值公式，当n取不同的值时，得到相应图象如下： n=4时，max[L(X)-g(X)]=0.4458； n=8时，max[L(X)-g(X)]=0.3240； N=10时，max[L(X)-g(X)]=0.8066。 图象分析： 从图象可以看出随着插值节点数的增加出现异常的摆动，中间能较好的接近原函数，随着插值节点增加有较好的插值重叠，但两边却出现很大的误差。 （3）. 以上实验是从对区间的等距划分而作出的拉格朗日多项式，也因此产生了龙格振荡现象，现在我们通过用切比雪夫点来对上述各函数的区间进行划分，以此来建立拉格朗日多项式，看这样建立的拉格朗日多项式是否还会出现龙格振荡现象。 还是考虑函数 ，为了与等距节点进行比较，我们仍然选取