单一附合导线测量控制网平差程序(共13页).docx
精选优质文档-倾情为你奉上单一符合导线测量控制网平差程序梁彩艳(西北石油局工程服务大队,新疆乌鲁木齐)摘要:由于全站仪(光电仪) 的不断普及,导线布网已在测量控制网中广泛应用。为了简化计算,提高平差精度,节约平差经费,结合实际工作中的一个实例,运用PC E500 计算器编制了单一符合导线高斯平差程序。关键词:单一;符合导线;平差程序1 差原理及程序框图(见图1)2 实例在单一符合导线(如图2) 中,已知各起始数据,观测了各折角方向值,测量了各边长(见表1) ,观测方向中误差r = ±5,边长丈量中误差si = ±0. 5 simm ,求各导线点的坐标平差值和最弱点4 号点的点位精度(可假设单位权中误差0 = ±5,则观测角权为0. 5 ,观测边权Psi = 100/ Si) 。图2 单一符合导线示意图表1 测量数据一览表点号已知坐标(m)X Y编号已知坐标方位角B 3020. 348 - 9049. 801 AB 226°4459C 3702. 437 - 10133. 399 CD 57°5931编号观测角度编号观测边长(m)1 230°3237 1 204. 9522 180°0042 2 200. 1303 170°3922 3 345. 1534 236°4837 4 278. 0595 192°1425 5 451. 6926 260°59013 平差程序10 :REM SUB 10 - 190 QIU DAO XIAN DIAN GAI LUEZUO BIAO20 :REM SHU ZU”A”WEI FANG WEI J IAO ,”B”WEIGUAN CE ZUO J IAO ,”S”WEI BIAN CHANG,”P”WEIQUAN ,”V”WEI GAI ZHENG SHU ,”MX、MY、M”WEIDAI QIU DIAN ZHONG WU CHA ,”X、Y”WEI ZUOBIAO30 : INPUT”GUAN CE J IAO GE SHU :N = ”;N40 :DIM A (N) ,B (N) ,S (N) ,X(2 3 N) , Y(2 3 N) , P (2 3N) ,V(2 3 N) ,MX(N) ,MY(N) ,M(N)45 :FOR I = 1 TO N50 :READ B( I) ,S( I)60 :B( I) = DEG (B( I) )70 :NEXT I80 : INPUT”QI SHI FANG WEI J IAO :A0 = ”;A(0)90 :A(0) = DEG (A(0) )100 :FOR I = 1 TO N110 :A( I) = A( I - 1) + B( I) + 180120 : IF A( I) > 360 THEN A( I) = A( I) - 360130 :NEXT I140 : INPUT”QI SHI DIAN ZUO BIAO : X0 = ”; X(1) ,”Y0= ”; Y(1)总第93 期2004 年第2 期西部探矿工程WEST - CHINA EXPLORATION ENGINEERINGseries No. 93Feb. 2004145 : PRINT ”DAO XIAN DIAN GAI LUE ZUO BIAO ZHI :”150 :FOR I = 2 TO N160 :X( I) = X( I - 1) + S( I - 1) 3 COS(A( I - 1) ) : Y( I) = Y( I- 1) + S( I - 1) 3 SIN(A( I - 1) )170 :PRINT ”X(”; I ;”) = ”;X( I) ,”Y( ;”I ;”) = ”; Y( I)180 :NEXT I185 :LOCATE 12190 :PRINT ”# # # # # ”200 :REM SUB 200 - 370 QIU TIAO J IAN FANG CHENGXI SHU ”A1 ,BI ,C1”SHU ZU J I QUAN HAN SHU XISHU ”FX ,FY”SHU ZU210 :DIM A1 (2 3 N - 1) ,B1 (2 3 N - 1) ,C1 (2 3 N - 1) ,FX(23 N - 1) ,FY(2 3 N - 1)220 :FOR I = 1 TO N230 :READ A1( I)240 :NEXT I250 :FOR I = 1 TO 2 3 N - 1260 :B1 ( I) = ( - Y(N) + Y( I) ) / 206. 26270 :C1 ( I) = (X(N) - X( I) ) / 206. 26280 : IF N < I THEN B1( I) = COS(A( I - N) :C1 ( I) = SIN (A( I - N) )290 :NEXT I300 : INPUT ”ZUI RUO DIAN HAO :FN = ”;FN310 :FOR I = 1 TO 2 3 N - 1320 :FX( I) = ( - Y(FN) + Y( I) ) / 206. 26330 :FY( I) = (X(FN) - X( I) ) / 206. 26340 : IF FN < I THEN FX( I) = 0 :FY( I) = 0350 : IF N < I THEN FX( I) = COS(A ( I - N) ) : FY( I) = SIN(A( I - N) )360 : IF I > N + FN - 1 THEN FX( I) = 0 :FY( I) = 0370 :NEXT I414 : INPUT ”FU HE FANG WEI J IAO :A2 = ”;A2 ,”FU HEDIAN ZUO BIAO :X2 = ”;X2 ,”Y2 = ”; Y2415 :WA = - DEG(A2) + A(N) :WX = - X2 + X(N) :WY =- Y2 + Y(N)416 :WA = DMS(WA) : IF WA > 360 THEN WA = WA - 360419 :WA = WA 3 10000 :WX = WX3 1000 :WY= WY3 1000420 :REM SUB 420 - 740 J IE SUAN FA FANG CHENG425 :REM SUB 430 - 470 J I SUAN GUAN CE J IAO、BIANDE QUAN430 : INPUT ”J IA SHE DAN WEI QUAN ZHONG WUCHA :D0 = ”;D0440 :FOR I = 1 TO 2 3 N - 1450 :P( I) = 0. 5460 : IF I > N THEN P( I) = (D0 3 D0) / (0. 25 3 S( I - N) )470 :NEXT I500 :REM ”E”WEI DAI QUAN QIU HE FU HAO520 : EAA = 0 : EAB = 0 : EAC = 0 : EAFX = 0 : EAFY= 0530 : EBB = 0 : EBC = 0 : EBFX = 0 : EBFY= 0540 : ECC = 0 : ECFX = 0 : ECFY= 0 : EFFX = 0 : EFFY= 0550 :FOR I = 1 TO 2 3 N - 1560 :EAA = A1 ( I) 3 A1 ( I) / P ( I) + EAA : EAB = A1 ( I) 3 B1( I) / P( I) + EAB570 :EAC = A1 ( I) 3 B1 ( I) / P( I) + EAC: EAFX = A1 ( I) 3 FX( I) / P( I) + EAFX: EAFY = A1 ( I) 3 FY( I) / P ( I) +EAFY580 :EBB = B1 ( I) 3 B1 ( I) / P ( I) + EBB : EBC = B1 ( I) 3 C1( I) / P( I) + EBC590 :EBFX = B1 ( I) 3 FX( I) / P ( I) + EBFX: EBFY= B1 ( I) 3FY( I) / P( I) + EBFY600 :ECC = C1 ( I) 3 C1 ( I) / P ( I) + ECC: ECFX = C1 ( I) 3 FX( I) / P( I) + ECFX: ECFY= C1( I) 3 FY( I) / P( I) + ECFY610 :EFFX = FX( I) 3 FX( I) / P( I) + EFFX: EFFY= FY( I) 3FY( I) / P( I) + EFFY620 :NEXT I630 :EBB1 = EBB - EAB 3 EAB/ EAA : EBC1 = EBC - EAB 3EAC/ EAA : EWX1 = WX - EAB 3 WA/ EAA640 :EBFX1 = EBFX - EAB 3 EAFX/ EAA : EBFY1 = EBFYEAB3 EAFY/ EAA650 :ECC2 = ECC - EAC3 EAC/ EAA - EBC1 3 EBC1/ EBB1660 :EWY2 = WY- EAC 3 WA/ EAA - EBC1 3 EWX1/ EBB1670 :ECFX2 = ECFX - EAC 3 EAFX/ EAA - EBC1 3 EBFX1/EBB1680 :ECFY2 = ECFY- EAC3 EAFY/ EAA - EBC1 3 EBFY1/EBB1690 : KC = - EWY2/ ECC2700 : KB = - EBC1 3 KC/ EBB1 - EWX1/ EBB1710 : KA = - EAB 3 KB/ EAA - EAC 3 KC/ EAA - WA/ EAA720 : EPVV = WA 3 WA/ EAA + EWX1 3 EWX1/ EBB1 +EWY2 3 EWY2/ ECC2730 : EFFX3 = EFFX - EAFX 3 EAFX/ EAA - EBFX1 3EBFX1/ EBB1 - ECFX2 3 ECFX2/ ECC2740 : EFFY3 = EFFY - EAFY 3 EAFY/ EAA - EBFY1 3EBFY1/ EBB1 - ECFY2 3 ECFY2/ ECC2750 :REM SUB 750 - 910 QIU GUAN CE BIAN、J IAODEGAI ZHENG SHU ,DAO XIAN DIAN ZUO BIAO PINGCHA ZHI J I ZUI RUO DIAN J ING DU PING DING748 :PRINT ”V (16) WEI GUAN CE J IAO GAI ZHENGSHU , V ( 7 11 ) WEI BIAN CHANG GAI ZHENGSHU :”750 :FOR I = 1 TO 2 3 N - 1760 :V( I) = (A1 ( I) 3 KA + B1 ( I) 3 KB + C1 ( I) 3 KC) / P( I)770 : IF I < = N THEN B ( I) = B ( I) + DEG(V ( I) / 10000)ELSE S( I - N) = S( I - N) + V( I) / 1000780 :PRINT ”V(”; I ;”) = ”;V( I) ,790 :NEXT I795 :LOCATE 12798 :PRINT ”# # # # # ”68 西部探矿工程Feb. 2004No. 2800 :FOR I = 1 TO N810 :A( I) = A( I - 1) + B( I) + 180820 : IF A( I) > 360 THEN A( I) = A( I) - 360830 :NEXT I835 :PRINT ”DAO XIAN DIAN PING CHA ZHI :”840 :FOR I = 2 TO N850 :X( I) = X( I - 1) + S( I - 1) 3 COS(A( I - 1) )860 : Y( I) = Y( I - 1) + S( I - 1) 3 SIN(A( I - 1) )870 :PRINT ”X(”; I ;”) = ”;X( I) ,”Y( ;”I ;”) = ”; Y( I)880 :STOP890 :NEXT I894 :LOCATE 12896 :PRINT ”# # # # # ”900 :M0 = SQR( EPVV/ 3) :MX( FN) = M0 3 SQR( EFFX3) :MY( FN) = M0 3 SQR ( EFFY3) :M ( FN) = SQR (MX(FN) 3 MX(FN) + MY(FN) 3 MY(FN)910 : PRINT ”J I SUAN DAN WEI QUAN ZHONG WUCHA :M0 = ”;M0 ,”ZUI RUO DIAN WU CHA :M( ”;FN ;”) = ”;M(FN) ;”mm”911 :LOCATE 12912 :PRINT ”# # # # # ”914 :LOCATE 15916 :PRINT ”END !”918 :REM 920 HANG FANG ZHI GUAN CE ZUO J IAO、BIAN CHANG ZHI ,930 HANG FANG ZHI FANG WEIJ IAO TIAO J IAN XI SHU920 : : DATA 230. 3237 , 204. 952 , 180. 0042 , 200. 130 ,170. 3922 , 345. 153 , 236. 4837 , 278. 059 , 192. 1425 ,451. 692 ,260. 5901 ,0930 :DATA 1 ,1 ,1 ,1 ,1 ,1 ,0 ,0 ,0 ,0 ,0 ,0940 : END4 运行结果RUNGUAN CE J IAO GE SHU :N = 6QI SHI FANG WEI J IAO :A0 = 226. 4459QI SHI DIAN ZUO BIAO :X0 = 3020. 348Y0 = - 9049. 801DAO XIAN DIAN GAI LUE ZUO BIAO ZHI :X(2) = 3046. 36649 Y(2) = - 9253. 09478X(3) = 3071. 81326 Y(3) = - 9451. 60039X(4) = 3059. 53346 Y(4) = - 9796. 53488X(5) = 3286. 66789 Y(5) = - 9956. 93052X(6) = 3702. 48850 Y(6) = - 10133. 33597# # # # #ZUI RUO DIAN HAO :FN = 4FU HE FANG WEI J IAO :A2 = 57. 5931FU HE DIAN ZUO BIAO :X2 = 3702. 437Y2 = - 10133. 399J IA SHE DAN WEI QUAN ZHONG WU CHA :D0 = 5V(16) WEI GUAN CE J IAO GAI ZHENG SHU ,V(711) WEI BIAN CHANG GAI ZHENG SHU :V(1) = - 4. 01 V(2) = - 3. 79V(3) = - 3. 57 V(4) = - 3. 93V(5) = - 1. 03 V(6) = 4. 33V(7) = 2. 77 V(8) = 2. 70V(9) = 4. 67 V(10) = 2. 31V(11) = 2. 64# # # # #DAO XIAN DIAN ZUO BIAO PING CHA ZHI :X(2) = 3046. Y(2) = - 9253. X(3) = 3071. Y(3) = - 9451. X(4) = 3059. Y(4) = - 9796. X(5) = 3286. 62793 Y(5) = - 9956. X(6) = 3702. Y(6) = - 10133. 399# # # # #J I SUAN DAN WEI QUAN ZHONG WU CHA : M0 =4. ZUI RUO DIAN WU CHA :M(4) = 13. mm# # # # #END !(上接第66 页)参考文献: 1 姜永智,廖雪珍,郭亮明,赵国育,刘宝山. 铜湿法冶金技术的研究于应用评述J . 甘肃有色金属,第17 卷,第1 期. 2 王成彦. 低品位铜湿法冶炼的现状及发展趋势J . 新疆地质,2001 ,12.3 向鹏成,谢英. 溶浸技术在低品位矿石回收利用中的应用与发展J . 中国资源综合利用,2002 ,3.4 刘大星. 我国铜湿法冶金技术的进展 J . 有色金属(矿山部分) ,2002 ,5.Leaching Technique and its Development Trend of Low- Grade Copper OresZHAO Guo2dong , WU Heng2shan , L IU Qing , L I Chao( School of A rchitect ural Engi neeri ng , Resources and Envi ronment ,Nanhua University , Hengyang , Hunan , , Chi na)Abstract : The presence of the hydro - metallurgy for low - grade copper oresis introduced. The mechanism and technique of leaching are analyzed. Opin2ions are put forward about the developing trend of leaching for low - gradecopper ores.Key words : low - grade copper ore ; leaching ; development69 2004 年第2 期梁彩艳:单一符合导线测量控制网平差程序专心-专注-专业
