《数学建模与数学实验》上机实验报告

  • 格式:doc
  • 大小:143.00 KB
  • 文档页数:4

1 / 4 成都信息工程大学

《数学建模与数学实验》上机实验报告

专业 信息与计算科学 班级 姓名 学号

实验日期 成绩等级 教师 评阅日期

[问题描述]

下表给出了某一海域以码为单位的直角坐标Oxy 上一点(x,y)(水面一点)以英尺为单位的水深z,水深数据是在低潮时测得的,船的吃水深为5英尺,问在矩形区域(75,200)x (-50,150)里那些地方船要避免进入。

低潮时测得的水深数据

[模型]

设水面一点的坐标为(x,y,z),用基点和插值函数在矩形区域(75,200)*(-50,150)内做二维插值、三次插值,然后在作出等高线图。

x 129 140 103.5 88 185.5 195 105.5

157.5 107.5 77 81 162 162 117.5

y 7.5 141.5 23 147 22.5 137.5 85.5

-6.5 -81 3 56.5 -66.5 84 -33.5

z 4 8 6 8 6 8 8

9 9 8 8 9 4 9 2 / 4 [求解方法]

使用matlab求解:

M文件:water.m

x=[129 140 103.5 88 185.5 195 105.5 157.5 107.5 77 81

162 162 117.5];

y=[7.5 141.5 23 147 22.5 137.5 85.5 -6.5 -81 3 56.5 -66.5

84 -33.5];

z=[-4 -8 -6 -8 -6 -8 -8 -9 -9 -8 -8 -9 -4 -9];

cx = 75:0.5:200;

cy = -50:0.5:150;

[cx,cy]=meshgrid(cx,cy);

作出曲面图:

代码如下:

>> water

>> cz=griddata(x,y,z,cx,cy,'cubic');

>> meshz(cx,cy,cz)

>> xlabel('X'),ylabel('Y'),zlabel('Z')

>>

作出等高线图:

代码如下:

>> water

>> cz=griddata(x,y,z,cx,cy,'cubic');

>> figure(2)

>> contour(cx,cy,cz,[-5,-5],'r')

>> hold on

>> plot(x,y,'*')

>> xlabel('X'),ylabel('Y')

[结果] 3 / 4 插值结果

等值图:

4 / 4

[结果分析及结论]

根据等值图可看出:红色区域为危险区域,所以船只要避免进入。