單一附合導(dǎo)線測量控制網(wǎng)平差程序

單一符合導(dǎo)線測量控制網(wǎng)平差程序梁彩艷(西北石油局工程服務(wù)大隊(duì),新疆烏魯木齊830011)摘要:由于全站儀(光電儀)的不斷普及,導(dǎo)線布網(wǎng)已在測量控制網(wǎng)中廣泛應(yīng)用。為了簡化計(jì)算,提高平差精度,節(jié)約平差經(jīng)費(fèi),結(jié)合實(shí)際工作中的一個(gè)實(shí)例,運(yùn)用PC—E500計(jì)算器編制了單一符合導(dǎo)線高斯平差程序。關(guān)鍵詞:單一;符合導(dǎo)線;平差程序1差原理及程序框圖(見圖1)2實(shí)例在單一符合導(dǎo)線(如圖2)中,已知各起始數(shù)據(jù),觀測了各折角方向值,測量了各邊長(見表1),觀測方向中誤差δr=±5″,邊長丈量中誤差δsi=±0.5simm,求各導(dǎo)線點(diǎn)的坐標(biāo)平差值和最弱點(diǎn)4號點(diǎn)的點(diǎn)位精度(可假設(shè)單位權(quán)中誤差δ0=±5″,則觀測角權(quán)為0.5,觀測邊權(quán)Psi=100/Si)。圖2單一符合導(dǎo)線示意圖表1測量數(shù)據(jù)一覽表點(diǎn)號已知坐標(biāo)(m)XY編號已知坐標(biāo)方位角B3020.348-9049.801αAB226°44′59″C3702.437-10133.399αCD57°59′31″編號觀測角度編號觀測邊長(m)1230°32′37″1204.9522180°00′42″2200.1303170°39′22″3345.1534236°48′37″4278.0595192°14′25″5451.6926260°59′01″3平差程序10:REMSUB10-190QIUDAOXIANDIANGAILUEZUOBIAO20:REMSHUZU”A”WEIFANGWEIJIAO,”B”WEIGUANCEZUOJIAO,”S”WEIBIANCHANG,”P”WEIQUAN,”V”WEIGAIZHENGSHU,”MX、MY、M”WEIDAIQIUDIANZHONGWUCHA,”X、Y”WEIZUOBIAO30:INPUT”GUANCEJIAOGESHU:N=”;N40:DIMA(N),B(N),S(N),X(23N),Y(23N),P(23N),V(23N),MX(N),MY(N),M(N)45:FORI=1TON50:READB(I),S(I)60:B(I)=DEG(B(I))70:NEXTI80:INPUT”QISHIFANGWEIJIAO:A0=”;A(0)90:A(0)=DEG(A(0))100:FORI=1TON110:A(I)=A(I-1)+B(I)+180120:IFA(I)>360THENA(I)=A(I)-360130:NEXTI140:INPUT”QISHIDIANZUOBIAO:X0=”;X(1),”Y0=”;Y(1)總第93期2004年第2期西部探礦工程WEST-CHINAEXPLORATIONENGINEERINGseriesNo.93Feb.2004145:PRINT”DAOXIANDIANGAILUEZUOBIAOZHI:”150:FORI=2TON160:X(I)=X(I-1)+S(I-1)3COS(A(I-1)):Y(I)=Y(I-1)+S(I-1)3SIN(A(I-1))170:PRINT”X(”;I;”)=”;X(I),”Y(;”I;”)=”;Y(I)180:NEXTI185:LOCATE12190:PRINT”#####”200:REMSUB200-370QIUTIAOJIANFANGCHENGXISHU”A1,BI,C1”SHUZUJIQUANHANSHUXISHU”FX,FY”SHUZU210:DIMA1(23N-1),B1(23N-1),C1(23N-1),FX(23N-1),FY(23N-1)220:FORI=1TON230:READA1(I)240:NEXTI250:FORI=1TO23N-1260:B1(I)=(-Y(N)+Y(I))/206.26270:C1(I)=(X(N)-X(I))/206.26280:IFN<ITHENB1(I)=COS(A(I-N):C1(I)=SIN(A(I-N))290:NEXTI300:INPUT”ZUIRUODIANHAO:FN=”;FN310:FORI=1TO23N-1320:FX(I)=(-Y(FN)+Y(I))/206.26330:FY(I)=(X(FN)-X(I))/206.26340:IFFN<ITHENFX(I)=0:FY(I)=0350:IFN<ITHENFX(I)=COS(A(I-N)):FY(I)=SIN(A(I-N))360:IFI>N+FN-1THENFX(I)=0:FY(I)=0370:NEXTI414:INPUT”FUHEFANGWEIJIAO:A2=”;A2,”FUHEDIANZUOBIAO:X2=”;X2,”Y2=”;Y2415:WA=-DEG(A2)+A(N):WX=-X2+X(N):WY=-Y2+Y(N)416:WA=DMS(WA):IFWA>360THENWA=WA-360419:WA=WA310000:WX=WX31000:WY=WY31000420:REMSUB420-740JIESUANFAFANGCHENG425:REMSUB430-470JISUANGUANCEJIAO、BIANDEQUAN430:INPUT”JIASHEDANWEIQUANZHONGWUCHA:D0=”;D0440:FORI=1TO23N-1450:P(I)=0.5460:IFI>NTHENP(I)=(D03D0)/(0.253S(I-N))470:NEXTI500:REM”E”WEIDAIQUANQIUHEFUHAO520: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:FORI=1TO23N-1560:EAA=A1(I)3A1(I)/P(I)+EAA:EAB=A1(I)3B1(I)/P(I)+EAB570:EAC=A1(I)3B1(I)/P(I)+EAC:EAFX=A1(I)3FX(I)/P(I)+EAFX:EAFY=A1(I)3FY(I)/P(I)+EAFY580:EBB=B1(I)3B1(I)/P(I)+EBB:EBC=B1(I)3C1(I)/P(I)+EBC590:EBFX=B1(I)3FX(I)/P(I)+EBFX:EBFY=B1(I)3FY(I)/P(I)+EBFY600:ECC=C1(I)3C1(I)/P(I)+ECC:ECFX=C1(I)3FX(I)/P(I)+ECFX:ECFY=C1(I)3FY(I)/P(I)+ECFY610:EFFX=FX(I)3FX(I)/P(I)+EFFX:EFFY=FY(I)3FY(I)/P(I)+EFFY620:NEXTI630:EBB1=EBB-EAB3EAB/EAA:EBC1=EBC-EAB3EAC/EAA:EWX1=WX-EAB3WA/EAA640:EBFX1=EBFX-EAB3EAFX/EAA:EBFY1=EBFYEAB3EAFY/EAA650:ECC2=ECC-EAC3EAC/EAA-EBC13EBC1/EBB1660:EWY2=WY-EAC3WA/EAA-EBC13EWX1/EBB1670:ECFX2=ECFX-EAC3EAFX/EAA-EBC13EBFX1/EBB1680:ECFY2=ECFY-EAC3EAFY/EAA-EBC13EBFY1/EBB1690:KC=-EWY2/ECC2700:KB=-EBC13KC/EBB1-EWX1/EBB1710:KA=-EAB3KB/EAA-EAC3KC/EAA-WA/EAA720:EPVV=WA3WA/EAA+EWX13EWX1/EBB1+EWY23EWY2/ECC2730:EFFX3=EFFX-EAFX3EAFX/EAA-EBFX13EBFX1/EBB1-ECFX23ECFX2/ECC2740:EFFY3=EFFY-EAFY3EAFY/EAA-EBFY13EBFY1/EBB1-ECFY23ECFY2/ECC2750:REMSUB750-910QIUGUANCEBIAN、JIAODEGAIZHENGSHU,DAOXIANDIANZUOBIAOPINGCHAZHIJIZUIRUODIANJINGDUPINGDING748:PRINT”V(1～6)WEIGUANCEJIAOGAIZHENGSHU,V(7～11)WEIBIANCHANGGAIZHENGSHU:”750:FORI=1TO23N-1760:V(I)=(A1(I)3KA+B1(I)3KB+C1(I)3KC)/P(I)770:IFI<=NTHENB(I)=B(I)+DEG(V(I)/10000)ELSES(I-N)=S(I-N)+V(I)/1000780:PRINT”V(”;I;”)=”;V(I),790:NEXTI795:LOCATE12798:PRINT”#####”68西部探礦工程Feb.2004No.2800:FORI=1TON810:A(I)=A(I-1)+B(I)+180820:IFA(I)>360THENA(I)=A(I)-360830:NEXTI835:PRINT”DAOXIANDIANPINGCHAZHI:”840:FORI=2TON850:X(I)=X(I-1)+S(I-1)3COS(A(I-1))860:Y(I)=Y(I-1)+S(I-1)3SIN(A(I-1))870:PRINT”X(”;I;”)=”;X(I),”Y(;”I;”)=”;Y(I)880:STOP890:NEXTI894:LOCATE12896:PRINT”#####”900:M0=SQR(EPVV/3):MX(FN)=M03SQR(EFFX3):MY(FN)=M03SQR(EFFY3):M(FN)=SQR(MX(FN)3MX(FN)+MY(FN)3MY(FN)910:PRINT”JISUANDANWEIQUANZHONGWUCHA:M0=”;M0,”ZUIRUODIANWUCHA:M(”;FN;”)=”;M(FN);”mm”911:LOCATE12912:PRINT”#####”914:LOCATE15916:PRINT”END!”918:REM920HANGFANGZHIGUANCEZUOJIAO、BIANCHANGZHI,930HANGFANGZHIFANGWEIJIAOTIAOJIANXISHU920::DATA230.3237,204.952,180.0042,200.130,170.3922,345.153,236.4837,278.059,192.1425,451.692,260.5901,0930:DATA1,1,1,1,1,1,0,0,0,0,0,0940:END4運(yùn)行結(jié)果RUNGUANCEJIAOGESHU:N=6QISHIFANGWEIJIAO:A0=226.4459QISHIDIANZUOBIAO:X0=3020.348Y0=-9049.801DAOXIANDIANGAILUEZUOBIAOZHI:X(2)=3046.36649Y(2)=-9253.09478X(3)=3071.81326Y(3)=-9451.60039X(4)=3059.53346Y(4)=-9796.53488X(5)=3286.66789Y(5)=-9956.93052X(6)=3702.48850Y(6)=-10133.33597#####ZUIRUODIANHAO:FN=4FUHEFANGWEIJIAO:A2=57.5931FUHEDIANZUOBIAO:X2=3702.437Y2=-10133.399JIASHEDANWEIQUANZHONGWUCHA:D0=5V(1～6)WEIGUANCEJIAOGAIZHENGSHU,V(7～11)WEIBIANCHANGGAIZHENGSHU:V(1)=-4.01V(2)=-3.79V(3)=-3.57V(4)=-3.93V(5)=-1.03V(6)=4.33V(7)=2.77V(8)=2.70V(9)=4.67V(10)=2.31V(11)=2.64#####DAOXIANDIANZUOBIAOPINGCHAZHI:X(2)=3046.362887Y(2)=-9253.098035X(3)=3071.802485Y(3)=-9451.607297X(4)=3059.503514Y(4)=-9796.545774X(5)=3286.62793Y(5)=-9956.959587X(6)=3702.436997Y(6)=-10133.399#####JISUANDANWEIQUANZHONGWUCHA:M0=4.336768793ZUIRUODIANWUCHA:M(4)=13.04462451mm#####END!(上接第66頁)參考文獻(xiàn):[1]姜永智,廖雪珍,郭亮明,趙國育,