误差理论与测量平差课程设计报告Word文档下载推荐.docx
《误差理论与测量平差课程设计报告Word文档下载推荐.docx》由会员分享,可在线阅读,更多相关《误差理论与测量平差课程设计报告Word文档下载推荐.docx(14页珍藏版)》请在冰豆网上搜索。
![误差理论与测量平差课程设计报告Word文档下载推荐.docx](https://file1.bdocx.com/fileroot1/2022-10/9/da941f43-be56-4668-b0ca-8a974f9b47f4/da941f43-be56-4668-b0ca-8a974f9b47f41.gif)
14
10.参考资料·
15
一.课程设计的目的
误差理论与测量平差是一门理论与实践并重的课程。
该课程设计是测量数据处理理论学习的一个重要的实践环节,它是在学生学习了专业基础课“误差理论与测量平差基础”课程后进行的一门实践课。
其目的是增强我对误差理论与测量平差基础理论的理解,牢固掌握测量平差的基本原理和基本公式,熟悉测量数据处理的基本技能和计算方法,灵活准确的应用于解决各类数据处理的实际问题,并能用所学的计算机理论知识,编制简单的计算程序。
这次课程设计中所用的数学模型和计算方法是我们在误差理论与测量平差课程中所学的内容,所使用的C程序语言使我们在计算机基础课程中所学知识。
误差理论与测量平差课程设计是测量平差和计算机程序设计等课程的综合实践与应用,同时也为我们今后步入工作岗位打下了一定基础。
二.课程设计题目及相关要求
本题目为水准网严密平差及精度评定。
要求正确运用平差模型列出观测值条件方程、误差方程和法方程,并解算法方程,得出平差后的平差值及各待定点的高程平差值和各高程平差值的精度。
具体题目如下:
如图所示水准网,有2个已知点,3个未知点,7个测段。
各已知数据及观测值见下表
(1)已知点高程H1=5.016mH2=6.016m
(2)高差观测值(m)
(3)求各待定点的高程;
3-4点的高差中误差;
3号点、4号点的高程中误差。
(提示,本网可采用以测段的高差为平差元素,采用间接平差法编写程序计算。
)
端点号
高差观测值(m)
测段距离
序号
1—3
1.359
1.1
1
1—4
2.009
1.7
2
2—3
0.363
2.3
3
2—4
1.012
2.7
4
3—4
0.657
2.4
5
3—5
0.238
1.4
6
5—2
-0.595
2.6
7
相关要求:
本次课程设计重点是培养我们正确应用公式、综合分析和解决问题的能力,以及计算机编程能力。
另外它要求我们完成1-2个综合性的结合生产实践的题目。
如目前生产实践中经常用到的水准网严密平差及精度评定,边角网(导线)严密平差及精度评定等。
本题目需正确应用平差模型列出观测值条件方程、误差方程和法方程;
解算法方程,得出平差后的平差值及各待定点的高程平差值;
评定各平差值的精度和各高程平差值的精度。
三.设计思路
(1).根据网型确定相关的数据。
包括,已知水准点数2个,未知水准点数u=3,总点数n=7,总的观测高差段数,必要观测数t=3,多余观测数r=n-t=4。
因此先选定三个参数,即3、4、5点的最或然高程X3、X4、X5(X=X0+x,X30=6.375、X40=7.025、X50=6.611;
其中X0为参数的近似值,x为其改正值)为参数。
(2).列立条件方程.左侧为观测值(系数为1),右侧为参数和常数项,并进一步改化成误差方程,最终写成矩阵形式,即将每一个观测量的平差值分别表达成所选参数的函数,H1+h1=X3、H1+h2=X4、H2+h3=X3、H2+h4=X4、X3+h5=X4、X3+h6=X5、X5+h7=H2;
整理后得出误差方程,v1=x3、v2=x4、v3=x3-4、v4=x4-3、v5=-x3+x4-7、v6=-x3+x5-2、v7=-x5,即v=Bx-l的形式。
得到系数矩阵A和常数项。
(3)定权.令每千米的观测高差为单位权观测,即Pi=1/Si,从而可写出权阵P;
根据误差方程式又可得其系数矩阵B和自由项l,并由它们组成法方程NBBx-W=0(其中NBB=BTPB,W=BTPl),法方程的个数等于所选参数的个数。
(4)
.列立法方程,并解求法方程.由于该水准网间接平差误差方程个数为7个而未知数个数为10个,所列的误差方程是一组相容方程,有无数组解,所以必须在最小二乘原则(VTPV=min)的基础上利用拉格朗日乘数法求解.令F=VTPV-2KT(V-A
+L),分别对V和
求导,并令其导数为零,得到2VTP-2KT=0,ATK=0,将二式合并即得法方程:
ATPV=ATPA
-ATPL=0。
求出Naa=ATPA,W=ATPL,即得到相应的法方程。
求解法方程,得到
=N-1aaW加上Xi即可得到待定点的高程平差值,将
代入误差方程得到相应的V值,hi+Vi得到各段高差的平差值。
(5).精度评定.
计算单位权中误差的估值:
评定各待定点的高程中误差:
各待定点的精度即为:
评定高程平差值的精度:
四、程序流程图
根据题目列出条件方程并写成误差方程的形式V=Bx-l
↓
根据误差方程得到矩阵B、l进而写出BT
运用MatLab程序语言定义出权P
用MatLab程序求出参数的改正数x=NBB-1W,其中NBB=BTPB、W=BTPl
根据MatLab程序语言求Bx,进而由V=Bx-l写出各观测值的改正数
根据L=L+V求出各观测值的平差值
检验所求各值是否正确,若无误则往下进行,反之检查各步骤查出错误并改正
由程序计算VTP进而求出VTPV,求单位权中误差,再根据权函数式、协因数传播定律评定各观测值及所求高程的精度
五.程序源代码及说明
%--12-12-18上午9:
22--%
26--%
AJAHD
52--%
disp
Ha=5.016
Hb=6.016
H=[1.359;
2.009;
0.363;
1.012;
0.657;
0.238;
-0.595]
B=[100;
010;
100;
-110;
-101;
00-1]
l=[0;
0;
4;
3;
7;
2;
0.]
S=[1.1,1.7,2.3,2.7,2.4,1.4,2.6]
P=inv(diag(S'
))
x0=inv(B'
*P*B)*B'
*P*l
V=B*x0-l
H0=H+V/1000
n=6
t=2
del=sqrt(V'
*P*V/(n-t))
BB=B'
*P*B
Q34=[1-10]*inv(BB)*[1;
-1;
0]
3=[100]*inv(BB)*[1;
]
Q33=[100]*inv(BB)*[1;
Q44=[010]*inv(BB)*[0;
1;
m4=del*sqrt(Q44)
m3=del*sqrt(Q33)
m34=del*sqrt(Q34)
六.执行调试,得出计算结果
>
disp('
水准网间接平差事例'
水准网间接平差事例
Ha=5.0160%已知点高程,单位m
Ha=
5.0160
Hb=6.016%已知点高程,单位m
Hb=
6.0160
观测高差,单位m'
观测高差,单位m
H=[1.359;
H=
1.3590
2.0090
0.3630
1.0120
0.6570
0.2380
-0.5950
系数矩阵B'
系数矩阵B
B=[100;
B=
100
010
-110
-101
00-1
l=[0;
0]%常数项
l=
0
3
S=[1.1,1.7,2.3,2.7,2.4,1.4,2.6]%是线路长度
S=
1.10001.70002.30002.70002.40001.40002.6000
P=inv(diag(S'
))%定义观测值的权
P=
0.9091000000
00.588200000
000.43480000
0000.3704000
00000.416700
000000.71430
0000000.3846
参数的解'
参数的解
x0=inv(B'
x0=
-0.2427
2.8552
1.1423
误差V(mm),各待定点的高程平差值L1(m)'
误差V(mm),各待定点的高程平差值L1(m)
V=B*x-l%误差方程(mm)
V=
-4.2427
-0.1448
-3.9021
-0.6151
-1.1423
H0=H+V/1000%观测值的平差值
H0=
1.3588
2.0119
0.3588
1.0119
0.6531
0.2374
-0.5961
精度评定'
精度评定
n=6;
%观测值的个数
t=2;
%必要观测数
del=sqrt(V'
*P*V/(n-t))
del=
2.2248
NBB=B'
BB=
2.4748-0.4167-0.7143
-0.41671.37530
-0.714301.0989
Q34=[1-10]*inv(BB)*[1;
0]%h34=X1-X2
Q34=
0.9850
Q33=[100]*inv(BB)*[1;
]%H3=X1
Q33=
0.5307
Q44=[010]*inv(BB)*[0;
]%H4=X2
Q44=
0.7758
m4=del*sqrt(Q44)
m4=
1.9597
m3=del*sqrt(Q33)
m3=
1.6208
m34=del*sqrt(Q34)
m34=
2.2080
七.题目计算及精度评定过程
已知数据及观测值见下表:
(1)已