控制测量学编程实验报告

控制测量学

实验报告专业：测绘工程班级学号：编写人：指导老师：日期：实验四一、实验名称三差改正与电测波测距归算二、实验目的与要求为了强化对三差改正与电测波测距归算的学习，通过已经掌握的VB语言的知识，利用VB编写一个简单的程序，用于将地面观测的水平方向归算至椭球面（包括垂线偏差改正，标高差改正，截面差改正）和将地面观测的长度归算至椭球面（电磁波测距的归算）。三、实验原理和过程标高差改正：'cosA)cotZumm1垂线偏差改正公式为：6''=H(1)cos2Bsin2Ah22221e2截面差改正：6”二— S2(2)2cos2Bsin2Ag 12p” i i i(H—H¥电磁波测距的归算：D224R3AD丿+

(HF

1电磁波测距的归算：D224R3A四、程序代码PrivateSubCommand1_Click()'计算垂线偏差改正Text5.Enabled=TrueText5.Text=(Text2.Text*Cos(Text3.Text)-Text1.Text*Sin(Text3.Text))*Tan(Text4.Text)EndSubPrivateSubCommand2_Click()'求标高差改正Text11.Enabled=TrueText11.Text=0.5*Text9.Text*Text9.Text*Text8.Text*206265/Text10.Text*Sin(2*Text7.Text)*Cos(Text6.Text)*Cos(Text6.Text)EndSubPrivateSubCommand3_Click()Text16.Enabled=True'计算截面差改正Text16.Text=Text12.Text*Text12.Text/12*206265*Text13.Text*Text13.Text/Textl4.Text/Textl4.Text*Cos(Textl5.Text)*Cos(Textl5.Text)EndSubPrivateSubCommand4_Click()Text21.Enabled=TrueDima#,b#,c#,d#'电磁波测距的归算a=Text17.Text:b=Text18.Text:c=Text19.Text:d=Text20.TextText21.Text=a*Sqr((1-(c-b)*(c-b)/a/a)/(1+b/e)/(1+c/e))+a*a*a/24/e/eEndSubPrivateSubForm_Load()'把所有的text清空Text1.Text=0:Text2.Text=0:Text3.Text=0:Text4.Text=0:Text5.Text=0Text6.Text=0:Text7.Text=0:Text8.Text=0:Text9.Text=0:Text10.Text=0Text11.Text=0:Text12.Text=0:Text13.Text=0:Text14.Text=0:Text15.Text=0Text16.Text=0:Text17.Text=0:Text18.Text=0:Text19.Text=0:Text20.Text=0Text21.Text=0Text21.Enabled=False:Text11.Enabled=False:Text5.Enabled=False:Text16.Enabled=False:EndSub五、实验结果六、实验体会通过本次实验，我们进一步的将vb程序设计与测量专业紧密的联系在了一起，让我们进一步学会怎样用编程的方法将测量变得简单化，力求将所有复杂的计算通过编程方法快速的计算出来。控制测量学实验报告专业：测绘工程班级学号：编写人：指导老师：日期：

实验五实验五、实验名称高斯坐标正算、实验要求与目的要求在VB环境中设计一个界面，在输入大地坐标的经度和纬度后能够计算出高斯平面内的x坐标和y坐标，并把计算结果输出显示在界面上。三、实验安排实验设计界面如下:四、实验原理N Nx四、实验原理N Nx=X+ tcos2Bl2+ t(5-12+9n2+4n4)cos4Bl4+2 24720t⑹-58t2+14)COS6B16N=N=NCOSBl+ (1-12+n2)COS3Bl3+6120(5-18t2+14+14n2-58n2t2)COS5Bl5五、程序代码PrivateSubForm_Load()Textl.Text="":Text2.Text="":Text3.Text=""Text4.Text="":Text5.Text="":Text6.Text=""Text7.Text="":Text8.Text=""EndSub'高斯坐标正算PrivateSubDadiZs()DimtAsDouble,ItpAsDouble,X0AsDouble,NAsDouble,L0AsDoubleDimvAsDouble,llAsDouble,WAsDouble,MAsDoubleConstPi=3.14159265358979,p=206264.806Lat=Radian(Lat)Lon=Radian(Lon)L0=Radian(Lo)a=6378140 '75椭球参数b=6356755.28815753

