使用夏普计算机导线平差程序的编写及应用.doc

SHARPPC-E500S袖珍计算机编写导线平差程序及应用沈重庆江苏东南交通工程咨询监理有限公司，南京，210018 摘要 目前进行导线平差计算的手段较多，但笔者通过工程实践认为使用PC-E500S袖珍计算机编程进行导线平差计算，既快捷准确，又便于携带，尤其适合野外作业，可以随时计算出结果，并可直接打印。关键词 SHARPPC-E500S袖珍计算机编程 快捷 便携 适合野外作业1 序言在公路等建设领域，无论是勘察设计阶段还是施工阶段，测量工作自始至终起到非常重要的作用，而导线测量成果，是整个测量工作成败的关键。笔者曾参加过高速公路、一级公路、二级公路及市政道路桥梁的建设工作，并通过工程实践的摸索，使用SHARPPC-E500S袖珍计算机，采用QUIKBASIC语言，编写了一套测量程序，可进行复杂的匝道中、边线坐标计算、高程计算、导线平差计算等。通过实践证明，该套程序切实可靠，方便使用，输出的结果清晰易懂。尤其在蚌埠南京（蚌明段）高速公路建设中得到了成功运用。下面具体谈一下导线平差程序的编程及应用。2 导线平差程序的编程2.1第4046为数据输入区40：DATA 8，263.0941，263.094142： DATA 295.102460，262.4814，167.49475，163.54235，233.05536，313.39126，184.08035，179.2411544：DATA 406.5335，572.3271，244.8648，308.8288，278.9507，592.6529，422.7168，046： DATA 52938.656，72822.028，52992.937，73273.9642.2以下为主程序区55：PRINT “DAO XIAN PING CHA”60：READ N，V，Z64：DIM A（N），D（N），D（X），D（Y），DD（X），DD（Y），FX（N），FY（N）66：SUMA=0，AA=DEGV68：FOR I=1 TO N70：READ A（I）72：A（I）=DEG A（I）74：AL=AL+A（I）76：IF AL360 THEN AL=AL-36078：SUMA=SUMA+A（I）80：NEXT I88：BB=AL-DEG Z90：BF=-BB/N92：BF=-BB/N94：PRINT“fa=”;INT(360000*BB)/100; “s” 96： PRINT ABS（360000*BB）10*SQR N THEN PRINT “fa10*SQR N=”；INT（100*SQR N）/10；“ GOOD”98： PRINT ABS（360000*BB）10*SQR N THEN PRINT “fa10*SQR N=”；INT（100*SQR N）/10；“ BAD”100：PRINT “SUM A=”；DMS SUMA；“（DMS）”；：PRINT“ Bf=”；INT（10*DMS（BF*3600）/10；“s”102：FOR I=1 TO N104：A（I）=A（I）+BF106：AA=AA+A（I）-180108：IF AA360 THEN AA=AA-360110：IF AA0 THEN AA=AA+360112：PRINT“AA（”；I；“）=”；INT（10000*DMS AA+0.5）/10000：“（DMS）”114：READ D（I）122：DX（I）=D（I）*COS AA124：DY（I）=D（I）*SIN AA130：NEXT I132：DD=0134：FOR I=1 TO N136：DD=D（I）+DD138：NEXT I140：PRINT“SUM D=”；DD；“mm”142：AVER=DD/（N-1）144：SDX=0；SDY=0146：READ O，P，Q，R148：FOR I=1 TO N150：SDX=SDX+DX（I）152：SDY=SDY+DY（I）154：NEXT I156：FDX=SDX+（Q-O）158：FDY=SDY+（R-P）160：DS=SQR（FDX2+ FDY2）162：KK=DD/DS：DS=INT（10000*DS）/10164：PRINT “fx=；”INT（10000*FDX）/10；“mm”；：PRINT “fy=”；INT（10000*FDY）/10；“mm”166：PRINT“K=1/”；INT（10*KK+0.5）/10168：PRINT“AVER D=”；AVER;“m”170：XX=O，YY=P172：FOR I=1 TO N-1174：FX（I）=-（FDX/DD）*D（I）176：FY（I）=-（FDY/DD）*D（I）178：XX=XX+DX（I）+FX（I）180：YY=YY+DY（I）+FY（I）182：XX= INT（10000*XX+0.5）/1000184：YY= INT（10000*YY+0.5）/1000186：PRINT “X（”；I；“）=”；XX；：PRINT “Y（”；I“）=”；YY190：NEXT I200：END3 程序的注释该程序只需一次性依次输入导线测量的有关参数：测站数 、起始边方位角、终边方位角、各个测角（左角）、各个测边距离、及起点、终点坐标，即可自动计算出有关结果：角度闭和差（fa）、按一级导线精度10N相比较的偏差结果（GOOD或BAD）、测角总和（SUM A）和角度分配值（Bf）、经分配后的各边方位角（A（I）、纵横坐标增量闭和差（fx、fy）、导线增量闭和差（Ds）、相对闭和差（K）、测边平均值（AVER D），最终输出平差后的坐标值（X（I）、Y（I）。如需打印结果，只需输入有关命令连接热敏打印机或EPSON1600-K等打印机即可输出结果。4 程序的应用实例4.1 数据整理现以蚌明高速公路十三标为例，介绍该程序的使用情况，该标段地处蚌埠市北大门，是一座大型互通立交，工程规模1.04亿元人民币，主线路基长1.349，主线桥上跨津浦铁路货运线，长1.436,有十一条匝道（四座匝道桥及七条匝道路基组成）。闭合导线以BNCBN1为基准边，按四等导线精度控制，使用TOPCON711全站仪观测，采用六测回数进行读数，数据经整理如下：基准点数据：基准点坐标值方位角XY经计算：起始边：=263-09-04.1终止边：=263-09-04.1因为是闭合导线，所以方位角相同BNC52992.937073273.9640BN152938.656072822.0280测量数据：测站名点号夹角（左角）测距（米）BN1BNC295-10-24.6455.1880AN8406.5335AN8BN1262-48-14AN6572.3271AN6AN8167-49-47.5AN5244.8648AN5AN6163-54-23.5AN4308.8288AN4AN5233-05-536AN3278.9507AN3AN4313-39-126BNB592.6529BNBAN3184-08-035BNC422.7168BNCBNB179-24-115BN1455.1880导线示意图如下：4.2数据的输入：第4046行为数据输入行，在BASIC状态下，第40行依次输入测站数、起始方位角、终止边方位角，第42行依次输入各测站夹角，第44行依次输入各测边的距离，最后一条边输0，第46行依次输入起点坐标及终点坐标，这样就完成了数据的输入过程。如果需要计算另一条导线，只需将相关数据从第40行重新输入一次即可。4.3程序的执行：将计算机设置在RUN状态下，输入RUN55，然后连续按下回车键，就可以依次计算所需要的结果。4.4计算成果汇总：十三标段导线计算结果汇总表点号实测坐标（m）原坐标（m）差值（mm）XYXYXYBN152938.6560 72822.0280 52938.6560 72822.0280 00AN853324.5821 72949.8362 53324.5816 72949.8450 0.5-8.8AN653214.1341 73511.3979 53214.1300 73511.3960 4.11.9AN5 53218.594473756.219953218.5920 73756.2070 2.412.9AN453309.596374051.3355 53309.5850 74051.3430 11.3-7.5AN353145.7906 74277.1196 临时点BNB53038.9620 73694.1709 53038.9550 73694.1680 7.02.9BNC52992.937073273.9640 52992.9370 73273.9640 00Fa=10.858=14.1 ，Fx=-48.5mm ，Fy=25.7mm，Ds=54.8mm ， K=1/515101/35000。 精度达到四等导线要求（闭和导线） 5 结束语 使用该程序进行计算，操作简单快捷，从数据的输入到结果的输出，仅需要约五分钟时间；如果采用手算，需要约六十分钟，且容易出错；如果采用其他电脑软件计算，就必须在办公室进行，不能即时在现场得出结果。而使用SHARPPC-E500S袖珍计算机编写的程序，