实验二 拉格朗日插值龙格现象
数值分析实验报告--实验2--插值法
1 / 21数值分析实验二：插值法1 多项式插值的震荡现象1.1 问题描述考虑一个固定的区间上用插值逼近一个函数。
显然拉格朗日插值中使用的节点越多，插值多项式的次数就越高。
我们自然关心插值多项式的次数增加时， 是否也更加靠近被逼近的函数。
龙格（Runge ）给出一个例子是极著名并富有启发性的。
设区间[-1,1]上函数21()125f x x=+ (1)考虑区间[-1,1]的一个等距划分，分点为n i nix i ,,2,1,0,21 =+-= 则拉格朗日插值多项式为201()()125nn ii iL x l x x ==+∑(2)其中的(),0,1,2,,i l x i n =是n 次拉格朗日插值基函数。
实验要求：(1) 选择不断增大的分点数目n=2, 3 …. ，画出原函数f(x)及插值多项式函数()n L x 在[-1，1]上的图像，比较并分析实验结果。
(2) 选择其他的函数，例如定义在区间[-5，5]上的函数x x g xxx h arctan )(,1)(4=+=重复上述的实验看其结果如何。
(3) 区间[a,b]上切比雪夫点的定义为 (21)cos ,1,2,,1222(1)k b a b ak x k n n π⎛⎫+--=+=+ ⎪+⎝⎭(3)以121,,n x x x +为插值节点构造上述各函数的拉格朗日插值多项式，比较其结果,试分析2 / 21原因。
1.2 算法设计使用Matlab 函数进行实验, 在理解了插值法的基础上，根据拉格朗日插值多项式编写Matlab 脚本，其中把拉格朗日插值部分单独编写为f_lagrange.m 函数，方便调用。
1.3 实验结果1.3.1 f(x)在[-1,1]上的拉格朗日插值函数依次取n=2、3、4、5、6、7、10、15、20，画出原函数和拉格朗日插值函数的图像，如图1所示。
Matlab 脚本文件为Experiment2_1_1fx.m 。
可以看出，当n 较小时，拉格朗日多项式插值的函数图像随着次数n 的增加而更加接近于f(x)，即插值效果越来越好。
实验二 插值法
实验二 插值法实验2.1（多项式插值的振荡现象）问题提出：考虑一个固定的区间上用插值逼近一个函数。
显然拉格朗日插值中使用的节点越多，插值多项式的次数就越高。
我们自然关心插值多项式的次数增加时，)(x L n 是否也更加靠近被逼近的函数。
龙格（Runge ）给出一个例子是极著名并富有启发性的。
设区间[-1,1]上函数22511)(xx f +=实验内容：考虑区间[-1,1]的一个等距划分，分点为 n i nix i ,,2,1,0,21 =+-= 则拉格朗日插值多项式为 201()()125nn ii iL x l x x ==+∑其中的n i x l i ,,2,1,0),( =是n 次拉格朗日插值基函数。
实验要求：（1） 选择不断增大的分点数目n=2,3….，画出原函数f(x)及插值多项式函数)(x L n 在[-1，1]上的图像，比较并分析实验结果。
（2）选择其他的函数，例如定义在区间[-5，5]上的函数x x g x xx h arctan )(,1)(4=+=重复上述的实验看其结果如何。
（3）区间[a,b]上切比雪夫点的定义为 1,,2,1,)1(2)12(cos 22+=⎪⎪⎭⎫⎝⎛+--++=n k n k a b a b x k π 以121,,+n x x x 为插值节点构造上述各函数的拉格朗日插值多项式，比较其结果,试分析原因。
实验2.2（样条插值的收敛性）问题提出：多项式插值是不收敛的，即插值的节点多，效果不一定就好。
对样条函数插值又如何呢？理论上证明样条插值的收敛性是比较困难的，但通过本实验可以验证这一理论结果。
实验内容：请按一定的规则分别选择等距或者非等距的插值节点，并不断增加插值节点的个数。
考虑实验2.1中的函数或选择其他你有兴趣的函数。
实验要求：（1）随节点个数增加，比较被逼近函数和样条插值函数误差的变化情况。
分析所得结果并与拉格朗日多项式插值比较(可以用MATLAB 的函数“spline”作此函数的三次样条插值，取n=10、20，分别画出插值函数及原函数的图形)。
实验二、拉格朗日插值法
‘linear’ ： 线性插值； ‘spline’ ： 三次样条插值； ‘cubic’ ： 立方插值。 缺省时： 分段线性插值。
9
例：在1-12的11小时内，每隔1小时测量一次温度，测 得的温度依次为：5，8，9，15，25，29，31，30，22，25， 27，24。试估计每隔1/10小时的温度值。
15
例：测得平板表面3*5网格点处的温度分别为： 82 81 80 82 84 79 63 61 65 81 84 84 82 85 86 试作出平板表面的温度分布曲面z=f(x,y)的图形。 1.先在三维坐标画出原始数据，画出粗糙的温度分布曲图. 输入以下命令： x=1:5; y=1:3; [x,y]=meshgrid(x,y); t=[82 81 80 82 84;79 63 61 65 81;84 84 82 85 86]; mesh(x,y,t) 2．以平滑数据,在x、y方向上每隔0.2个单位的地方进行插值.
课堂任务：
1、用MATLAB实现拉格朗日插值法； 2、用MATLAB分析龙格现象，绘制教材40页图2-5； 3、用自己编写的程序做课后题: 48页：2题. 50页：3题 (1).
用MATLAB作插值计算
一维插值函数：
yi=interp1(x，y，xi，'method')
xi处的插 值结果 插值节点 被插值点 插值方法
h=1:12; t=[5 8 9 15 25 29 31 30 22 25 27 24]; h1=1:0.1:12;
t1=interp1(h,t,h1,'spline')
plot(h1,t1)
10
例 已知飞机机翼下轮廓线上数据如下，求x每改变0.1时的y 值。
计算方法拉格朗日插值
2012~ 2013学年第 2学期 计算方法 教案 计1101/02，1181 开课时间：2012-02第二章 插值法知识点：拉格朗日插值法，牛顿插值法，误差，龙格现象，分段插值。
1.背景实践活动中，表现事务变化的信息往往只是一些离散点值，例如 每个6小时记录一次温度，以此反映一天的气温变化状况，如下表图能从已知这些离散点值信息知道10时的气温是多少吗？如果能通过这些离散点值找到气温变化的规律，也就是说能找到一个反映气温变化规律的“原”函数，就可以知道10时的气温是多少。
但我们能采集到的信息只有这些离散点值，时常给不出反映气温变化规律“原”函数的解析表达式，怎么办？通常可以用近似的办法解决这个问题，办法是构造一个通过所有离散点值的“近似”函数，用这个“近似”函数逼近“原”函数。
如图构造这个“近似”函数的方法称为插值方法。
时间/时温度。
C24.3623.81231.51831.024253432302826242220时间（时）温度（。
C ）3432302826242220温度（。
C ）2012~ 2013学年第 2学期计算方法教案计1101/02，1181 开课时间：2012-022.概念实际问题中，能采集到的信息只是一些离散点值{x i,f(x i)}(i=0,1,2,…n)，时常给不出一个函数f（x）的解析表达式，因之，转而考虑选择一个简单的函数ϕ（x）近似替代（原来）f（x）。
定义：设f（x）为定义在区间[a，b]上的函数，x0,x1,…,x n为[a，b]上的互异点，y i=f（x i）。
若存在一个简单函数ϕ（x），满足（插值条件）ϕ（x i）=f（x i），i=0，1，…，n。
则称 ϕ（x）为f（x）插值函数，f（x）为被插函数,点x0,x1,…,x n为插值节点,点{x i,f(x i)},i=0,1,2,…n为插值点。
若用ϕ（x）≈f（x），则计算f（x）就转换为计算 ϕ（x）。
插值需要解决：插值函数是否存在唯一；插值函数如何构造；插值函数与被插函数的误差估计和收敛性。
数值分析实验报告2——Runge现象
数值分析课程实验报告——插值逼近题目一.Runge 函数的插值1. Runge 函数Runge 函数的表达式为：21()125R x x =+ 其在[-1,1]区间上的函数图像如图1.1。
在课程学习中我们知道，对Runge 函数进行高次插值时有可能在两端出现不收敛的情况，即Runge 现象。
下面将分别用四种不同的插值方法在[-1,1]区间上对Runge 函数进行插值，并分析是否产生Runge 现象，比较插值效果。
图1.1.Runge 函数在[-1,1]区间的函数图像2.Newton 插值首先根据课本上的Newton 插值算法进行编程（代码略）。
核心思想就是用符号变量进行中间运算，以便将最终的插值函数用符号表达式表示出来，并进一步生成图像。
此处插值节点选择为等距插值节点，即：0.1(0,1,2,,)i x ih i =-+= (20)其中h=0.1。
插值曲线与原曲线的对比如图1.2（蓝色为原曲线，红色为插值曲线）。
从图中看出，在区间中部，二者吻合较好；但在区间两端二者则产生了明显偏差，甚至可以达到一个非常大的数值（e20量级）。
因此，在等距节点的20次Newton 插值下，产生了明显的Runge 现象。
图1.2.Newton 插值曲线与原曲线对比3. Lagrange 插值此处同样是根据Lagrange 插值的具体算法进行编程。
但插值节点不再是等距分布，而是如下形式：21cos()(0,1,2,,)42i i x i π+==…20 插值曲线与原曲线的对比如图1.3（蓝色为原曲线，红色为插值曲线）。
从图中看出，插值曲线与原曲线吻合的很好，没有产生明显的Runge 现象。
对比产生了明显Runge 现象的20次Newton 插值，Lagrange 插值的最高次数虽然也是20，但由于此处的插值节点不是等距分布的（事实上，此处采用的插值节点正是Chebyshev 多项式的零点），而是中间疏两边密，因此两侧较密的节点很好地抑制了Runge 现象。
拉格朗日龙格现象
实 验 报 告一．实验名称：拉格朗日插值的龙格现象二．实验目的：理解高阶插值的病态性，观察拉格朗日插值的龙格现象。
三．实验内容：在区间[－5，5]上取节点数n=11，等距离h=1的节点为插值点，对于函数25()1f x x =+ 进行拉格朗日插值，把f(x)与插值多项式的曲线花在同一张图上。
四． 实验基础知识及原理：1）拉格朗日插值函数定义:对某个多项式函数，已知有给定的k + 1个取值点：其中对应著自变数的位置，而对应著函数在这个位置的取值。
假设任意两个不同的x j 都互不相同，那麼应用拉格朗日插值公式所得到的拉格朗日插值多项式为：其中每个为拉格朗日基本多项式（或称插值基函数），其表达式为：[3] 拉格朗日基本多项式的特点是在 上取值为1，在其它的点 上取值为0。
2）龙格现象：在计算方法中，有利用多项式对某一函数的近似逼近，这样，利用多项式就可以计算相应的函数值。
一般情况下，多项式的次数越多，需要的数据就越多，而预测也就越准确。
但是在有的情况下，并非取节点（日期数）越多多项式就越精确。
3）matlab:MATLAB 是矩阵实验室（Matrix Laboratory ）的简称，是美国MathWorks 公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB 和Simulink 两大部分。
该实验中可利用matlab 的绘图功能将拉格朗日函数跟原函数显示出来，通过对比得出实验结论。
五. 具体实验过程1. 用matlab的editor将老师提供的拉格朗日插值的调用函数编辑为一个m文件lagrange.m2.在脚本窗口中输入以下命令:执行上面的命令可以得到输出的曲线如下：在以下图中可以看到相应变量的数据值：通过龙格现象可知，并不是插值多项式的次数越高（即插值节点越多）精度就越高，在这个实验中插值函数在两个端点处发生剧烈的波动，造成较大的误差，从数值计算上可解释为高次插值多项式的计算会带来舍入误差的增大，从而引起计算失真，因此，实际应用时，我们一般只用一次、二次最多是三次插值多项式，若想提高插值精度，我们可以采用分段插值。
拉格朗日多项式插值
拉格朗日多项式插值法浅析摘要拉格朗日插值多项式是一种最常见的多项式插值法，也是一种最常用的逼近工具。
“学以致用 ”是每一门学科都致力追求的境界，数学自然也不例外。
下面，探讨拉格朗日插值法的基本原理、如何构造拉格朗日多项式、拉格朗日多项式的误差界，并用 MATLAB 程序来实现这一数学算法的自动化，为复杂的分析研究提供了一条数学算法的捷径。
【关键词】：拉格朗日多项式 算法实现 MATLAB在科学研究和实际的工程设计中，几乎所有的问题都可以用)(x f y =来表示其某种内在规律的数量关系。
但理想化的函数关系在实际工程应用中是很难寻找 的，对于那些没有明显解析式的函数关系表达式则只能通过实验观察的数据，利用多项式对某一函数的进行逼近，使得这个逼近函数能够反映)(x f 的特性，而且利用多项式就可以简便的计算相应的函数值。
例如我们不知道气温随日期变化的具体函数关系，但是我们可以测量一些孤立的日期的气温值，并假定此气温随日期变化的函数满足某一多项式。
这样，利用已经测的数据，应用待定系数法便可以求得一个多项式函数f （x ）。
应用此函数就可以计算或者说预测其他日期的气温值。
一般情况下，多项式的次数越多，需要的数据就越多，而预测也就越 准确。
当然，构造组合多项式方法比较多，如线性方程求解、拉格朗日系数多项式以及构造牛顿多项式的分段差分和系数表等等，这里只对拉格朗日多项式插值法进行深入探讨。
一、拉格朗日多项式插值算法基本原理函数)(x f y =在区间[a,b]上有定义,在是[ a,b]上取定的 N + 1个互异节点, 且在这些点处的函数值)(0x f , )(1x f ,…,)(n x f 为已知, 即 yi =f (xi ) , (N i ...1,0=)，若存在一个和)(x f 近似的函数)(x P N ，满足)()(i i N x f x P = (N i ...1,0=) （1）则称 φ(x) 为 f (x) 的一个插值函数, 点i x 为插值节点，（1）称为插值条件, 区间[a,b]称为插值区间, 而误差函数)()(x P x f E N N -=称为插值余项。
拉格朗日插值公式的证明及其应用
拉格朗日插值公式的证明及其应用一、拉格朗日插值公式的证明：假设给定n+1个不同的数据点(x0, y0), (x1, y1), ..., (xn, yn)，其中xi不等于xj，i≠j。
我们要找到一个满足这些数据点的多项式函数P(x)，使得P(xi) = yi，i = 0, 1, ..., n。
设P(x)的表达式为P(x) = a0 + a1x + a2x^2 + ... + anxn，其中a0, a1, ..., an是待确定的系数。
由于希望P(x)满足P(xi) = yi，可以得到以下等式：a0 + a1x0 + a2x0^2 + ... + anx0^n = y0a0 + a1x1 + a2x1^2 + ... + anx1^n = y1...a0 + a1xn + a2xn^2 + ... + anx0n = yn将这些等式进行展开，可以得到如下的一个线性方程组：a0 + a1x0 + a2x0^2 + ... + anx0^n = y0a0 + a1x1 + a2x1^2 + ... + anx1^n = y1...a0 + a1xn + a2xn^2 + ... + anx0n = yn我们令Li(x)表示一个满足Li(xi) = 1，Li(xj) = 0 (j≠i)的多项式函数。
将P(x)的表达式代入上述的线性方程组，我们可以得到以下等式：y0Li(x) + y1Li(x) + ... + ynLi(x) = P(x)将P(x)表示成等于数据点的线性组合的形式，即拉格朗日插值公式的形式。
二、拉格朗日插值公式的应用：1.数据拟合：拉格朗日插值公式可以通过已知的数据点，得到一个满足数据点的多项式函数。
通过拟合已有数据，可以进行数据的预测和预估。
2.函数逼近：对于已知的函数，可以通过拉格朗日插值公式插值得到一系列的数据点。
这样可以将原函数进行逼近，并在所插入的数据点上进行具体的计算。
3.误差估计：通过拉格朗日插值公式得到的多项式函数可以作为原函数的近似函数。
计算方法 实验报告 拉格朗日 龙贝格 龙格库塔
主界面：
/*lagrange.c*/
float real_value(float x) /*由被插值函数计算真实值*/
c=getchar();
if(c=='N'||c=='n') break;
}
}
/*romberg.c*/
double function(double x) /*被积函数*/
{
return 4.0/(1+(x)*(x));
}
double t(double a,double b,int m) /*计算T1*/
实验二（龙贝格公式）
§公式
§算法描述
§流程图
§运行结果
§结果分析：Romberg积分法是在积分步长逐步折半的过程中，用低精度求积公式的组合得到更高精度求积公式的一种方法，它算法简单，且收敛加速效果极其显著。
实验三（四阶龙格库塔）
§公式
k1=h*f(xn,yn);
k2=h*f(xn+h/2,yn+k1/2);
T1=t(a,b,0);
T2=T1/2.0+t(a,b,1);
S1=(4*T2-T1)/3.0;
T1=T2;
T2=T1/2.0+t(a,b,2);
S2=(4*T2-T1)/3.0;
C1=(16*S2-S1)/15.0;
T1=T2;
T2=T1/2.0+t(a,b,3);
S1=S2;
S2=(4*T2-T1)/3.0;
数值分析实验报告--实验2--插值法
1 / 21数值分析实验二：插值法1 多项式插值的震荡现象1.1 问题描述考虑一个固定的区间上用插值逼近一个函数。
显然拉格朗日插值中使用的节点越多，插值多项式的次数就越高。
我们自然关心插值多项式的次数增加时， 是否也更加靠近被逼近的函数。
龙格（Runge ）给出一个例子是极著名并富有启发性的。
设区间[-1,1]上函数21()125f x x=+ (1)考虑区间[-1,1]的一个等距划分，分点为n i nix i ,,2,1,0,21 =+-= 则拉格朗日插值多项式为201()()125nn ii iL x l x x ==+∑(2)其中的(),0,1,2,,i l x i n =是n 次拉格朗日插值基函数。
实验要求：(1) 选择不断增大的分点数目n=2, 3 …. ，画出原函数f(x)及插值多项式函数()n L x 在[-1，1]上的图像，比较并分析实验结果。
(2) 选择其他的函数，例如定义在区间[-5，5]上的函数x x g xxx h arctan )(,1)(4=+=重复上述的实验看其结果如何。
(3) 区间[a,b]上切比雪夫点的定义为 (21)cos ,1,2,,1222(1)k b a b ak x k n n π⎛⎫+--=+=+ ⎪+⎝⎭(3)以121,,n x x x +为插值节点构造上述各函数的拉格朗日插值多项式，比较其结果,试分析2 / 21原因。
1.2 算法设计使用Matlab 函数进行实验, 在理解了插值法的基础上，根据拉格朗日插值多项式编写Matlab 脚本，其中把拉格朗日插值部分单独编写为f_lagrange.m 函数，方便调用。
1.3 实验结果1.3.1 f(x)在[-1,1]上的拉格朗日插值函数依次取n=2、3、4、5、6、7、10、15、20，画出原函数和拉格朗日插值函数的图像，如图1所示。
Matlab 脚本文件为Experiment2_1_1fx.m 。
可以看出，当n 较小时，拉格朗日多项式插值的函数图像随着次数n 的增加而更加接近于f(x)，即插值效果越来越好。
计算方法实验二
《计算方法》实验报告实验二插值法二级学院：计算机学院专业：计算机科学与技术指导教师：爨莹班级学号：姓名：实验二插值法1、实验目的：1、掌握直接利用拉格郎日插值多项式计算函数在已知点的函数值；观察拉格郎日插值的龙格现象。
2、了解Hermite插值法、三次样条插值法原理，结合计算公式，确定函数值。
2、实验要求:1、认真分析题目的条件和要求，复习相关的理论知识，选择适当的解决方案和算法；2、编写上机实验程序，作好上机前的准备工作；3、上机调试程序，并试算各种方案，记录计算的结果（包括必要的中间结果）；4、分析和解释计算结果；5、按照要求书写实验报告；3、实验内容：1、用拉格郎日插值公式确定函数值；对函数f(x)进行拉格郎日插值，并对f(x)与插值多项式的曲线作比较。
2、已知函数表：（0.56160,0.82741）、（0.56280,0.82659）、（0.56401,0.82577）、（0.56521,0.82495）用三次拉格朗日插值多项式求x=0.5635时函数近似值。
3、P129,(12)4、题目：插值法5、原理：拉格朗日插值法是以法国十八世纪数学家约瑟夫·路易斯·拉格朗日命名的一种多项式插值方法。
许多实际问题中都用函数来表示某种内在联系或规律，而不少函数都只能通过实验和观测来了解。
如对实践中的某个物理量进行观测，在若干个不同的地方得到相应的观测值，拉格朗日插值法可以找到一个多项式，其恰好在各个观测的点取到观测到的值。
这样的多项式称为拉格朗日（插值）多项式。
数学上来说，拉格朗日插值法可以给出一个恰好穿过二维平面上若干个已知点的多项式函数。
6、设计思想：通过拉格朗日插值法找到拉格朗日多项式，计算某个x对应的y值。
7、对应程序：//Lagrange.cpp#include <stdio.h>#include <conio.h>#define N 4int checkvalid(double x[], int n);void printLag (double x[], double y[], double varx, int n);double Lagrange(double x[], double y[], double varx, int n);void main (){double x[N+1] = {0.56160, 0.56280, 0.56401, 0.56521};double y[N+1] = {0.82741, 0.82659, 0.82577, 0.82459};double varx = 0.5635;if (checkvalid(x, N) == 1){printf("\n\n插值结果: P(%f)=%f\n", varx, Lagrange(x, y, varx, N));}else{printf("结点必须互异");}getch();}int checkvalid (double x[], int n){int i,j;for (i = 0; i < n; i++){for (j = i + 1; j < n+1; j++){if (x[i] == x[j])//若出现两个相同的结点，返回-1{return -1;}}}return 1;}double Lagrange (double x[], double y[], double varx, int n){double fenmu;double fenzi;double result = 0;int i,j;printf("Ln(x) =\n");for (i = 0; i < n+1; i++){fenmu = 1;for (j = 0; j < n+1; j++){if (i != j){fenmu = fenmu * (x[i] - x[j]);}}printf("\t%f", y[i] / fenmu);fenzi = 1;for (j = 0; j < n+1; j++){if (i != j){printf("*(x-%f)", x[j]);fenzi = fenzi * (varx - x[j]);}}if (i != n){printf("+\n");}result += y[i] / fenmu * fenzi;}return result;}8、实验结果：9、实验体会：拉格朗日插值法能很好的解决生产和科研中的一些问题，计算机编程的实现有利于其更好的应用。
拉格朗日插值
拉格朗日插值绘制龙格现象一、问题叙述龙格反例1/(1+x^2)说明高次代数插值会导致误差很大。
在区间[-5,5]上取等距结点构造10次拉格朗日插值多项式用计算机绘制图形显示龙格现象。
二、理论分析1. 拉格朗日插值：假设有（n+1）个拉格朗日插值结点01n x x x <<< ，已知函数值0011(x ),y (),,()n n y f f x y f x ===求n 次多项式()n L x 使其满足插值条件()(j 0,1,,)j j f x y n ==类似于二次插值方法，根据插值结点构造(n+1)个拉格朗日插值基函数)())(()()())(()()(110110n k k k k k k n k k k x x x x x x x x x x x x x x x x x l --------=+-+- 每一个基函数都是零点多项式1(),(0,1)0k j j k l x j n j k =⎧==⎨≠⎩()n L x 满足插值条件 L ()()(j 0,1,,)n j j x f x n == 拉格朗日插值基函数：0()()()n j k j k j j k x x l x x x =≠-=-∏拉格朗日插值多项式：0()y nn j j j L l x ==∑2. 切比雪夫插值：n 阶切比雪夫多项式定义为()cos(arccos ),(0,1,),x [1,1]n T x n x n ==∈-若令arccos x θ= ，则有cos n T n θ= 。
由余弦函数性质，有cos(1)cos(1)2cos cos n n n θθθθ++-=所以有递推关系 1()1()2(),(1,2)n n x n T x xT T x n +-=-=，又有011,T T x ==，所以n 阶切比雪夫多项式零点为21cos,(0,1,,1)2k k x k n n π+==- 。
3. Hermite 插值如果f(x)在区间[a,b]上连续可导，01,,x x L [,]n x a b ∈ 是互异的，那么存在唯一的多项式21()n H x +满足多项式在这些点上的值与函数f(x)的值相等、多项式在这些点的一阶导数值与函数的一阶导数值相等。
拉格朗日插值讲解
特征提取
在计算机视觉中，拉格朗日插值可以用于提 取图像中的特征点，为后续的图像识别和分 析提供基础。
THANKS FOR WATCHING
感谢您的观看
02
它是由意大利数学家约瑟夫·拉格 朗日于18世纪提出的一种数学工 具，广泛应用于科学、工程和经 济学等领域。
拉格朗日插值的重要性
拉格朗日插值方法为数据分析和预测提供了一种重要的工具 ，特别是在数据量较小或数据分布不均匀的情况下，可以通 过插值方法来填补数据空白或提高数据精度。
它可以帮助我们更好地理解数据的内在规律和趋势，为决策 提供科学依据。
基于拉格朗日插值拟合出的多项 式，可以进一步预测未来数据点 的趋势和走向，为决策提供依据 。
工程计算与设计
工程建模
在工程计算中，拉格朗日插值可以用 于建立数学模型，模拟复杂系统的行 为和性能。
优化设计
通过拉格朗日插值，工程师可以对设 计方案进行优化，提高产品的性能和 效率。
图像处理与计算机视觉
图像修复
多项式插值的精度较高，适用于数据 点之间变化较大的情况，但构造多项 式的过程较为复杂，需要选择合适的 基函数和节点。
拉格朗日插值公式
拉格朗日插值公式是利用拉格朗日多 项式进行插值的方法，通过已知数据 点构造拉格朗日多项式，然后利用这 个多项式计算出需要插值的点的值。
拉格朗日插值公式的优点是构造简单、 精度较高，适用于任意数据点的情况， 但当数据点较多时，计算量较大，可 能会出现龙格现象。
拉格朗日插值的历史背景
拉格朗日插值方法的发展经历了漫长的历史过程。最早的插值方法可以追溯到古 希腊时期，而现代的插值方法则是在17世纪和18世纪随着数学的发展而逐步完善 的。
拉格朗日插值法
Lagrange polynomial 维基百科，自由的百科全书
定义
对某个多项式函数，已知有给定的 个取值点：
其中 对应着自变量的位置，而 对应着函数在这个位置(每一个y值都不等于0)的取 值。 假设任意两个不同的 都互不相同，那么应用拉格朗日插值公式所得到的拉格朗 日插值多项式为：
其中每个
几何性质
拉格朗日插值法中用到的拉格朗日基本多项式 （由某一组 确定）可以看做是由次数不超过 n 的多项式所组成的线性空 间： 的一组基底。首先，如果存在一组系数： 使得， ， 那么，一方面多项式P是满足 的拉格 朗日插值多项式，另一方面P是零多项式，所以取值永远是0。所以 。 这证明了 是线性无关的。同时它一共包含 个多项式，恰好等于 的维数。所以 构成了 的一组基底。 拉格朗日基本多项式作为基底的好处是所有的多项式都是齐次的（都是 n 次多项 式）。
重心拉格朗日插值法
重心拉格朗日插值法是拉格朗日插值法的一种改进。在拉格朗日插值法中，运用 多项式 可以将拉格朗日基本多项式重新写为：
定义重心权
上面的表达式可以简化为：
于是拉格朗日插值多项式变为：
即所谓的重心拉格朗日插值公式（第一型）或改进拉格朗日插值公式。它的优点 是当插值点的个数增加一个时，将每个 都除以 ，就可以得到新的重 心权 ，计算复杂度为 ，比重新计算每个基本多项式所需要的复杂度 降了一个量级。 将以上的拉格朗日插值多项式用来对函数 插值，可以得到：
优点与缺点
拉格朗日插值法的公式结构整齐紧凑，在理论分析中十分方便，然而在计算中， 当插值点增加或减少一个时，所对应的基本多项式就需要全部重新计算，于是整 个公式都会变化，非常繁琐。这时可以用重心拉格朗日插值法或牛顿插值法来代 替。此外，当插值点比较多的时候，拉格朗日插值多项式的次数可能会很高，因 此具有数值不稳定的特点，也就是说尽管在已知的几个点取到给定的数值，但在 附近却会和“实际上”的值之间有很大的偏差（如右下图）。这类现象也被称为龙 格现象，解决的办法是分段用较低次数的插值多项式。
《拉格朗日插值》课件
通过对比已知数据点和插值函数的结果，评估拟合效果和预测准确度。
总结和展望
本课件介绍了拉格朗日插值的背景、原理、计算方法和应用，并总结了其优缺点。这门课程为数值分析和数据 插值领域提供了基础知识，并为未来的研究和应用提供了展望。
用于数据拟合、函数逼近和误差 估计。
图像处理
用于图像重建、图像修复和图像 插值。
信号处理
用于信号重构、信号滤波和信号 插值。
拉格朗日插值多项式的优缺点
拉格朗日插值多项式具有以下优点和缺点。 优点 简单易懂 计算效率高 适用于小范围插值
缺点 对于大量数据点计算复杂度高 对离散数据敏感 数值误差较大
示例演示：使用拉格朗日插值进行数据的 拟合
2 定义
插值问题是通过插值函数来逼近离散数据点 的函数曲线，以实现数据的拟合和预测。
拉格朗日插值多项式的原理
拉格朗日插值多项式是通过已知数据点构建的一个多项式函数，使得该函数在给定区间内的所有数据点上与已 知数据点完全一致。
插值多项式
拉格朗日插值多项式是通过拉 格朗日基函数构造出的一个多 项式函数。
《拉格朗日插值》PPT课 件
本课件旨在介绍拉格朗日插值的原理、计算方法、应用以及优缺点，帮助大 家更好地理解插值问题的背景和定义。
插值问题的背景和定义
插值问题是数值分析中的基础概念，指的是通过已知数据点构建一个函数，使得该函数在给定区间内的所有数 据点上与已知数据点完全一致。
1 背景
插值问题源于天文学和地理测量等领域中的 数据分析需求。
1
数据点坐日基函数计算
根据已知数据点的坐标值计算拉格朗日基函数。
3
插值多项式计算
将拉格朗日基函数与对应系数相乘并求和得到插值多项式。
实验二_拉格朗日插值龙格现象
汕 头 大 学 实 验 报 告学院: 工学院 系: 计算机系 专业: 计算机科学与技术 年级: 2010 姓名: 林金正 学号: 2010101032 完成实验时间: 5月24日一．实验名称：拉格朗日插值的龙格现象二．实验目的：通过matlab 处理,观察拉格朗日插值的龙格现象.三．实验内容：（1）学习matlab 的使用（2）以实验的方式，理解高阶插值的病态性，观察拉格朗日插值的龙格现象。
四．实验时间、地点，设备：实验时间：5月24日实验地点： 宿舍 实验设备：笔记本电脑五,实验任务在区间[－5，5]上取节点数n=11，等距离h=1的节点为插值点，对于函数25()1f x x =+进行拉格朗日插值，把f(x)与插值多项式的曲线花在同一张图上。
六.实验过程拉格朗日插值函数定义:对某个多项式函数，已知有给定的k + 1个取值点：其中对应著自变数的位置，而对应著函数在这个位置的取值。
假设任意两个不同的xj 都互不相同，那麼应用拉格朗日插值公式所得到的拉格朗日插值多项式为：其中每个为拉格朗日基本多项式（或称插值基函数），其表达式为：[3] 拉格朗日基本多项式的特点是在 上取值为1，在其它的点 上取值为0。
1.使用matlab,新建function.m 文件,使用老师所给代码,构建拉格朗日函数.%lagrange.mfunction y=lagrange(x0,y0,x)n=length(x0);m=length(x);for i=1:mz=x(i);s=0;for k=1:nL=1;for j=1:nif j~=kL=L*(z-x0(j))/(x0(k)-x0(j));endends=s+L*y0(k);endy(i)=s;endy;程序解释:(x0,y0):已知点坐标x:所求点的横坐标,y:由(x0,y0)所产生的插值函数,以x 为参数,所的到的值2.再一次新建function.m 文件.构建自定义函数: 25()1f x x=+ %f.mfunction y = f(x)y = 5/(1+x*x);end3.在脚本窗口中输入:>>a = [-10:0.2:10]>>for I = 1:length(a)b(i) = f(a(i))end ;%画出原函数(a,b)>>>>for i = 1:length( c)d(i) = f(c(i))end ;%获取插值坐标(c,d)>>e = [-5:0.2:5]>>z = largange(c,d,e);%获取插值坐标函数(e,z) >>plot(a,b,’r-‘,e,z);%画图过程及插图七:实验所得:这次实验是我初步学会Matlab的使用，学会新建function函数，在matlab命令窗口敲入一些基础的命令，同时更深刻地了解了拉格朗日插值的龙格现象。
《拉格朗日插值》课件
计算结果即为待插值点的 预测值。
根据需要，可以对插值结 果进行优化和调整，以提 高精度和适应性。
04
拉格朗日插值法的优缺点
Chaptபைடு நூலகம்r
优点
精确度高
拉格朗日插值法能够提供非常精 确的插值结果，尤其是在数据点 较多的情况下。
易于理解
拉格朗日插值法的原理相对简单 ，容易理解，方便教师和学生进 行教学和交流。
在数据点较少或数据变化较大的情况下， 拉格朗日插值法的插值结果可能会不稳定 ，甚至出现较大的误差。
对异常值敏感
对初值敏感
如果数据集中存在异常值，拉格朗日插值 法的插值结果可能会受到较大影响，导致 误差增大。
该方法对初值的选择较为敏感，如果初值 选择不当，可能会导致迭代不收敛或得到 不正确的插值结果。
VS
高精度计算
通过拉格朗日插值法，可以构造高精度的 数值积分公式，提高计算结果的精度。
在微分方程数值解中的应用
初值问题
在求解微分方程的初值问题时，可以使用拉格朗日插值法对解进行近似，得到数值解。
边值问题
在求解微分方程的边值问题时，可以利用拉格朗日插值法构造数值逼近方案，得到近似 解。
06
总结与展望
《拉格朗日插值》PPT课件
目录
• 引言 • 拉格朗日插值法的基本原理 • 拉格朗日插值法的实现步骤 • 拉格朗日插值法的优缺点 • 拉格朗日插值法的应用实例 • 总结与展望
01
引言
Chapter
拉格朗日插值法的背景
拉格朗日插值法是一种数学方法，用于通过已知的离散 数据点来构造一个多项式，该多项式可以用来估计未知 的数据点。 该方法由意大利数学家约瑟夫·拉格朗日于18世纪提出， 是数值分析中的重要工具之一。
龙格现象实验报告
y(i)=L;
end
a=1时，
1）取n=4
x0=[-5:2:5];
y0=5./(1+x0.^2);
x=[-5:0.1:5];
y=lagrange(x0,y0,x);
y1=5./(1+x.^2);
plot(x,y,'--r')
hold on
plot(x,y1,'-b')
hold off
2）取n=10
y1=5./(1+x.^2);
plot(x,y,'9;-b')
hold off
实
验
结
果
分
析
a=1时，
1）取n=4
a=1时，
2）取n=10
a=0.25时，
1）取n=4
a=0.25时，
2）取n=10
上述现象和定理告诉我们，并不是插值多项式的次数越高(即插值节点越多)精度越高，从数值计算上可解释为高次插值多项式的计算会带来舍入误差的增大，从而引起计算失真。因此，实际应用做插值时一般只用一次、二次最多用三次插值多项式。
x0=[-5:1:5];
y0=5./(1+x0.^2);
x=[-5:0.1:5];
y=lagrange(x0,y0,x);
y1=5./(1+x.^2);
plot(x,y,'--r')
hold on
plot(x,y1,'-b')
hold off
2、a=0.25时，
1）取n=4
x0=[-5:2:5];
y0=5./(0.25*0.25+x0.^2);
数值计算实验报告（四）
数值分析实验报告--Runge现象的产生和克服
数值分析实验报告（四）题目：Runge现象的产生和克服学院：机电工程学院（二专业）专业：机械设计制造及其自动化班级：1008108班姓名：***学号：**********号Runge现象的产生和克服摘要：对于多项式插值运算，随着插值阶数的逐渐增多，如果带入离散点过于密集，使得定义域中的“边缘区域”，没有有效的点，将导致插值函数的边缘区域大幅度的偏离函数的真值，该现象称之为“Runge现象”。
0 前言（目的与意义）：了解Runge现象，体会插值运算的不准确性，以及其差值带来的误差甚至是错误。
1 数学背景：插值运算的误差公式：|w n (x)||R n (x)|<M n+1(n+1)!M n+1=max{f(n+1)(x i)}于是，如果函数的n+1阶导数一旦很大，则会出现函数的误差很大的情况。
2 程序及代码：（1）lagrange多项式插值函数syms f x p dp lx L;f=1/(1+25*x^2);N=input('请输入插值节点数N=');xx=-1:2/N:1;p=1; L=0;ff=zeros(1,length(xx));for i=1:(N+1)x=xx(i);ff(i)=eval(f);syms x;p=p*(x-xx(i));enddp=diff(p);for j=1:(N+1)x=xx(j);k=eval(dp);syms x;lx=p/((x-xx(j))*k);L=L+lx*ff(j);endaa=[-0.96:0.1:-0.06,0,0.06:0.1:0.96];for i=1:length(aa)x=aa(i);S(i)=eval(L);fff(i)=eval(f);ende=0;for i=1:length(aa)e=e+(S(i)-fff(i))^2;ende=sqrt(e/(20*21));fprintf('插值偏差为e=%.6f\n',e)ezplot(f,[-1,1])hold onezplot(L,[-1,1])hold onplot(xx,ff,'*')hold onplot(aa,S,'o')hold off（2）分段线性插值函数syms f x p lx;f=1/(1+25*x^2);N=input('请输入插值节点数N=');xx=-1:2/N:1;p=1; L=0;ff=zeros(1,length(xx));for i=1:(N+1)x=xx(i);ff(i)=eval(f);endsyms xfor i=1:Nfor j=1:(N+1)if j==ilx(i,j)=(x-xx(i+1))/(xx(i)-xx(i+1)); else if j==i+1lx(i,j)=(x-xx(i))/(xx(i+1)-xx(i)); elselx(i,j)=0;endendendendp=lx*ff';aa=[-0.96:0.1:-0.06,0,0.06:0.1:0.96];for i=1:length(aa)x=aa(i);for j=1:N+1if x<xx(j)breakendendS(i)=eval(p(j-1));fff(i)=eval(f);ende=0;for i=1:length(aa)e=e+(S(i)-fff(i))^2;ende=sqrt(e/(20*21));fprintf('插值偏差为e=%.6f\n',e)ezplot(f,[-1,1])hold onxxx=(-1:0.01:1);for i=1:length(xxx)x=xxx(i);for j=1:N+1if x<xx(j)breakendendSS(i)=eval(p(j-1));endplot(xxx,SS,'r')hold onplot(xx,ff,'*')hold onplot(aa,S,'o')hold off（3）：三转角插值法函数syms f x df s s1 s2 s3 s4;f=1/(1+25*x^2);df=diff(f);N=input('请输入插值节点数N=');h=2/N;xx=-1:2/N:1;p=1; L=0;ff=zeros(1,length(xx));for i=1:(N+1)x=xx(i);ff(i)=eval(f);dff(i)=eval(df);endsyms xfor i=1:Ns1=(x-xx(i+1))^2*(h+2*(x-xx(i)))*ff(i)/h^3; s2=(x-xx(i))^2*(h+2*(xx(i+1)-x))*ff(i+1)/h^3; s3=(x-xx(i+1))^2*(x-xx(i))*dff(i)/h^2;s4=(x-xx(i))^2*(x-xx(i+1))*dff(i+1)/h^2;s(i)=s1+s2+s3+s4;endaa=[-0.96:0.1:-0.06,0,0.06:0.1:0.96]; for i=1:length(aa)x=aa(i);for j=1:N+1if x<xx(j)breakendendS(i)=eval(s(j-1));fff(i)=eval(f);ende=0;for i=1:length(aa)e=e+(S(i)-fff(i))^2;ende=sqrt(e/(20*21));fprintf('插值偏差为e=%.6f\n',e) ezplot(f,[-1,1])hold onxxx=(-1:0.01:1);for i=1:length(xxx)x=xxx(i);for j=1:N+1if x<xx(j)breakendendSS(i)=eval(s(j-1));endplot(xxx,SS,'r')hold onplot(xx,ff,'*')hold onplot(aa,S,'o')hold off（4）.三弯矩插值法函数：syms f x ddf s;f=1/(1+25*x^2);ddf=diff(diff(f));N=input('请输入插值节点数N=');h=2/N;xx=-1:2/N:1;p=1; L=0;ff=zeros(1,length(xx));for i=1:(N+1)x=xx(i);ff(i)=eval(f);ddff(i)=eval(ddf);endsyms xfor i=1:NA=(ff(i+1)-ff(i))/h-h*(ddff(i+1)-ddff(i))/6;B=ff(i)-h^2*ddff(i)/6;s(i)=(xx(i+1)-x)^3*ddff(i)/(6*h)+(x-xx(i))^3*ddff(i+1)/(6*h)+A*(x-xx(i))+B; endaa=[-0.96:0.1:-0.06,0,0.06:0.1:0.96];for i=1:length(aa)x=aa(i);for j=1:N+1if x<xx(j)breakendendS(i)=eval(s(j-1));fff(i)=eval(f);ende=0;for i=1:length(aa)e=e+(S(i)-fff(i))^2;ende=sqrt(e/(20*21));fprintf('插值偏差为e=%.6f\n',e)ezplot(f,[-1,1])hold onxxx=(-1:0.01:1);for i=1:length(xxx)x=xxx(i);for j=1:N+1if x<xx(j)breakendendSS(i)=eval(s(j-1));endplot(xxx,SS,'r')hold onplot(xx,ff,'*')hold onplot(aa,S,'o')hold off3 总结与评价：函数的Runge现象可以通过三转角插值和三弯矩插值来解决，而且对于三转角和三弯矩插值来说，带入的数据越多，其插值效果越好4 实验结果：图1：观察Runge现象图2：分段线性插值图3：三转角插值：图4：三弯矩插值：。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
汕 头 大 学 实 验 报 告
学院: 工学院系: 计算机系专业: 计算机科学与技术年级:2010
姓名: 林金正学号:2010101032完成实验时间: 5月24日
一．实验名称：拉格朗日插值的龙格现象
二．实验目的：
通过matlab 处理,观察拉格朗日插值的龙格现象.
三．实验内容：
（1）学习matlab 的使用
（2）以实验的方式，理解高阶插值的病态性，观察拉格朗日插值的龙格现象。
四．实验时间、地点，设备：
实验时间：5月24日
实验地点：宿舍 实验设备：笔记本电脑
五,实验任务
在区间[－5，5]上取节点数n=11，等距离h=1的节点为插值点，对于函数2
5()1f x x =+进行拉格朗日插值，把f(x)与插值多项式的曲线花在同一张图上。
六.实验过程
拉格朗日插值函数定义:
对某个多项式函数，已知有给定的k + 1个取值点：
其中对应著自变数的位置，而对应著函数在这个位置的取值。
假设任意两个不同的xj 都互不相同，那麼应用拉格朗日插值公式所得到的拉格朗日插值多项式为：
其中每个为拉格朗日基本多项式（或称插值基函数），其表达式为：
[3]
拉格朗日基本多项式
的特点是在上取值为1，在其它的点上取值为0。
1.使用matlab,新建function.m 文件,使用老师所给代码,构建拉格朗日函数.
%lagrange.m
function y=lagrange(x0,y0,x)
n=length(x0);
m=length(x);
fori=1:m
z=x(i);s=0;
for k=1:n
L=1;
for j=1:n
if j~=k
L=L*(z-x0(j))/(x0(k)-x0(j));
end
end
s=s+L*y0(k);
end
y(i)=s;
end
y;
程序解释:
(x0,y0):已知点坐标
x:所求点的横坐标,
y:由(x0,y0)所产生的插值函数,以x 为参数,所的到的值
2.再一次新建function.m 文件.
构建自定义函数:25()1f x x
=
+ %f.m
function y = f(x)
y = 5/(1+x*x);
end
3.在脚本窗口中输入:
>>a = [-10:0.2:10]
>>for I = 1:length(a)
b(i) = f(a(i))
end ;%画出原函数(a,b)
>>c = [-5:1:5]
>>for i = 1:length( c)
d(i) = f(c(i))
end;%获取插值坐标(c,d)
>>e = [-5:0.2:5]
>>z = largange(c,d,e);%获取插值坐标函数(e,z) >>plot(a,b,’r-‘,e,z);%画图
过程及插图
七:实验所得:
这次实验是我初步学会Matlab的使用，学会新建function函数，在matlab命令窗口敲入一些基础的命令，同时更深刻地了解了拉格朗日插值的龙格现象。