拉格朗日插值算法及应用实验报告

合集下载

数学建模实验报告(3)插值

数学建模实验报告(3)插值

数学模型实验报告——插值专业:姓名:李学号:姓名:刘学号:姓名:汪学号:数学模型实验报告(插值)一、 实验目的:1、了解插值的基本内容。

2、掌握用数学软件包求解插值问题。

二、实验内容:(一)一维插值一、插值的定义 已知n+1个节点,,1,0(),(n j y x j j =其中 j x 互不相同,不妨设),10b x x x a n =<<<= 求任一插值点 )(*j x x ≠处的插值.*y构造一个(相对简单的)函数),(x f y =通过全部节点, 即 ),1,0()(n j y x f j j ==再用)(x f 计算插值,即).(**x f y =二、插值的方法拉格朗日(Lagrange)插值已知函数f (x )在n +1个点x 0,x 1,…,xn 处的函数值为 y 0,y 1,…,yn 。

求一n 次多项式函数Pn (x ),使其满足:Pn (xi )=yi ,i =0,1,…,n .解决此问题的拉格朗日插值多项式公式如下∑=⋅=ni i i n y x L x P 0)()(其中Li (x ) 为n 次多项式:)())(())(()())(())(()(11101110n i i i i i i i n i i i x x x x x x x x x x x x x x x x x x x x x L ----------=+-+-称为拉格朗日插值基函数。

特别地:两点一次(线性)插值多项式:()101001011y x x x x y x x x x x L --+--=三点二次(抛物)插值多项式:()()()()()()()()()()()()()2120210121012002010212y x x x x x x x x y x x x x x x x x y x x x x x x x x x L -⋅--⋅-+-⋅--⋅-+-⋅--⋅-=().,满足插值条件直接验证可知x L n例55,11)(2≤≤-+=x xx g 采用拉格朗日多项式插值:选取不同插值节点个数n +1,其中n 为插值多项式的次数,当n 分别取2,4,6,8,10时,绘出插值结果图形.拉格朗日多项式插值的这种振荡现象叫 Runge 现象 解:编写M 文件程序如下: m=101;x=-5:10/(m-1):5; y=1./(1+x.^2);z=0*x;plot(x,z,'r',x,y,'LineWidth',1.5), gtext('y=1/(1+x^2)'),pause n=3; x0=-5:10/(n-1):5; y0=1./(1+x0.^2); y1=lagr1(x0,y0,x); hold on ,plot(x,y1,'b'),gtext('n=2'),pause,hold off n=5; x0=-5:10/(n-1):5; y0=1./(1+x0.^2); y2=lagr1(x0,y0,x); hold on ,plot(x,y2,'b:'),gtext('n=4'),pause,hold offn=7;x0=-5:10/(n-1):5; y0=1./(1+x0.^2);y3=lagr1(x0,y0,x);hold on , plot(x,y3,'r'),gtext('n=6'), pause,hold off n=9; x0=-5:10/(n-1):5; y0=1./(1+x0.^2); y4=lagr1(x0,y0,x);hold on ,plot(x,y4,'r:'),gtext('n=8'),pause,hold off n=11; x0=-5:10/(n-1):5; y0=1./(1+x0.^2); y5=lagr1(x0,y0,x);hold on , plot(x,y5,'m'),gtext('n=10')分段线性插值计算量与n 无关; n 越大,误差越小.n n n x x x x g x L ≤≤=∞→0),()(lim例66,11)(2≤≤-+=x xx g 用分段线性插值法求插值,并观察插值误差. 1. 在[-6,6]中平均选取5个点作插值 2. 在[-6,6]中平均选取11个点作插值 3. 在[-6,6]中平均选取21个点作插值 4. 在[-6,6]中平均选取41个点作插值 解:编写M 文件程序如下:x=linspace(-6,6,100); y=1./(x.^2+1);x1=linspace(-6,6,5);%第三个参数表示插值点的个数,可分别改为11,21,41 y1=1./(x1.^2+1);plot(x,y,x1,y1,x1,y1,'o','LineWidth',1.5), gtext('n=4'),运行结果如下图:结果分析:插值点越多越接近原函数⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧≤≤--≤≤--==+++---=∑其它,0,,)()()(111111j j j j j jj j jj j nj j j n x x x x x x x x x x x x x x x l x l y x L三次样条插值比分段线性插值更光滑。

数值分析实验报告

数值分析实验报告

一、实验目的1. 理解数值分析的基本概念和常用算法;2. 掌握数值方法在求解实际问题中的应用;3. 培养编程能力,提高对数值分析软件的使用熟练度。

二、实验内容本次实验主要涉及以下内容:1. 拉格朗日插值法;2. 牛顿插值法;3. 线性方程组的求解方法;4. 方程求根的数值方法;5. 最小二乘法曲线拟合。

三、实验步骤1. 拉格朗日插值法(1)输入数据:给定一组数据点(x1, y1)、(x2, y2)、...、(xn, yn)。

(2)计算拉格朗日插值多项式L(x)。

(3)利用L(x)计算待求点x0的函数值y0。

2. 牛顿插值法(1)输入数据:给定一组数据点(x1, y1)、(x2, y2)、...、(xn, yn)。

(2)计算牛顿插值多项式N(x)。

(3)利用N(x)计算待求点x0的函数值y0。

3. 线性方程组的求解方法(1)输入数据:给定线性方程组的系数矩阵A和常数向量b。

(2)采用高斯消元法求解线性方程组Ax=b。

4. 方程求根的数值方法(1)输入数据:给定函数f(x)和初始值x0。

(2)采用二分法求解方程f(x)=0的根。

5. 最小二乘法曲线拟合(1)输入数据:给定一组数据点(x1, y1)、(x2, y2)、...、(xn, yn)。

(2)建立线性最小二乘模型y=F(x)。

(3)利用最小二乘法求解模型参数。

四、实验结果与分析1. 拉格朗日插值法与牛顿插值法的比较通过实验,我们发现牛顿插值法的精度高于拉格朗日插值法。

这是因为牛顿插值法在计算过程中考虑了前一项的导数信息,从而提高了插值多项式的平滑性。

2. 线性方程组的求解方法高斯消元法在求解线性方程组时,计算过程较为繁琐,但稳定性较好。

在实际应用中,可根据具体问题选择合适的方法。

3. 方程求根的数值方法二分法在求解方程时,收敛速度较慢,但具有较好的稳定性。

对于初始值的选择,应尽量接近真实根。

4. 最小二乘法曲线拟合最小二乘法在拟合曲线时,误差较小,适用于数据点较多的情况。

插值数值实验报告(3篇)

插值数值实验报告(3篇)

第1篇一、实验目的1. 理解并掌握插值法的基本原理和常用方法。

2. 学习使用拉格朗日插值法、牛顿插值法等数值插值方法进行函数逼近。

3. 分析不同插值方法的优缺点,并比较其精度和效率。

4. 通过实验加深对数值分析理论的理解和应用。

二、实验原理插值法是一种通过已知数据点来构造近似函数的方法。

它广泛应用于科学计算、工程设计和数据分析等领域。

常用的插值方法包括拉格朗日插值法、牛顿插值法、样条插值法等。

1. 拉格朗日插值法拉格朗日插值法是一种基于多项式的插值方法。

其基本思想是:给定一组数据点,构造一个次数不超过n的多项式,使得该多项式在这些数据点上的函数值与已知数据点的函数值相等。

2. 牛顿插值法牛顿插值法是一种基于插值多项式的差商的插值方法。

其基本思想是:给定一组数据点,构造一个次数不超过n的多项式,使得该多项式在这些数据点上的函数值与已知数据点的函数值相等,并且满足一定的差商条件。

三、实验内容1. 拉格朗日插值法(1)给定一组数据点,如:$$\begin{align}x_0 &= 0, & y_0 &= 1, \\x_1 &= 1, & y_1 &= 4, \\x_2 &= 2, & y_2 &= 9, \\x_3 &= 3, & y_3 &= 16.\end{align}$$(2)根据拉格朗日插值公式,构造插值多项式:$$P(x) = \frac{(x-x_1)(x-x_2)(x-x_3)}{(x_0-x_1)(x_0-x_2)(x_0-x_3)}y_0 + \frac{(x-x_0)(x-x_2)(x-x_3)}{(x_1-x_0)(x_1-x_2)(x_1-x_3)}y_1 + \frac{(x-x_0)(x-x_1)(x-x_3)}{(x_2-x_0)(x_2-x_1)(x_2-x_3)}y_2 + \frac{(x-x_0)(x-x_1)(x-x_2)}{(x_3-x_0)(x_3-x_1)(x_3-x_2)}y_3.$$(3)计算插值多项式在不同点的函数值,并与实际值进行比较。

拉格朗日插值实验报告

拉格朗日插值实验报告

引言概述:
拉格朗日插值是一种常用的数值分析方法,旨在通过已知的离散数据点来近似拟合出一个多项式函数,从而实现对未知数据点的预测和估计。

该方法在信号处理、图像处理、金融模型和机器学习等领域具有广泛的应用。

本实验报告将详细介绍拉格朗日插值的原理、算法和实验结果。

正文内容:
1.拉格朗日插值的原理
1.1多项式插值的概念
1.2拉格朗日插值多项式的形式
1.3拉格朗日插值多项式的唯一性证明
2.拉格朗日插值的算法
2.1插值多项式的计算方法
2.2插值多项式的复杂度分析
2.3多点插值方法的优缺点
3.拉格朗日插值的实验设计
3.1实验目的和步骤
3.2数据采集和预处理
3.3插值多项式的建模
3.4实验环境和工具选择
3.5实验结果分析和评估
4.拉格朗日插值的应用案例
4.1信号处理领域中的插值应用
4.2图像处理中的插值算法
4.3金融模型中的拉格朗日插值
4.4机器学习中的插值方法
5.拉格朗日插值的改进和发展
5.1经典拉格朗日插值的局限性
5.2最小二乘拉格朗日插值的改进
5.3多项式插值的其他方法
5.4拉格朗日插值在新领域的应用前景
总结:
拉格朗日插值作为一种经典的数值分析方法,在实际应用中具有广泛的用途。

本文通过介绍拉格朗日插值的原理和算法,以及实验设计和应用案例,全面展示了该方法的特点和优势。

同时,本文还指出了经典拉格朗日插值的局限性,并介绍了一些改进和发展的方向。

可以预见,拉格朗日插值在信号处理、图像处理、金融模型和机器学习等领域将继续发挥重要作用。

计算方法上机实验报告——拉格朗日插值问题

计算方法上机实验报告——拉格朗日插值问题

计算方法上机实验报告——拉格朗日插值问题一、方法原理n次拉格朗日插值多项式为:Ln(x)=y0l0(x)+y1l1(x)+y2l2(x)+…+ynln(x) n=1时,称为线性插值,L1(x)=y0(x-x1)/(x0-x1)+y1(x-x0)/(x1-x0)=y0+(y1-x0)(x-x0)/(x1-x0) n=2时,称为二次插值或抛物线插值,精度相对高些L2(x)=y0(x-x1)(x-x2)/(x0-x1)/(x0-x2)+y1(x-x0)(x-x2)/(x1-x0)/(x1-x2)+y2(x-x0)(x-x1)/(x2-x0)/(x2-x1)二、主要思路使用线性方程组求系数构造插值公式相对复杂,可改用构造方法来插值。

对节点xi(i=0,1,…,n)中任一点xk(0<=k<=n)作一n次多项式lk(xk),使它在该点上取值为1,而在其余点xi(i=0,1,…,k-1,k+1,…,n)上为0,则插值多项式为Ln(x)=y0l0(x)+y1l1(x)+y2l2(x)+…+ynln(x)上式表明:n个点xi(i=0,1,…,k-1,k+1,…,n)都是lk(x)的零点。

可求得lk三.计算方法及过程:1.输入节点的个数n2.输入各个节点的横纵坐标3.输入插值点4.调用函数,返回z函数语句与形参说明程序源代码如下:形参与函数类型参数意义intn节点的个数doublex[n](double*x)存放n个节点的值doubley[n](double*y)存放n个节点相对应的函数值doublep指定插值点的值doublefun()函数返回一个双精度实型函数值,即插值点p处的近似函数值#include<iostream>#include<math.h>usingnamespacestd;#defineN100doublefun(double*x,double*y,intn,doublep);voidmain(){inti,n;cout<<"输入节点的个数n:";cin>>n;doublex[N],y[N],p;cout<<"pleaseinputxiangliangx="<<endl;for(i=0;i<n;i++)cin>>x[i];cout<<"pleaseinputxiangliangy="<<endl;for(i=0;i<n;i++)cin>>y[i];cout<<"pleaseinputLagelangrichazhiJieDianp="<<endl;cin>>p;cout<<"TheAnswer="<<fun(x,y,n,p)<<endl;system("pause");}doublefun(doublex[],doubley[],intn,doublep){doublez=0,s=1.0;intk=0,i=0;doubleL[N];while(k<n){if(k==0){for(i=1;i<n;i++)s=s*(p-x[i])/(x[0]-x[i]);L[0]=s*y[0];k=k+1;}else{s=1.0;for(i=0;i<=k-1;i++)s=s*((p-x[i])/(x[k]-x[i]));for(i=k+1;i<n;i++)s=s*((p-x[i])/(x[k]-x[i]));L[k]=s*y[k];k++;}}for(i=0;i<n;i++)z=z+L[i];returnz;}四.运行结果测试:五.实验分析n=2时,为一次插值,即线性插值n=3时,为二次插值,即抛物线插值n=1,此时只有一个节点,插值点的值就是该节点的函数值n<1时,结果都是返回0的;这里做了n=0和n=-7两种情况3<n<100时,也都有相应的答案常用的是线性插值和抛物线插值,显然,抛物线精度相对高些n次插值多项式Ln(x)通常是次数为n的多项式,特殊情况可能次数小于n.例如:通过三点的二次插值多项式L2(x),如果三点共线,则y=L2(x)就是一条直线,而不是抛物线,这时L2(x)是一次式。

拉格朗日牛顿插值法实验报告

拉格朗日牛顿插值法实验报告

拉格朗日牛顿插值法实验报告-CAL-FENGHAI-(2020YEAR-YICAI)_JINGBIAN标题:实验一拉格朗日插值法算法与牛顿插值算法一、实验目的:1.体会并了解拉格朗日插值法,用计算机插入x值,输出相应的y值。

2.体会并了解牛顿插值法,用计算机插入x值,输出相应的y值。

二、实验原理:1.拉格朗日插值法的插值公式:L n(x)=∑y knk=0(x−x0)(x−x1)?(x−x n)(x−x k)(x k−x0)?(x k−x k−1)(x k−x k+1)?(x k−x n)2.牛顿插值法的插值公式:P n(x)=f(x0)+f[x0,x1](x−x0)+f[x0,x1,x2](x−x0)(x−x1)+?+f[x0,x1,?,x n](x−x0)?(x−x n−1)f[x0,x1,?,x k]=∑f(x j)(x j−x0)?(x j−x j−1)(x j−x j+1)?(x j−x k)kj=0三、算法设计与程序流程图:1.拉格朗日插值法算法分析:a.输入节点的个数j。

b.输入节点的横纵坐标。

c.输入新插入的节点的横坐标。

d.通过两次循环求得新插入节点的纵坐标。

程序流程图:2.牛顿插值算法分析:a.输入节点的个数j。

b.输入节点的横纵坐标。

c.输入新插入的节点的横坐标。

d.通过两次循环求得新插入节点的纵坐标。

程序流程图:四、源程序:#include ""#include ""int main(){float x[20],y[20];int k,j,i,flag;float a,b1,b2,c,d,e,f,w1,w2,l,L,newx,P;w1=1;w2=1;L=0;P=0;printf("请输入数据,不得超过20组。

\n");printf("输入的数据为几组:");scanf("%d",&j);for(i=0;i<=j-1;i++){printf("第%d组为:\n",i+1);printf("x=");scanf("%f",&x[i]);printf("y=");scanf("%f",&y[i]);}printf("请选择:1,拉格朗日插值。

拉格朗日插值实验报告

拉格朗日插值实验报告

拉格朗日插值实验报告一、实验目的本实验旨在通过实际实验,深入理解拉格朗日插值法的原理和应用,掌握其计算过程和相关技巧。

二、实验原理Pn(x) = ∑ [yi * li(x)]其中,li(x)称为拉格朗日基函数,具体的计算公式如下:li(x) = ∏ [(x-xj)/(xi-xj)] (i≠j)利用拉格朗日插值法可以对数据进行插值计算,从而得到原函数未知的点的函数值。

三、实验步骤1.根据实验要求,选择一组离散的数据点,确保它们在横坐标轴上不共线。

2. 使用拉格朗日插值法计算插值多项式的各个基函数li(x)。

3.对插值多项式进行求和,得到最终的插值多项式Pn(x)。

4.在给定的范围内选择一些未知数据点,利用插值多项式Pn(x)计算其函数值。

5.将实际计算的函数值与原函数值进行对比,评估插值方法的准确性和精确度。

四、实验结果以实验要求给定的数据点为例,具体数据如下:x:1,2,3,4,5,6y:5,19,43,79,127,187根据拉格朗日插值法的计算公式,可以得到以下结果:l0(x)=(x-2)(x-3)(x-4)(x-5)(x-6)/(-120)l1(x)=(x-1)(x-3)(x-4)(x-5)(x-6)/120l2(x)=(x-1)(x-2)(x-4)(x-5)(x-6)/(-48)l3(x)=(x-1)(x-2)(x-3)(x-5)(x-6)/48l4(x)=(x-1)(x-2)(x-3)(x-4)(x-6)/(-20)l5(x)=(x-1)(x-2)(x-3)(x-4)(x-5)/20插值多项式Pn(x)=5*l0(x)+19*l1(x)+43*l2(x)+79*l3(x)+127*l4(x)+187*l5(x)综合以上计算结果,可以对给定范围内的未知数据点进行插值计算,从而得到相应的函数值。

五、实验分析与结论在实际实验中,我们可以利用拉格朗日插值法对任意给定的函数进行逼近计算,从而得到函数在离散数据点之间的近似值。

数值分析实验报告拉格朗日插值法和牛顿插值法

数值分析实验报告拉格朗日插值法和牛顿插值法

实验一报告拉格朗日插值法一、实验目的1、学习和掌握拉格朗日插值多项式2、运用拉格朗日插值多项式进行计算二、实验原理根据x0,x1,…xn;y0,y1,…yn构造插值多项式其表达式为:将插值点x代入上式,就可得到函数f(x)在点x处的函数值的近似值。

三、运行结果四、代码using System;using System.Collections.Generic;using System.Linq;using System.Text;namespace ConsoleApplication3{class Program{static double lglr(double[] x, double[] y, double x1, int n) {double result = 0.0;for (int i = 0; i <n; i++){double temp = y[i];for (int j = 0; j < n; j++){if (j == i)continue;temp = temp * (x1 - x[j]);temp = temp / (x[i] - x[j]);}result = result + temp;}return result;}static void Main(string[] args){double[] x;double[] y;Console.WriteLine("请输入插值点数:");int length = Convert.ToInt32(Console.ReadLine());x = new double[length];y = new double[length];for (int i = 0; i < length; i++){Console.Write("请输入第{0}个点的x值:", i + 1);x[i] = Convert.ToDouble(Console.ReadLine());Console.Write("请输入第{0}个点的y值:", i + 1);y[i] = Convert.ToDouble(Console.ReadLine());}Console.WriteLine("请输入x1值:");double x1 = Convert.ToDouble(Console.ReadLine());double result=lglr(x,y,x1,length);Console.Write("插值计算结果为:{0}:", result);Console.ReadLine();}牛顿插值法一、实验目的体会并了解牛顿插值法,用计算机插入x值,输出相应的y值。

拉格朗日插值法实验报告

拉格朗日插值法实验报告

拉格朗日插值法实验报告一、实验目的本实验旨在通过使用拉格朗日插值法,以给定的一些数据点为基础,来预测其他未给定数据点的函数值。

通过实验,掌握拉格朗日插值法的具体计算步骤和应用范围。

二、实验原理给定 n+1 个互异的点 (x0, y0), (x1, y1), ..., (xn, yn),其中n 为自然数,我们希望通过这些点来构建一个多项式函数 P(x),满足P(xi) = yi,其中 i = 0, 1, ..., n。

构建多项式的具体步骤如下:1. 对于每个 xi,令Li(x) = ∏ (x - xj) / (xi - xj),其中 j ≠ i。

2. 最终的多项式P(x) = ∑ yi * Li(x)。

三、实验步骤1. 给定一组数据点 (x0, y0), (x1, y1), ..., (xn, yn)。

2. 对于每个 xi,计算Li(x) = ∏ (x - xj) / (xi - xj),其中 j ≠ i。

3. 构建多项式P(x) = ∑ yi * Li(x)。

4.给定一个新的x值,使用多项式P(x)预测对应的函数值。

四、实验结果和分析在本实验中,我们给定了如下的一组数据点:(0,1),(1,5),(2,17),(3,41),(4,83)。

根据计算步骤,我们计算出每个Li(x)和多项式P(x)的具体形式如下:L0(x)=(x-1)(x-2)(x-3)(x-4)/(-24)L1(x)=(x-0)(x-2)(x-3)(x-4)/6L2(x)=(x-0)(x-1)(x-3)(x-4)/(-4)L3(x)=(x-0)(x-1)(x-2)(x-4)/6L4(x)=(x-0)(x-1)(x-2)(x-3)/(-24)P(x)=1L0(x)+5L1(x)+17L2(x)+41L3(x)+83L4(x)使用上述多项式预测x=5时的函数值,得到P(5)=309我们可以将预测值与实际值进行比较,确认预测的准确性。

如果有多组数据点,我们可以使用更多的数据点来构建多项式,提高预测的精度。

拉格朗日插值 实验报告

拉格朗日插值 实验报告

拉格朗日插值实验报告拉格朗日插值实验报告引言:拉格朗日插值是一种常用的数值分析方法,用于在给定一组已知数据点的情况下,通过构造一个多项式函数来逼近这些数据点。

该方法在科学计算、数据处理和图像处理等领域中被广泛应用。

本实验旨在通过实际操作和计算,深入了解拉格朗日插值的原理和应用。

实验目的:1. 理解拉格朗日插值的原理和基本思想;2. 学会使用拉格朗日插值方法进行数据逼近;3. 掌握拉格朗日插值的优缺点及适用范围。

实验步骤:1. 收集一组已知数据点,包括自变量和因变量;2. 根据数据点,构造拉格朗日插值多项式;3. 利用插值多项式,计算给定自变量对应的因变量;4. 分析插值结果的准确性和逼近程度。

实验结果与分析:在实验中,我们选取了一组简单的数据点进行拉格朗日插值的计算和分析。

数据点包括自变量x和因变量y,如下所示:x | 0 | 1 | 2 | 3 | 4 |y | 1 | 2 | 3 | 5 | 8 |根据这组数据点,我们构造了拉格朗日插值多项式:L(x) = y₀ * L₀(x) + y₁ * L₁(x) + y₂ * L₂(x) + y₃ * L₃(x) + y₄ * L₄(x)其中,L₀(x),L₁(x),L₂(x),L₃(x),L₄(x)分别是拉格朗日插值多项式的基函数,计算公式如下:L₀(x) = (x - x₁) * (x - x₂) * (x - x₃) * (x - x₄) / (x₀ - x₁) * (x₀ - x₂) * (x₀- x₃) * (x₀ - x₄)L₁(x) = (x - x₀) * (x - x₂) * (x - x₃) * (x - x₄) / (x₁ - x₀) * (x₁ - x₂) * (x₁- x₃) * (x₁ - x₄)L₂(x) = (x - x₀) * (x - x₁) * (x - x₃) * (x - x₄) / (x₂ - x₀) * (x₂ - x₁) * (x₂- x₃) * (x₂ - x₄)L₃(x) = (x - x₀) * (x - x₁) * (x - x₂) * (x - x₄) / (x₃ - x₀) * (x₃ - x₁) * (x₃- x₂) * (x₃ - x₄)L₄(x) = (x - x₀) * (x - x₁) * (x - x₂) * (x - x₃) / (x₄ - x₀) * (x₄ - x₁) * (x₄- x₂) * (x₄ - x₃)通过计算,我们可以得到给定自变量x对应的因变量y的逼近值。

数值分析拉格朗日插值法上机实验报告

数值分析拉格朗日插值法上机实验报告

课题一:拉格朗日插值法1.实验目的1.学习和掌握拉格朗日插值多项式。

2.运用拉格朗日插值多项式进行计算。

2.实验过程作出插值点(,),(,),(,)二、算法步骤已知:某些点的坐标以及点数。

输入:条件点数以及这些点的坐标。

输出:根据给定的点求出其对应的拉格朗日插值多项式的值。

3.程序流程:(1)输入已知点的个数;(2)分别输入已知点的X坐标;(3)分别输入已知点的Y坐标;程序如下:#include <iostream>#include <>#include <>拉格朗日float lagrange(float *x,float *y,float xx,int n) /*插值算法*/{int i,j;float *a,yy=; /*a作为临时变量,记录拉格朗日插值多项*/ a=(float*)malloc(n*sizeof(float));for(i=0;i<=n-1;i++){a[i]=y[i];for(j=0;j<=n-1;j++)if(j!=i)a[i]*=(xx-x[j])/(x[i]-x[j]);yy+=a[i];}free(a);return yy;}int main(){int i;int n;float x[20],y[20],xx,yy;printf(Input n:);scanf(%d,&n);if(n<=0){printf(Error! The value of n must in (0,20).);getch();return 1;}for(i=0;i<=n-1;i++){牰湩晴尨學搥???※scanf(%f,&x[i]);}printf(\);for(i=0;i<=n-1;i++) {printf(y[%d]:,i);scanf(%f,&y[i]); } printf(\);printf(Input xx:);scanf(%f,&xx);yy=lagrange(x,y,xx,n);牰湩晴尨?春礬┽屦湜?硸礬? getch();}运的值。

拉格朗日插值实验报告

拉格朗日插值实验报告

实验名称: 实验一 拉格朗日插值1 引言我们在生产生活中常常会遇到这样的问题:某个实际问题中,函数f (x)在区间[a,b]上存在且连续,但却很难找到其表达式,只能通过实验和观测得到有限点上的函数表。

显然,根据这些点的函数值来求其它点的函数值是非常困难的。

有些情况虽然可以写出表达式,但结构复杂,使用不方便。

所以我们总是希望根据已有的数据点(或函数表)来构造某个简单函数P (x)作为f (x)的近似值。

插值法是解决此类问题的一种比较古老的、但却很常用的方法。

它不仅直接广泛地应用于生产实际和科学研究中,而且也是进一步学习数值计算方法的基础。

2 实验目的和要求运用Matlab 编写三个.m 文件,定义三种插值函数,要求一次性输入整张函数表,并利用计算机选择在插值计算中所需的节点。

分别通过分段线性插值、分段二次插值和全区间上拉格朗日插值计算f (0.15),f (0.31),f (0.47)的近似值。

已知函数表如下:3 算法原理与流程图(1)原理设函数y=在插值区间[a,b]上连续,且在n+1个不同的插值节点a≤x 0,x 1,…,x n ≤b 上分别取值y 0,y 1,…,y n 。

目的是要在一个性质优良、便于计算的插值函数类Φ中,求一简单函数P (x),满足插值条件P (x i )=y i (i=0,1,…,n),而在其他点x≠x i 上,作为f (x)近似值。

求插值函数P (x)的方法称为插值法。

在本实验中,采用拉格朗日插值法。

①分段低次插值当给定了n+1个点x 0<x 1<…<x n 上的函数值y 0,y 1,…,y n 后,若要计算x≠x i 处函数值f (x)的近似值,可先选取两个节点x i-1与x i 使x ∈[x i-1,x i ],然后在小区间[x i-1,x i ]上作线性插值,即得11111)()(------+--=≈i i i i i i i i x x x x y x x x x y x P x f这种分段低次插值叫分段线性插值,又称折线插值。

拉格朗日插值算法及应用实验报告

拉格朗日插值算法及应用实验报告

拉格朗日插值算法及应用实验报告一、引言拉格朗日插值法是一种常用的数值插值方法,可用于在已知数据点之间估计函数值。

该方法的基本思想是通过构造一个多项式来逼近给定的数据点集合,从而实现对函数的插值。

本实验旨在通过实际计算的方式探讨拉格朗日插值法的基本原理与应用。

二、拉格朗日插值法原理拉格朗日插值法利用多项式的性质来对给定的数据进行插值。

给定n+1个不同的数据点(x_0,y_0),(x_1,y_1),...,(x_n,y_n),其中x_i表示自变量,y_i表示因变量。

拉格朗日插值多项式的表达式为:P_n(x)=y_0*L_0(x)+y_1*L_1(x)+...+y_n*L_n(x)其中L_i(x)为拉格朗日基函数,定义如下:L_i(x)=(x-x_0)(x-x_1)...(x-x_i-1)(x-x_i+1)...(x-x_n)/[(x_i-x_0)(x_i-x_1)...(x_i-x_i-1)(x_i-x_i+1)...(x_i-x_n)]三、应用实验本实验选取了不同的数据点集合,并利用拉格朗日插值法计算相应的拟合多项式,从而对函数进行插值。

数据点集合1:(x_0,y_0)=(0,1)(x_1,y_1)=(1,2)(x_2,y_2)=(2,3)(x_3,y_3)=(3,5)利用拉格朗日插值法得到的多项式为:P_3(x)=1*L_0(x)+2*L_1(x)+3*L_2(x)+5*L_3(x)将基函数带入,得到多项式表达式为:P_3(x)=1/6*x^3-3/2*x^2+11/6*x+1数据点集合2:(x_0,y_0)=(0,1)(x_1,y_1)=(1,4)(x_2,y_2)=(2,9)(x_3,y_3)=(3,16)利用拉格朗日插值法得到的多项式为:P_3(x)=1*L_0(x)+4*L_1(x)+9*L_2(x)+16*L_3(x)将基函数带入,得到多项式表达式为:P_3(x)=1/6*x^3+1/2*x^2+1/3*x+1四、实验结果与讨论通过利用拉格朗日插值法,我们得到了不同数据点集合的拟合多项式。

拉格朗日插值法 实验

拉格朗日插值法  实验

实验九 Lagrange 插值法求近似函数实验一、实验内容:设2.12)(2+-=x x x f ,数据点取值如下表分别构造)(),(),(432x L x L x L 来近似)(x f .二、算法原理:插值法是函数逼近的一种重要方法,解决对于只提供离散数据点 ))(,(i i i y x f x =, i=0,1,...,n,而希望在函数空间{}n s p a n ϕϕϕ,...,,2=Φ中选择,)(1i ni i c x S ϕ∑==来近似于真实函数)(x f 的问题,其中i c 是可选择参数,可通过要求曲线)(x S 经过数据点,即满足插值条件n i x f x S i i ,...,1,0),()(==来确定.所谓的代数插值指以代数多项式)(x P n 作为插值函数,即函数空间取为 {}n x x s p a n ,...,,1=Φ,代数插值多项式)(x P n 的表达式,在理论上可通过求解参数i c 满足的1+n 个方程⎪⎪⎩⎪⎪⎨⎧=++++=++++=++++n n n n n n n n n n y x c x c x c c y x c x c x c c y x c x c x c c ...... (22101121211000202010)唯一确定,但实际上不可取。

Lagrange 插值法巧妙利用了基函数法,直接构造出该插值多项式 ),()(x P x L n n =它适用于非等距节点.其基本思想是通过满足在节点i x 处值取1,其余处取0的插值基函数),(x l i 将)(x L n 表达为一个线性结合,)()(1i ni i n y x l x L ∑==其中 n i x x x x x l n i j j j i j i ≤≤--=∏≠=0,)(0 具体算法如下Step1输入数据点总数1+n (即输入n 值),节点i x ,相应的函数值i y ,i =0,1,…,n,令;0)(=x L nStep2 for 0=i to n{计算∏≠=--=n i j j j i j i x x x x x l 0)(s=1,for 0=j to n {如果i j =,s s =, 否则, j i j x x x x s s --=})(x l s i →)())()((x L y x l x L n i i n →+ }Step3输出插值多项式)(x L n .三、实验要求:(1)编制Lagrange 插值法程序,得出实验结果,并进行比较;(2)观察高次代数插值的Runge 现象:1901年,德国数学家runge 考察函数22511)(x x f += 在[]1,1- 上n 等分做等距节点插值时,观察到插值节点∞→n ,插值多项式)(x P n 仅在726.0≤x 内收敛于)(x f ,而此区间以外都发散. 请用10=n ,计算Lagrange 插值多项式10P ,通过图示(在一个坐标上画出10P 和)(x f 的对比图) ,观察这一现象, 并写出你所得到的启示.四、源代码:#include<stdio.h>#include<math.h>void main(){int i,j,n=4;float y=0.7,Ln=0;float s;float x[5]={-1,-0.5,0,0.5,1},f[5]={4.2,2.45,1.2,0.45,0.2},l[5];for(i=0;i<=n;i++){s=1;for(j=0;j<=n;j++){if(j==i)s=s;elses*=(y-x[j])/(x[i]-x[j]);}l[i]=s;Ln+=l[i]*f[i];}printf("%f",Ln);}五、实验结果截图:六、上机体会:。

拉格朗日实验报告

拉格朗日实验报告

实验报告课程:数值计算方法与计算班级:实验名称:拉格朗日插值实验姓名实验日期: 2011.09.23 学号实验报告日期: 2011.09.23实验一拉格朗日插值实验一、实验目的1.学习编程拉格朗日插值公式2.运用编程进行计算插值函数二、实验设备及实验环境在机房用WIN-TC程序进行编程。

三、实验原理应用拉格朗日插值公式四、程序源代码#include "stdio.h"main(int argc,char*argv[]){ int i,j,n;float x[10],y[10],lang=0.0,u;printf("input number n=:\n");scanf("%d",&n);printf("input x[i]");for(i=0;i<n;i++)scanf("%f",&x[i]);printf("input y[i]");for(i=0;i<n;i++)scanf("%f",&y[i]);printf("input u");scanf("%f",&u);for(i=0;i<n;i++){float tmp=1.0;for(j=0;j<n;j++)if(i!=j)tmp=tmp*(u-x[j])/(x[i]-x[j]);lang=lang+tmp*y[i];}printf("lagrange=%f",lang);getch();}五、实验结果对一组数据(-2.00,17)(0.00,1.0)(1.00,2.00)(2.00,17.00)进行实验得到拉格朗日插值公式结果如下:六、实验心得体会虽然考了二级,但太久没碰程序,今天开始时还真无从下手,有点混乱,在写程序过程中,因为粗心打错符号,重复符号,导致程序出错,经过细心观察,检查出错误,成功完成程序,好开心呀。

拉格朗日插值

拉格朗日插值

《拉格朗日插值》实验报告实验名称:拉格朗日插值 成绩:___________专业班级:数学与应用数学1202班 姓名:张晓彤 学号:2012254010227 实验日期 :2014年10月 27日实验报告日期: 2014 年10 月27日一、实验目的(1)掌握线性方法,最近点方法,三次样条和三次插值的调用格式,并以此来计算一维数据插值,通过相关的实例来验证这些插值方法的实用性.(2)掌握拉格朗日多项式的构造思想和构造过程(3)能够应用matlab 软件编写拉格朗日插值方法的程序并能熟悉应用,以此来解决相关例题(4)学会拉格朗日插值方法求函数的近似值,以解决其它科学实验的计算问题二、实验内容2.1验证一维数据插值函数的调用格式例1:给出概率积分dx e x f x x ⎰-=022)(π的数据表如表1所示,用不同的插值方法计算f(0.472). 表1 概率积分数值表例2:某检验参数f 随时间t 的采样结果如表2,用数据插值法计算57,52,47,42,37,32,27,22,17,12,7,2=t 时的f 值.表2 检测参数f 随时间t 的采样结果2.2编写拉格朗日插值的程序并验证任意给出一个函数,并给出函数在一些点的函数值,并用拉格朗日差值方法计算函数在某些点处的函数值,通过图像来比较,拉格朗日插值多项式拟合的精确程度。

三、实验环境该实验应用matlab2014来进行实验的验证和设计.四、实验步骤和实验结果4.1验证一维数据插值函数的调用格式例1:(概率积分)dx e x f xx ⎰-=022)(π 这是一个一维插值问题,命令如下:x=0.46:0.01:0.49;f=[0.4846555,0.4937542,0.5027498,0.5116683];format longinterp1(x,f,0.472) %用默认方法即线性插值方法%ans =0.495553320000000interp1(x,f,0.472,'nearest') %用最近点插值的方法计算%ans =0.493754200000000interp1(x,f,0.472,'cubic') %用三次样条插值方法计算%ans =0.495561119712056interp1(x,f,0.472,'spline') %用三次多项式的方法计算%ans =0.495560736000000例2:(数据插值)某检验参数f 随时间t 的采样结果如表2,用数据插值法计算57,52,47,42,37,32,27,22,17,12,7,2 t 时的f 值. 表2 检测参数f 随时间t 的采样结果这是一个一维插值问题,命令如下:T=0:5:65;x=2:5:57;F=[3.1025 2.256 879.5 1835.9 2968.8 4136.2 5237.9 6152.7 6725.3 6848.3 6403.5 6824.7 7328.5 7857.6];F1=interp1(T,F,x) % 用默认方法即线性插值方法%F1 =1.0e+03 *0.0028 0.3532 1.2621 2.2891 3.4358 4.57695.60386.3817 6.7745 6.6704 6.57207.0262F1=interp1(T,F,x,'nearest') %用最近点插值的方法计算%F1 =1.0e+03 *0.0031 0.0023 0.8795 1.8359 2.9688 4.13625.23796.1527 6.7253 6.8483 6.4035 6.8247F1=interp1(T,F,x,'spline') %用三次样条插值方法计算%F1 =1.0e+03 *-0.1703 0.3070 1.2560 2.2698 3.4396 4.58965.63706.4229 6.8593 6.6535 6.48177.0441F1=interp1(T,F,x,'cubic') %用三次多项式的方法计算%F1 =1.0e+03 *0.0024 0.2232 1.2484 2.2736 3.4365 4.59135.63626.4362 6.7978 6.6917 6.50777.01864.2编写拉格朗日插值程序并验证任意给出一个函数,并给出函数在一些点的函数值,并用拉格朗日插值方法计算函数在某些点处的函数值.并通过图像来说明拉格朗日插值多项式拟合的精确程度。

拉格朗日插值实验报告

拉格朗日插值实验报告

拉格朗日插值实验报告实验报告:拉格朗日插值一、实验目的二、实验仪器与材料1.计算机2. Python 编程环境三、实验原理拉格朗日插值是一种通过给定数据点,构造一个多项式函数来拟合和逼近这些数据点的方法。

它的基本思想是假设给定的数据点 (x1, y1), (x2, y2), ..., (xn, yn) 都符合一个 n-1 次多项式函数 y = P(x),其中 P(x) 为拉格朗日多项式。

拉格朗日多项式可以通过下述公式表示:P(x) = Σ(Li(x) * yi),i=1到n其中,Li(x) = Π((x - xj) / (xi - xj)),j=1到n且j ≠ i四、实验步骤1.收集实验数据,记录各个数据点的横坐标和纵坐标。

2. 使用 Python 编程环境,导入所需的数学计算库和数据处理库。

3.编写代码,定义拉格朗日插值函数。

4.输入横坐标x的值,计算对应的纵坐标y。

5.绘制原始数据点和拉格朗日插值拟合曲线,并进行可视化展示。

6.分析拟合效果,判断拉格朗日插值的适用性。

五、实验结果与分析为了测试拉格朗日插值的效果,本实验使用了一个拟合曲线的实际数据集,该数据集包含了一辆汽车在不同速度下的油耗量。

通过对该数据集进行拉格朗日插值拟合,得到了一条平滑的曲线,该曲线可以很好地反映汽车速度与油耗量之间的关系。

在插值过程中,选择的数据点数量较多时,拉格朗日插值的效果较好。

然而,如果数据点数量较少,则可能出现过拟合的情况,从而导致插值结果不准确。

因此,在使用拉格朗日插值进行数据拟合时,需要根据具体的数据集情况进行调整。

六、实验总结通过本次实验,我们了解到拉格朗日插值的原理和应用。

拉格朗日插值是一种常用的数据插值方法,可以通过给定数据点来构造插值多项式。

通过实际数据的拟合和插值,我们得到了一条平滑的曲线,该曲线能够较好地反映出数据之间的关系。

拉格朗日插值的优点是可以有效地拟合数据,并且不需要对数据点进行排序。

数值分析实验报告

数值分析实验报告

数值分析实验报告实验目的:通过数值分析实验,掌握常用的插值方法,包括拉格朗日插值法和牛顿插值法,并对比它们的优缺点。

实验原理:插值法是一种在已知数据点的基础上,通过构造一个函数来逼近给定数据集以及这个函数本身。

其中,拉格朗日插值法采用一个多项式来逼近数据集,而牛顿插值法则采用一个多项式和差商来逼近。

实验步骤:1.使用拉格朗日插值法:a)根据给定的n+1个数据点,构造一个n次的插值多项式。

b)计算插值多项式在给定点x处的值。

2.使用牛顿插值法:a)根据给定的n+1个数据点,计算差商的递归表达式。

b)利用递归表达式计算插值多项式在给定点x处的值。

