GPS网平差系列之求误差椭圆代码
误差椭圆
2 E[∆x ] = E[( x− x)2 ] = E[( x − E(x)) 2 ] = σ x 2 2 E[∆y ] = E[( y− y)2 ] = E[( y − E( y)) 2 ] = σ y 2 ~
~
σ = σ +σ
p
ϕ
p′′
p′′′
∆ϕ
y
由广义误差传播律: 由广义误差传播律
Qϕϕ = Qxx cos2 ϕ + Qyy sin2 ϕ + Qxy sin 2ϕ
2 2 2 σϕ = σ 0 Qϕϕ = σ 0 (Qxx cos2 ϕ + Qyy sin2 ϕ + Qxy sin(2ϕ)
三、位差的极大值 E和极小值 F
上
E = σ QEE =
2 2 0
σ02
2
(Qxx + Qyy + K),
∆ψ = cosψ∆E + sinψ∆F Q = QEE cos2 ψ + QFF sin2 ψ + QEF sin 2 ψ ψψ
QEF = 0
Qψψ = QEE cos2 ψ + QFF sin 2 ψ
2 2 2 σψ = σ 0 Qψψ = σ 0 (QEE cos2 ψ + QFF sin 2 ψ ),
∧
Q∧ Q∧ Q∧ Q∧ Q∧ Q∧
∧
X1 X i
∧
X1Y i
∧
L Q∧ L Q∧ L
∧
X1 X u
∧
Y1 Y1
Y1 X i
Y1 Y i
Y1 X u
误差椭圆
Qxx Qxx
Qyy Qyy
K
其中 K Qxx Qyy 2 4Q2xy
22
7.5 点位落入误差椭圆内的概率
二维正态分布的联合概率密度函数：
f
x,
y
2
1
x y
1
2
exp
1
21 2
PP PP PP
x cos y sin cos
sin
x y
6
Q cos
s
in
Qxx Qxy
Qxy Qyy
c s
os in
Qxx cos2 Qyy sin 2 Qxy sin 2
• 误差椭圆反映的是待定点的点位分布情况，是一 个概率形象。当坐标轴旋转E角后，可得标准椭 圆方程。
17
• 在误差椭圆上量取任意方向的位差
X´
P
• 方法：过椭圆作方向的正交切线PF
F
• P－切点，F－垂足，则
OF
O Y´
• 中误差曲线在测量中的应用举例
18
•误差椭圆绘制： • 1) E角无负值； • 2) E>F； • 3) 测量坐标系绘制；
用求
特征
向量
的方
法求极
大值
和极
小值
的方
向
E、
。
F
Qxx Qxy
1
Qxy Qyy
1
c s
osE in E
0
,
得
tan E
QEE Qxx Qxy
误差椭圆
第十章——误差椭圆
2 令： K (Qxx Q yy ) 2 4Qxy K为算术平方根，恒大于零。 1 则有： Q Q xx Q yy K 2 用E表示位差的极大值，F表示位差的极小值，则有： 1 2 2 2 E 0 Q E E 0 Q xx Q yy K 2 （5） 1 2 2 2 F 0 Q F F 0 Q xx Q yy K 2 （5）式就是计算位差极大值与极小值的实用公式。
第十章——误差椭圆
（7）式和（8）式就是用极值E、F计算纵横坐标中误差 的公式。 若规定任何方向都由E 轴起算，则纵坐标轴X相对于E轴 的方位角为 360 E （如图）。故（7）式可写为：
2 x E 2 cos2 (360 E ) F 2 sin 2 (360 E )
第十章——误差椭圆
GPS 网 三 维 无 约 束 平 差 误 差 椭 圆
第十章——误差椭圆
GPS 网 三 维 无 约 束 平 差 误 差 椭 圆
第十章——误差椭圆
§10-5 相对误差椭圆
在平面控制网中，绘出各待定点的位误差椭圆后，就可应用点位误 差椭圆图解各待定点与已知点之间的边长中误差与方位角中误差。 但不能用同样的方法图解待定点与待定点之间的边长中误差与方位 角中误差。而在实际工作中，重要的却是任意两个待定点之间的相 对精度。为此，有必要研究任意两个待定点之间的相对精度问题。 设有任意两个待定点 为：
第十章——误差椭圆
极值方向
当
当
tan 2 0
2Qxy Qxx Q yy
Qxy 0 时，极大值在一、三象限；
Qxy
极小值在二、四象限。 时，极大值在二、四象限； 0 极小值在一、三象限。
测量平差基础课件——误差椭圆
tg2 0
2Q xy (Qxx Q yy )
2020/6/11
确定极值方 向的公式
两个根：2 0 和 20 180 即，使Q 取得极值的方 向值为 0 和 0 90 ，其 中一个为极大值方向， 另一个为极小值方向。 那么，哪个是极大值方 向？哪个是极小值方向 呢？下面作进一步的推 证：
9
§6-2 点位ห้องสมุดไป่ตู้差
三、位差的极大值 E 和极小值 F
1.极值方向的确定
总之：
tg2 0
2Q xy (Qxx Q yy )
两个根： 2 0
20 180
两个极值方向：0 0 90
当 当Qxy 0 ： 极大值方向 E 在一、
三象限，极小值在二、四象限。
当 当Qxy 0 ： 极大值方向 F 在一、
三象限，极小值在二、四象限。
2 P
02(Qxx
Qyy)
02(
1 Px
1 )
Py
P 0
Qxx Qyy 0
11 Px Py
Qx1 y1
Qx1 x2
Qx1 y2
Qx1 xk
Q x1
yk
Q y1 y1 Qx2 y1 Q y2 y1
Q y1 x2 Qx2 x2 Q y2 x2
Q y1 y2 Qx2 y2 Q y2 y2
三、位差的极大值 E 和极小值 F
1.极值方向的确定
cos2
0
1
cos 20
2
,
sin2
0
1
cos 20
2
或
Q
(Qxx
1 cos 20 2
Qyy
1 cos 20 2
Qxy sin20 )
误差理论与平差基础_第10章_误差椭圆
一．点位中误差 二．点位误差的计算 三．误差曲线 四．误差椭圆
一、点位中误差
控制点的平面位置是用一对平 面直角坐标来确定的。坐标是 由观测值的平差值计算所得的， 因此不可避免地带有误差。
x
A
O
Dy P¢(x, y) Du
Dx DP Ds
P(x, y)
y
在平面控制网的平差计算中，往往要评定待定点的点位精度； 待定点的点位精度通常用点位中误差简称“点位误差”的大 小来评定； 经过平差后的坐标（坐标的平差值）是估值，而不是真值！
Qxy Qyy
c s
os in
Qxx cos2 Qyy sin 2 Qxy sin 2
Dx
j
P
P¢¢ Dj
s
2 j
=
s
2 0
éëQxx
cos2
j
+
Qyy
sin2
j
+
Qxy
sin
2jùû
坐标方位角
P¢
P¢¢¢
y
二、点位任意方向的位差
x
与 j 垂直方向的位差如何求？
Q
=
æ èç
2 0.5
0.5 3
öø÷(dm2 / ('')2 )
单位权方差 0 0.5''
待定点P点到已知点A的距离为6.45km，方位角为45°，求P 点在AP方向的纵向误差和横向误差及AP边的边长相对中误 差。
s
2 j
=
s
2 0
éëQxx
cos2
j
+
Qyy
sin2
j
误差椭圆
Qmin 0
QFF 0
1 2
Qxx Qyy K
或写成
E2
1 2
2 x
2 y
2 x
2 y
2
4
2 xy
F 2
1 2
2 x
2 y
2 x
2 y
2
4
2 xy
其中：E 2
F2
2 x
2 y
2 P
11
第七章 误差椭圆
1
• 7.1 概述 • 一、点位误差 • A－－已知点 • P－－待放样点P的真位置 • P´－－P点平差后的平差值点位
P点在x、y方向的真误差：X
x y
~xP ~yP
xˆP yˆ P
x
P点点位真误差： A
y P´ u
s P
P x2 y2
QXX I C 0
就是QXX的两个特征根。
9
由特征方程：
QXX
I
Qxx
Qxy
Qxy 0
Qyy
,
得
1 2
Qxx Qyy
1 2
Qxx Qyy 2 4 QxxQyy Qx2y
即
1
Qmax QEE
• 4）误差椭圆绘制的比例尺比地形图的比例尺要 大，一般为1：10或1：1等。
19
• 7.4 相对误差椭圆
• 点位误差椭圆－待定点对已知点的点位精度情况
科傻GPS平差软件说明书
科傻系统(COSA)系列软件GPS工程测量网通用平差软件包（CosaGPS V5.1）使用说明书2007年11月CosaGPS版权所有不得翻录Tel: 131******** 136******** 133******** 133******** Email: **************************CosaGPS目录目录 (1)１.简介 (3)1.1 功能全面 (3)1.2 整体性好 (3)1.3 解算容量大，运算速度快 (3)1.4 操作简明，使用方便 (4)2.“文件”下拉菜单 (6)2.1 工程与文件 (6)2.2 “文件”菜单项 (8)2.2.1新建 (8)2.2.2打开 (9)2.2.3关闭 (9)2.2.4保存 (9)2.2.5另存为 (9)2.2.6新建工程 (9)2.2.7 打开工程 (15)2.2.8 打印 (16)2.2.9 打印预览 (16)2.2.10 打印设置 (16)2.2.11 退出 (16)3.“GPS数据处理”下拉菜单 (17)3.1 已知数据 (17)3.1.1 三维已知坐标 (18)3.1.2 二维已知坐标 (19)3.1.3 一维高程点 (19)3.1.4 输入地面边长 (19)3.1.5 输入地面方位 (20)3.2 基线数据 (20)3.3 GPS三维向量网平差（无约束平差或约束平差） (21)3.4 二维网联合/约束平差 (22)CosaGPS 13.4.1 联合/约束平差 (22)3.4.2 输出用户自定义任意两点相对精度 (23)3.5 椭球面上三维平差 (23)3.6 工程网(一点一方向)平差 (24)3.7 GPS高程拟合 (25)3.8 GPS三维秩亏自由网平差 (26)3.9 稳定性分析 (27)3.10 设置 (28)4.“查看”下拉菜单 (28)5.“工具”下拉菜单 (29)5.1 闭合差计算 (30)5.2 重复基线差 (30)5.3 网图显绘 (31)5.4 贯通误差影响值计算 (31)5.5 GPS网设计 (32)5.6 输出AutoCAD格式的GPS网图 (33)6.“坐标转换”下拉菜单 (34)6.1 XYZ-〉BLH (34)6.2 BLH->XYZ (35)6.3 BL->XY (36)6.4 XY->BL (36)6.5 XY1->XY2 (37)6.6 XYZ1->XYZ2 (39)6.7 高程面坐标变换 (41)7.“帮助”下拉菜单 (42)附录1. 功能菜单框图 (43)附录2. 算例及说明 (44)附录3. 基线解文件格式说明 (46)附录4. 方向及经纬度的角度格式说明 (54)附录5. 简要操作步骤 (55)CosaGPS 2１.简介基于全球卫星定位系统(GPS)的现代测量理论和技术改变了传统的测量模式，使工程测量行业发生了革命性变化，测量外业工作自动化程度大大提高，测量内业软件的作用更加重要。
二维网约束平差报告.doc
----------------------------------------------------------------------------------5工区补点二维网平差结果多余观测数= 20已知点数= 5总点数= 6GPS基线向量数= 12地面边长数= 0地面方位角数= 0旋转角(GPS->地面)= -0.000016006(dms)尺度(GPS->地面)= 22.9071(ppm)中央子午线= 130.060000000(dms)椭球长轴= 6378137.000(m)椭球扁率分母= 298.257222101PVV= 2.554(cm^2)M0= 0.357(cm)----------------------------------------------------------------------------------已知坐标(X,Y)----------------------------------------------------------------------------------No. Name X(m) Y(m)----------------------------------------------------------------------------------1 CPI1056-1 4755885.4720 510504.41312 CPI1057 4756731.7460 510540.34803 CPI1060 4755303.6060 516849.12564 CPI1061 4754298.2900 517281.99255 JM7502 4755340.5495 512924.4067----------------------------------------------------------------------------------GPS二维基线向量----------------------------------------------------------------------------------No. From To VectorDX(m) VectorDY(m) Distance(m)----------------------------------------------------------------------------------1 7502-1 CPI1056-1 -21.743 -2655.885 2655.9742 7502-1 CPI1057 824.514 -2619.9552746.6323 7502-1 JM7502 -566.646 -235.947 613.8064 7502-1 JM7502 -566.644 -235.945 613.8055 CPI1057 CPI1056-1 -846.255 -35.930 847.0186 CPI1060 7502-1 603.591 -3688.679 3737.7377 CPI1060 JM7502 36.945 -3924.627 3924.8018 CPI1061 7502-1 1608.884 -4121.537 4424.4309 CPI1061 CPI1060 1005.291 -432.859 1094.52210 CPI1061 JM7502 1042.238 -4357.483 4480.39311 JM7502 CPI1056-1 544.904 -2419.940 2480.53012 JM7502 CPI1057 1391.159 -2384.010 2760.222--------------------------------------------------------------------------------------平差坐标(X,Y)--------------------------------------------------------------------------------------No. Name X(m) Y(m) Mx(cm) My(cm) Mp(cm) --------------------------------------------------------------------------------------1 CPI1056-1 4755885.4720 510504.41312 CPI1057 4756731.7460 510540.34803 CPI1060 4755303.6060 516849.12564 CPI1061 4754298.2900 517281.99255 JM7502 4755340.5495 512924.40676 7502-1 4755907.2116 513160.3612 0.14 0.14 0.20 --------------------------------------------------------------------------------------最弱点--------------------------------------------------------------------------------------No. Name MX(cm) MY(cm) MP(cm)--------------------------------------------------------------------------------------6 7502-1 0.14 0.14 0.20--------------------------------------------------------------------------------------误差椭圆参数(E,F,T)--------------------------------------------------------------------------------------No. Name E(cm) F(cm) T(dms)--------------------------------------------------------------------------------------1 CPI1056-1 0.00 0.00 0.00002 CPI1057 0.00 0.00 0.00003 CPI1060 0.00 0.00 0.00004 CPI1061 0.00 0.00 0.00005 JM7502 0.00 0.00 0.00006 7502-1 0.15 0.13 46.3509--------------------------------------------------------------------------------------二维基线向量残差--------------------------------------------------------------------------------------No. From To V_DX(cm) V_DY(cm)--------------------------------------------------------------------------------------1 7502-1 CPI1056-1 0.19 -0.222 7502-1 CPI1057 -0.02 0.103 7502-1 JM7502 -0.36 -0.154 7502-1 JM7502 -0.48 -0.365 CPI1057 CPI1056-1 0.07 -0.386 CPI1060 7502-1 -0.18 -0.147 CPI1060 JM7502 -0.58 -0.218 CPI1061 7502-1 -0.28 -0.099 CPI1061 CPI1060 0.12 0.0910 CPI1061 JM7502 -0.62 -0.3611 JM7502 CPI1056-1 0.39 0.1212 JM7502 CPI1057 0.42 0.46--------------------------------------------------------------------------------------二维基线向量可靠性--------------------------------------------------------------------------------------No. From To 内部可靠性DX DY--------------------------------------------------------------------------------------1 7502-1 CPI1056-1 0.77 0.722 7502-1 CPI1057 0.81 0.743 7502-1 JM7502 0.81 0.794 7502-1 JM7502 0.81 0.895 CPI1057 CPI1056-1 0.99 0.996 CPI1060 7502-1 0.68 0.627 CPI1060 JM7502 0.86 0.918 CPI1061 7502-1 0.70 0.679 CPI1061 CPI1060 0.97 0.9410 CPI1061 JM7502 0.85 0.8811 JM7502 CPI1056-1 0.90 0.8712 JM7502 CPI1057 0.92 0.90累计内部可靠性 10.08 9.92平均内部可靠性 0.84 0.83内部可靠性总和 20.00内部可靠性均值 0.83--------------------------------------------------------------------------------------------------平差后方位角、边长及精度--------------------------------------------------------------------------------------------------No. FROM TO A(dms) MA(s) S(m) MS(cm) MS:S ppm--------------------------------------------------------------------------------------------------1 7502-1 CPI1056-1 269.31517 0.11 2656.0371 0.14 1/1887000 0.532 7502-1 CPI1057 287.28088 0.11 2746.6937 0.13 1/2036000 0.493 7502-1 JM7502 202.36235 0.45 613.8244 0.15 1/419000 2.384 7502-1 JM7502 202.36235 0.45 613.8244 0.15 1/419000 2.385 CPI1057 CPI1056-1 182.25533 0.00 847.0366 0.00 0 0.006 CPI1060 7502-1 279.17353 0.08 3737.82320.14 1/2719000 0.377 CPI1060 JM7502 270.32215 0.00 3924.89280.00 0 0.008 CPI1061 7502-1 291.19253 0.07 4424.53090.13 1/3308000 0.309 CPI1061 CPI1060 336.42157 0.00 1094.54740.00 0 0.0010 CPI1061 JM7502 283.27052 0.00 4480.49760.00 0 0.0011 JM7502 CPI1056-1 282.41238 0.00 2480.58650.00 0 0.0012 JM7502 CPI1057 300.15552 0.00 2760.28320.00 0 0.00--------------------------------------------------------------------------------------------------最弱边--------------------------------------------------------------------------------------------------No. FROM TO A(dms) MA(s) S(m)MS(cm) MS:S ppm--------------------------------------------------------------------------------------------------3 7502-1 JM7502 202.36235 0.45 613.82440.15 1/419000 2.38--------------------------------------------------------------------------------------------------相邻点坐标分量及其相对中误差--------------------------------------------------------------------------------------------------No. FROM TO dX(m) dY(m) mdX(cm) mdY(cm)--------------------------------------------------------------------------------------------------1 7502-1 CPI1056-1 -21.7396 -2655.9481 0.14 0.14 合格2 7502-1 CPI1057 824.5344 -2620.0132 0.14 0.14 合格3 7502-1 JM7502 -566.6621 -235.9545 0.14 0.14 合格4 7502-1 JM7502 -566.6621 -235.9545 0.14 0.14 合格5 CPI1057 CPI1056-1 -846.2740 -35.9349 0.00 0.00 合格6 CPI1060 7502-1 603.6056 -3688.7644 0.14 0.14 合格7 CPI1060 JM7502 36.9435 -3924.7189 0.00 0.00 合格8 CPI1061 7502-1 1608.9216 -4121.6313 0.14 0.14 合格9 CPI1061 CPI1060 1005.3160 -432.8669 0.00 0.00 合格10 CPI1061 JM7502 1042.2595 -4357.5858 0.00 0.00 合格11 JM7502 CPI1056-1 544.9225 -2419.9936 0.00 0.00 合格12 JM7502 CPI1057 1391.1965 -2384.0587 0.00 0.00 合格--------------------------------------------------------------------------------------------------。
误差椭圆
§6-1 概 论在测量中，点P 的平面位置常用平面直角坐标P P y x ,来确定。
为了确定待定点的平面直角坐标，通常由已知点与待定点构成平面控制网，并对构成控制网的元素（角度、边长等）进行一系列观测，进而通过已知点的平面直角坐标和观测值，用一定的数学方法（平差方法）求出待定点的平面直角坐标。
由于观测条件的存在，观测值总是带有观测误差，因而根据观测值通过平差计算所获得的待定点的平面直角坐P P y x ~,~面位置并不是 P 点的真位置，而是最或然点位， 记为 P '，在 P 和 P '对应的这两对坐标之间 存在着坐标真误差 x∆和 y∆。
由图6-1知⎭⎬⎫-=∆-=∆P P y P P x y y x x ˆ~ˆ~ (6-l-1) 由于x ∆和y ∆的存在而产生的距离P ∆称为 P 点的点位真误差，简称真位差。
由图6-1知222yxP∆+∆=∆222y xPσσσ+=(6-1-2)2．点位真误差的随机性P 点的最或然坐标Px ˆ和P yˆ是由一组带有观测误差的观测值通过平差所求得的结果，因此，它们是观测值的函数。
设P xˆ和P y ˆ与观测值向量L 之间的线性函数关系为 ⎭⎬⎫++=++=00ˆˆββααL y y L x xA P A P（6-1-3）设有两组不同的观测值向量1L 、2L ，分别代入式（6-1-3）可得010111ˆˆββαα++=++=L y yL x xA P A P 和020222ˆˆββαα++=++=L y yL x xA P A P对于同一控制网而言，如果观测量相同（如同样的角度、边长等），采取同样的平差方法，则式中的00βαβα、、、是不变量，但观测值向量1L 、2L 不会相等，因此21ˆˆP P x x ≠、21ˆˆP P y y ≠。
可见，随着观测值L 的不同，P x ˆ和P y ˆ也将取得不同的数值。
但P 点的真坐标P x ~和P y ~是唯一的，由式(6-l-1)、(6-l-2)知，就会出现不同的x ∆和y∆值以及P∆，所以说点位真误差随观测值不同而变化，即点位真误差具有随机性。
测量平差学习情境6 误差椭圆
教学内容 主要介绍点位真误差和点位误差、任意方向上 的位差、待定点的误差曲线与误差椭圆以及点与点
能正确陈述点位真误差和点位误差及其计算方 法；能正确陈述任意方向上的位差及其位差的极值； 能正确陈述误差曲线和误差椭圆；能基本正确陈述 相对误差椭圆。
1
技能目标 能正确计算点在任意方向上的位差，能正确计 算位差的极值方向和极值；能正确绘制误差曲线和 误差椭圆，能根据误差椭圆求点在任意方向上的位 差，能计算点与点之间的相对误差椭圆参数并绘制 相对误差椭圆。
13
2. 点位任意方向上的位差σψ：从椭圆的中心作 方向线，然后再作该方向线的垂线（要与椭圆上一 点相切），则垂足到椭圆中心的长度便是点位在该 方向上的位差（见图6-8），图中线段OD的长度就 等于该方向上的位差，即σψ=OD
图6-8
14
子情境3 相对误差椭圆 一、两点之间相对位置的精度 1.相对位置的表示 两点相对位置可用其两点的坐标差来表示,即 用矩阵表达为
11
②确定点位中误差。 ③待定点P至任意已知三角点（视其无误差）的边 长中误差。 ④待定点P至任意已知三角点（视其无误差）的方 位角中误差。
12
误差曲线不是一种典型曲线，作图也不方便， 因此降低了它的实用价值。但其形状与以E、F为 长、短半轴的椭圆很相似。在以xe、ye为坐标轴 的坐标系中，该椭圆的方程为 1.误差椭圆的绘制 误差椭圆是一种规则图形,作图比较容易。因 此，实际应用中常以E、F为长、短半轴来绘制标 准的椭圆来代替相应的误差曲线，用来计算待定点 在各方向上的位差，故称该椭圆为误差椭圆。
3
子情境1 点位真误差及点位误差 一、点位真误差 在测量工作中，通过野外所进行的一系列的观 测，然后对观测数据进行平差处理便可得到点的平 面坐标平差值(x^，y^）。 但是，观测值总是带有观测误差的，而由观测 值所计算的平差值虽然较观测值更合理、可靠，但 是，它是不可能消除误差的，即待定点坐标的平差 值(x^,y^），不是待定点坐标的真值(x，y)，这两 者之间是有差异的。
科傻GPS平差软件说明书
科傻系统(COSA)系列软件GPS工程测量网通用平差软件包（CosaGPS V5.1）版权所有不得翻录Email:目录2.2.1新建 ..............................2.2.2打开 ..............................3.“GPS数据处理”下拉菜单.......................3.1 已知数据................................3.4.2 输出用户自定义任意两点相对精度 ...3.5 椭球面上三维平差........................5.3 网图显绘................................5.4 贯通误差影响值计算......................6.7 高程面坐标变换..........................7.“帮助”下拉菜单..............................１.简介种坐标转换。
1.2 整体性好式或文本方式进行数据录入，大部分操作采用“傻瓜“式选项。
对于输入量较少的已知数据和参数，采用表格方式输入；对于大批量的数据，则采用文件方式输入。
表格方式输入时，屏幕上显示格式如图1.1。
图1.1 输入数据表格表格中各列的宽度可以改变，将鼠标移到表格中各列标题的结合处，按下左键拖动，调节到合适的宽度即可。
表格中的行数是不受限制的，输满后将向上滚动，底下弹出新的空白行。
表格式输入的数据被保存到文本文件（文件名参见2.1），用户也可直接对相应的文本文件（例如：工程名. GPS3dKnownXYZ）进行修改，重新进入表格后，表格中的数据将自动进行更新。
表格中各列的宽度可以改变，将鼠标移到表格中各列标题的结合处，按下左键拖动，调节到合适的宽度即可。
表格中的行数是不受限制的，输满后将向上滚动，底下弹出新的空白行。
表格式输入的数据被保存到文本文件（文件名参见2.1），用户也可直接对相应的文本文件（例如：工程名. GPS3dKnownXYZ）进行修改，重新进入表格后，表格中的数据将自动进行更新。
TGO的使用方法
a． 确认选中了b2点。
b． 在左侧小窗口中，点击b2点名前的+号，显示出该点的所有观测量。
c． 选（b2—a2）基线。
另外一种直接方法是：在屏幕上直接用鼠标选中b2—a2基线
2、显示GPS基线处理报告。选报告。选 报告 / 基线处理报告。在报告的后面是卫星残差图，本例中，显示的是卫星20（SV 20）的数据残差，卫星的残差不宜大于0.15周，可以使用时序器来检查这些卫星的数据。
现在可以进行平差。
6、最小约束网平差
最小约束网平差是整个网中之有一个固定控制点的平差，这一部分讲述如下：
u 显示点的误差椭圆：
u 在网中固定一个控制点：
u 实施最小约束网平差
u 查阅网平差结果
6.1、显示误差椭圆
要显示误差椭圆开关按钮，做下面其中之一：
4.2评估解算结果
要查看a2点：
1、双击在屏幕上的点，出现下面的属性窗口：
在此窗口中可以查看该点的坐标，测量时间，天线高度，基线解算的结果及质量指标等。
4.2、查阅GPS 基线处理报告
要查阅从b2点到a2点的基线处理报告：
1、选择从b2点到a2点的基线。在属性窗口中的方法如下：
4、 点 确认。
5、 要进行重新平差，选 平差 / 平差。
6、 查看网平差报告，选 报告 / 网平差报告。
查看报告，可以看到网参考因子在1左右，差方检验通过，坐标变化量均为0.000
7、 关闭报告。
8、 在工具栏，点 误差椭圆 来关闭椭圆显示。
显示GPS环的闭合差报告，请：
1、选 报告 / GPS闭合环报告。
显示GPS环的闭合差报告。在 总结 中，显示闭合环总数，通过的和超限的闭合环个数。
误差理论与平差基础_第10章_误差椭圆
Q
=
æ èç
2 0.5
0.5 3
öø÷(dm2 / ('')2 )
单位权方差 0 0.5''
待定点P点到已知点A的距离为6.45km，方位角为45°，求P 点在AP方向的纵向误差和横向误差及AP边的边长相对中误 差。
s
2 j
=
s
2 0
éëQxx
cos2
j
+
Qyy
sin2
j
+
Qxy
sin
2jùû
Qxn y1 Qyn y1
Qx1xs Qy1xs
Qxn xs Qyn xs
Qx1ys
ö ÷
Qy1ys
÷ ÷
Qxn ys
÷ ÷
Qyn ys
÷ ø
二、点位任意方向的位差
Dj = pp¢¢ + p¢¢p¢¢¢ = Dx cosj + Dysinj
x
Dy
j
Q cos
s
in
Qxx Qyx
QFF
=
1 2
(Qxx
+ Qyy
-
K)
= 100.015
tan j E
=
QEE - Qxx Q
= -0.276
jE =164°33¢ / 344°33¢
xy
tan j F
=
Q FF
-
Q xx
Qxy
= 3.618
jF = 74°33¢ / 254°33¢
B P
b
A xA = 4578.67 yA = 3956.74 aAB = 345°18¢12¢¢
误差椭圆.
仿式(3.5-3)可得
2 P
(x23 .5-4)y2
这说明，尽管点位真误差△P
在不同坐标系的两个坐标轴上的投
影长度不等，但点位方差 总P2 是等 于两个相互垂直的方向上的坐标方
差之和，即它与坐标系的选择无关。
图3.5-2
如果再将点P的真位差△P投影于AP方向和垂直于AP的
方向上，则得 s和 （见u 图3.5-1）， 、s 为点u 的纵向误差和 横向误差，此时有
(2)计算P2点的误差椭圆的元素
由
tan 20
2Qxˆ2 yˆ2 Qxˆ2 - Qyˆ2
2 0.2106 -1.1353 0.4912 - 0.8624
得
= E624.33
14
14
误差椭圆
K2
(Qxˆ2
- Qyˆ2
)2
4Q 2 xˆ2 yˆ2
0.561
E2
1 2
给出后，可根据这个图得到坐标平差值在任一方向的位差大
小。如图3.5-6为控制网中P点的点位误差曲线，A、B、C为已
知点。由图3.5-6可知，
，
，
xP
Pa
，
yP Pb
E Pc E
F Pd F
由图还可得到坐标平差值函数的中
误差。例如要想得到平差后方位角
垂直P的A 于中P误A方差向上，的可P位A 先差从Pg图，中这量是出 PA
2
2
误差椭圆
知识准备
1.点位真误差 在测量中，为了确定待定点的平面直角坐标，通常需进
行一系列观测。由于观测值总是带有观测误差，因而根据观
测 而值不，是通 待过 定平点差 坐计 标算 的所 真获 值得~x，的是~y。待定点坐标的平差值 xˆ ， yˆ，
误差理论与平差基础-误差椭圆.
误差理论与测量平差基础
—误差椭圆
本章教学内容
7.1 点位误差 7.2 误差曲线与误差椭圆 7.3 相对误差椭圆
1）直观：把各方向的位差清楚地图解出来了； 2）任意方向ψ上的向径0P就是该方向的位差σψ。 3）图形是关于E轴和F轴对称的。
2、误差曲线图的用途
B
A 图解点位点位中误差、任意方向上的位差等。
7.2.2 误差椭圆
误差曲线优点：
能直观地反映点位在任意方向上的位差； 能根据图找出点位在各个方向上的位差。
极大、极小方向的计算公式：
tan E
QEE Qxx Qxy
tan F
QFF Qxx Qxy
7.1.4 以E、F表示任意方向ψ上的位差
说明：任意方向ψ指以E轴为起算的方向！（与φ不同。）
E
∆F
P’
ψ
∆E
∆P
P P’’ ∆ψ P’’’
F
由上图，可得：
cos E sin F
表示的？
两点的坐标差： 坐标差的协因数：
xik xk xi yik yk yi
Qxx Qxk xk Qxixi 2Qxk xi
Qyy Qyk yk Qyi yi 2Qyk yi
Q Q Q Q Q xy
xk yk
xk yi
(Qxx Qyy ) sin 2 2Qxy cos 2
0
设φ 0为位差的极值方向，则有：
第五章GPS信号的误差
式中： fs ---卫星时钟的频率； f ---同类静止的时钟频率 ；
Vs ---卫星的运行速度； C0---真空光速。
若用GPS卫星的运行速度Vs=3874m/s，而 C0=299792458m/s，则可算得GPS卫星时 钟相对于地面同类时钟的频率之差是
测距码观测值:
d ion
C
40 .28 f2
s Ne ds
载波相位观测值:
dion
C
40 .28 f2
s Ne ds
将二者结合处理可基本消除电离层折射误
差的影响,使单频GPS接收机的测程扩大到
200Km左右。
（5）选择有利观测时段
时 延
0 4 8 12 16 20 24 地方时
对流层折射误差
3.减弱电离层影响的措施
（1）利用双频观测（对于双频接收机） 现令: dion= A / f2 对于双频接收机,可以同时接收两个载波信号。 则有: S = ρ1 + A / f12 ；
S = ρ2 + A / f22 ；
式上中的：PS码为信星号站进的行理测论量距分离别, ρ获1 得和的ρ伪2 为距对观两测个值载。波
GPS信号电离层折射率为： nGPS= 1+40.28Nef -2
GPS信号在电离层中传播速度为： Vg=C0/nGPS=C0(1-40.28Nef-2)
若伪距测量中信号的传播时间为Δt,那么 S=vg Δt= C0(1-40.28Nef-2) Δt = C0 Δt- C0 40.28Nef-2 Δt
但由于其能反映全球的平均状况,与各地的 实际情况必然会有一定的差异,所以其改正 效果仅能改正电离层折射误差的 75﹪左右。
GPS基线向量网平差VB程序设计
GPS基线向量网平差程序设计前言GPS技术以其定位精度高，观测自动化,不需测站间通视及网型与精度关系不大的优势,已成为建立城市及工程控制网的主要技术手段之一。
而与常规地面网相比,GPS控制网的数据处理有其自身的特点，由于基线向量是不可独立于坐标系而存在的特殊观测值，所以在平差时或平差后必须转入测区所在的坐标系统。
本论文讨论了GPS基线向量的转换和平差问题及工程控制测量实用的方法，并运用VB程序设计语言完成了大地空间直角坐标向大地坐标的转换、大地坐标向高斯平面坐标的转换、二维基线向量网平差的功能。
1GPS原理1.1 GPS的简介全球定位系统（全局位置系统，GPS）是美国从上世纪70年代开始研制，历时20年，耗资200亿美元，于1994年全面建成的利用导航卫星进行测时和测距，具有在海、陆、空进行全方位实时三维导航与定位能力的新一代卫星导航与定位系统。
它是继阿波罗登月计划、航天飞机后的美国第三大航天工程。
如今，GPS已经成为当今世界上最实用，也是应用最广泛的全球精密导航、指挥和调度系统。
它主要由三大子系统构成：空间卫星系统、地面监控系统、用户接收系统。
1.2 GPS定位原理GPS系统采用高轨测距体制，以观测站至GPS卫星之间的距离作为基本观测量。
为了获得距离观测量，主要采用两种方法：一是测量GPS卫星发射的测距码信号到达用户接收机的传播时间，即伪距测量；一是测量具有载波多普勒频移的GPS卫星载波信号与接收机产生的参考载波信号之间的相位差，即载波相位测量。
采用伪距观测量定位速度最快，而采用载波相位观测量定位精度最高。
通过对４颗或４颗以上的卫星同时进行伪距或相位的测量即可推算出接收机的三维位置。
按定位方式，GPS 定位分为单点定位和相对定位（差分定位）。
单点定位就是根据一台接收机的观测数据来确定接收机位置的方式，它只能采用伪距观测量。
相对定位（差分定位）是根据两台以上接收机的观测数据来确定观测点之间的相对位置的方法，它既可采用伪距观测量也可采用相位观测量。
各种位置椭圆编程
各种位置椭圆编程1、毛坯Ф30，加工右半个椭圆（Z向有偏心）。
主程序：WGF1.MPFG95G23G90G71T1D1M03S500F0.3G00X32Z2R20＝14.2（取最大切削余量30的一半为15，考虑到每次单边1mm的切削深度及预留0.2mm 的单边精加工余量，则取14.2）MA1：G158 X＝R20WGF2R20＝R20－1IF R20>=0.2 GOTOB MA1G00X32Z2M03S800F0.15G158R20=0WGF2G00X60Z80M05M02子程序：WGF2.SPFR1=20(椭圆长半轴) R2=12（短半轴）R3=20(以椭圆中心为坐标原点的坐标系O1X1Z1中的坐标，是所加工椭圆轮廓起始点的Z坐标)MA2：R4=R2*SQRT(R1*R1-R3*R3)/R1 (通过本公式算出对应的椭圆坐标系O1X1Z1中的X 值)G01 X=2*R4 Z=R3-20 (将上述O1X1Z1坐标系中的X、Z转换到工件坐标系OXZ中，进行直线插补)R3=R3-0.05IF R3>=0GOTOB MA2 (R3=0是在椭圆坐标系O1X1Z1中的坐标，是所加工椭圆轮廓终点的Z坐标)G91X2G90Z2M022、毛坯Ф30，加工右小半个椭圆（Z向有偏心）。
主程序：WGF1.MPFG95G23G90G71T1D1M03S500F0.3G00X32Z2R20＝6.2（取最大切削余量的一半7.06，考虑到每次单边1mm的切削深度及预留0.2mm的单边精加工余量，则取6.2）MA1：G158 X＝R20WGF2R20＝R20－1IF R20>=0.2 GOTOB MA1G00X32Z2M03S800F0.15G158R20=0WGF2G00X60Z80M05M02子程序：WGF2.SPFR1=20(椭圆长半轴) R2=12（短半轴）R3=15(以椭圆中心为坐标原点时椭圆轮廓的Z坐标起始点)MA2：R4=R2*SQRT(R1*R1-R3*R3)/R1G01X=2*R4 Z=R3-15R3=R3-0.05IF R3>=0GOTOB MA2G91X2G90Z2M023、毛坯Ф30，加工右半个椭圆（X、Z向都有偏心）。