Excel环境中控制测量三角网平差及坐标计算

  • 格式:pdf
  • 大小:290.63 KB
  • 文档页数:6

收稿日期:2003󰀁03󰀁25作者简介:余梁蜀(1964󰀁),女,四川重庆人,副教授,长期从事工程建设管理与工程测量的教学与科研工作.第37卷第1期2004年2月武汉大学学报(工学版)EngineeringJournalofWuhanUniversityVol.37No.1Feb.2004

文章编号:1671󰀁8844(2004)01󰀁063󰀁05

Excel环境中控制测量三角网平差及坐标计算

余梁蜀1,张宗孝2,马󰀁斌2,孟文杰2

(1.西安交通大学,陕西西安󰀁710049;2.西安理工大学,陕西西安󰀁710048)

摘要:将Excel的表格特性和测量计算传统的列表格式结合起来进行控制测量三角网的平差和坐标计算,有效地克服了Fortran或Basic语言编程后台计算、外业资料和已知条件的输入不直观、出现错误比较难检查的缺陷.充分地利用Excel强大的功能和良好的人机交互对话界面,方便地进行数据处理和函数运算,迅速的制作复杂图表的特性,将整个计算过程展示在前台,易操作、易修改、透明度高,是控制测量计算的有效方法之一.关键词:电子表格;控制测量;三角点;坐标中图分类号:TB22󰀁󰀁󰀁文献标识码:A

Areatriangulationadjustmentcorrectionandcoordinatecalculation

ofcontrolsurveyingbasedonExcelplatform

YULiang󰀁shu1,ZHANGZong󰀁xiao2,MABin2,MENGWen󰀁jie2

(1.Xi󰀂anJiaotongUniversity,Xi󰀂an710049,China;2.Xi󰀂anUniversityofTechnology,Xi󰀂an710048,China)

Abstract:Amethodforareatriangulationadjustmentcorrectionandcoordinatecalculationofcontrolsurvey󰀁ingwithcombiningelectronicformofExcelandcalculationformoftraditionisputforward.Thismethod

overcomestheshortcomingsofcalculationprocedureusingtheFortranandBasic,suchascomputationin

background,inputtingfilesindirectly,findingmistakesdifficultly,etc..ThismethodbasedonExcelhasthe

followingadvantages:thecomputationprocessisofcompleteman󰀁machineinteractionanddialogueprocess;convenientindataprocessing;quicklytomakesophisticateddiagramandform;computationonthestage;

convenientinoperation;convenientlytomodify;higherintransparency.Sothemethodisoneofefficient

methodsofcalculationoftriangulationcontrolsurveying.Keywords:electronicform;controlsurveying;trigpoint;coordinate

󰀁󰀁当测量区域为带状时,平面控制测量通常采用

三角网控制.三角网的平差及坐标计算,工作量

大,程序繁琐,计算精度要求高;特别是当控制网中

三角点数较多时,计算更加复杂,如果借助计算器

手算,则计算时间长,且容易出错;如果用高级语言

如Fortran或Basic编程计算,虽然计算比手算快,

但由于计算是在后台进行,外业资料和已知条件的

输入不直观,出现错误比较难检查.MicrosoftEx󰀁

cel是微软公司在Windows操作系统中开发的一

种通用软件[1],它具有强大的功能和良好的人机交互对话界面,可以方便地进行数据处理和函数运

算,而且能够迅速地制作复杂的图表,计算过程展

示在前台,易操作、易修改、透明度高,在工程测量

中利用它进行三角形网的计算非常方便[2].

1󰀁三角网计算模型和程序步骤

本文针对两端有基线的单三角锁的平差及坐

标计算编制Excel程序,中点多边形、线形三角锁

等其他形式的三角网可仿照编制.三角网计算的

目的是根据多余条件进行平差计算、推算出各三角点的坐标值.

外业测量工作完成后,获得三角网中各角的角

值ai,bi,ci,两端基线的边长d0,dn,一条边的方位角󰀁AB和一个点的坐标XA,YA,并绘制计算简图

(计算简图、外业资料和起算数据列在界面1中).

两端有基线的单三角锁的平差及坐标计算程

序为(1)按三角形图形条件对角值ai,bi,ci进行第

一次平差改正,求第一次改正后的角值a i,b i,c i.

a i=ai-f󰀂i/3,b i=bi-f󰀂i/3,c i=ci-f󰀂i/3

(1)

式中:f󰀂i为三角形角度闭合差,f󰀂i=ai+bi+ci-

180!.(2)按三角锁基线条件对角值ai,bi,ci进行第

二次平差改正,求第二次改正后的角值a∀i,b∀i,c∀i.

a∀i=a i-v a,b∀i=b i+v b,c∀i=c i(2)式中:v a,v b为角度的第二次改正数,v a=-v b=