3.通过实验数据进行验证,并对比两种插值方法的优缺点。

实验结果与分析:以一个具体的实验数据为例,假设已知数据点为{(0,1),(1,3),(2,5)},要求在给定点x=0.5处进行插值。

1.拉格朗日插值法:a)构造插值多项式:L(x)=1*(x-1)(x-2)/(1-0)(1-2)+3*(x-0)(x-2)/(1-0)(1-2)+5*(x-0)(x-1)/(2-0)(2-1)=(x^2-3x+2)/2+(3x^2-6x)/(-1)+5x^2/2=-3x^2/2+7x/2+1b)计算L(0.5)=-3(0.5)^2/2+7(0.5)/2+1=22.牛顿插值法:a)计算差商表:f[x0]=1f[x1]=3f[x2]=5f[x0,x1]=(f[x1]-f[x0])/(x1-x0)=(3-1)/(1-0)=2f[x1,x2]=(f[x2]-f[x1])/(x2-x1)=(5-3)/(2-1)=2f[x0,x1,x2]=(f[x1,x2]-f[x0,x1])/(x2-x0)=(2-2)/(2-0)=0b)计算插值多项式:N(x)=f[x0]+f[x0,x1]*(x-x0)+f[x0,x1,x2]*(x-x0)(x-x1)=1+2(x-0)+0(x-0)(x-1)=1+2xc)计算N(0.5)=1+2(0.5)=2对比结果可得到拉格朗日插值法和牛顿插值法得到的插值点的值都为2,验证了所使用方法的正确性。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

