用于函数近似表达、求导和积.PDFVIP

南京林产工业学院学报 一九八三年 No.2 JOURNAL OF NANJING TECHNOLQGICAL 第二期 COLLEGE OF FOREST PRODUCTS 198 3 用于函数近似表达、求导和积 分的一组计算公式 潘平奇 (基础课吾~) 一、引言 在实际中，往往需要把从实验而来的列表函数一点列〈冉、 y，) , ì = 0 , 1 ,…, n , 用较简单的分析式近似表达，求其近似导函数或积分，这实际上是一个平面点列的插值或拟 合问题。流行的方法常不能保证对曲线的凸性要求，产生多余拐点或波动，尤其当点子很 多，或点列蕴含曲率剧变、平直微曲、曲直相接时更是难以适应。 从处理方式讲，流行方法大都可归于整体性方法，在作插值或拟合时着眼于点列整 体。例如多项式插值，一般是寻求一个通过全体点子的多项式，因而有n+l个点就得到n次多 项式。当n很大时该多项式一般拐点很多就很难保凸。另外，再以三次样条函数括值为例z 连 续性方程〈如所谓三弯距方程〉加上适当的边界条件构成一个n+l阶线性代数方程组，从中 解出各点的一阶或二阶导数，由此确定一条曲线。由于如此而来的各点的导数值显然同全体 点子有关，故{王一点的改变必然会引起远离该点处曲线的改变，这对于实际计算是不利的。 实际上，绘曲线总是一小段一小段绘的。绘每一小段时，我们也不去注意全体点子，而 只着眼于邻近几个点，因为它们可以确定这段曲线的性态，这就启示我们采用局部性方 法。本文提出的就是这种方法，它克服了整体性方法所固有的种种缺点而具有以下长 处z 1 .对各类点列适应性好，描述能力强并能保凸， 2 .一点的改变只影响邻近，对曲线的其余部分无影响F 3. 点子可以很多 (n很大) ，计算误差不随点子数目增加而增加， 4 .使用计算机绘图时，可边计算边绘图，而不必象整体性方法那样必须待全部计 算好才能绘图F 5 .编制程序简单。 -141 一 - 多重节点插值 -、 i支XOX …X ，称Xl 为节点，被逼近函数f(x)在X\ 处的 0 至l 阶导数YI 他)为巳知， n 1 k =0， 1 ，… ， 1 。 我们在每个小区间 [XI ， Xl+1J 上，用一个21 + 1 次多项式P21+1 ， 1(XH乍f( 川的近似，要求 它满足21 + 2个条件， {pzftiJ1 =yrtJ • k =O. 1. …• 1 ( 1 ) P2βid 〈 Xi+1 〉 =Yt(JT 这是关于P l+ l(X)的系数的一个正规的21 + 2阶线性代数方程组，可以唯一地确定这组系 2 1 数，由此得到3\J是一条逐段21+1次多项式曲线。显然每段由该段两端点处的条件独立确 定，与其他节点处的条件无关，然而整条曲线有直至l 阶的连续导数，即 l 阶平滑。 当然没有必要从( 1 )去解出P2+1 ， t( 川的21 +