ep=0.006694384999588epl=0.006739501819473f=(a-b)/ac=a人2/bd=b人2/aX0=111133.0047*(Lat*180/Pi)-(32009.8575*Sin(Lat)+133.9602*(Sin(Lat))人3+0.6976*(Sin(Lat))人5+0.0039*(Sin(Lat))人7)*Cos(Lat)EndIf11=Lon-L0t=Tan(Lat)Itp=ep1*Cos(Lat)人2W=Sqr(1-ep*Sin(Lat)人2)v=Sqr(1+ep1*Cos(Lat)人2)M=c/v人3N=a/Wx=X0+N*t*(Cos(Lat))人2*11人2/2+N*t*(5 -t*t+9*Itp人2+4*Itp人4) *(Cos(Lat))人4*11人4/24+N*t*(61-58*t人2+t人4+270*Itp人2-330*t人2*Itp人2)*(Cos(Lat))人6*11人6/720+N*t*(1385-3111*t人2+543*t人4-t人6)*Cos(Lat)人8*11人8/40320y=N*Cos(Lat)*11+N*(1-t*t+Itp)*(Cos(Lat))人3*11人3/6+N*(5-18*t*t+t人4+14*Itp-58*Itp*t*t)*(Cos(Lat))人5*11人5/120+N*(61-479*t人2+179*t人4-t人6)*Cos(Lat)人7*11人7/5040r=Sin(Lat)*11+Sin(Lat)*(Cos(Lat))人2*11人3*(1+3*Itp+2*Itp人2)/3+Sin(Lat)*(Cos(Lat))人4*11人5*(2-t*t)/15r=Degree(r)y=y+500000#Text7.Text=X:Text8.Text=YEndSub六、实验结果以教材上的例题为例：输入数据 计算结果输出:以教材上的例题为例：输入数据 计算结果输出:七、实验体会本次高斯坐标正算编程并不是一件简单的事，并不是因为计算复杂，而是公式分散得太散，但通过看书终于找到了一些相关的资料和例题，帮助我完成了本次实验任务，希望下次能做得更好。实验六、实验名称高斯坐标反算、实验要求与目的要求在VB环境中设计一个界面，在输入咼斯平面内的x坐标和y坐标后能够计算出大地坐标的经度和纬度，并把计算结果输出显示在界面上。三、实验原理B=B-1v211(、y21 (5+3t2+耳2—9耳212)(、y+丄(61+9012+4512)( \6yf2ff1N< f丿12 f f ffN< f丿360 f fN< f丿180冗y-—(1+212+n2)(\y+1(5+2812+2412+6n2+8耳212)(\yINJ6 f f120 f f f ff:N丿_1_cosB180冗四、 程序代码PrivateSubDadiFs()DimtAsDouble,ItpAsDouble,X0AsDouble,BfAsDouble,NAsDoubleDimvAsDouble,llAsDouble,WAsDouble,MAsDouble,L0AsDoubleDima,b,c,d,e,fAsDoubleX=Text1.Text:Y=Text2.TextL0=Radian(Lo)X0=X*0.000001Y=Y-500000#IfTq=0Thena=6378245b=6356863.01877305ep=0.006693421622966ep1=0.006738525414683f=(a-b)/ac=a人2/bd=b人2/aIfX0<3ThenBf=9.04353301294*X0-0.00000049604*X0人2-0.00075310733*X0人3-0.00000084307*X0人4-0.00000426055*X0人5-0.00000010148*X0人6ElseIfX0<6ThenBf=27.11115372595+9.02468257083*(X0-3)-0.00579740442*(X0-3)人2-0.00043532572*(X0-3)人3+0.00004857285*(X0-3)人4+0.00000215727*(X0-3)人5-0.00000019399*(X0-3)人6EndIfElsea=6378140 '75椭球参数b=6356755.28815753ep=0.006694384999588ep1=0.006739501819473f=(a-b)/ac=a人2/bd=b人2/aIfX0<3ThenBf=9.04369066313*X0-0.00000049618*X0人2-0.00075325505*X0人3-0.0000008433*X0人4-0.00000426157*X0人5-0.0000001015*X0人6ElseIfX0<6ThenBf=27.11162289465+9.02483657729*(X0-3)-0.00579850656*(X0-3)人2-0.00043540029*(X0-3)人3+0.00004858357*(X0-3)人4+0.00000215769*(X0-3)人5-0.00000019404*(X0-3)人6EndIfEndIfBf=Bf*Pi/180#t=Tan(Bf)Itp=ep1*Cos(Bf)人2W=Sqr(1-ep*Sin(Bf)人2)v=Sqr(1+ep1*Cos(Bf)人2)M=c/v人3N=a/WLat=Bf-0.5*v人2*t*((Y/N)人2-(5+3*t*t+Itp-9*Itp*t*t)*(Y/N)人4/12+(61+90*t*t+45*t人4)*(Y/N)人6/360)11=((Y/N)-(1+2*t*t+Itp)*(Y/N)人3/6+(5+28*t*t+24*t人4+6*Itp+8*Itp*t*t)*(Y/N)人5/120)/Cos(Bf)r=Y*t/N-Y人3*t*(1+t*t-Itp)/(3*N人3)+Y人5*t*(2+5*t*t+3*t人4)/(15*N人5)Lat=Degree(Lat)Lon=Degree(L0+11)a=Lat*180/Pi:Text3.Text=Int(a)b=a-Int(a):b