实验六 测量数据可视化实验
一、实验内容（练习以下题目，并写出程序运行结果）
要求：读懂每一程序，掌握有关命令或函数的用法，运行这些程序，观察运行的结果，最后给出实验小结或心得。
1．测量平差中偶然误差分布图形绘制
根据偶然误差Δ服从Δ～ N (0 ,σ2)，可以应用MATLAB绘制出均方差为σ= 1 ,σ= 2 的正态分布概率密度函数的误差分布曲线，具体为 x = - 4∶0.1∶4;
y1 = normpdf ( x ,0 ,1) ; plot（x , y1 , ′r′） hold on
y2 = normpdf ( x ,0 ,2) ; plot（x , y2 , ′b′） hold off 2.
测量平差中误差椭圆的绘制
平差后任一待定点坐标平差值可以计算出，而且还可计算出误差椭圆的三要素：长半轴E,短半轴F,极大值方向。
function ellipse(x,y,A,B,angle)
plot(y+A*sin(angle)*cos(0:pi/360:2*pi)+B*cos(angle)*sin(0:pi/360:2*pi),… x+A*cos(angle)*cos(0:pi/360:2*pi)-B*sin(angle)*sin(0:pi/360:2*pi));
axis('equal'); xlabel('x'); ylabel('y'); title('eclipse');
运行ellipse(400,400,100,50,pi/4) 3.
对某点进行位移观测，得到8个不同时间的位移值，见下表，试进行傅立叶变换，
并作出频谱图。
1
时间k 0 1 2 3 4 5 6 7 位移f(k) 2.3 3.0 3.8 4.2 5.1 4.5 3.2 1.8 f=[2.3 3 3.8 4.2 5.1 4.5 3.2 1.8]' F=fft(f) FM=abs(f) t1=(0:1:7) for i=1:8 t=0:0.001:FM(i) plot(t1(i),t,'b-') hold on end
title('频域分析') xlabel('频率值') ylabel('能量')
4. 要在某山区方圆大约27平方公里范围内修建一条公路，从山脚出发经过一
个居民区，再到达一个矿区。横向纵向分别每隔400米测量一次，得到一些点的高程，下表为数据表，请利用matlab绘制该山区地貌图、等高线图和立体等高线图。
提示：meshz(x,y,z)绘制地貌图
contour(x,y,z)绘制等高线图
contour3(x,y,z)绘制立体等高线图 y 1200 1600 2000 2400 2800 3200 3600 y 1200 1600 x 1200 1600 2000 2400 2800 3200 3600 4000 1130 1250 1280 1230 1040 900 500 700 1320 1490 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 1540 1450 1300 1200 980 x 1200 1600 2000 2400 2800 3200 3600 4000 1130 1250 1280 1230 1040 900 500 700 1320 1490 1420 1400 1300 700 900 850 2
2000 2400 2800 3200 3600 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 1540 1450 1300 1200 980 x=1200:400:4000; y=1200:400:3600;
[xx,yy]=meshgrid(x,y); %xx,yy均为7行8列的矩阵 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); %地貌图 hold on;
contour(xx,yy,zz); %等高线图 hold off
contour3(xx,yy,zz); %立体等高线图
[c,h]=contour(xx,yy,zz); clabel(c,h);
5.已知测量碎部点三维坐标如下表所示，根据下述程序经过内插绘制等高线图
表 观测数据表
序号 1 2 3 4 5 6 7 8 9 10 11
X(10M) 1 1 1 2 2 2 3 3 3 4 4
Y(10M) 1 2 3 1 2 3 1 2 3 2 3
H(M) 5.14 9.56 8.66 12.55 10.51 2.70 22.03 22.97 32.94 55.58 24.55
3