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

测绘程序设计(六) 2.大地四边形条件平差 共有4个点，其中2个为起算点，2个未知点，应有4个必要观测（t=4），共有8个观测值(n=8)，条件方程的个数为： r = n - t = 4 对于大地四边形， 可以列出7个图形条件，但是只有 3 个是相互独立的，其余几个可以由这 3 个方程推导出来： 其中图形条件3个： 极条件1个: 以A点为极有: 作 业 　　 1 实验一的实验报告? 2 间接平差计算程序编制? 3 条件平差计算程序编制? 4 水准网通用平差计算程序设计? 5 方向误差方程和距离误差方程的线性化形式? 1、通用水准网间接平差程序的设计: (1) 按文件方式进行输入, 约定数据格式 (2) 高程近似值计算 (3) 列立观测值的误差方程 (4) 解误差方程并求高程平差值 (5) 精度估计 (6) 计算结果输出. 三、水准网的间接平差 2、水准路线误差方程的组成 i j Xi Xj hij 当i点已知时： 当j点已知时： 平面控制网可分为测边网、测角网、边角网等。 （1）数据输入 （2）坐标近似值计算 （3）列立误差方程 （4）平差计算 线性化的误差方程是近似式，略去了参数和观测值二次以上的各项的影响，只考虑到一次项的影响。如果近似值计算精确度很差，就要把第一次的平差结果作为未知数的近似值再进行一次平差。就是迭代平差。 四、平面控制网的间接平差 测角网间接平差算例： A B D C 1 2 3 4 5 6 7 8 9 12 11 10 13 14 15 16 17 18 P2 P1 设有一测角三角网，A、B、C、D为已知点，P1、P2为待定点，同精度观测了18个角度，按间接平差求平差后P1、P2点的坐标及精度。已知数据见下表。 平面控制网平差示例 10156.11 12168.60 10232.16 11879.60 边长 216°49’06.5 A 95°53’29.1 49923.19 17814.63 D 34°40’56.3 37818.86 19063.66 C 274°39’38.4 31996.50 10649.55 B 43836.82 9684.28 A Y（m) X（m) 方位角 坐标（m） 点名 29°52’55.4 18 47°02’11.4 12 31°10’22.6 6 120°08’31.1 17 67°02’49.4 11 31°26’50.0 5 29°58’35.5 16 65°55’00.8 10 117°22’46.2 4 66°46’08.2 15 27°53’59.3 9 30°05’46.7 3 66°34’54.7 14 130°03’14.2 8 23°39’46.9 2 46°38’56.4 13 22°02’43.0 7 126°14’24.1 1 观测值 角度 编号 观测值 角度 编号 观测值 角度 编号 解：n=18, t=2*6-4-4=4, r=18-4=14 设P1、P2点的坐标作为未知参数X1、Y1、X2、Y2，根据前方交会可以求出P1、P2的近似坐标： 根据角度的误差方程： V B x l 定权，P为单位阵，形成法方程为： 精度评定： * * 第四讲　测量平差编程计算(四) ----测量平差程序设计 主要内容 一、大地四边形的条件平差 二、单导线的间接平差 三、水准网的间接平差 四、平面控制网的间接平差 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 系数矩阵