《计算方法上机答案》由会员分享，可在线阅读，更多相关《计算方法上机答案（2页珍藏版）》请在金锄头文库上搜索。
1、计算方法上机实习一实习内容一 拉格朗日插值多项式的构造1、方法简介：设y=f(x)是区间a,b上的连续函数，已知它在a,b上的n+1个互不相同的点 x0,X,.xn上取值y0,y,.yn希望构造一个函数L(x)作为fx)的近似表达式，且满 足：L(xi)=f(xi)=y. (1=0,1,n),这类问题即为插值问题。拉格朗日插值多项式的公式为：y = Ln (X)=工(已亠) yix .-x.ii=0j=0 i j i丰j2、使用说明输入参数：N 整变量，插值节点个数。X0 N 个元素的一维实数组，存放给定的插值节点。Y0 N 个元素的一维实数组，存放给定的插值节点上相应的函数值。X 实变量，插
2、值点。输出参数：Y 实变量，插值结果。3、实例：(1)据探测资料已知：各标准等压面(300hPa, 500hPa,700hPa, 850hPa, 1000hPa) 上各个网格点(100 X 100)的位势高度值(8280gpm, 5800gpm, 3140gpm, 1520gpm, 770gpm),用拉格朗日插值多项式求400hPa, 600hPa, 800hPa个层网 格上的位势高度值。Subroutine LAGR (X0,Y0,N,X,Y)Dimension X0(N),Y0(n)Y=0.0Do 30 I=1,nP=1,0Do 20 J=1,nIf(i.ne.j) p=p*(x-x0(j)/(x0(i)-x0(j)20 continue30y=y+p*y0(i)ReturnEnd4、流程图：5 实习要求及实习报告要求按以上过程完成实习内容，完成实习报告。实习报告包括 实习内容，要求，基本原理，运行结果，结果分析和总结。
