资源描述:
《插值法实验报告》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、1.(1).x=linspace(-5,5,11);y=1./(1+x.^2);x0=linspace(-5,5,101);y0=1./(1+x.^2);y1=interp1(x,y,x0,'spline')y2=interp1(x,y,x0);A=[ones(11,1)x'(x.^2)'(x.^3)'(x.^4)'(x.^5)'(x.^6)'(x.^7)'(x.^8)'(x.^9)'(x.^10)']a=Ay';y3=a(1)+a(2).*x0+a(3).*x0.^2+a(4).*x0.^3+a(5).*x0.^4+a(6).*x0.^5+a(7).*x0.^6+a(8).*x0.^7+a
2、(9).*x0.^8+a(10).*x0.^9+a(11).*x0.^10;plot(x0,y3,'r'),gtext('Lagr.'),holdon,plot(x0,y2,'b'),gtext('Pies.Lin.'),holdon,plot(x0,y1,'m'),gtext('Spline')holdoff(2).x=linspace(0,2*pi,11);y=cos(x);x0=linspace(0,pi,101);y0=cos(x0);y1=interp1(x,y,x0,'spline')y2=interp1(x,y,x0);A=[ones(11,1)x'(x.^2)'(x.^3)'(
3、x.^4)'(x.^5)'(x.^6)'(x.^7)'(x.^8)'(x.^9)'(x.^10)']a=Ay';y3=a(1)+a(2).*x0+a(3).*x0.^2+a(4).*x0.^3+a(5).*x0.^4+a(6).*x0.^5+a(7).*x0.^6+a(8).*x0.^7+a(9).*x0.^8+a(10).*x0.^9+a(11).*x0.^10;plot(x0,y3,'r'),gtext('Lagr.'),holdon,plot(x0,y2,'b'),gtext('Pies.Lin.'),holdon,plot(x0,y1,'m'),gtext('Spline')hold
4、off(3)x=linspace(0,pi,11);y=cos(x).^10;x0=linspace(0,pi,101);y0=cos(x0).^10;y1=interp1(x,y,x0,'spline')y2=interp1(x,y,x0);A=[ones(11,1)x'(x.^2)'(x.^3)'(x.^4)'(x.^5)'(x.^6)'(x.^7)'(x.^8)'(x.^9)'(x.^10)']a=Ay';y3=a(1)+a(2).*x0+a(3).*x0.^2+a(4).*x0.^3+a(5).*x0.^4+a(6).*x0.^5+a(7).*x0.^6+a(8).*x0.^7+a(
5、9).*x0.^8+a(10).*x0.^9+a(11).*x0.^10;plot(x0,y3,'r'),gtext('Lagr.'),holdon,plot(x0,y2,'b'),gtext('Pies.Lin.'),holdon,plot(x0,y1,'m'),gtext('Spline')holdoff2.(1).t=linspace(0,4,10);x=linspace(5,25,10);[t,x]=meshgrid(t,x);z=sin((t-1/200.*x)*2);[t1,x1]=meshgrid(linspace(0,4,100),linspace(5,25,100));z1=
6、interp2(t,x,z,t1,x1,'cubic');z2=interp2(t,x,z,t1,x1,'nearest');z3=interp2(t,x,z,t1,x1,'linear');mesh(t1,x1,z1);holdonmesh(t1,x1,z2);holdonmesh(t1,x1,z3);holdoff三次插值,在这三种插值方法中效果最好最邻近插值,效果不是很好线性插值,效果一般(2)取e=1.5M文件:x=linspace(-1,1,10);y=linspace(-1,1,10);[x,y]=meshgrid(x,y);z=3/10+sin(16/15.*x-1.5)+(si
7、n(16/15.*x-1.5)).^2+sin(16/15.*y-1.5)+(sin(16/15.*y-1.5)).^2[x1,y1]=meshgrid(linspace(-1,1,100),linspace(-1,1,100));z1=interp2(x,y,z,x1,y1,'cubic');z2=interp2(x,y,z,x1,y1,'nearest');z3=interp2(x,y,z,x