- ∀wdn(#ctga+#ctgb);w为基线闭合差,w

=d0∃sina i∃sinb i-dn; ∀为弧度与角度换

算系数206265∀.(3)计算三角形各边的边长dij.自d0开始用a∀i,

b∀i,c∀i角值,及正弦定律推算各三角形边长,推算出终

了基线的边长应与dn值相等,进行计算检核.

(4)推算各边的方位角󰀁ij.利用起算基线的

方位角󰀁0和各三角形的角值a∀i,b∀i,c∀i依次推算各

边的方位角.(5)计算各边的坐标增量!Xij,!Yij和三角点

的坐标Xi,Yi.计算式为!Xij=dijcos󰀁ij,!Yij=dijsin󰀁ij(3)

Xj=Xi+!Xij,󰀁Yj=Yi+!Yij(4)󰀁󰀁三角锁平差及各三角点坐标计算程序见图1.

图1󰀁三角锁平差及坐标计算程序

󰀁󰀁以上计算相当繁杂,通过Excel列表和编写语

句,可使计算过程简单明了,计算线路清晰.

2󰀁三角网Excel计算程序编制

本文结合实例介绍三角网平差坐标计算Excel

程序编制方法.(1)外业资料和起算数据的输入.

打开计算机,进入Excel界面(见界面1)Sheet1工作簿中(重命名为%平差表%),在相应单元

格(Excel地址)输入各三角形编号、各三角点点名

和各三角形角号、外业资料(ai,bi,ci,d0,dn)和起算数据(󰀁AB,XA,YA)等.例如可分别在D6,E6,F6单元格中输入第&个三角形a1角值用度、分、

秒值,同理输入各三角形各角值.测量外业观测时

各内角值用度、分、秒表示,角值输入时也以度、分、秒分别占用不同单元格,而计算机计算时角度是以

度运算.因此,首先将各角值化为度的形式.击活G6单元格,在公式编辑栏编辑语句G6=D6+E6/60+F6/3600,确定后即完成将度分秒形式表示的

角值a1化为度,并存放入G6单元格中.随后击活G6,将鼠标移到其右下角%自动填充柄%黑十处,向

下拖动鼠标复制G6语句,则将第&个三角形中c1,b1的角值也化为度的形式表达,并分别存入

G7、G8单元格中.利用相对地址的原理[2],给G10单元格复制G6语句,将第∋个三角形中b2的角值化为度的形式,再击活G10,利用%自动填充柄%给G11、G12复制G10语句,将第∋个三角形中c2,b2的角值也化为度的形式.其他各三角形角度的转化与此相同.

(2)计算三角形的内角和及第一次改正数f󰀂1.

利用Excel的求和函数%Sum%就能方便地求出三角形内角和.在G9单元格中编写语句Sum(G6(G8),即完成第&个三角形内角求和.该语句的含

义是G6,G8单元格中的角度值求和,对于其他各

三角形内角求和采用复制G9语句完成.在单元格

H6中,编辑公式:-($G$9-180)/3,即对第&64武汉大学学报(工学版)200465󰀁第1期余梁蜀等:Excel环境中控制测量三角网平差及坐标计算个三角形角度闭合差(f󰀂i=ai+bi+ci-180!)按反

号平均平差的原则求得角值的改正数∀11(∀ij为三

角形内角的第一次改正数,i为三角形编号、j为各

三角形内角编号,例∀13表示第&个三角形、第3个

角的改正数).然后击活H6,用%自动填充柄%复制H6语句,求得∀12,∀13.击活H6,单击%复制%按

钮;击活H10,单击%粘贴%按钮,则得到语句-

($G$9-180)/3,在公式编辑栏,移动光标改$G$9中行号9为13,即可得到∀21,用%自动填充柄%

复制H10语句,求得∀22,∀23.其余的∀ij求法操作

同上.

(3)计算第一次改正后的角值a i,b i,c i.将原角值加上∀ij即为第一次改正后角值,在

I6中,编辑公式I6=G6+H6,即完成计算第一次

改正后的角值,其各角值的计算仍然采用复制I6

语句的操作.(4)计算第二次改正数v a,v b.

计算模型为式(2),首先需计算a i,b i,c i角的

正弦值和余切值.在J6中编制语句%sin(I6)#/

180)%.在K6中编制语句%(I6)#/180)/tan%.完成第&个三角形中角a 1正弦值和余切值的计算.

通过复制J6、K6的语句,完成其他各角正弦值和

余切值的计算.在K26中编辑语句:%G26)J6)J10)~)J22/J8/J12/~J24-I26%计算基线闭合

差;在K27中编辑语句:%-K26/J26/Sum(K6,

K8,~,K24))180/#%计算角值第二次改正值v a=-v b.并分别在L6、L8中编辑引用语句%=K27%、%=-K27%,将第二次改正值存放入L6,L8,

L10∗∗

(5)第二次改正后的角值a∀i,b∀i,c∀i.

由式(2)知,第二次改正后的角值a∀等于第一

次改正后角值与第二次改正数之和.在M6中编辑语句%I6+L6%,即完成计算,其余各角只要采用

复制M6语句即可.同时复制G9的求和函数,计

算各三角形内角和应等于理论值180!,进行计算检核.

(6)各边边长.

运用正弦定理,计算三角形各边边长.首先运

用%复制%、%粘贴%和%自动填充柄%给N6,N10∗∗及N7,N8,N11,N12∗∗单元格复制J6语句,计

算第二次改正后角值的正弦函数值.在O8中编

辑%=G26%语句,即连通引用存放在G26中的起算边边长值,在O6,O7中编辑公式%N6/N8)O8%、%

N7/N8)O8%,即完成第一个三角形待求边d12和

间隔边s13长度的计算.在O12中编辑%=O6%语句,即连通引用已计算出存放在O6中的第一个三

角形待求边d12值(等于第二个三角形的已知边

d21),在O10,O11中复制O6,O7语句,即完成第二