数学模型第四次实验报告.docx
《数学模型第四次实验报告.docx》由会员分享，可在线阅读，更多相关《数学模型第四次实验报告.docx（10页珍藏版）》请在冰豆网上搜索。
数学模型第四次实验报告
实验报告
（二）
院（系）理学院课程名称：
数学建模日期2012.03.01
班级
A0921
学号
01
实验室
文理楼209
专业
信息与计算科学
姓名
计算机号
12
实验
名称
Matlab三维图形绘法
成绩评定
所用
软件
Matlab
教师签名
实
验
目
的
或
要
求
学会Matlab软件中三维绘图的方法
掌握空间曲线图、曲面图。
等高线图等图形的绘制
实
验
步
骤
、
心
得
体
会
实验题目：
1.画出由参数方程
确定的空间曲线图形
（1）新建一个M-file
（2）输入命令：
t=0:
pi/50:
10*pi;
x=sin（t）;
y=cos（t）;
z=t;
plot3（x,y,z）
（3）结果为：
2.绘出旋转抛物面
的图形
（1）新建一个M-file
（2）输入命令：
x=-4:
0.1:
4;
y=-4:
0.1:
4;
[X,Y]=meshgrid（x,y）;
Z=X.^2+Y.^2;
surf（X,Y,Z）;
（3）结果为：
3.用ezsurf命令画出
输入命令：
ezsurf（'u*cos（v）','u*sin（v）','u',[0402*pi]）
结果为：
4.在范围
绘出曲面
的等高线
输入：
clear
x=-2:
0.1:
2;
y=-2:
0.1:
3;
[X,Y]=meshgrid（x,y）;
Z=X.*exp（-X.^2-Y.^2）;
[c1,h1]=contour（X,Y,Z）;
clabel（c1,h1）
输出：
5.（1994年全国大学生数学建模竞赛A题）要在一山区修建公路，首先测得一些地点的高程，数据如表所列：
4800
1350
1370
1390
1400
1410
960
940
880
800
690
570
430
290
210
150
4400
1370
1390
1410
1430
1440
1140
1110
1050
950
820
690
540
380
300
210
4000
1380
1410
1430
1450
1470
1320
1280
1200
1080
940
780
620
460
370
350
3600
1420
1430
1450
1480
1500
1550
1510
1430
1300
1200
980
850
750
550
500
3200
1430
1450
1460
1500
1550
1600
1550
1600
1600
1600
1550
1500
1500
1550
1500
2800
950
1190
1370
1500
1200
1100
1550
1600
1550
1380
1070
900
1050
1150
1200
2400
910
1090
1270
1500
1200
1100
1350
1450
1200
1150
1010
880
1000
1050
1100
2000
880
1060
1230
1390
1500
1500
1400
900
1100
1060
950
870
900
930
950
1600
830
980
1180
1320
1450
1420
1400
1300
700
900
850
840
380
780
750
1200
740
880
1080
1130
1250
1280
1230
1040
900
500
700
780
750
650
550
800
650
760
880
970
1020
1050
1020
830
800
700
300
500
550
480
350
400
510
620
730
800
850
870
850
780
720
650
500
200
300
350
320
0
370
470
550
600
670
690
670
620
580
450
400
300
100
150
250
y/x
0
400
800
1200
1600
2000
2400
2800
3200
3600
4000
4400
4800
5200
5600
试利用表中数据，绘制这一山区的地貌网格图、平滑地貌图、等高线图。
输入：
x=0:
400:
5600;%给出x轴的坐标
y=4800:
-400:
0;%给出y轴的坐标
%——————————录入每行（x，y）的高程————————————————
z1=[13501370139014001410960940880800690570430290210150];
z2=[13701390141014301440114011101050950820690540380300210];
z3=[138014101430145014701320128012001080940780620460370350];
z4=[1420143014501480150015501510143013001200980850750550500];
z5=[143014501460150015501600155016001600160015501500150015501500];
z6=[9501190137015001200110015501600155013801070900105011501200];
z7=[9101090127015001200110013501450120011501010880100010501100];
z8=[88010601230139015001500140090011001060950870900930950];
z9=[830980118013201450142014001300700900850840380780750];
z10=[740880108011301250128012301040900500700780750650550];
z11=[650760880970102010501020830800700300500550480350];
z12=[510620730800850870850780720650500200300350320];
z13=[370470550600670690670620580450400300100150250];
Z=[z1;z2;z3;z4;z5;z6;z7;z8;z9;z10;z11;z12;z13]
[X,Y]=meshgrid（x,y）;%生成（x，y）坐标的网格
surf（X,Y,Z）;%网格阴影图
figure
（2）;%新开一个窗口
contour（X,Y,Z,20）;%画出平面等高线
figure（3）;%再新开一个窗口
contour3（X,Y,Z,20）;%画三维等高线
xi=linspace（0,5600,50）;yi=linspace（0,4800,50）;%给出新的插值坐标
[XI,YI]=meshgrid（xi,yi）;
ZI=interp2（X,Y,Z,XI,YI,'*cubic'）;%对数据（xi,yi,zi）使用样条在网格{X,Y}上插值
figure（4）;
surf（XI,YI,ZI）;
shadinginterp
输出：
备注：
本实验报告用于各学科与计算机应用相关课程的实验，务必按时完成。
不交此报告者，本次实验为“不合格”。
九江学院
升级会员 