用MATLAB进行数据插值探讨.ppt

例 [x,y]=meshgrid(-3:.6:3,-2:.4:2); z=(x.^2-2*x).*exp… (-x.^2-y.^2-x.*y); surf(x,y,z), 选较密的插值点，用默认的线性插值算法进行插值 [x1,y1]=meshgrid(-3:.2:3,-2:.2:2); z0=interp2(x,y,z,x1,y1); surf(x1,y1,z0) 立方和样条插值： z1=interp2(x,y,z,x1,y1,cubic); z2=interp2(x,y,z,x1,y1,spline); surf(x1,y1,z1),axis([-3,3,-2,2,-0.7,1.5]) figure;surf(x1,y1,z2),axis([-3,3,-2,2,-0.7,1.5]) 算法误差比较 z=(x1.^2-2*x1).*exp(-x1.^2-y1.^2-x1.*y1); surf(x1,y1,abs(z-z1)) figure;surf(x1,y1,abs(z-z2)) figure;surf(x1,y1,abs(z-z0)) 山区地形地貌图 已知某处山区地形选点测量坐标数据为： x=0? 0.5? 1? 1.5? 2? 2.5? 3? 3.5? 4? 4.5? 5 y=0? 0.5? 1? 1.5? 2? 2.5? 3