数值分析实验(四)插值法实验.docVIP

插值法实验 班别： 学生姓名： 学号： 实验目的 1.通过进行不同类型的插值，比较各种插值的效果，明确各种插值的优越性； 2.通过比较不同次数的多项式拟合效果，了解多项式拟合的原理； 3.利用matlab编程，学会matlab命令； 4.掌握拉格朗日插值法和牛顿插值法； 实验题目 (1) 调用拉格朗日插值程序求拉格朗日插值4次多项式在0.45 0.5 0.6 0.8上的值，并画出原函数与拉格朗日插值4次多项式的图像进行比较。 (2) 调用牛顿插值程序求牛顿插值4次多项式， 并求其在0.45 0.5 0.6 0.8上的值。 （3）选做：调用分段线性插值程序求其在0.45 0.5 0.6 0.8上的值，并画出原函数与分段线性插值多项式的图像进行比较。 实验原理 1. 拉格朗日插值 2. 牛顿插值的原理 实验内容与结果 1. 拉格朗日插值法 （1）相关程序 function yy=mlagr(x,y,xx) %用途：拉格朗日插值法求解 %格式：yy=mlagr(x,y,xx), x是节点向量, y是节点对应的函 %数值向量, xx是插值点(可以是多个), yp返回插值结果 n=length(x); m=length(xx); yy=zeros(1,m); c1=ones(n-1,1); c2=ones(1,m); for i=1:n xp=x([1:i-1,i+1:n]); yy=yy+y(i)*prod((c1*xx-xp*c2)./(x(i)-xp*c2)); End (原函数与拉格朗日插值4次多项式的图像进行比较的程序） x=[0.4 0.55 0.65 0.8 0.9];y=[0.41075 0.57815 0.69675 0.88811 1.02652];plot(x,y,+); x=0.4:0.01:0.6;y0=sin(x); plot(x,y0,*) 实验结果 x=[0.4 0.55 0.65 0.8 0.9];y=[0.41075 0.57815 0.69675 0.88811 1.02652];xx=[0.45 0.5 0.6 0.8]; yy=mlagr(x,y,xx) yy = 0.4653 0.5211 0.6367 0.8881 (原函数与拉格朗日插值4次多项式的图像进行比较) x=[0.4 0.55 0.65 0.8 0.9];y=[0.41075 0.57815 0.69675 0.88811 1.02652];plot(x,y,+); x=0.4:0.01:0.6;y0=sin(x); plot(x,y0,*) （3）实验结果分析 2. 牛顿插值法相关程序 （1）相关程序 %程序5.2--mnewp.m function yy=mnewp(x,y,xx) %用途：牛顿差值 %格式：yy=mnewp(x,y,xx), x是节点向量, y是节点对应的函 %数值向量, xx是插值点(可以是多个), yy返回插值结果 n=length(x); syms t; yy=y(1); y1=0; lx=1; for i=1:n-1 for j=i+1:n y1(j)=(y(j-1)-y(j))/(x(j-i)-x(j)); %计算差商 end c(i)=y1(i+1); lx=lx*(t-x(i)); yy=yy+c(i)*lx; %计算牛顿插值多项式的值 y=y1; end if nargin==3 yy=subs(yy,t,xx); else yy=collect(yy); yy=vpa(yy,6); end （2）实验结果 x=[0.4 0.55 0.65 0.8 0.9];y=[0.41075 0.57815 0.69675 0.88811 1.02652];xx=[0.45 0.5 0.6 0.8]; yy=mnewp(x,y,xx) yy = 0.4653 0.5211 0.6367 0.8881 （3）实验结果分析