多项式与插值.ppt

对n＋1个数据点，设，则得到n+1个线性方程，可以表示为矩阵形式求解该方程组可确定系数（或用polyfit(x,y,n)确定）第29页，共60页，星期日，2025年，2月5日例求用下列数据点拟合出的多项式的系数，并求当x取2.101、4.234时y的值，并画出已知数据点和拟合曲线。clear,clf,holdoffx=[1.1,2.3,3.9,5.1];y=[3.887,4.276,4.651,2.117];n=length(x)-1;c=polyfit(x,y,n);xi=[2.101,4.234];yi=polyval(c,xi)plot(x,y,ro)holdonxp=1.1:0.05:5.1;yp=polyval(c,xp);pauseplot(xp,yp)yi=4.14574.3007第30页，共60页，星期日，2025年，2月5日第31页，共60页，星期日，2025年，2月5日三、Lagrange插值多项式1、给定n+1个数据点：插值基函数：满足插值条件：的次数不超过n的lagrange多项式g(x)为：第32页，共60页，星期日，2025年，2月5日functiony=lagrange1(x0,y0,x)np=length(x0);y=zeros(size(x));fori=1:npz=ones(size(x));forj=1:npifj~=iz=z.*(x-x0(j))/(x0(i)-x0(j));endendy=y+y0(i)*zend2、MATLAB程序（法1）:调用格式：y=lagrange1(x0,y0,x)第33页，共60页，星期日，2025年，2月5日clearx0=[1.1,2.3,3.9,5.1];y0=[3.887,4.276,4.651,2.117];x=[2.101,4.234];y=lagrange1(x0,y0,x)例：写出拟合下面四个数据点的Lagrange插值公式，并计算x＝2.101、4.234时y的值。x0y01.13.8872.34.2763.94.6515.12.117结果为y=4.14574.3007第34页，共60页，星期日，2025年，2月5日由于lagrange插值基函数li(x)在节点xk(k=1,2,…,n)处满足又由于过n+1个数据点的n次插值多项式是唯一的,故li(x)也可看做是拟合下列数据点的n次多项式可编写程序求出2、MATLAB程序实现（法2）:进而求出lagrange插值公式第35页，共60页，星期日，2025年，2月5日functionp=lag_base(x)np=length(x);forj=1:npy=zeros(1,np);y(j)=1;p(j,:)=polyfit(x,y,np-1);end其调用格式为：p=lag_base(x)其中:x是数据点的横坐标数组，p是一个矩阵，它的第i行为li(x)的幂系数。先编写程序求出第36页，共60页，星期日，2025年，2月5日如前例也可求解如下：x=[1.1,2.3,3.9,5.1];y=[3.887,4.276,4.651,2.117];xi=[2.101,4.234];p=lag_base(x);np=length(x);s=0;fori=1:nps=s+p(i,:).*y(i);endsyi=polyval(s,xi)结果为：yi＝4.14574.3007第37页，共60页，星期日，2025年，2月5日3、Lagrange插值公式的微分插值公式微分为计算Lagrange插值多项式的一阶导数，可用polyder函数将p的每一行转换为一阶导数的系数数组。第38页，共60页，星期日，2025年，2月5日例：求出拟合下面四个数据点的Lagrange插值多项式在x＝2.101、4.234处的一阶导数值。x0y01.13.8872.34.2763.94.6515.