函数插值.docVIP

插值在数学发展史上是个古老问题.它是和拉格朗日（Lagrange）,牛顿（Newton），高斯（Gauss）等著名数学家的名字连在一起的.它最初来源于天体计算——由若干观测值（即节点）计算任意时刻星球的位置（即插值点和插值）——的需要.现在，虽然人们已很少需要用它从函数表计算函数值了，但是插值仍然在诸如机械加工等工程技术和数据处理等科学研究中有着许多直接的应用.另一方面，插值又是数值微分、数值积分、常微分方程数值解等数值计算的基础.本章介绍插值问题的提法和常用的插值.例如，拉格朗日多项式插值、牛顿（Newton）插值、埃尔米特（Hermite）插值、高元插值、分段插值和三次样条插值及其MATLAB程序等. 6.1 插值问题及其误差 本节介绍插值问题的有关概念，求解的基本思路，插值多项式存在和唯一性，插值余项和误差限.最后介绍MATLAB命令POLY(A), CONV(A, B) 和DECONV(A, B)等的功能和调用方法. 6.1.2 与插值有关的MATLAB 函数 在这里我们把多项式的按降幂排列的系数所构成的向量V称为多项式的系数向量.下面介绍MATLAB函数POLY2SYM,POLY, CONV 和DECONV等的功能和调用方法. (一) POLY2SYM 函数 POLY2SYM 的功能是把多项式的系数向量转化为符号多项式，其调用格式如下： 调用格式一：poly2sym (C) 如果输入维向量C，其中C是多项式的系数向量，则运行后输出是以为自变量的次多项式. 调用格式二：f1=poly2sym(C,V) 或 f2=poly2sym(C, sym (V) ), POLY2SYM(C,V) 和 POLY2SYM(C,SYM(V) ) 两者都用于返回以V为自变量，系数向量C的符号多项式. (二) POLYVAL 函数 POLYVAL函数的主要功能是估计多项式的值. 调用格式：Y = polyval(P,X) （1）如果输入数X 和维向量P，其中P是次多项式的系数向量，则运行后输出在的函数值，即 . (三) POLY 函数 POLY命令的主要功能是把根转换成多项式的系数向量. 调用格式：Y = poly (V) （1）如果输入的V是维向量，V的元素是方程的根，则运行后输出维的行向量P，其中 . （2）如果输入的V是阶方阵，运行后输出维的行向量P，P是方阵V的特征多项式的系数向量. (四) CONV 函数 CONV 的功能是卷积和多项式的积. 调用格式：C =conv (A, B) 如果输入维向量A 和维向量B，其中A是多项式的系数向量，B是多项式的系数向量，则运行后输出维的行向量C ，其中C是多项式与的积的系数向量. 例 6.1.2 求三个一次多项式、和的积.它们的零点分别依次为0.4,0.8,1.2. 解 我们可以用两种MATLAB程序求之. 方法1 如输入MATLAB程序 X1=[0.4,0.8,1.2]; l1=poly(X1), L1=poly2sym (l1) 运行后输出结果为 l1 = 1.0000 -2.4000 1.7600 -0.3840 L1 = x^3-12/5*x^2+44/25*x-48/125 方法2 如输入MATLAB程序 P1=poly(0.4);P2=poly(0.8);P3=poly(1.2); C =conv (conv (P1, P2), P3) , L1=poly2sym (C) 运行后输出的结果与方法1相同. (五) DECONV 函数 DECONV 的功能是逆卷积和多项式的除法.商和余式. 调用格式：[Q,R] =deconv (B,A) 如果输入维向量A 和维向量B，其中A是多项式的系数向量，B是多项式的系数向量，则运行后输出 被除所得的商的系数向量Q和余式的系数向量R，使得B = conv(A,Q) + R. (六) roots(poly(1:n))命令 调用格式：roots(poly(1:n)) 如果输入的1: n是1到的等差数列，且此等差数列的元素是方程的全部根，则运行后输出维的行向量P，P元素是多项式的零点，即方程的全部近似根. (七) det(a*eye(size (A)) - A)命令 调用格式：b=det(a*eye(size (A)) - A) 如果输入的是阶方阵和数a，运行后输出的值b.如果b是零，则数a是的特征值；如果b不是零，则数a不是的特征值. 6.2 拉格朗日（Lagrange）插值及其MATLAB程序 先考察几种简单的情形，然后推广到一般地情形. 6.2.1 线性插值及其MATLAB程序 例6.2.1 已知函数在上具有二阶连续导数，，且满足条件.求线性插值多项式和函数值，并估计其误差. 解