導(dǎo)線網(wǎng)平差資料

'***********************************************************************************'本程序用于單一附（閉）合導(dǎo)線嚴(yán)密平差計(jì)算，采用按角度條件平差法。以方向觀測值中'誤差的先驗(yàn)值作為單位權(quán)中誤差。計(jì)算結(jié)果可求得各待定點(diǎn)的坐標(biāo)平差值及其點(diǎn)位精度Mx，My，'及M，并計(jì)算出各待定點(diǎn)誤差橢圓元素E，F(xiàn)，Z''參考文獻(xiàn)：郭久訓(xùn)．《控制網(wǎng)平差程序設(shè)計(jì)》北京：原子能出版社，2004．8''平差數(shù)據(jù)來源：潘正風(fēng)等．《數(shù)字測圖原理與方法》武漢大學(xué)出版社186頁表6-5''等級：所用平差數(shù)據(jù)為首級圖根導(dǎo)線數(shù)據(jù)（精度很不高）。本程序中寫出了方位角和導(dǎo)線全長相對閉合差的判'別，但考慮到程序的通用性，將這些限差判斷當(dāng)作了注釋處理，而不實(shí)際運(yùn)行。''程序不足：沒有導(dǎo)線網(wǎng)的圖形表達(dá)。'***********************************************************************************Privatei%,j%,n%,tc#,tb#,B_x!,B_y!,e1!,e2!,m!,m0#,z#,aaAsBoolean,bbAsBoolean,ccAsBoolean',ddAsBooleanPrivateNaa#(2,2),Naa逆#(2,2),W#(2),K#(2),qq#(2),fx#(2),fy#(2)PrivateA#(),Q#(),V#(),C#(),mx#(),my#(),mk#(),e#(),f#(),zz#()PrivatePo()AsPoint'文件格式說明：'文件格式詳見文件"平差數(shù)據(jù).txt"'PrivateSub打開文件_Click()'打開文件Dimff$,temp$,A_name$,A_x!,A_y!,A_l#,A_s!,B_name$,B_l#,C_name$,C_x!,C_y!,D_name$,D_x!,D_y!Form1.Cls'清屏CommonDialog1.DialogTitle="打開數(shù)據(jù)文件"CommonDialog1.FileName=""CommonDialog1.ShowOpen'出錯(cuò)處理OnErrorGoToFileErrff=CommonDialog1.FileName'ff是文件路徑名OpenffForInputAs#1'以順序文件方式打開文件，使用inputLineInput#1,temp'讀取文件中的說明語句LineInput#1,temp'讀取文件中的說明語句LineInput#1,temp'讀取文件中的說明語句Input#1,n'讀取n，n為（測站數(shù)-1）ReDimPo(n)AsPoint'定義Po(n)，其中Po(0)存A點(diǎn)數(shù)據(jù)，Po(n)存B點(diǎn)數(shù)據(jù)，Po(1)到Po(n-1)存n-1個(gè)未知點(diǎn)數(shù)據(jù)。LineInput#1,temp'讀取文件中的說明語句LineInput#1,temp'讀取文件中的說明語句Input#1,m,e1,e2'讀取先驗(yàn)方向觀測值中誤差m，測距儀固定誤差e1，比例誤差e2LineInput#1,temp'讀取文件中的說明語句LineInput#1,temp'讀取文件中的說明語句Input#1,C_name,C_x,C_y'讀取已知點(diǎn)CInput#1,A_name,A_x,A_y,A_l,A_s'讀取已知點(diǎn)AFori=1Ton-1'讀取n-1個(gè)未知點(diǎn)Input#1,Po(i).name,Po(i).l,Po(i).sNextiInput#1,B_name,B_x,B_y,B_l'讀取已知點(diǎn)BInput#1,D_name,D_x,D_y'讀取已知點(diǎn)DClose#1tc=ZBiaoFSuan(C_x,C_y,A_x,A_y)'坐標(biāo)反算,求點(diǎn)C到A的坐標(biāo)方位角，并記作tc，單位是度tb=ZBiaoFSuan(B_x,B_y,D_x,D_y)'坐標(biāo)反算,求點(diǎn)B到D的坐標(biāo)方位角，并記作tb，單位是度Po(0).name=A_name:Po(0).x=A_x:Po(0).y=A_y:Po(0).l=A_l:Po(0).s=A_s'Po(0)存A點(diǎn)數(shù)據(jù)Po(n).name=B_name:Po(n).x=B_x:Po(n).y=B_y:Po(n).l=B_l'Po(n)存B點(diǎn)數(shù)據(jù)Fori=0TonPo(i).l=deg(Po(i).l)'將觀測方向左角的單位度分秒化作度NextiReDimQ#(2*n),A#(2,2*n),V#(2*n),C#(1,2*n-2),mx#(n-1),my#(n-1),_mk#(n-1),e#(n-1),f#(n-1),zz#(n-1)'變量重新定義bb=TrueMsgBox"文件已成功打開",,"提示"顯示平差數(shù)據(jù)ffA(1,i)=Cos(Po(i).t/p0)A(2,i)=Sin(Po(i).t/p0)Q(i)=(e1+e2*Po(i).s*0.0001)^2/m^2'求距離觀測值的權(quán)倒數(shù)Nexti'求WW(0)=(Po(n).t-tb)*3600'單位是sW(1)=(Po(n).x-B_x)*100'單位是cmW(2)=(Po(n).y-B_y)*100'單位是cmFori=0Ton-1sums=sums+Po(i).sNexti'限差判斷'這是首級圖根導(dǎo)線的精度要求，其中n+1為測站數(shù)'IfAbs(W(0))>40*Sqr(n+1)OrSqr((W(1)/100)^2+(W(2)/100)^2)/sums>0.00025Then'方位角閉合差和導(dǎo)線全長相對閉合差'ExitSub'EndIf'dd=TrueEndSub'''PrivateSub組法方程式()Dimg%Fori=0To2Forj=iTo2'由于Naa是對稱的，只求其上三角即可。Forg=0To2*nNaa(i,j)=Naa(i,j)+A(i,g)*A(j,g)*Q(g)NextgNextjNextiEndSub'''對于3*3的矩陣來說，利用公式Naa逆=Naa*/|Naa|去求Naa的逆即可，其中Naa*是伴隨矩陣。PrivateSub求逆求K()Dimee#,ii#,jj#,kk#,det#'det是|Naa|的值ee=1Fori=0To2ee=ee*Naa(i,i)Nextiee=ee+2*Naa(0,1)*Naa(1,2)*Naa(0,2)ii=Naa(0,1)^2jj=Naa(0,2)^2kk=Naa(1,2)^2det=ee-(ii*Naa(2,2)+jj*Naa(1,1)+kk*Naa(0,0))'求伴隨矩陣Naa*，用Naa逆存儲Naa逆(0,0)=Naa(1,1)*Naa(2,2)-kkNaa逆(0,1)=Naa(0,2)*Naa(1,2)-Naa(0,1)*Naa(2,2)Naa逆(1,0)=Naa逆(0,1)Naa逆(1,1)=Naa(0,0)*Naa(2,2)-jjNaa逆(0,2)=Naa(0,1)*Naa(1,2)-Naa(0,2)*Naa(1,1)Naa逆(2,0)=Naa逆(0,2)Naa逆(2,2)=Naa(0,0)*Naa(1,1)-iiNaa逆(1,2)=Naa(0,2)*Naa(0,1)-Naa(0,0)*Naa(1,2)Naa逆(2,1)=Naa逆(1,2)'求Naa逆Fori=0To2Forj=0To2Naa逆(i,j)=Naa逆(i,j)/detNextjNexti'求KFori=0To2Forj=0To2K(i)=K(i)-Naa逆(i,j)*W(j)NextjNextiEndSub''''PrivateSub求改正數(shù)和平差值()Dimpvv#pvv=0Fori=0To2*nForj=0To2V(i)=V(i)+K(j)*A(j,i)*Q(i)'求改正數(shù)VNextjpvv=pvv+V(i)*V(i)/Q(i)'求pvvNextim0=Sqr(pvv/3)'求單位權(quán)中誤差m0Fori=0Ton-1Po(i).s=Po(i).s+V(i)/100'求導(dǎo)線邊長的平差值，其中邊長改正數(shù)的單位是厘米，邊長的單位是米NextiFori=0TonPo(i).l=Po(i).l+V(n+i)/3600'求導(dǎo)線左角的平差值，其中角度改正數(shù)的單位是秒，角度的單位是度NextiEndSub''列出待定點(diǎn)坐標(biāo)的權(quán)函數(shù)式，據(jù)此進(jìn)行精度評定。'PrivateSub精度評定()Dimuu#,vv#,ww#,ss#,r#,kk!Forj=1Ton-1uu=0:vv=0:ww=0Fori=0Toj-1C(0,i)=Cos(Po(i).t/p0)C(1,i)=Sin(Po(i).t/p0)C(0,n+i)=(Po(i).y-Po(j).y)/2062.65C(1,n+i)=(Po(j).x-Po(i).x)/2062.65uu=uu+C(0,i)^2*Q(i)+C(0,n+i)^2*Q(n+i)vv=vv+C(1,i)^2*Q(i)+C(1,n+i)^2*Q(n+i)ww=ww+C(0,i)*C(1,i)*Q(i)+C(0,n+i)*C(1,n+i)*Q(n+i)NextiForkk=0To2ss=0:r=0Fori=0Toj-1ss=ss+A(kk,i)*C(0,i)*Q(i)+A(kk,n+i)*C(0,n+i)*Q(n+i)r=r+A(kk,i)*C(1,i)*Q(i)+A(kk,i+n)*C(1,n+i)*Q(n+i)Nextifx(kk)=ss:fy(kk)=rNextkkForkk=0To2ss=0:r=0Fori=0To2ss=ss-Naa逆(kk,i)*fx(i)r=r-Naa逆(kk,i)*fy(i)NextiK(kk)=ss:qq(kk)=rNextkkFori=0To2uu=uu+fx(i)*K(i)vv=vv+fy(i)*qq(i)ww=ww+fx(i)*qq(i)Nextimx(j)=m0*Sqr(Abs(uu))my(j)=m0*Sqr(Abs(vv))mk(j)=Sqr((mx(j))^2+(my(j))^2)'每個(gè)點(diǎn)的點(diǎn)位中誤差，單位cmr=uu+vvss=Sqr((uu-vv)^2+4*ww^2)e(j)=m0*Sqr((r+ss)/2)'點(diǎn)位誤差橢圓元素E，單位cmf(j)=m0*Sqr(Abs(r-ss)/2)'點(diǎn)位誤差橢圓元素F，單位cmz=Atn(2*ww/(uu-vv))z=z/2*p0Ifz<0Thenz=z+90Ifww<0Thenz=z+180zz(j)=dms(z)'點(diǎn)位誤差橢圓元素ZNextjEndSub'''PrivateFunction顯示平差數(shù)據(jù)(ByValfname$)DimlinesfromFile$,Nextline$OpenfnameForInputAs#1DoUntilEOF(1)LineInput#1,NextlinelinesfromFile=linesfromFile+Nextline+Chr(13)LoopClose#1PrintlinesfromFileEndFunction'''PrivateSub顯示平差結(jié)果()Print"平差結(jié)果："Print"---------------------------------