误差理论与测量平差基础（MATLAB）
误差理论与测量平差基础
（MATLAB） 实习报告
学 号： 姓 名： 班 级： 1420501Z 专 业： 测 绘 工 程 课程名称： 误差理论与测量平差基础 任课老师：
2016 年 5 月
1
误差理论与测量平差基础（MATLAB）
一.水准网间接平差
1.列出误差方程
设P1和P2点高程平差值为 和 ，相应的近似值取为
按已知数据及观测数据表列出观测方程后，将有关观测数据带入即得误差方程
V1= V2= V3= V4= V5= V6=
式中常数项以m为单位。 2.列出权函数式
P1至P2间高差平差值的权函数式为 3.组成法方程
以1km水准测量的观测高差为单位权观测值，各观测值互相独立，定权式为Pi=1/Si ， 得权阵为
p = 0.9091 0 0 0 0 0
0 0.5882 0 0 0 0 0 0 0.4348 0 0 0 0 0 0 0.3704 0 0 0 0 0 0 0.4167 0
2
误差理论与测量平差基础（MATLAB）
0 0 0 0 0 0.2500
由此组成法方程为
NBB-1= 0.5307 0.1608 0.1608 0.7758 4.计算V
V= 5.精度评定
单位权中误差 P1，P2点高程中误差 P1至P2点高差平差值中误差
6.MATLAB解算过程
function [v,ch,cx]=szw(s,h,B,x0,d,n,t,f)%改正数，高差中误差，高程中误差 p=diag(1./s);%定义权阵 l=h-B*x0-d; W=B'*p*l; Nbb=B'*p*B; x=(inv(Nbb))*W; disp('改正数') v=B*x-l;%改正数 disp(v);
c0=sqrt((v'*p*v)/(n-t));%单位中误差 Nbb=B'*p*B;
Qh=f*(inv(Nbb))*f';%h5的协因数阵 disp('高差平差值中误差') ch=c0*sqrt(Qh);%高差平差值中误差 disp(ch);
3
误差理论与测量平差基础（MATLAB）
Nbbn=inv(Nbb);%求逆矩阵 disp('高程平差值中误差')
cx=c0*sqrt(diag(Nbbn));%高程平差值中误差 disp(cx); return
load B.txt; load d.txt; load s.txt; load x0.txt; load h.txt; load n.txt; load t.txt load f.txt
[v,ch,cx]=szw(s,h,B,x0,d,n,t,f);
计算结果展示
二.导线网间接平差
本题n=7，既有7个误差方程，其中有4个角度误差方程，3个边长误差方程。必要观测数t=4。现选取待定点坐标平差值为参数，即
4
误差理论与测量平差基础（MATLAB）
1.计算待定点近似坐标
根据坐标正算公式，算得E,F俩点的近似坐标 XE0=203046.366 YE0= -59253.0948 XF0=203071.813 YF0= -59451.6004 2列出误差方程
3.确定边和角的权
已知测角中误差 ，则角度观测值的权为 各导线边的权为
P= 0.4879 0 0 0 0 0 0 0 0.4997 0 0 0 0 0 0 0 0.2897 0 0 0 0 0 0 0 1.0000 0 0 0 0 0 0 0 1.0000 0 0 0 0 0 0 0 1.0000 0 0 0 0 0 0 0 1.0000
4.计算边长和角度误差方程系数和常数项，V列为边长和角度改正数，在解出坐标改正数后给出
5
2016测量平差实习报告