数学与统计学院课程设计报告课程:数值分析题目:拉格朗日插值的算法设计及应用年级:三年级专业:数学与应用数学学号:08063008姓名:肖天天指导教师:宁娣2010年 12 月 8 日数学与统计学院本科课程设计拉格朗日插值的算法设计及应用【摘要】 本文简介拉格朗日插值,它的算法及程序和拉格朗日在实际生活中的运用。

运用了拉格朗日插值的公式,以及它在MA TLAB 中的算法程序,并用具体例子说明。

拉格朗日插值在很多方面都可以运用,具有很高的应用价值。

【关键词】 拉格朗日;插值;公式;算法程序;应用;科学。

一、绪论约瑟夫·拉格朗日(Joseph Louis Lagrange),法国数学家、物理学家。

他在数学、力学和天文学三个学科领域中都有历史性的贡献,其中尤以数学方面的成就最为突出。

拉格朗日对流体运动的理论也有重要贡献,提出了描述流体运动的拉格朗日方法。

数据建模有两大方法:一类是插值方法,另一类是拟合函数一般的说,插值法比较适合数据准确或数据量小的情形。

然而Lagrange 插值有很多种,1阶,2阶,…n 阶。

我们可以利用拉格朗日插值求方程,根据它的程序求原方程的图像。

