山区城区地貌matlab.docx

《山区城区地貌matlab.docx》由会员分享，可在线阅读，更多相关《山区城区地貌matlab.docx（7页珍藏版）》请在第一文库网上搜索。

1、数学模型与数学实验实验报告教学班：，姓名：学号：实险:茨:目名和山区地貌与城区地貌绘制问题实脸地点实验日期实验环境实聆目的及意义目的及意义：1、掌握二维插值方法2、掌握一些基础命令的使用方法3、培养利用MAT1AB画出地貌图以及画出等高线的能力4、养成将算法转化为程序的思维能力实验yT务任务：(问题：、掌握插值基点为网格结点和散乱节点的二维插值方法入利用函数画出所确定的函数图像以及其等高线1、山区地貌：在某山区测得一些地点的高题表1所示。平面区域为：1200=x=4000,1200=yyz):contour(X,Y,Z)指定Z中各值的X和y坐标。7、c1abe1:c1abe1(C,h)为当前等

2、高线图添加标签，将旋转文本插入每条等高线。等高线必须足够长以容纳标签，否则c1abe1无法插入标签。如果没有等高线矩阵C,则将C替换为o8、SUbP1Ot:画多幅图4、算法代码：算法：StePI:提取文档中的数据Step2:网格化x,yStep3:二维插值Step4:绘图代码：山区：%输入已知量x=12001600200024002800320036004000;y=1200160020002400280032003600;Xi=I200:4000;yi=1200:3600;%将xi,yi网格化并进行插值X,Y=meshgrid(xizyi);z=interp2(x,y,height,X,Y,

3、cubic,);%画地貌图subp1ot(1,2,1);mesh(X,Y,z);tit1e山区地貌图,);%画等高线subp1ot(1,2,2);C=contour(X,Y,z);c1abe1(C);tit1e(,山区等高线,);城区：%提取原始数据X=S1(:,1)1；Y=S1(:,2),;Z=S1(:,3),;x=0:40:28700;y=0:40:18500;%网格化并进行散点插图xi,yi=meshgrid(x,y);z=griddata(X,Y,Z,xi,yi,cubic,);%画地貌图subp1ot(1,2,1);mesh(x,y,z);tit1e(城区地貌图D;%画等高线subp

4、1ot(1,2,2);C=contour(x,y,z);c1abe1(C);tit1e城区等高线D;山区地貌图扩展z=1130,1250,1280,1230,1040,900,500,700;1320,1450,1420,1400,1300,700,900,85x=1200:400:4000;y=1200:400:3600;xi=1200:10:4000;%X的范困yi=1200:10:3600;%y的范围XI,Y1=meshgrid(xi,yi);%建*.】维网格Z=griddata(x,y,z,XI,YI,cubic*);%“方插值SUbP1ot(2,2,DQb在乎铺位置创建两行两列的坐标

5、区.位广一行例mesh(X1Y1Z)%网格图tit1e(山区地貌图(：次多项式)用标题subp1ot(2,2,2)%位于第：行第一列contour(XI,YI,Z,Showtext*,On)%山区地貌等高线绘制tit1e(等高线)%标题Z=griddata(x,y,z,XI,YI,nearest)%最邻近插值SUbPIot(2,2,3)%在平轴位置创建两行两列的坐标区.(；/J-行：列IneSh(X1,YI,Z)%网格图tit1e(山区地貌图(临近插值法)设标题Z=griddata(x,y,z,XI,YI,1inear*);%最邻近插值SUbPIOt(2,2,4)(WFr铺位置创建两行两列的坐标区，位行：列IneSh(X1Y1Z)%网格图tit1e(山区地貌图(线性插值法)%标题本次实验模型以两种不同数据类型为基础，建立起网格结点和散乱结点的插值模型的过程较为轻松，算法设计过程中，先要建立网格结点，再确定插值方法，再进行后面的求值和画图任务。此次软件应用较不熟练，许多函数也未能灵活使用。但此次掌握许多函数用法建立起构造数学模型的思维方式，应继续发扬认真钻研的态度，为以后的数学建模打下良好基础。