插值实验讲义 计算方法课件及实验 教学课件.docVIP

与插值有关的MATLAB 函数 预备知识：matlab多项式计算 在matlab里，多项式由一个行向量表示，它的系数是按降序排列。 例如，输入多项式x4－12x3＋0x2＋25x＋116 ? p=[1 -12 0 25 116] 注意，必须包括具有零系数的项。除非特别地辨认，MATLAB无法知道哪一项为零。给出这种形式，用函数roots找出一个多项式的根。 r=roots(p) 结果： r =11.7473 2.7028 -1.2251 + 1.4672i -1.2251 - 1.4672i 在MATLAB中，无论是一个多项式，还是它的根，都是向量，MATLAB按惯例规定，多项式是行向量，根是列向量。 (一) poly函数 求多项式的系数（由已知根求多项式的系数） 调用格式：Y = poly (V) X1=[1,2]; l1=poly(X1), l1=[ 1 -3 2] (二) poly2sym函数 将多项式系数向量转为解析式形式 调用格式：poly2sym (C) poly2sym(l1) ans = x^2 - 3*x + 2 (三) polyval函数 给定多项式求值Y = polyval(P,X) polyval(l1,2) ans = 0 (四) conv函数 求多项式相乘后的多项式系数。(执行两个数组的卷积) 调用格式：poly2sym (C) 例如 求(x^2 - 3*x + 2)* (x-3) l1=[ 1 -3 2]; l2=[ 1 -3]; l=conv(l1,l2) 结果l = 1 -6 11 -6 poly2sym(l) 结果ans =x^3 - 6*x^2 + 11*x – 6 poly([1,2,3]) 结果ans = 1 -6 11 -6 (五) 其它多项式函数 deconv(a,b)—多项式除法 polyder(p)—多项式求导 1.2拉格朗日（Lagrange）插值及其MATLAB程序 1.2.1 线性插值及其MATLAB程序 例6.2.1 已知函数在上具有二阶连续导数，，且满足条件.求线性插值多项式和函数值，并估计其误差. 解 输入程序 X=[1,3];Y=[1,2]; l01= poly(X(2))/( X(1)- X(2)); l11= poly(X(1))/( X(2)- X(1)); l0=poly2sym (l01), l1=poly2sym (l11) P = l01* Y(1)+ l11* Y(2), L=poly2sym (P), x=1.5; Y = polyval(P,x) 运行后输出基函数l0和l1及其插值多项式的系数向量P（略）、插值多项式L和插值Y为 l0 = l1 = L = Y = -1/2*x+3/2 1/2*x-1/2 1/2*x+1/2 1.2500 输入程序 M=5;R1=M*abs((x-X(1))* (x-X(2)))/2 运行后输出误差限为 R1 = 1.8750 1.2.2 拉格朗日多项式和基函数的MATLAB程序 求拉格朗日插值多项式和基函数的MATLAB主程序 function [C, L,L1,l]=lagran1(X,Y) m=length(X); L=ones(m,m); for k=1: m V=1; for i=1:m if k~=i V=conv(V,poly(X(i)))/(X(k)-X(i)); end end L1(k,:)=V; l(k,:)=poly2sym (V) end C=Y*L1;L=Y*l 例1.2 给出节点数据，，，， ，，作五次拉格朗日插值多项式和基函数，并写出估计其误差的公式. 解 在MATLAB工作窗口输入程序 X=[-2.15 -1.00 0.01 1.02 2.03 3.25]; Y=[17.03 7.24 1.05 2.03 17.06 23.05]; [C, L ,L1,l]= lagran1(X,Y) 运行后输出五次拉格朗日插值多项式L及其系数向量C，基函数l及其系数矩阵L1如下 C = -0.2169 0.0648 2.1076 3.3960 -4.5745 1.0954 L = 1.0954-4.5745*x+3.3960*x^2+2.1076*x^3+0.0648*x^4-0.2169*x^5 L1 = -0.0056 0.0299 -0.