Lab02_多项式插值计算及其收敛性实验.docVIP

实用标准文案 PAGE 文档 Lab02．多项式插值计算及其收敛性实验 【实验目的和要求】 1．使学生深入理解Langrage插值法和Newton插值法以两者之间的异同，能用Matlab语言编写按Langrage插值法和Newton插值法计算插值的程序； 2．用所编写的程序进行插值计算、验证Runge现象、分析插值多项式的收敛性； 3．使学生深入理解教材介绍的两种分段低次插值法，熟悉掌握函数interp1的使用； 4．使用函数interp1用不同方法进行插值计算，对教材介绍的几种分段低次插值法进行分析比较。 【实验内容】 1．根据Matlab语言特点，描述Langrage插值法和Newton插值法。 2．用Matlab语言编写按Langrage插值法和Newton插值法计算插值的程序。 3．对，分别取3个，5个、9个、11个等距节点，用所编写的程序进行插值计算并画图，以验证Runge现象、分析插值多项式的收敛性。 4．用函数interp1，对，用n=11个节点(等分)作分段线性插值、分段Hermit插值和三次样条插值，用m=101个插值点(等分)作图，比较结果。 【实验仪器与软件】 1．CPU主频在1GHz以上，内存在128Mb以上的PC； 2．Matlab 6.0及以上版本。 实验讲评： 实验成绩： 评阅教师： 200 年 月 日 Lab02．多项式插值计算及其收敛性实验 一、算法描述 1、lagrange插值描述： 构造通过n+1个节点的n次插值多项式，假定它满足条件=，j=0,1,…,n.定义n次插值基函数.若n次多项式（j=0,1,…,n）在n+1个节点上满足条件，…,n,则根据基函数可以确定 由此则可以求出lagrange插值多项式 2、Newton插值描述： f在插值点上的值为，要求n次插值多项式满足条件则插值函数为 通过线性插值、均差可确定，即：由此下去.最后得到Newton插值公式为： 其中根据均差表可以把均差表示为：表示第i行第j列的均差值。例如、 二、算法程序 lagrange程序： function y=xm_lagrange(x0,y0,x) if length(x0)~=length(y0) error(x0y0的维数不相等); else n=length(x0); m=length(x); for i=1:m z=x(i); s=0.0; for k=1:n p=1.0; for j=1:n; if j~=k p=p*(z-x0(j))/(x0(k)-x0(j)); end end s=s+p*y0(k); end y(i)=s; end end end newton插值公式： function y=xm_newdon(x0,y0,x) if length(x0)~=length(y0) error(x0y0的维数不相等); else n=length(x0); d=y0; for j=2:n for i=j:n y0(i)=(d(i)-d(i-1))/(x0(i)-x0(i-j+1)); end d=y0; end m=length(x); for k=1:m s=0;p(1)=1.0;z=x(k); for A=1:n p(A+1)=p(A)*(z-x0(A)); s=s+d(A)*p(A); end y(k)=s; end end end 三、插值计算 1、lagrange插值和newton插值分别取3个等距节点进行插值计算 。 clear clc x0=linspace(-5,5,3); y0=1./(1+x0.^2); x=-5:0.1:5; y=1./(1+x.^2); y1=xm_lagrange(x0,y0,x); y2=xm_newton(x0,y0,x); plot(x,y,b,x,y1,r) 红色曲线为lagrange插值函数图。 plot(x,y,b,x,y2,g) 浅绿色曲线为Newton插值函数图。 2、1、lagrange插值和newt