控制测量学编程实验报告

控制测量学实验报告专 业：测绘工程 班级学号：编 写 人： 指导老师：日 期：实验四一、实验名称三差改正与电测波测距归算二、实验目的与要求为了强化对三差改正与电测波测距归算的学习，通过已经掌握的VB语言的知识，利用 VB 编写一个简单的程序，用于将地面观测的水平方向归算至椭球面（包括垂线偏差改正， 标高差改正，截面差改正）和将地面观测的长度归算至椭球面（电磁波测距的归算）。三、实验原理和过程标高差改正：cos A )cot Zu m m 1垂线偏差改正公式为：6 = H (1) cos2 B sin 2A h 2 2 2 2 1e2截面差改正：6”二S2(2)2cos2B sin2Ag12p”iii(H H 电磁波测距的归算：D 224R 3AD丿+(H F1 + 2Ra丿四、程序代码Private Sub Command1_Click()计算垂线偏差改正Text5.Enabled = TrueText5.Text = (Text2.Text * Cos(Text3.Text) - Text1.Text * Sin(Text3.Text) * Tan(Text4.Text)End SubPrivate Sub Command2_Click()求标高差改正Text11.Enabled = TrueText1 1.Text = 0.5 * Text9.Text * Text9.Text * Text8.Text * 206265 / Text10.Text* Sin(2 *Text7.Text) * Cos(Text6.Text) * Cos(Text6.Text)End SubPrivate Sub Command3_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)End SubPrivate Sub Command4_Click()Text21.Enabled = TrueDim a#, 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 / eEnd SubPrivate Sub Form_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:End Sub五、实验结果六、实验体会通过本次实验，我们进一步的将vb程序设计与测量专业紧密的联系在了一起，让我们 进一步学会怎样用编程的方法将测量变得简单化，力求将所有复杂的计算通过编程方法快速 的计算出来。控制测量学实验报告专 业：测绘工程 班级学号：编 写 人： 指导老师：日 期：实验五、实验名称高斯坐标正算、实验要求与目的要求在VB环境中设计一个界面，在输入大地坐标的经度和纬度后能够计算 出高斯平面内的x坐标和y坐标，并把计算结果输出显示在界面上。三、实验安排实验设计界面如下:四、实验原理NNx = X +t cos2 Bl 2 +t (5 -12 + 9n 2 + 4n 4)cos4 Bl 4 +224720t -58t 2 +14)COS6 B16N=N COS Bl +(1- 12 +n 2)COS3 Bl 3 +6 120(5 - 18t2 + 14 + 14n 2 - 58n 2t2)COS5 Bl5五、程序代码Private Sub Form_Load()Textl.Text = : Text2.Text = : Text3.Text =Text4.Text = : Text5.Text = : Text6.Text =Text7.Text = : Text8.Text =End Sub高斯坐标正算Private Sub DadiZs()Dim t As Double, Itp As Double, X0 As Double, N As Double, L0 As DoubleDim v As Double, ll As Double, W As Double, M As DoubleConst Pi = 3.14159265358979, p = 206264.806Lat = Radian(Lat)Lon = Radian(Lon)L0 = Radian(Lo)a = 637814075椭球参数b = 6356755.28815753ep = 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)End If11 = 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 = YEnd Sub六、实验结果以教材上的例题为例：输入数据计算结果输出:七、实验体会本次高斯坐标正算编程并不是一件简单的事，并不是因为计算复杂，而是公式分散得 太散，但通过看书终于找到了一些相关的资料和例题，帮助我完成了本次实验任务，希望下 次能做得更好。实验六、 实验名称高斯坐标反算、 实验要求与目的要求在VB环境中设计一个界面，在输入咼斯平面内的x坐标和y坐标后能够计算出大地坐标的经度和纬度，并把计算结果输出显示在界面上。三、实验原理B = B - 1 v 21 1( 、 y21(5 + 3t 2 +耳 2 9 耳 2 12)( 、 y+ 丄(61 + 9012 + 4512)( 6yf 2 f f 1Nf丿12fff fNf丿360ffNf丿180冗y- (1 + 212 + n 2 )( y+ 1 (5 + 2812 + 2412 + 6n 2 + 8 耳 212)( yI N J6ff120ffff f:N丿_1_cos B180冗四、程序代码Private Sub DadiFs()Dim t As Double, Itp As Double, X0 As Double, Bf As Double, N As DoubleDim v As Double, ll As Double, W As Double, M As Double, L0 As DoubleDim a, b, c, d, e, f As DoubleX = Text1.Text: Y = Text2.TextL0 = Radian(Lo)X0 = X * 0.000001Y = Y - 500000#If Tq = 0 Thena = 6378245b = 6356863.01877305ep = 0.006693421622966ep1 = 0.006738525414683f = (a - b) / ac = a 人 2 / bd = b 人 2 / aIf X0 3 ThenBf = 9.04353301294 * X0 - 0.00000049604 * X0 人 2 - 0.00075310733 * X0 人 3 - 0.00000084307 * X0 人 4 - 0.00000426055 * X0 人 5 - 0.00000010148 * X0 人 6ElseIf X0 6 ThenBf = 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)人 6End IfElsea = 637814075椭球参数b = 6356755.28815753 ep = 0.006694384999588 ep1 = 0.006739501819473f = (a - b) / ac = a 人 2 / bd = b 人 2 / aIf X0 3 ThenBf = 9.04369066313 * X0 - 0.00000049618 * X0 人 2 - 0.00075325505 * X0 人 3 - 0.0000008433 * X0 人 4 - 0.00000426157 * X0 人 5 - 0.0000001015 * X0 人 6ElseIf X0 6 ThenBf = 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)人 6End IfEnd IfBf = Bf * Pi / 180# t = Tan(Bf)Itp = ep1 * Cos(Bf)人 2 W = 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 = b * 60: Text5.Text = Int(b) c = b - Int(b): Text7.Text = c * 60 d = Lat * 180 / Pi: Text4.Text = Int(d) e = d - Int(d): e = e * 60: Text6.Text = Int(e) f = e - Int(e): Text8.Text = f * 60End SubPrivate Sub Form_Load()Textl.Text = : Text2.Text = : Text3.Text =Text4.Text = : Text5.Text = : Text6.Text =Text7.Text = : Text8.Text =End Sub五、实验结果数据输入：奩择反算口回区I输代尘标黑 |2E34561. 381弗茂&了 36高斯坐标反算E |匿分抄L |度廿秒数据输出：奩择反算匚亘凶输入坐标黑 |2E3456L. 381|5弗茂& 了 36高斯坐标反算B |25 匿 分t-L |1H 度 | 廿 |OO7S秒六、实验体会本次高斯坐标反算编程跟上次正算一样，所以这次比较顺利地完成了此次实验任务， 经过多次程序编程，切实的感受到VB程序设计给测量工作内业计算带来的方便，公式简单 化，只需要输入数据，那些复杂的计算就能由计算机来完成，节约时间，方便快捷，在今后 的学习工作中将会有更广的应用空间。