莫言颁奖词简写50字:拉格朗日插值算法

来源:百度文库 编辑:中科新闻网 时间:2024/05/04 19:22:02

全区间拉格朗日插值

功用 
本程序用拉格朗日插值公式对一元不等距观测数据进行程组插值 。

方法概要
对给定的n个插值节点x1,x2,…,xn及对应的函数值y1,y2,…,yn,计算给定点x的函数值y(x)。

本程序可以在插值区间内对给定的NJ个插值点进行插值。

程序说明
(1)、程序名

LGRCZ1.FOR (此文件为免费软件,你可以从本站点 下载。)

(2)、子程序语句

SUBROUTINE LGRCZ1(Z,N,CZ,NJ)

(3)、形参说明

Z  输入参数,两维实数组,容量为N×2,存放给定的插值节点数列及对应的观测值。
N  输入参数,整变量,插值节点数。
NJ  输入参数,整变量,插值节点数。
CZ  输入输出参数,二维实数组,容量NJ×2,第一列输入插值点数列,第二列输出插值

结果。

计算实例
(1)、例题

设已知下表观测数据

x
1
2
3
4
5
6
7

Y1
1
4
9
16
25
36
49

要在给出2,2.5,3,3.5,4,4.5,5,5.5,6,6.5上进行插值。这里N=7,NJ=10。

(2)、主程序设计

N、NJ用参数说明语句定义,观测数据及插值点数用DATA语句输入,屏幕显示插值结果。

(3)、计算结果
 

程序清单

http://ftp.iap.ac.cn/cal/html/basic/interpolate2.html