实验六测量数据可视化Word文档下载推荐.docx
《实验六测量数据可视化Word文档下载推荐.docx》由会员分享，可在线阅读，更多相关《实验六测量数据可视化Word文档下载推荐.docx（17页珍藏版）》请在冰豆网上搜索。
3、图形加注功能
例：
t=0:
0.1:
10；
y1=sin（t）;
y2=cos（t）;
plot（t,y1,'
r'
t,y2,'
b--'
x=[1.7*pi;
1.6*pi];
y=[-0.3;
0.8];
s=[‘sin（t）’;
‘cos（t）’];
text（x,y,s）;
%指定位置加标注
正弦和余弦曲线'
legend（'
正弦'
'
余弦'
）
label（'
时间t'
ylabel（'
正弦、余弦'
%自动打开
grid
axissquare
gridoff
4.再看：
下面表格是某城市的采样点的具体位置及海拔，取前50个样本点（共319个取样点）。
附件1.取样点位置及其所属功能区
编号
x（m）
y（m）
海拔（m）
功能区
1
74
781
5
4
26
5635
7965
29
2
1373
731
11
27
5394
8631
12
3
1321
1791
28
5291
7349
10
1787
4742
7293
9
1049
2127
30
4948
6
1647
2728
31
5567
6782
7
2883
3617
15
32
7004
6226
8
2383
3692
33
7304
5230
2708
2295
22
34
7048
4600
24
2933
1767
35
8180
4496
4233
895
36
9328
4311
4043
1895
14
37
9090
5365
20
13
2427
3971
38
8049
5439
18
3526
4357
39
8077
6401
5062
4339
40
8017
7210
16
4777
4897
41
6869
7286
17
5868
4904
42
7056
8348
6534
5641
43
7747
8260
49
19
5481
6004
44
8457
8991
21
4592
4603
45
9460
8311
2486
5999
46
9062
7639
3299
6018
47
9319
6799
23
3573
6213
48
10631
6472
57
4741
6434
10685
5528
25
5375
8643
50
10643
4472
生活区
工业区
山区
交通区
公园绿地区
下表是测量所得50个样本点的Cd浓度
表二为前50个样本点的Cd的浓度：
附件2.8种主要重金属元素的浓度
Cd（ng/g）
153.80
583.40
146.20
366.40
439.20
323.90
223.90
424.50
525.20
630.00
1092.90
635.30
269.80
463.40
1066.20
532.00
1123.90
778.70
267.10
754.80
201.40
396.30
287.00
687.80
193.70
526.00
359.50
449.10
516.40
852.70
1044.50
459.00
445.40
337.30
347.90
568.10
345.70
599.00
614.00
635.50
257.20
600.70
1213.50
567.60
325.80
228.50
212.10
568.60
90.50
214.70
根据数据可知，取样点之间的地理距离是无关的，在考虑重金属污染元素的分布式，只需把海拔换成浓度。
在编写程序中发现这不是一般的函数曲面，而是不规则采样点的集合。
调用linespace函数在采样点范围内生成非均匀分布的数
。
联合调用griddata函数将数据内插到图中。
程序清单
我们给出地形图和Cd元素的分析程序代码：
程序1：
A=xlsread（'
附件_数据.xls'
1,'
B4:
E53'
%先将数据放在电子表格中，再读取把第1页中数据读取,必须会!
!
x1=A（:
1）;
y1=A（:
2）;
z1=A（:
3）;
[X1,Y1,Z1]=griddata（x1,y1,z1,linspace（min（x1）,max（x1）,200）'
linspace（min（y1）,max（y1）,200）,'
v4'
%插值，v4是一种插值算法
surf（X1,Y1,Z1）;
shadinginterp;
%在flat的基础上进行色彩的插值处理，使色彩平滑过渡
含量分布三维曲面'
程序2：
B=xlsread（'
2,'
H4:
I53'
D=[AB];
x1=D（:
y1=D（:
z1=D（:
6）;
Cd含量分布三维曲面'
实现图形分析
第一个图是城市地形图即地形地貌，海拔越高越趋向于红。
在分析重金属污染源时，研究人员可以根据流体相关知识结合本图及当地降水量等信息，进行数据分析，以促进研究的准确性。
第二个图是Cd元素在这个城市分布的三位立体图，视觉越高颜色越深则此重金属元素含量越高。
第三个图是和第二个图相对性的二维，颜色越趋近于红色则所含浓度越高。
6.要在某山区方圆大约27平方公里范围内修建一条公路，从山脚出发经过一个居民区，再到达一个矿区。
横向纵向分别每隔400米测量一次，得到一些点的高程，下表为数据表，请利用matlab绘制该山区地貌图、等高线图和立体等高线图。
提示：
meshz（x,y,z）绘制地貌图
contour（x,y,z）绘制等高线图
contour3（x,y,z）绘制立体等高线图
y
x
12001600200024002800320036004000
1200
11301250128012301040900500700
13201490142014001300700900850
139015001500140090011001060950
15001200110013501450120011501010
15001200110015501600155013801070
15001550160015501600160016001550
1480150015501540145013001200980
1600
2000
2400
2800
3200
3600
x=1200:
400:
4000;
y=1200:
3600;
[xx,yy]=meshgrid（x,y）;
zz=[
1130,1250,1280,1230,1040,900,500,700;
1320,1450,1420,1400,1300,700,900,850;
1390,1500,1500,1400,900,1100,1060,950;
1500,1200,1100,1350,1450,1200,1150,1010;
1500,1200,1100,1550,1600,1550,1380,1070;
1500,1550,1600,1550,1600,1600,1600,1550;
1480,1500,1550,1510,1430,1300,1200,980;
];
surf（xx,yy,zz）;
holdon;
7.已知测量碎部点三维坐标如下表所示，根据下述程序经过内插绘制等高线图
表观测数据表
序号
X（10M）
Y（10M）
H（M）
5.14
9.56
8.66
12.55
10.51
2.70
22.03
22.97
32.94
55.58
24.55
参考程序：
fid=fopen（'
D:
\MyDocuments\MATLAB\等值线.txt'
data=fscanf（fid,'
%f%f%f'
[3,inf]）;
data=data'
;
fclose（fid）;
scale=1e+8;
x=data（:
y=data（:
z=data（:
scale;
xb=min（x）;
xe=max（x）;
%得到最小值
yb=min（y）;
ye=max（y）;
Xl=linspace（xb,xe,200）;
%沿x
Yl=linspace（yb,ye,200）;
%得到y点
[X,Y]=meshgrid（Xl,Yl）;
%得到值的点坐标矩阵
method='
cubic'
%插值方法
Z=griddata（x,y,z,X,Y,method）;
%曲面
figure;
%准备图形窗口
[C,h]=contour（X,Y,Z）;
%画等高
clabel（C,h）;
colormapcool;
%标示等
boxon;
axisequal;
等值线图:
（单位:
scale）'
%设置坐标轴和标题
holdoff;
%停止在当前图形窗口将网格
请绘制等高线图形：
（此题目写到报告中）
四、实验注意事项
（1）对没有程序的题目，分析实验内容，写出程序大致框架或完整的程序代码。
并写到实验报告中。
（2）进入MATLAB7.0集成环境。
（3）编辑或验证程序，并进行保存。
（4）运行程序，若有错误，修改错误后再次运行，如此反复进行到不显示出错为止。
（5）检查程序输出结果。
五、实验心得
本次上机实验是测量数据可视化，这很复杂其中包含了许多其他方面的东西。
自己在课前做了预习，课堂上听取了老师所讲的知识。
然后做实验报告时，一些简单的基础操作没有什么问题。
但是其中有一部分拟合与差值的部分不是很懂，加上实验中这块的内容很少，将数据放在电子表格中，再数据读取这一块也不是很懂，需要多加的练习，很快就做出来了。
实验成绩__________________________
实验时间和地点____________________
实验指导老师______________________
升级会员 