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

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

ID:38698416

大小:267.00 KB

页数:8页

时间:2019-06-17

实验八  测量程序设计应用(最新)_第1页
实验八  测量程序设计应用(最新)_第2页
实验八  测量程序设计应用(最新)_第3页
实验八  测量程序设计应用(最新)_第4页
实验八  测量程序设计应用(最新)_第5页
资源描述:

《实验八 测量程序设计应用(最新)》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库

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,

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。