资源描述：
《插值与多项式逼近的数组计算方法实验.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、插值与多项式逼近的数组计算方法实验郑发进22【摘要】计算机软件中经常要用到库函数，如，，，它们是用多项式逼近来计算的。虽然目前最先进的逼近方法是有理函数（即多项式的商），但多项式逼近理论更适于作为数值分析的入门课程。在已知数据具有高精度的情况下，通常用组合多项式来构造过给定数据点的多项式。构造组合多项式的方法有许多种，如线性方程求解、拉格朗日系数多项式以及构造牛顿多项式的方分和系数表。关键字泰勒级数、拉格朗日插值法、牛顿插值法、帕德逼近一、实验目的1.通过具体实验，掌握泰勒级数、拉格朗日插值法、牛顿插值法、帕
2、德逼近的编程技巧。2.比较各插值方法的优劣并掌握。二、实验原理1.泰勒级数在数学中，泰勒级数（英语：Taylorseries）用无限项连加式——级数来表示一个函数，这些相加的项由函数在某一点的导数求得。如果 在点x=x0具有任意阶导数，则幂级数称为 在点x0处的泰勒级数。在泰勒公式中，取x0=0，得到的级数称为麦克劳林级数。函数的麦克劳林级数是x的幂级数，那么这种展开是唯一的，且必然与的麦克劳林级数一致。2.拉格朗日插值法如对实践中的某个物理量进行观测，在若干个不同的地方得到相应的观测值，拉格朗日插值法可
3、以找到一个多项式，其恰好在各个观测的点取到观测到的值。这样的多项式称为拉格朗日（插值）多项式。数学上来说，拉格朗日插值法可以给出一个恰好穿过二维平面上若干个已知点的多项式函数。在平面上有(x1,y1)(x2,y2)...(xn,yn)共n个点，现作一条函数f（x）使其图像经过这n个点。作n个多项式pi(x),i=1,2,3...,n,使得最后可得 3.牛顿插值法插值法利用函数f(x)在某区间中若干点的函数值，作出适当的特定函数，在这些点上取已知值，在区间的其他点上用这特定函数的值作为函数f(x)的近似值。如果
4、这特定函数是多项式，就称它为插值多项式。利用插值基函数很容易得到拉格朗日插值多项式，公式结构紧凑，在理论分析中甚为方便，但当插值节点增减时全部插值基函数均要随之变化，整个公式也将发生变化，这在实际计算中是很不方便的，为了克服这一缺点，提出了牛顿插值。牛顿插值通过求各阶差商，递推得到的一个公式：牛顿插值与拉格朗日插值具有唯一性。4.帕德逼近它不仅与逼近论中其他许多方法有着密切的关系，而且在实际问题特别是许多物理问题中有着广泛的应用。设是在原点某邻域内收敛的、具有复系数的麦克劳林级数。欲确定一个有理函数，式中,使
5、得前次方的系数为0，即使得此处约定qk＝0（k>n）。虽然所求得的Pm(z)和Qn(z)不惟一，但是比式却总是惟一的。有理函数称为F(z)的(m,n)级帕德逼近,记为(m/n)。由(m/n)所形成的阵列称为帕德表。三、实验内容1.P154.1:用plot命令，在同一幅图中绘制区间-1≤x≤1上的sin(x)，以及P5(x),P7(x)和P9(x)。其中：·········································式（1）··································
6、··式（2）····························式（3）2.P171.2:下表（表一）给出了11月8号美国洛杉矶的一个郊区在5小时内的测量温度。（a）对表一中的数据构造一个拉格朗日插值多项式。（b）估计这5小时内的平均温度。（c）在同一坐标系中画出表中的数据和由（a）得到的多项式。讨论用（a）中的多项式计算平均温度可能产生的误差。时间（下午）华氏度166266365464563663表13.P178.1:用牛顿插值多项式计算实验P171.2的内容。4.P194.1:比较对于函数的逼近：泰勒
7、多项式逼近：····························式（4）帕德逼近：································式（5）（a）在同一坐标系中画出f(x)，T4(x)，R2,2(x)的曲线。（b）分别求出在区间[-1，1]上用T4(x)和R2,2(x)逼近f(x)的最大误差。1.P194.3:比较对于函数f(x)=tan(x)的逼近：泰勒多项式逼近：·······················式（6）帕德逼近：·····························式（
8、7）（a）在同一坐标系中画出f(x)，T9(x)，R5,4(x)的曲线。（b）分别求出在区间[-1，1]上用T9(x)和R5,4(x)逼近f(x)的最大误差。四、实验结果及分析1.P154.1：实验描述：（1）plot绘图的原理为连续点绘图，只需输入一组等间距的坐标点即可完成；(2)坐标点的计算使用C++完成，计算完成后输入文件中；(3)绘图使用matlab的plot函数完成，具体方法为从文件中读取
