资源描述:
《数学实验报告二》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、贵州师范大学数学与计算机科学学院学生实验报告课程名称:数学建模班级:08数本(2)班实验H期:2011年6月25日学号:080702010069姓名:田双标指导教师:周晓军老师实验成绩:一、实验名称用插值法绘制地貌图和等高线图二、实验目的及要求1.了解插值的基本内容。2.掌握用数学软件包求解插值问题。3•比较不同插值方法作山区的地貌图和等高线图。三、实验环境每人一台计算机,要求安装WindowsXP操作系统,Microsoftoffice2003>MATLAB6.5四、实验内容山区地貌:在某山区测得一些地点的高程如下表:(平面区域1200<=x<=
2、4000,1200<=y<=3600),试作出该山区的地貌图和等高线图,并对儿种插值方法进行比较。120016002000240028003200360040001200113012501280123010409005007001600132014501420140013007009008502000139015001500140090011001060950240015001200U0013501450120011501010280015001200noo1550160015501380107032001500155016001550160016
3、001600155036001480150015501510143013001200980通过此例对最近邻点插值、双线性插值方法和双三次插值方法的插值效果进行比较。五、实验步骤及算法描述1.先在三维坐标画出原始数据,画出粗糙的温度分布曲图.2.以平滑数据,在x、y方向上每隔50个单位的地方进行插值3.画出不同插值法的山区等高线图和地貌曲面图.六、调试过程及实验结果1.作该山区的地貌图(1)作粗糙的地貌图»x二1200:400:4000;7=1200:400:3600;沪[11301250128012301040900500700;...132014
4、50142014001300700900850;•••139015001500140090011001060!)50;...15001200110013501450120011501010;...15001200110015501600155013801070;...15001550160015501600160016001550;...1480150015501510143013001200980];>>figure(1);moshz(x,y,z)xlabel('X'),ylabel('Y'),zlabel('Z')(2)最邻近插值的地貌图»xi
5、=1200:50:4000;yi=1200:50:3600;»figure(2)z1i=interp2(x,y,z,xi,yi1/nearest);surfc(xi,yi,zli)xlabel(X),ylabel('Y‘)忆label(Z)(2)双线性插值作地貌图>>figure(3)z2i=interp2(x,y,z,xi,yi‘);surfc(xi,yi,z2i)xlabel('X’),ylabelCY‘),zlabel('Z‘)(2)双三次插值作地貌图>>figure(4)z3i=interp2(x,y,z,xi,yi,'cubic');su
6、rfc(xi,yi,z3i)xlabel('X’),ylabelCY‘),zlabelCT)20002000Y060004000400020006000比较不同插值方法的等高线图:»figure(5)subplot(1,3,l),contour(xi,yi,z1i,10/r1);4500400035003000250015001OOO500O七•总结:清楚插值数学模型的建立,掌握插值的基本方法,并会应用简单的插值方法,掌握实际生活中有哪些应用,在插值的时候要注意维数的变化和数据变化,还要注意程序的编写。八、附录(源程序清单)x=1200:400:4
7、000;y=1200:400:3600;z=[11301250128012301040900500700;...13201450142014001300700900850.139015001500140090011001060950;・・15001200110013501450120011501010;15001200110015501600155013801070;15001550160015501600160016001550;1480150015501510143013001200980];figure(1);meshz(x,y,z)xlab
8、el('X'),ylabel('Y')fzlabel('Z')xi=1200:50:4000;st');surfc(xi,