原创力文档- 1、本文档共107页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
例 [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? 3.5? 4? 4.5? 5? 5.5? 6 海拔高度数据为： z=89 90 87 85 92 91 96 93 90 87 82 ?? 92 96 98 99 95 91 89 86 84 82 84 ?? 96 98 95 92 90 88 85 84 83 81 85 ?? 80 81 82 89 95 96 93 92 89 86 86 ?? 82 85 87 98 99 96 97 88 85 82 83 ?? 82 85 89 94 95 93 92 91 86 84 88 ?? 88 92 93 94 95 89 87 86 83 81 92 ?? 92 96 97 98 96 93 95 84 82 81 84 ?? 85 85 81 82 80 80 81 85 90 93 95 ?? 84 86 81 98 99 98 97 96 95 84 87 ?? 80 81 85 82 83 84 87 90 95 86 88 ?? 80 82 81 84 85 86 83 82 81 80 82 ?? 87 88 89 98 99 97 96 98 94 92 87 山区地形地貌图 程序 原始地貌图程序： x=0:.5:5; y=0:.5:6; [xx,yy]=meshgrid(x,y); z=[89 90 87 85 92 91 96 93 90 87 82 92 96 98 99 95 91 89 86 84 82 84 96 98 95 92 90 88 85 84 83 81 85 80 81 82 89 95 96 93 92 89 86 86 82 85 87 98 99 96 97 88 85 82 83 82 85 89 94 95 93 92 91 86 84 88 88 92 93 94 95 89 87 86 83 81 92 92 96 97 98 96 93 95 84 82 81 84 85 85 81 82 80 80 81 85 90 93 95 84 86 81 98 99 98 97 96 95 84 87 80 81 85 82 83 84 87 90 95 86 88 80 82 81 84 85 86 83 82 81 80 82 87 88 89 98 99 97 96 98 94 92 87]; mesh(xx,yy,z) 山区地形地貌图 结果 x=-3+6*rand(200,1);y=-2+4*rand(200,1); z=(x.^2-2*x).*exp(-x.^2-y.^2-x.*y); [x1,y1]=meshgrid(-3:.2:3,-2:.2:2); z1=griddata(x,y,z,x1,y1,cubic); surf(x1,y1,z1),axis([-3,3,-2,2,-0.7,1.5]) z2=griddata(x,y,z,x1,y1,v4); figure;surf(x1,y1,z2),axis([-3,3,-2,2,-0.7,1.5]) 误差分析 z0=(x1.^2-2*x1).*exp(-x1.^2-y1.^2-x1.*y1); surf(x1,y1,abs
文档评论（0）