实验八 测量程序设计应用(最新).docVIP

实验八 测量程序设计应用初步 实验目的和要求: 初步掌握利用matlab软件进行测量程序设计的思路和技巧; 初步掌握matlab软件在数字测量技术计算中应用; 初步掌握利用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=input(请输入AB的坐标方位角秒:) a=ad+am/60+as/3600; a=a*pi/180; bx = ax + s * cos(a) by = ay + s * sin(a) 结果：bx = -576.1604 by = -13.4568 (2) 已知导线网中五个点的坐标 X X 2507.69 2299.82 2192.43 2179.72 2166.72 Y 1215.63 1303.80 1383.98 1556.41 1757.29 编程计算每条导线的长度、坐标方位角，并绘出导线图形（要求导线点用红色圆圈表示） 程序： x=[2507.69 2299.82 2192.43 2179.72 2166.72] y=[1215.63 1303.80 1383.98 1556.41 1757.29] plot(x,y) hold on for i=1:5 plot(x(i),y(i),o) end for i=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);end for i=1:4 if a(i)~=0??? m(i)=atand(abs(b(i)/a(i))); ??? k(i)=m(i); ??? if a(i)0 b(i)0??????? c(i)=k(i);??? elseif a(i)0 b(i)0??????? c(i)=180-k(i);??? elseif a(i)0 b(i)0??????? c(i)=180+k(i);??? el c(i)=360-k(i);??? end se??????? elseif b(i)0 c(i)=90; else c(i)=270; end end d=fix(c) e=fix((c-d)*60) f=((c-d)*60-e)*60 g=[d e f] 运行结果： 2简单的地形图绘制：地形图绘制（有11个碎部点坐标） 1 1 5.14 1 2 9.56 1 3 8.66 2 1 12.55 2 2 10.51 2 3 2.70 3 1 22.03 3 2 22.97 3 3 32.94 4 2 55.58 4 3 24.55 要求：1）建立.txt文件，用于存储数据； 2）利用函数fopen( )和fscanf( )对数据读取 程序： fid=fopen(F:\work\dengzhixian.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); figure; hold on; contour(X,Y,Z); colormap cool; box on; xlabel(x);ylabel(y); axis equal; title(等值线图:(单位:scale)); hold off 绘图结果： 3.进行精密平差计算： A和B是已知高程的水准点，C、D、E是待定点，A、B高程、观测高差和相应