文档介绍：测绘工程matlab论文课程编号:课程性质:MATLAB及其应用课程论文学院:测绘学院专业:测绘工程姓名:学号:目录一、题目内容 3二、程序编写思路 41、数学模型 42、程序解析 4三、程序运行结果 5四、程序源代码 6注:这是《误差理论与测量平差基础》这一课程在学****其基本原理之后为解决实际问题而编写的一个程序。一、题目内容要求:对于给出的导线控制网图2,根据已知条件和观测数据,设计平差方案,编制平差程序,给出平差结果,评定精度。题目:A、B为已知点,已知点坐标已给出如下。为待定点。同精度观测了24个角度,测角中误差。测量了17条边长,其观测结果及中误差列于表4中,试按间接平差法求:(1)平差后单位权中误差;(2)14个待定点的的坐标平差值及中误差;(3)平差后边的相对中误差。已知点坐标:=,=,,表一::导线网边长观测值及精度编号边观测值S(m)中误差(mm)编号边长观测值S(m)中误差(mm)编号边长观测值S(m)中误差(mm) 图一:二、程序编写思路1、数学模型◆函数模型观测方程误差方程◆随机模型◆方程的解与精度评定令法方程即为可得解为:故可得平差结果:从而解得单位权中误差为2、程序解析◆计算框图将导线网的角度观测数据和边长观测数据按照自己编程的思路以一定的格式存入两个S1、S2中编写程序,首先由已知点计算出各个待定点的近似坐标由计算出的近似坐标在计算出相应的方位角由导线网的形状设计编程思路计算出系数矩阵B和常数项矩阵l由于导线网有两种不同类型的数据,故求P的时候尤其要注意由B、P、l求出、V、、、Q等等三、程序运行结果◆平差后单位权中误差:=(由于程序取的先验单位权中误差为10,两者相比较,故可以认为这个结果是可以接受的)◆14个待定点的坐标平差值及中误差见下表:表三:点号坐标平差值(m)中误差(mm)点号坐标平差值(m)中误差(mm):点号中误差(mm)点号中误差(mm)点号中误差(mm)◆平差后边的相对中误差为:四、程序源代码这是要输入的文件◆%%从文件中输入数据fid1=fopen('E:\学****平差\','rt');[S1,count]=fscanf(fid1,'%f',[724]);%#ok<NASGU>fclose(f