测量平差实习程序设计报告

X=-inv(B'*P*B)*B'*P*l;
V=B*X+l;
[row,col]=size(B);
mo=sqrt(V'*P*V)/(row-col);
M=mo^2*inv(B'*P*B);
msgbox('计算成功了!')
function [ok,ptTab,lineTab]=getH0(ptTab,lineTab)
计算改正数
计算平差值
，得：
又
可得：
精度评定
其中 ,
又
且
故有：
平差模型的正确性检验
原假设和备选假设为
又
计算统计量
以自由度 =4， =0.05查得 的分布表：
可见， 在 内，该平差模型正确，平差结果可用。
3.2
边角网平差资料如下：
上图为一平面控制网，A、B、C、D、E为已知点，P1P2为待定点，同精度观测了9个角度，L1– L2，测角中误差为2.5″；观测了5条边长，L10– L14，观测结果及中误差列于表中，按间接平差法对该控制网进行平差。
近似方位角、近似边长
测站
照准点
近似方位角
°′″
近似边长（m）
A
B
35 00 17.11
350 54 27.0
2185.070
1484.781
B
C
A
0 52 06.0
77 43 46.73
170 54 27.0
3048.650
1522.827
1484.781
C
B
180 52 06.0
152 06 45.1
1009.021
E
D
289 31 44.
确定角和边的权
取单位权中误差 ，则角度观测值的权位
各导线的权为： 。由此可得出：
误差方程系数、常数项、权值表
P
角
1
2
3
4
5
6
7
8
9
-0.5415
1.3235
-0.7820
-1.3235
0.3130
1.0105
1.6691
-1.3561
0.7732
-0.2879
break;
end
end
end
if beginIndexh~=0 &&endIndexh==0%如果起点高程非零，终点高程为0
lineIndex=i;
return;
else
%if not(isContr(beginIndex,ptTab))&&isContr(endIndex,ptTab)
ifbeginIndexh==0 &&endIndexh~=0
msgbox,'由起算边计算失败!!'
end
if not(hOk(ptTab))
[ok,ptTab,lineTab]=getH0(ptTab,lineTab);
end
ok=true;
function [lineIndex]=getBegin(ptTab,lineTab);
[rowNo1,colNo1]=size(lineTab);
对平差模型进行正确性检验
已知测角中误差：
后验方差检验假设为：
已算得
计算统计量
可见 不在区域 内，应拒绝H。
该平差模型不正确。
4
4.1
clear;
load('ptTabl','ptTabl');
load('lineTabl','lineTabl');
ptTab=ptTabl;
lineTab=lineTabl;
（3）p2和p3点之间平差后高差值 的中误差；
（4）平差模型正确性检验（四等水准测量每公里高差观测中误差为±5毫米）。
3.1.1
求平差值
列条件方程：
由题意可知：n=7，t=4，r=n-t=3.观测方程为
线性化得条件方程：
其中系数阵A= ，W=
组建法方程：
令1km的观测高差为单位权观测，即 ，又 。即
°′″
编号
观测值L
°′″
编号
观测值/m
中误差/cm
1
44 05 44.8
6
74 22 55.1
10
2185.070
3.3
2
93 10 43.1
7
127 25 56.1
11
1522.853
2.3
3
42 43 27.2
8
201 57 34.0
12
3082.621
4.6
4
76 51 40.7
9
168 01 45.2
3048.650
3082.631
C
B
A
332 06 45.1
257 43 46.73
215 0017.1
99 32 25.42
3082.631
1522.827
2185.016
1499.933
D
279 32 25.2
121 29 59.7
1499.931
1009.021
D
E
301 29 59.7
109 31 44.9
13
1500.017
2.2
5
28 45 20.9
14
1009.021
1.5
表4
平差如下：
求待定点近似坐标
由题意知：n=14，必要观测数t=4，可设 、 点坐标为参数，记为 、 、 、 。
由前方交会得P1点坐标：
故P1点近似坐标为
P2点的近似坐标为：
列误差方程
由题意可知列如下方程：
由 、 点近似坐标以及观测数据可计算近似方位角近似边长，现列表如下
总体思路：现有等级导线网的全部观测数据及网型、起算数据。要求对该导线网，用间接进行严密平差，并进行精度评定及平差模型的正确性检验。
边角网的间接平差：
①列观测值平差值方程、误差方程、法方程；
②解算基础方程，求出观测值的平差值、待定点的坐标平差值；
③评定观测值平差值的精度，坐标平差值的精度，求点位误差。
已知数据
表3
点
坐标/m
至点
边长/m
坐标方位角
X
Y
A
3143.237
5260.334
B
1484.781
350 54 27.0
B
4609.361
5025.696
C
3048.650
0 52 06.0
C
7657.661
5071.897
D
D
4157.197
8853.254
E
109 31 44.9
角
边
编号
观测值L
高差观测值(m)
对应线路长度(km)
已知点高程（m）
h1= 1.359
h2= 2.009
h3= 0.363
h4=-0.640
h5= 0.657
h6= 1.000
h7= 1.650
1
1
2
2
1
1
2
H1= 35.000
H2= 36.000
要求：按条件以及间接平差法分别求：
（1）待定点高程平差值；
（2）待定点高程中误差；
0.1658
0.5763
0.9772
0.4677
-0.9862
0
0
0
-0.1658
0.5225
0
0
0
0.9862
-0.8526
3.8
-0.5
-2.9
11.5
0
0.57
1.18
0.30
1.29
2.78
组成的法方程：
。法方程为
由矩阵运算程序可得：
又 并由矩阵运算程序可得
其中
精度计算：
单位权中误差，即测角中误差
《测量平差课程设计》报告
设计名称：测量平差课程设计
学生学号：*********
学生班级：0614121
******
专 业：测绘工程
指导教师：李军杰马晓玉
时 间：2014.12.22—2014.12.26
2014年 12 月 19 日
1
《测量平差》是一门理论与实践并重的课程，测量平差课程设计是测量数据处理理论学习的一个重要实践环节，是在学生学习了专业基础理论课《误差理论与测量平差基础》课程后进行的一门实践课程，其目的是增强学生对测量平差基础理论的理解，牢固掌握测量平差的基本原理和公式，熟悉测量数据处理的基本原理和方法，灵活准确地应用于解决各类数据处理的实际问题，并能用所学的计算机基础知识，编制简单的计算程序。
[isok,ptTab,lineTab]=getH0(ptTab,lineTab);
%[ok,ptTab,lineTab]=getH0(ptTab,lineTab);
if not(isok)
%if not(ok)
msgbox('计算概略高程失败!')
return;
end
[P,B,l,ptTab,lineTab]=getV(ptTab,lineTab);
要求结算：
⑴ 待定点坐标平差值，点位中误差；
⑵ 最弱边边长中误差，边长相对中误差；
⑶ 待定点误差椭圆参数、相对误差椭圆参数、绘出误差椭圆及相对误差椭圆、图解求出P1 P2点点位中误差、边长相对中误差（与计算比较）、最弱边方位角中误差。
⑷ 对平差模型进行正确性检验；