资源描述:
《实验八 测量程序设计应用(最新)》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、实验八测量程序设计应用初步一、实验目的和要求:1.初步掌握利用matlab软件进行测量程序设计的思路和技巧;2.初步掌握matlab软件在数字测量技术计算中应用;3.初步掌握利用matlab软件进行测量数据处理的思想;二、实验内容:1.进行坐标正、反算:(1)A(308.675,456.98),S=1002.12,delta=207°59′53″程序:ax=input('请输入A的纵坐标:')ay=input('请输入A点的横坐标:')s=input('请输入AB长度:')ad=input('请输入AB的坐标方位角度:')am=input('请输入AB的坐标方位角分:')as=inp
2、ut('请输入AB的坐标方位角秒:')a=ad+am/60+as/3600;a=a*pi/180;bx=ax+s*cos(a)by=ay+s*sin(a)结果:bx=-576.1604by=-13.4568(2)已知导线网中五个点的坐标XX2507.692299.822192.432179.722166.72Y1215.631303.801383.981556.411757.29编程计算每条导线的长度、坐标方位角,并绘出导线图形(要求导线点用红色圆圈表示)程序:x=[2507.692299.822192.432179.722166.72]y=[1215.631303.801383.9
3、81556.411757.29]plot(x,y)holdonfori=1:5plot(x(i),y(i),'o')endfori=1:4 a(i)=x(i+1)-x(i); b(i)=y(i+1)-y(i); s(i)=sqrt(a(i)^2+b(i)^2);endfori=1:4ifa(i)~=0 m(i)=atand(abs(b(i)/a(i))); k(i)=m(i); ifa(i)>0&b(i)>0 c(i)=k(i); elseifa(i)<0&b(i)>0 c(i)=180-k(i); elseifa(i)<0&b
4、(i)<0 c(i)=180+k(i); elc(i)=360-k(i); endse elseifb(i)>0c(i)=90;elsec(i)=270;endendd=fix(c)e=fix((c-d)*60)f=((c-d)*60-e)*60g=[d'e'f']运行结果:2简单的地形图绘制:地形图绘制(有11个碎部点坐标)115.14129.56138.662112.552210.51232.703122.033222.973332.944255.584324.55要求:1)建立.txt文件,用于存储数据;2)利用函数fopen()和fscanf(
5、)对数据读取程序:fid=fopen('F:workdengzhixian.txt','r');data=fscanf(fid,'%f%f%f',[3,inf]);data=data';fclose(fid);x=data(:,1);y=data(:,2);z=data(:,3);xb=min(x);xe=max(x);yb=min(y);ye=max(y);Xl=linspace(xb,xe,200);Yl=linspace(yb,ye,200);[X,Y]=meshgrid(Xl,Yl);method='cubic';Z=griddata(x,y,z,X,Y,method);
6、figure;holdon;contour(X,Y,Z);colormapcool;boxon;xlabel('x');ylabel('y');axisequal;title('等值线图:(单位:scale)');holdoff绘图结果:3.进行精密平差计算:A和B是已知高程的水准点,C、D、E是待定点,A、B高程、观测高差和相应的水准路线长度见下表,试按间接平差求各待定点的高程平差值。线路编号观测高差(m)线路长度(km)已知高程(m)1234567+1.359+2.009+0.363+1.012+0.657+0.238-0.5951.11.72.32.72.41.42.6HA=5
7、.016HB=6.016要求:利用函数式M文件进行编程。disp('水准网间接平差示例')disp('已知高程')Ha=5.015Hb=6.016disp('观测高差,单位m')L=[1.359;2.009;0.363;1.012;0.657;-0.357]disp('系数矩阵B')B=[10;01;10;01;-11;-10]l=[0;0;4;3;7;2];disp('C是单位权观测高差的线路公里数,S是线路长度')C=ones(1,6);S=[1.1,