内容发布更新时间 : 2024/10/20 3:27:02星期一 下面是文章的全部内容请认真阅读。

《数值方法》实验报告

1

插值与多项式逼近的数组计算方法实验

郑发进 2012042020022

【摘要】计算机软件中经常要用到库函数，如sin（x），cos（x），ex，它们是用多项式逼近来计算的。虽然目前最先进的逼近方法是有理函数（即多项式的商），但多项式逼近理论更适于作为数值分析的入门课程。在已知数据具有高精度的情况下，通常用组合多项式来构造过给定数据点的多项式。构造组合多项式的方法有许多种，如线性方程求解、拉格朗日系数多项式以及构造牛顿多项式的方分和系数表。

关 键 字 泰勒级数、拉格朗日插值法、牛顿插值法、帕德逼近

一、实验目的

1.通过具体实验，掌握泰勒级数、拉格朗日插值法、牛顿插值法、帕德逼近的编 程技巧。

2.比较各插值方法的优劣并掌握。

二、实验原理

1.泰勒级数

在数学中，泰勒级数（英语：Taylor series）用无限项连加式——级数来表示一个函数，这些相加的项由函数在某一点的导数求得。

如果

在点x=x0具有任意阶导数，则幂级数

称为

在点x0处的泰勒级数。

在泰勒公式中，取x0=0，得到的级数

称为麦克劳林级数。函数

的麦克劳林级数是x的幂级数，那么这种展开

《数值方法》实验报告

2

是唯一的，且必然与

2.拉格朗日插值法

的麦克劳林级数一致。

如对实践中的某个物理量进行观测，在若干个不同的地方得到相应的观测值，拉格朗日插值法可以找到一个多项式，其恰好在各个观测的点取到观测到的值。这样的多项式称为拉格朗日（插值）多项式。数学上来说，拉格朗日插值法可以给出一个恰好穿过二维平面上若干个已知点的多项式函数。

在平面上有(x1,y1)(x2,y2)...(xn,yn)共n个点，现作一条函数f（x）使其图像经过这n个点。

作n个多项式pi(x),i=1,2,3...,n,使得

最后可得

3.牛顿插值法

插值法利用函数f (x)在某区间中若干点的函数值，作出适当的特定函数，在这些点上取已知值，在区间的其他点上用这特定函数的值作为函数f (x)的近似值。如果这特定函数是多项式，就称它为插值多项式。利用插值基函数很容易得到拉格朗日插值多项式，公式结构紧凑，在理论分析中甚为方便，但当插值节点增减时全部插值基函数均要随之变化，整个公式也将发生变化， 这在实际计算中是很不方便的，为了克服这一缺点，提出了牛顿插值。

牛顿插值通过求各阶差商，递推得到的一个公式：

PN(x)?PN?1(x)?aN(x?x0)(x?x1)(x?x2)牛顿插值与拉格朗日插值具有唯一性。

(x?xN?1)

4.帕德逼近

它不仅与逼近论中其他许多方法有着密切的关系，而且在实际问题特别是许多物理问题中有着广泛的应用。设是在原点某邻域内收敛的、具有复系数的麦克劳林级数。欲确定一个有理函数，式中,使得前次方的系数为0，即使得 此处约定qk＝0（k>n）。虽然所求得的Pm(z)和Qn(z)不惟一，但是比式却总是惟一的。有理函数称为F(z)的(m,n)级帕德逼近,记为(m/n)。由(m/n)所形成的阵列称为帕德表。

《数值方法》实验报告

3

三、实验内容

1.P154.1:

用plot命令，在同一幅图中绘制区间-1≤x≤1上的sin(x)，以及P5(x),P7(x)和P9(x)。 其中：

x3x5P5(x)?x?? ·········································式（1） 3!5!x3x5x7P7(x)?x??? ····································式（2） 3!5!7!x3x5x7x9P9(x)?x???? ····························式（3） 3!5!7!9!

2.P171.2: 下表（表一）给出了11月8号美国洛杉矶的一个郊区在5小时内的测量温度。 （a）对表一中的数据构造一个拉格朗日插值多项式。 （b）估计这5小时内的平均温度。 （c）在同一坐标系中画出表中的数据和由（a）得到的多项式。讨论用（a）中的多项式计算平均温度可能产生的误差。 时间（下午） 华氏度 1 66 2 66 3 65 4 64 5 63 6 63 表1 3.P178.1:

用牛顿插值多项式计算实验P171.2的内容。

4.P194.1:

比较对于函数f(x)?ex的逼近：