前言总体感觉测量平差这门课还是不是很容易,需要比较深刻的理解概念还要会相关的推导,比如精密度精确度准确度粗值观测值平差值等。主体内容就是间接平差,条件平差。在这两者基础之上,间接平差参数选择比较多,就出现附有限制条件的间接平差。如果条件平差又有参数,就称作附有参数的条件平差。然后是误差椭圆,分析误差分布规律的,哪里误差大,哪里误差小。最后为了评定平差结果或者精度的好坏,又有一部分参数检验和假设检验的内容,不过基本都是概率论的内容,比如U检验,T检验等,所以还是需要熟悉一些参数的构造。第一章观测误差的分类及其处理给出误差分类的表达式,粗差、系统误差和偶然误差的定义。系统误差:在相同的观测条件下作
题目如下: 如下图所示为一简单GPS网,用两台GPS接收机观测,测得5条基线向量,每一条基线向量中三个坐标差观测值相关,由于只用两台GPS接收机观测,所以各观测基线向量互相独立。观测基线向量信息见表1。假定1号点为起算点坐标信息表2。 表1 GPS网平差观测数据及已知方差阵表2 GPS网平差起算数据点号XYZLC01-1974638.73404590014.81903953144.9235要求:1)基于Matlab或其他编程语言(如C++等)编程实现该GPS网间接平差过程通用程序,包括误差方程、法方程的组成与解算。得出平差后各基线向量观测值的平差值及各待定点的坐标平差值;评定各待定点坐标平差值
条件平差P=input('请输入权阵');A=input('请输入矩阵A');w=input('请输入矩阵w');n=A*inv(P)*A.';k=-inv(n)*w;v=inv(P)*A.'*kdisp('单位权中误差的估值为');sqrt(v.'*P*v/size(w,1))disp('观测值平差值的协因数阵为')Qh=inv(P)-inv(P)*A.'*inv(n)*A*inv(P)f=input('请输入平差值函数的权函数')disp('平差值的权倒数为');q=f*Qh*f.'有参数的条件平差P=input('请输入权阵');A=input('请输入矩阵A');B=input('请输