资源描述:
《《数学建模与数学实验》上机实验报告》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、-----用心做文档值得下载值得拥有-------《数学建模与数学实验》上机实验报告专业班级姓名学号实验日期成绩等级教师评阅日期[问题描述]1某海域测得一些点(x,y)处的水深z处由下表给出,船的吃水深度为5英尺在矩形区域(75,200)×(-50,150)里的哪些地方船要避免进入X129140103.588185.5195105.5157.5107.57781162162117.5Y7.5141.52314722.5137.585.5-6.5-81356.5-66.584-33.5Z48686889988949[模型]模型假设:1表中的统计数据均真实可
2、靠2矩形区域外的海域不对矩形海域造成影响3短期内海水不涨潮4无人为干涉导致z发生变化符号说明:x表示海域的横向位置,y表示海域的纵向位置,z表示海域的深度建立微分模型:a输入插值基点数据b在矩形区域(75,200)×(-50,150)作二维插值,运用三次插值法c作海底曲面图d作出水深小于5的海域范围,即z=5的等高线[求解方法]建立M文件hd1.m如下%数据输入x=[129140103.588185.5195105.5157.5107.57781162162117.5];y=[7.5141.52314722.5137.585.5-6.5-81356.5-
3、66.584-33.5];z=[-4-8-6-8-6-8-8-9-9-8-8-9-4-9];%插值cx=75:0.3:200;-----用心做文档值得下载值得拥有-------cy=-50:0.3:150;cz=griddata(x,y,z,cx,cy','cubic');%作表面图meshz(cx,cy,cz);xlabel('x-axis'),ylabel('y-axis'),zlabel(‘z-axis’)colorbarfigure(2)[C,h]=contour(cx,cy,cz,[-5,5])grid,holdonplot(x,y,'*');
4、[cx,cy]=meshgrid(cx,cy);xlabel('x-axis'),ylabel('y-axis')[结果]-----用心做文档值得下载值得拥有-------图1图2C=1.0e+03*1至10列-0.00500.12590.12600.12630.12660.12690.12720.12750.12760.12781.1930-0.0026-0.0026-0.0027-0.0027-0.0027-0.0027-0.0026-0.0026-0.002611至20列0.12810.12840.12870.12870.12900.12930.1
5、294-----用心做文档值得下载值得拥有-------0.12960.12990.1301-0.0025-0.0024-0.0023-0.0023-0.0022-0.0021-0.0020-0.0019-0.0018-0.001721至30列0.13020.13050.13070.13080.13110.13120.13140.13170.13170.1320-0.0016-0.0015-0.0014-0.0013-0.0011-0.0011-0.0010-0.0008-0.0008-0.000631至40列0.13210.13230.13260.13
6、260.13290.13300.13320.13340.13350.1338-0.0005-0.0004-0.0002-0.00020.00000.00010.00020.00040.00050.000741至50列0.13380.13410.13420.13440.13450.13470.13490.13500.13520.13530.00070.00090.00100.00120.00130.00140.00160.00170.00190.002051至60列0.13560.13560.13590.13590.13620.13620.13650.13
7、650.13670.13680.00220.00220.00250.00250.00280.00280.00310.00310.00340.003561至70列0.13700.13710.13730.13740.13770.13770.13800.13810.13830.13860.00370.00380.00400.00410.00430.00430.00450.00460.00470.004971至80列0.13860.13890.13920.13920.13950.13980.13980.14010.14040.14040.00490.00510.
8、00520.00520.00540.00550.0055-----用心做文档值得