湖北国土资源资源职业学院
测
量
平
差
实
习
报
告
姓名:刘彬
班级:测量1002班
时间:2012年4月27
实习老师:毕靖
目录
实训一、误差传播应用示例………………………………….
实训二、条件平差方法应用………………………………….
2.1有关距离的平差:…………………………………………………2.2有关水准的平差……………………………………………………2.3有关角度的平差……………………………………………………
实训三、间接平差方法应用……………………………………3.1有关距离的平差…………………………………………………….
3.2有关水准的平差…………………………………………………….
3.3有关角度的平差…………………………………………………….. 实训四、平面控制网平差………………………………………
4.1导线网平差………………………………………………………….
4.2三角网平差……………………………………………………………实训五、高程控制网平差……………………………………….
5.1水准网平差…………………………………………………………5.2三角高程网平差……………………………………………………实训六、平差易软件的应用……………………………………………
6.1平面控制网平差………………………………………………..
6.2高程控制网平差………………………………………………..
实训一误差传播应用示例
1、解:D=[2 -1;-1 3]
D =
2 -1
-1 3
>> B=2/0.4
B =
5
>> Q=D/B
Q =
0.4000 -0.2000
-0.2000 0.6000
>> P=inv(Q)
P =
3 1
1 2
>> p1=1/0.4
p1 =
2.5000
>> p2=1/0.6
p2 =
1.6667
2、解:有已知条件可知:01
m m''
==p1=p2=p3=1
列线性函数:L1^+L3^=L2^
V=-(L1+L3-L2)/3 K=[1 -1 1]
K =
1 -1 1
>> P=[1 0 0;0 1 0;0 0 1]
P =
1 0 0
0 1 0
0 0 1
>> Q=inv(P)
Q =
1 0 0
0 1 0
0 0 1
>> QF=K*Q*K'
QF =
3
3、P =
1 0 0 0 0 0 0 0
0 1 0 0 0 0 0 0
0 0 1 0 0 0 0 0
0 0 0 1 0 0 0 0
0 0 0 0 1 0 0 0
0 0 0 0 0 1 0 0
0 0 0 0 0 0 1 0
0 0 0 0 0 0 0 1
>> Q=inv(P)
Q =
1 0 0 0 0 0 0 0
0 1 0 0 0 0 0 0
0 0 1 0 0 0 0 0
0 0 0 1 0 0 0 0
0 0 0 0 1 0 0 0
0 0 0 0 0 1 0 0
0 0 0 0 0 0 1 0
0 0 0 0 0 0 0 1 A=[1 -1 1 0 0 0 0 0;0 0 -1 1 0 0 0 1;0 0 0 0 1 1 -1 -1;]
A =
1 -1 1 0 0 0 0 0
0 0 -1 1 0 0 0 1
0 0 0 0 1 1 -1 -1 >> QH=A*Q*A'
QH =
3 -1 0
-1 3 -1
0 -1 4
4、解:S=∏/4X3*X3+X1*X3+1/2*X=2506.858m2取微分Ds=x3*dx1+1/2*x3*dx2+(∏/2*x3+x1+x2/2)dx3
D=[0.0016 0 0;0 0.004 0;0 0 0.009]
D =
0.0016 0 0
0 0.0040 0
0 0 0.0090
>> A=[30 15 107]
A =
30 15 107
>> MS=sqrt(A*D*A')
MS =
10.2655
实训二条件平差方法应用
2.1有关距离的平差
1、解:1.r=2
2.P=C/S
P的矩阵=Q
3.列出2个条件方程
L1+L2-L4=0
L2-L3=0
W1=L1+L2-L4
W2=L2-L3
4.组成法方程
K=-inv(N)*W
把K代入改正数方程求改正数
V=Q*A'*K
5.求平差值
>> Q=[200.010 0 0 0;0 300.050 0 0;0 0 300.070 0;0 0 0 500.090]
Q =
200.0100 0 0 0
0 300.0500 0 0
0 0 300.0700 0
0 0 0 500.0900
>> A=[1 1 0 -1;0 1 -1 0]
A =
1 1 0 -1
0 1 -1 0
>> W=[200.010+300.050-500.090;300.050-300.070]
W =
-0.0300
-0.0200
>> L=[200.010;300.050;300.070;500.090]
L =
200.0100
300.0500
300.0700
500.0900
>> N=A*Q*A'
N =
1.0e+003 *
1.0001 0.3001
0.3001 0.6001
>> K=-inv(N)*W
K =
1.0e-004 *
0.2353
0.2156
>> V=Q*A'*K
V =
0.0047
0.0135
-0.0065
-0.0118
>> LL=L+V
LL =
200.0147
300.0635
300.0635
500.0782
2.2有关水准的平差
2、B=[1 0;0 1;0 1;1 1]
B =
1 0
0 1
0 1
1 1
>>l=[0;0;-20;-30]
l=
-20
-30
>> Q=[200.010 0 0 0 ;0 300.050 0 0;0 0 300.070 0;0 0 0 500.090]/100 Q =
2.0001 0 0 0
0 3.0005 0 0
0 0 3.0007 0
0 0 0 5.0009
>> P=inv(Q)
P =
0.5000 0 0 0
0 0.3333 0 0
0 0 0.3333 0
0 0 0 0.2000
>> N=B'*P*B
N =
0.6999 0.2000
0.2000 0.8665
>> C=B'*P*l
C =
-5.9989
-12.6640
>> x=-inv(N)*C
x =
4.7055
13.5293
L=[200.010 ;300.050; 300.070; 500.090]
L =
200.0100
300.0500
300.0700
500.0900
V=(B*x+l)/1000
V =
0.0047
0.0135
-0.0065
-0.0118
>> LL=L+V
LL =
200.0147
300.0635
300.0635
500.0782
3、Q=[3.1 0 0 0 0 0 0 0;0 9.3 0 0 0 0 0 0;0 0 59.7 0 0 0 0 0;0 0 0 6.2 0 0 0 0;0 0 0 0 16.1 0 0 0;0 0 0 0 0 35.1 0 0;0 0 0 0 0 0 12.1 0;0 0 0 0 0 0 0 9.3]
Q =
3.1000 0 0 0 0 0 0 0
0 9.3000 0 0 0 0 0 0
0 0 59.7000 0 0 0 0 0
0 0 0 6.2000 0 0 0 0
0 0 0 0 16.1000 0 0 0
0 0 0 0 0 35.1000 0 0
0 0 0 0 0 0 12.1000 0
0 0 0 0 0 0 0 9.3000
A=[1 -1 0 1 0 0 0 0;0 0 0 1 -1 0 -1 0;0 0 0 0 1 -1 0 -1;1 0 -1 0 0 1 0 0]
A =
1 -1 0 1 0 0 0 0
0 0 0 1 -1 0 -1 0
0 0 0 0 1 -1 0 -1
1 0 -1 0 0 1 0 0
W=[189.404+547.576-736.977;547.576-273.528-274.082;273.528-187.274-86.261;189.404-376.6 07+187.274]
W =
0.0030
-0.0340
-0.0070
0.0710
>> L=[189.404;736.977;376.607;547.576;273.528;187.274;274.082;86.261]
L =
189.4040
736.9770
376.6070
547.5760
273.5280
187.2740
274.0820
86.2610
>> N=A*Q*A'
N =
18.6000 6.2000 0 3.1000
6.2000 34.4000 -16.1000 0
0 -16.1000 60.5000 -35.1000
3.1000 0 -35.1000 97.9000
>> K=-inv(N)*W
K =
-0.0004
0.0010
-0.0000
-0.0007
>> V=Q*A'*K
V =
-0.0034
0.0036
0.0431
0.0041
-0.0173
-0.0245
-0.0127
0.0002
>> LL=L+V
LL =
189.4006
736.9806
376.6501
547.5801
273.5107
187.2495
274.0693
86.2612
3.2有关水准的平差
(1)B=[1 0 0 0;1 1 0 0;1 0 0 1;0 1 0 0;0 0 1 0;0 0 0 1;0 1 -1 0;0 0 1 -1]
B =
1 0 0 0
1 1 0 0
1 0 0 1
0 1 0 0
0 0 1 0
0 0 0 1
0 1 -1 0
0 0 1 -1
>>
l=[0;189.404+547.576-736.977;189.404+187.274-376.607;0;0;0;547.576-273.528-274.082;273.5 28-187.274-86.261]
l =
0.0030
0.0710
-0.0340
-0.0070
Q=[3.1 0 0 0 0 0 0 0;0 9.3 0 0 0 0 0 0;0 0 59.7 0 0 0 0 0;0 0 0 6.2 0 0 0 0;0 0 0 0 16.1 0 0 0;0 0 0 0 0 35.1 0 0;0 0 0 0 0 0 12.1 0;0 0 0 0 0 0 0 9.3]
Q =
3.1000 0 0 0 0 0 0 0
0 9.3000 0 0 0 0 0 0
0 0 59.7000 0 0 0 0 0
0 0 0 6.2000 0 0 0 0
0 0 0 0 16.1000 0 0 0
0 0 0 0 0 35.1000 0 0
0 0 0 0 0 0 12.1000 0
0 0 0 0 0 0 0 9.3000
P=inv(Q)
P =
0.3226 0 0 0 0 0 0 0
0 0.1075 0 0 0 0 0 0
0 0 0.0168 0 0 0 0 0
0 0 0 0.1613 0 0 0 0
0 0 0 0 0.0621 0 0 0
0 0 0 0 0 0.0285 0 0
0 0 0 0 0 0 0.0826 0
0 0 0 0 0 0 0 0.1075 >> N=B'*P*B
N =
0.4469 0.1075 0 0.0168
0.1075 0.3515 -0.0826 0
0 -0.0826 0.2523 -0.1075
0.0168 0 -0.1075 0.1528
>> C=B'*P*l
C =
0.0015
-0.0025
0.0021
0.0019
>> x=-inv(N)*C
x =
-0.0034
0.0041
-0.0173
-0.0245
>> V=B*x+l
V =
-0.0034
0.0036
0.0431
0.0041
-0.0173
-0.0245
-0.0127
0.0002
>> L=[189.404;736.977;376.607;547.576;273.528;187.274;274.082;86.261]
L =
189.4040
736.9770
376.6070
547.5760
273.5280
187.2740
274.0820
86.2610
LL =
189.4006
736.9806
376.6501
547.5801
273.5107
187.2495
274.0693
86.2612
(2)f1=[1 0 0 1]
f1 =
1 0 0 1
>> M1=sqrt(V'*P*V/4)*sqrt(f1*inv(N)*f1')
M1 =
0.0155
(3)f2=[0 1 0 -1]
f2 =
0 1 0 -1
M2=sqrt(V'*P*V/4)*sqrt(f2*inv(N)*f2')
M2 =
0.0154
3.3有关角度的平差
(1)B=[1 0 0;0 1 0;1 0 0;0 0 1;-1 1 0;1 0 -1;0 1 -1 ]
B =
1 0 0
0 1 0
1 0 0
0 0 1
-1 1 0
1 0 -1
0 1 -1
>>
l=[0;0;50+10.356-40-20.360;0;15-10.360-4.651;50+10.356-40-14.501-5.856;50+15-40-14.501-10 .5]
l =
-0.0040
-0.0110
-0.0010
>> P=[2 0 0 0 0 0 0;0 2 0 0 0 0 0;0 0 1 0 0 0 0;0 0 0 1 0 0 0;0 0 0 0 2 0 0;0 0 0 0 0 2 0;0 0 0 0 0 0 1]
P =
2 0 0 0 0 0 0
0 2 0 0 0 0 0
0 0 1 0 0 0 0
0 0 0 1 0 0 0
0 0 0 0 2 0 0
0 0 0 0 0 2 0
0 0 0 0 0 0 1
>> N=B'*P*B
N =
7 -2 -2
-2 5 -1
-2 -1 4
>> C=B'*P*l
C =
0.0160
-0.0230
0.0030
>> x=-inv(N)*C
x =
-0.0012
0.0040
-0.0004
>> X0=[50+10.356;50+15;40+14.501]
X0 =
60.3560
65.0000
54.5010
>> X=X0+x
X =
60.3548
65.0040
54.5006
(2)f=[-1 1 0]
f =
-1 1 0
V=B*x+l
V =
-0.0012
0.0040
-0.0052
-0.0004
-0.0057
-0.0019
0.0034
>> M0=sqrt(V'*P*V/4)
M0 =
0.0061
M=M0*sqrt(f*inv(N)*f')
M =
0.0031
实训四、平面控制网平差
4.1导线网平差
1、平面闭合差计算结果
导线线路号:1
线路点名: D C 5 4 3
2 B A
角度闭合差: -96845.3 (Sec)
X坐标闭合差: -272.4838 (M)
Y坐标闭合差: -48.6812 (M)
总长度: 1.4800 (KM)
相对精度: 1 : 5
平差后点位坐标
==============
--------------------------------------------------------------------- 点名 X坐标 Y坐标 H高程代码
A 203158.1560 441095.1440 0.0000 K
B 203020.3480 440950.1990 0.0000 K
C 203702.4370 439866.6010 0.0000 K
D 203756.0300 439951.0270 0.0000 K
2 203001.9171 440743.6601 0.0000 U
3 203020.4840 440540.8448 0.0000 U
4 203033.0389 440192.0917 0.0000 U
5 203274.5508 440039.5761 0.0000 U
平面控制网平差成果表
━━━━━━━━━━
平差类型:约束平差
网型:
总点数:8 已知点数:4
观测值总数:17 方向数:12 测边数:5 条件总数:0 方位角条件:0 边长条件:0 多余观测值总数:3
先验单位权中误差:3.54
后验单位权中误差:29518.17
2、平面闭合差计算结果
闭合环号:1
线路点名: A 2 3 4 5
6 7 8
角度闭合差: 11.0 (Sec) 超限!!!
限差: 5.7 (Sec)
边长闭合差: 0.0157 (M)
总长度: 2.4274 (KM)
平差后点位坐标
==============
--------------------------------------------------------------------- 点名 X坐标 Y坐标 H高程代码
A 2272.0450 5071.3300 0.0000 K
B 2487.6980 5279.8820 0.0000 K
2 2355.2080 5308.0136 0.0000 U
3 2382.0001 5566.078
4 0.0000 U
4 2028.5162 5607.4100 0.0000 U
5 1815.8455 5844.7170 0.0000 U
6 1684.405
7 5621.8020 0.0000 U
7 1701.3446 5352.8534 0.0000 U
8 1832.4959 5113.7452 0.0000 U
平面控制网平差成果表
━━━━━━━━━━
网型:
总点数:9 已知点数:2
观测值总数:25 方向数:17 测边数:8 条件总数:0 方位角条件:0 边长条件:0 多余观测值总数:3
先验单位权中误差:0.71
后验单位权中误差:1.94
3、平面闭合差计算结果
导线线路号:1
线路点名: E F 13 14 15
5 4 3 2 B
A
角度闭合差: 8.2 (Sec)
X坐标闭合差: 0.0538 (M)
Y坐标闭合差: 0.0385 (M)
总长度: 2.0814 (KM)
相对精度: 1 : 31476
==============================================
导线线路号:2
线路点名: C D 6 7 8
9 10 5 4 3
2 B A
角度闭合差: 20.0 (Sec)
X坐标闭合差: -0.0805 (M)
Y坐标闭合差: 0.1028 (M)
总长度: 2.4088 (KM)
相对精度: 1 : 18454
==============================================
导线线路号:3
线路点名: C D 6 7 8
9 10 5 15 14
13 F E
角度闭合差: 11.8 (Sec)
X坐标闭合差: -0.1707 (M)
Y坐标闭合差: 0.0692 (M)
总长度: 2.5994 (KM)
相对精度: 1 : 14112
===============================================
平差后点位坐标
==============
--------------------------------------------------------------------- 点名 X坐标 Y坐标 H高程代码
A 5256.9530 4520.0680 0.0000 K
B 5163.7520 4281.2770 0.0000 K
C 3659.3710 3621.2100 0.0000 K
D 4119.8790 3891.6070 0.0000 K
E 4581.1500 5345.2920 0.0000 K
F 4851.5540 5316.9530 0.0000 K
2 5018.2585 4092.138
3 0.0000 U
3 4939.6575 4243.7205 0.0000 U
4 4723.7606 4214.4844 0.0000 U