contour3(x,y,z)绘制立体等高线图
11
y 1200 1600 2000 2400 2800 3200 3600 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: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); hold on;
12
7.已知测量碎部点三维坐标如下表所示，根据下述程序经过内插绘制等高线图
表 观测数据表
序号 X(10M) Y(10M) 1 1 1 2 1 2 3 1 3 4 2 1 5 2 2 6 2 3 7 3 1 8 3 2 9 3 3 10 4 2 11
4
3
参考程序：
fid=fopen('D:\\My Documents\\MATLAB\\等值线.txt','r'); data = fscanf(fid,'%f %f %f',[3,inf]); data=data'; fclose(fid); scale=1e+8;
x =data(:,1);y =data(:,2);z =data(:,3); 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;hold on;%准备图形窗口 [C,h]=contour(X,Y,Z);%画等高
13
H(M) 5.14 9.56 8.66 12.55 10.51 2.70 22.03 22.97 32.94 55.58 24.55
clabel(C,h);colormap cool;%标示等 box on;xlabel('x');ylabel('y');
axis equal;title('等值线图:(单位:scale)'); %设置坐标轴和标题
hold off;%停止在当前图形窗口将网格
请绘制等高线图形：（此题目写到报告中）
四、实验注意事项
（1）对没有程序的题目，分析实验内容，写出程序大致框架或完整的程序代码。并写到实验报告中。
（2）进入MATLAB7.0集成环境。 （3）编辑或验证程序，并进行保存。
（4）运行程序，若有错误，修改错误后再次运行，如此反复进行到不显示出错为止。 （5）检查程序输出结果。
五、实验心得
本次上机实验是测量数据可视化，这很复杂其中包含了许多其他方面的东西。自己在
课前做了预习，课堂上听取了老师所讲的知识。然后做实验报告时，一些简单的基础操作没有什么问题。但是其中有一部分拟合与差值的部分不是很懂，加上实验中这块的内容很少，将数据放在电子表格中，再数据读取这一块也不是很懂，需要多加的练习，很快就做出来了。
实验成绩__________________________
实验时间和地点____________________
实验指导老师______________________
14