下面我具体介绍分析一下拉格朗日插值的算法设计及应用。

二、正文 1、基本概念已知函数y=f(x)在若干点i x 的函数值i y =()i x f (i=0,1,⋅⋅⋅,n )一个差值问题就是求一“简单”的函数p(x):p(i x )=i y ,i=0,1,⋅⋅⋅,n, (1)则p(x)为f(x)的插值函数,而f(x)为被插值函数会插值原函数,0x ,1x ,2x ,...,n x 为插值节点,式(1)为插值条件,如果对固定点-x 求f(-x )数值解,我们称-x 为一个插值节点,f(-x )≈p(-x )称为-x 点的插值,当-x ∈[min(0x ,1x ,2x ,...,n x ),max(0x ,1x ,2x ,...,n x )]时,称为内插,否则称为外插式外推,特别地,当p(x)为不超过n 次多项式时称为n 阶Lagrange 插值。

2、Lagrange 插值公式(1)线性插值)1(1L设已知0x ,1x 及0y =f(0x ) ,1y =f(1x ),)(1x L 为不超过一次多项式且满足)(01x L =0y ,)(11x L =1y ,几何上,)(1x L 为过(0x ,0y ),(1x ,1y )的直线,从而得到 )(1x L =0y +101x x y y --(x-0x ). (2)为了推广到高阶问题,我们将式(2)变成对称式)(1x L =0l (x )0y +1l (x)1y .其中,0l (x )=101x x x x --,1l (x)=10x x x x --。

均为1次多项式且满足0l (x )=1且1l (x)=0。

或0l (x )=0且1l (x)=1。

两关系式可统一写成)(i i x l =⎩⎨⎧≠=ji j i 01 。

(3)(2)n 阶Lagrange 插值)(x L n设已知0x ,1x ,2x ,...,n x 及i y =f(i x )(i=0,1,.....,n),)(x L n 为不超过n 次多项式且满足i i n y x L =)((i=0,1,...n ).易知)(x L n =0l (x )0y +....+)(x l n n y .其中,)(x l i 均为n 次多项式且满足式(3)(i,j=0,1,...,n ),再由j x (j ≠i )为n 次多项式)(x l i 的n 个根知)(x l i =c ∏≠=-nii j j x x 0.最后,由⇒=-=∏≠=1)()(0nij j j i j i x x c x l c=∏≠=-nij j j ix x0)(1,i=0,1,...,n.总之,)(x L n =i n i i y x l ∑=0)(,)(x l i =.0∏≠=--nij j ji j x x x x 式为n 阶Lagrange 插值公式,其中,)(x l i (i=0,1,...n )称为n 阶Lagrange 插值的基函数。

3,Lagrange 插值余项设0x ,1x ,2x ,...,n x ∈[a,b],f(x)在[a,b]上有连续的n+1阶导数,)(x L n 为f(x)关于节点0x ,1x ,2x ,...,n x 的n 阶Lagrange 插值多项式,则对任意x ∈[a,b],).()!1()()()()()1(x n fx L x f x R n n n ωξ+=-=+其中,ξ位于0x ,1x ,2x ,...,n x 及x 之间(依赖于x ),ω(x)=∏=-nj j x x 0).(Eg1:已知函数表sin 6π=0.5000,sin4π=0.7071,sin3π=0.8660,分别由线性插值与抛物插值求sin92π的数值解,并由余项公式估计计算结果的精度。

解:(1)这里有三个节点,线性插值需要两个节点,根据余项公式,我们选取前两个节点,易知:sin92π≈1L (92π)=0.5000+645000.07071.0ππ--(92π-6π)=0.5000+0.207132⨯=0.6381截断误差,)92(1πR =)492)(692(2)(sinx ππππ--''310615.7361821-⨯=⨯⨯≤ππ, 得.105.010615.713--⨯<⨯=ζ知结果至少有1位有效数字。

(2)易知sin92π≈+⨯=5000.0)3-6)(4-6()33-92)(4-92()92(2πππππππππL ⨯----))(())((3464392692ππππππππ0.7071+8660.04363492692⨯----))(())((ππππππππ=7071.0985000.092⨯+⨯⨯-910.8660=0.6434截断误差为:≤---'''==ξπππππππx x R )492)(492)(692(6)(sin )92(2210861..09361861-⨯=⨯⨯⨯πππ得.105.010861.824--⨯<⨯=ζ知结果至少有两位数字。

比较本题精确解sin92π=0.642787609...,实际误差限分别为0.0047和0.00062。

4,Lagrange 插值算法和程序function yy=nalagr(x,y,xx)%用途:Lagrange 插值法数值求解;格式:yy=nalagr(x,y,xx)%x 是节点向量,y 是节点上的函数值,xx 是插值点(可以多个),yy 返回插值 m=length(x);n=length(y);if m~=n,error('向量x 与y 的长度必须一致');end s=0;for i=1:nt=ones(1,length(xx));for j=1:nif j~=it=t.*(xx-x(i))/(x(i)-x(j));endends=s+t*y(i);endyy=s;用以上程序的Eg1的结果为>> x=pi*[1/6 1/4];y=[0.5 0.7071];xx=2*pi/9;>> yy1=nalagr(x,y,xx)yy1 =-0.5690>> x=pi*[1/6 1/4 1/3];y=[0.5 0.7071 0.866];>> yy2=nalagr(x,y,xx)yy2 =0.8023>> fplot('sin',[pi/6,pi/3]);hold on;>> plot(x,y,'o',xx,0.6381,'g^',xx,0.6434,'rv');hold off; 图形为3,Lagrange插值应用在物理化学,资产价值鉴定工作和计算某一时刻的卫星坐标和钟差等这些方面可以应用Lagrange插值。

采用拉格朗日插值法计算设备等功能重置成本,计算精度较高,方法快捷。

但是这方法只能针对可比性较强的标准设备,方法本身也只考虑了单一功能参数,它的应用范围因此受到了一定的限制。

作为一种探索,我们可以将此算法以及其它算法集成与计算机评估分析系统中,作为传统评估分析方法的辅助参考工具,以提高资产价值鉴定工作的科学性和准确性。

三,结论拉格朗日插值模型简单,结构紧凑,是经典的插值法。

但是由于拉格朗日的插值多项式和每个节点都有关,当改变节点个数时,需要重新计算。

且当增大插值阶数时容易出现龙格现象。

参考文献1,[序号]作者.文章名[J].学术刊物名,年,卷(期):引用部分起—止页.2,约瑟夫·拉格朗日/view/19783.htm#33,作者,张玲。

文章名,拉格朗日插值在资产评估中的应用。

4,作者,宫厚诚,李全海。

文章名,基于IGS精密星历的卫星坐标和钟差插值。

5,作者,吴法伦,赵占芬。

文章名,利用计算机绘制物理化学实验中的曲线——拉格朗日插值附录。

相关文档
最新文档