二元拉格朗日插值Fortran程序设计实验7.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(),为不超过一次多项式且满足=,=,几何上，为过（，），（，）的直线，从而得到 =+（x-）. （2） 为了推广到高阶问题，我们将式（2）变成对称式 =（x）+(x). （3） 其中， （x）=,(x)=。 均为1次多项式且满足 （）=1且()=0；（）=0且()=1。 两关系式可统一写成 2.2 n阶Lagrange插值 设已知，，，...,及=f()(i=0,1,.....,n),为不超过n次多项式且满足（i=0,1,...n）. 易知 其中，均为n次多项式且满足式（4）（i,j=0,1,...,n）,再由（ji）为n次多项式的n个根，知=c.最后，由 c=,i=0,1,...,n. 总之得到： （6）式为n阶Lagrange插值公式，其中，（i=0,1,…,n）称为n阶Lagrange插值的基函数。 3 二元拉格朗日插值方法 对于一元函数y=f(x)，得到n+1个数据点（,）(i=0,1,…,n),可由（6）、（7）式求得n阶Lagrange插值公式，然后求函数在y=f(x)在x点的函数值。 对于二元函数，若知道数据点（i=1,…,m，j=1,…,n）,可利用两次拉格朗日插值计算在点（x,y）的函数值，方法如下： （1）对每个( i=1,…,m),以（ j=1,…,n）为插值节点， （ j=1,…,n）为对应函数值，y为插值变量，作一元函数插值得( i=1,…,m); （2） 以( i=1,…,m)为插值节点，( i=1,…,m)为对应函数值，x为插值变量，作一元函数插值求得（x,y）点的值z。 FORTRAN编程 开发环境 使用Compaq Visual Fortran 6.6进行程序设计，编程实现二元拉格朗日插值算法。 使用说明 先编出一元拉格朗日差值算法子程序lagrange，然后编写二元拉格朗日插值算法程序lagrange2,其中两次调用lagrange子程序。 Lagrange(xa,ya,n,x,y) n 整型变量，输入参数，节点个数 xa n个元素的一维实数型数组，输入参数，存放自变量插值节点xi（i=1,…,n） ya n个元素的一维实数型数组，输入参数,存放函数值（y1,…,yn）T x 实型变量，输入参数，插值自变量 y 实型变量，输出参数，所求值 ****************************************************** Lagrange2(x1a, x2a,ya,m,n,x1, x2,y) m 整型变量，输入参数，x自变量节点个数 n 整型变量，输入参数，y自变量节点个数 x1a m个元素的一维实数型数组，输入参数，存放x自变量插值节点xi（i=1,…,m） x2a n个元素的一维实数型数组，输入参数，存放y自变量插值节点yj（j=1,…,n） x1 实型变量，输入参数，插值x自变量 x2 实型变量，输入参数，插值y自变量 ya m×n个元素的二维实数型数组，输入参数,存放（xi,yj）（i=1,…,m，j=1,…,n）函数值（y1,…,yn）T