PAGE
PAGE 1
坐标正反计算
实 验 报 告
所属课程名称 MATLAB编程与应用
实 验 地 点
实 验 日 期 2019.12.19
班 级
学 号
姓 名
指导老师
一、实验目的:
坐标计算主要包括坐标正算和坐标反算,坐标正反计算是测量计算的基本公式,通过编写相应函数来实现对标量和向量(或矩阵)输入参数的坐标正算及坐标反算。
二、实验内容:
【实验过程及成果】(程序说明、实验代码、实验数据、实验结果)
程序说明
坐标正算:用if循环判断strcat横向连接路径名pathname和文件名filename中文件是否为空函数isempty,uigetfile实现文件的打开对话框来调用需要的数据文件。Dlmread来读取已知点数据,读取已知点编号yzd_n、个数yzd_count、X、Y坐标yzd_X、yzd_Y、水平距离yzd_Dis、方位角yzd_fwj,最后根据正算公式计算出坐标X、Y。
坐标反算:x1,y1起点坐标,x2,y2终点坐标,dis距离,azi方位角,通过反算公式计算方位角和距离。
实验代码
>>function [X,Y]=ZBZS(pathname)%坐标正算代码
if(nargin<1)
[filename,pathname]=uigetfile(...
{'*.txt;*.dat;','data Files (*.txt,*.dat)';
'*.txt','yzd-files(*.txt)';'*.dat','survey data(*.dat)';
'*.*','All-file.s(*.*)';},'Pick a file');
if(isempty(pathname))
return;
end
pathname=strcat(pathname,filename);
end
yzd=dlmread(pathname);
yzd_n=yzd(:,1);
yzd_count=length(yzd_n);
yzd_X=yzd(:,2);
yzd_Y=yzd(:,3);
yzd_Dis=yzd(:,4);
yzd_fwj=yzd(:,5);
X=yzd_X+yzd_Dis.*cos(yzd_fwj);
Y=yzd_Y+yzd_Dis.*sin(yzd_fwj);
end
>>function [dis,azi]=ZBFS(x1,y1,x2,y2)%坐标反算函数
dx=x2-x1;
dy=y2-y1;
tx=atan2(dy,dx);
dis=sqrt(dx.^2+dy.^2);
if(tx<0)
azi=tx+2*pi;
else
azi=tx;
end
end
实验数据
>>正算:
1 3628.022 6183.864 371.61 81.5725
2 4354.543 3215.654 432.12 120.34532
>>反算:
x1=[1 1 1 1 1];y1=[1 2 3 4 5]
x2=[2 2 2 2 2];y2=[5 4 3 2 1]
实验结果
【实验小结】(收获体会)
在平面测量坐标系中,已知一点a的坐标(Xa,Ya)及该点到另一点b的水平距离(Dab)和方位角(Aab),求另一点的坐标(Xb,Yb)是坐标正算,已知两点a和b的坐标(Xa,Ya)、(Xb,Yb)反求这两点之间的方位角Aab和距离Dab称为坐标反算。学习此次实验学会了实现坐标正反计算,了解了坐标基础计算方法和相应的程序编写。
三、指导教师评语及成绩:
评 语
评语等级
优
良
中
及格
不及格
1.实验报告按时完成,字迹清楚,文字叙述流畅,逻辑性强
2.实验方案设计合理
3.实验过程(实验步骤详细,记录完整,数据合理,分析透彻)
4实验结论正确.
成 绩:
指导教师签名:
批阅日期:
- 下载文档
- 收藏
- 0