二元拉格朗日插值Fortran程序设计实验.docVIP

二元拉格朗日插值 实验目的-程序功能 利用FORTRAN编程实现二元拉格朗日插值求解函数在给定点的函数值。 设已知插值节点（xi,yi）(i=1,…,m,j=1,…,n)及对应函数值zij=f(xi,yi) (i=1,…,m,j=1,…,n),用拉格朗日插值法求函数在给定点（x,y）处的对应函数值z。 实验内容 了解和学习FORTRAN程序语言，会编写一些小程序； 学习和理解拉格朗日插值的原理及方法，并拓展至二元拉格朗日插值方法； 利用FORTRAN编程实现二元拉格朗日插值法； 举例进行求解，并对结果进行分析。 实验原理及方法 1、基本概念 已知函数y=f(x)在若干点的函数值=（i=0,1,,n）一个差值问题就是求一“简单”的函数p(x)： p()=,i=0,1,,n, (1) 则p(x)为f(x)的插值函数，而f(x)为被插值函数会插值原函数，，，，...,为插值节点，式（1）为插值条件，如果对固定点求f()数值解，我们称为一个插值节点，f()p()称为点的插值，当[min(，，，...,)，max(，，，...,)]时，称为内插，否则称为外插式外推，特别地，当p(x)为不超过n次多项式时称为n阶Lagrange插值。 2、Lagrange插值公式 2.1 线性插值 设已知 ，及=f() ,=f(),为不超过一次多项式且满