工学b课程ppt课件测绘程序设计6六.pptVIP

测绘程序设计(六) ;第四讲　测量平差编程计算(四) ----测量平差程序设计; 1.条件平差解算程序 ‘通用的条件平差解算过程：输入系数矩阵A、权矩阵P、常数向量W和解向量V，求出V，并通过参数传出去 Public Sub CondiAdjust(b, P, W, V) Dim b1%, b2%, p1%, p2%, w1%, v1% 输入矩阵或向量的大小 Dim Q#(), Bt#(), QBt#(), Nbb#(), K#(), i% 几个中间矩阵 计算并检查输入矩阵或向量的大小 On Error Resume Next b1 = UBound(b, 1) - LBound(b, 1) + 1 If Err Then MsgBox 系数矩阵B大小错误！ Exit Sub End If ; On Error Resume Next b2 = UBound(b, 2) - LBound(b, 2) + 1 If Err Then MsgBox 系数矩阵B大小错误！ Exit Sub End If On Error Resume Next w1 = UBound(W) - LBound(W) + 1 If Err Then MsgBox 常数向量W大小错误！ Exit Sub End If On Error Resume Next v1 = UBound(V) - LBound(V) + 1 If Err Then MsgBox 改正数向量V大小错误！ Exit Sub End If On Error Resume Next p1 = UBound(P, 1) - LBound(P, 1) + 1 ; If Err Then MsgBox 权矩阵P大小错误！ Exit Sub End If On Error Resume Next p2 = UBound(P, 2) - LBound(P, 2) + 1 If Err Then MsgBox 权矩阵P大小错误！ Exit Sub End If If p1 p2 Then MsgBox 权矩阵P不是方阵！ Exit Sub End If If p1 b2 Then MsgBox 权矩阵P与系数矩阵A大小不符！ Exit Sub End If If b2 v1 Then MsgBox 系数矩阵B大小与解向量V大小不符！ Exit Sub End If ;If b1 w1 Then MsgBox 系数矩阵B大小与常数向量W大小不符！ Exit Sub End If 定义中间矩阵的大小 ReDim Bt(1 To b2, 1 To b1), QBt(1 To b2, 1 To b1) ReDim Nbb(1 To b1, 1 To b1), K(1 To b1), Q(1 To p1, 1 To p2) 组成法方程并计算 For i = 1 To p1 求Q矩阵 Q(i, i) = 1 / P(i, i) Next i MatrixTrans B, Bt Matrix_Multy QBt, Q, Bt Matrix_Multy Nbb, B, QBt 法方程系数矩阵 ShowMatrix Nbb kk = MRinv(Nbb) Matrix_Multy K, Nbb, W Matrix_Multy V, QBt, K End Sub;2.大地四边形条件平差;对于大地四边形，;其中图形条件3个：;极条件1个:;1. 间接平差解算过程 通用的间接平差解算过程：输入系数矩阵A、权矩阵P、常数向量L和解向量X，求出X，并通过参数传出去 Public Sub InAdjust(A, P, L, X) Dim a1%