实验六测量数据可视化.doc

实验六 测量数据可视化 一、实验目的 1、掌握用plot函数和fplot函数绘制曲线的方法 2、通过练习熟悉三维曲线和曲面图的绘制方法 3、掌握测量误差曲线和二维地形图等绘制方法 二、实验环境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个不同时间的位移值，见下表，试进行傅立叶变换，并作出频谱图。 时间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(能量) 5. 要在某山区方圆大约27平方公里范围内修建一条公路，从山脚出发经过一个居民区，再到达一个矿区。横向纵向分别每隔400米测量一次，得到一些点的高程，下表为数据表，请利用matlab绘制该山区地貌图、等高线图和立体等高线图。 提示：meshz(x,y,z)绘制地貌图 contour(x,y,z)绘制等高线图 contour3(x,y,z)绘制立体等高线图 y x 1200 1600 2000 2400 2800 3200 3600 4000 1200 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 1600 2000 2400 2800 3200 3600 x=1200:400:4000;y=1200:400: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