拉格朗日插值、牛顿插值的matlab代码.doc

实验五 多项式插值逼近 信息与计算科学金融 崔振威 201002034031 实验目的： 拉格朗日插值和牛顿插值的数值实现 实验内容：p171.1、p178.1、龙格现象数值实现 实验要求： 根据所给题目构造相应的插值多项式， 编程实现两类插值多项式的计算 试分析多项式插值造成龙格现象的原因 主程序 1、拉格朗日 function [c,l]=lagran(x,y) %c为多项式函数输出的系数 %l为矩阵的系数多项式 %x为横坐标上的坐标向量 %y为纵坐标上的坐标向量 w=length(x); n=w-1; l=zeros(w,w); for k=1:n+1 v=1; for j=1:n+1 if k~=j v=conv(v,poly(x(j)))/(x(k)-x(j)) %对多项式做卷积运算 end end l(k,:)=v; end c=y*l; 牛顿插值多项式主程序 function [p2,z]=newTon(x,y,t) %输入参数中x,y为元素个数相等的向量 %t为插入的定点 %p2为所求得的牛顿插值多项式 %z为利用多项式所得的t的函数值。 n=length(x); chaS(1)=y(1); for i=2:n x1=x;y1=y; x1(i+1:n)=[]; y1(i+1:n)=[]; n1=length(x1); s1=0; for j=1:n1 t1=1; for k=1:n1 if k==j %如果相等则跳出循环 continue; else t1=t1*(x1(j)-x1(k)); end end s1=s1+y1(j)/t1; end chaS(i)=s1; end b(1,:)=[zeros(1,n-1) chaS(1)]; cl=cell(1,n-1); %cell定义了一个矩阵 for i=2:n u1=1; for j=1:i-1 u1=conv(u1,[1 -x(j)]); %conv()用于多项式乘法、矩阵乘法 cl{i-1}=u1; end cl{i-1}=chaS(i)*cl{i-1}; b(i,:)=[zeros(1,n-i),cl{i-1}]; end p2=b(1,:); for j=2:n p2=p2+b(j,:); end if length(t)==1 rm=0; for i=1:n rm=rm+p2(i)*t^(n-i); end z=rm; else k1=length(t); rm=zeros(1,k1); for j=1:k1 for i=1:n rm(j)=rm(j)+p2(i)*t(j)^(n-i); end z=rm; end end plot(t,z,'y',x,y,'*r') %输出牛顿插值多项式的函数图 p171.1 、f(x)=ex 解：在matlab窗口中输入： >> x=[0 0.2 0.4 0.6 0.8 1]; >> y=[exp(0) exp(0.2) exp(0.4) exp(0.6) exp(0.8) exp(1)] y = 1.0000 1.2214 1.4918 1.8221 2.2255 2.7183 >> [c,l]=lagran(x,y) 可以得出输出结果为： c = 0.0139 0.0349 0.1704 0.4991 1.0001 1.0000 l = -26.0417 78.1250 -88.5417 46.8750 -11.4167 1.0000 130.2083 -364.5833 369.7917 -160.4167 25.0000 0 -260.4167 677.0833 -614.5833 222.9167 -25.0000 0 260.4167 -625.0000 510.4167 -162.5000 16.6667 0 -130.2083 286