高程网平差程序设计测绘软件设计学习报告.doc

文档介绍：《测绘软件设计》
综合测试作业

高程网平差程序设计
一、题目
请根据间接平差模型,编写一个水准网平差的程序,并以下面的水准网观测数据为算例进行验算(A,B为已知高程点,P1,P2,P3为待求高程点)。
路线号
观测高差(m)
水准路线长(km)
已知高程(m)
1
+

HA=
HB=
2
+

3
+

4
+

5
+

6
+

7
+

h6
h5
h3
h1
h7
h2
h4
A
B
P2
P1
P3
程序设计过程
程序的详细步骤
(1)建立系数阵B和常数阵L
(2)建立权矩阵P和法方程系数矩阵Naa
(3)建立法方程常系数矩阵W和求取观测值改正数Vi
(4)调试程序,解算高程平差值
2、所有相关数据
路线号
观测高差(m)
水准路线长(km)
已知高程(m)
1
+

HA=
HB=
2
+

3
+

4
+

5
+

6
+

7
+

程序清单及程序说明
'已知点名
Dim strKnowPName() As String
'已知点高程
Dim dblKnowH() As Double
'未知点名
Dim strUnKnowPName() As String
'已知点的数量
Dim intPntNum As Integer
'未知点的数量
Dim t As Integer
'观测个数
Dim n As Integer
'申明水准路线观测值
Dim strBegin() As String
Dim strEnd() As String
Dim diff_H() As Single
Dim sur_D() As Single
Dim LL() As Double
Dim BB() As Double
'定义变量
Dim Q() As Double
Dim BTPB() As Double
Dim BTPL() As Double
Dim NN() As Double
Dim XX() As Double
Dim v() As Double
Private Sub readdata_Click()

For Input As #1

Do While Not EOF(1)
Line Input #1, strData
strData
Loop
Close #1
'读取观测路线数据

'组建误差方程
BuildError
'组建法方程
BuildFaFangCheng
'BTPT求逆阵
Matrix_InverseOpp BTPB, NN
'求出未知数
MatrixMultiply_nx1 NN, BTPL, XX
Dim vtpv As Double
'精度评定
vtpv = Precision_Eval(n)
'输出结果
Dim strtemp As String
Dim strtemp2 As String
Dim i, j As Integer
vbCrLf
vbCrLf
"平差高程值为:" & vbCrLf
For i = 1 To t
strtemp = ""
For j = 1 To t
strtemp = strUnKnowPName(i) & "的高程为: "
strtemp2 = Format(XX(i), "")
strtemp = strtemp & strtemp2 & " "
Next j
strtemp
Next i
vbCrLf
"单位权中误差:" & Format(vtpv, "")
vbCrLf
End Sub
'*************************