2016测量平差实习报告

来源：用户分享时间：2025/8/5 9:02:40 本文由

loading 分享下载这篇文档手机版

说明：文章内容仅供预览，部分内容可能不全，需要完整文档或者需要复制内容，请下载word后使用。下载word有问题请添加微信号:xxxxxxx或QQ：xxxxxx 处理（尽可能给您提供完整文档），感谢您的支持与谅解。

误差理论与测量平差基础（MATLAB）

误差理论与测量平差基础

（MATLAB）实习报告

学号：姓名：班级： 1420501Z 专业：测绘工程课程名称：误差理论与测量平差基础任课老师：

2016 年 5 月

误差理论与测量平差基础（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

误差理论与测量平差基础（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);

误差理论与测量平差基础（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。现选取待定点坐标平差值为参数，即

误差理论与测量平差基础（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列为边长和角度改正数，在解出坐标改正数后给出

搜索更多关于： 2016测量平差实习报告的文档

2016测量平差实习报告.doc 将本文的Word文档下载到电脑，方便复制、编辑、收藏和打印

下载这篇word文档

本文链接：https://www.diyifanwen.net/c4xid376ios3ibqx7simj_1.html（转载请注明文章来源）