水准间接平差程序设计说明书

可编辑版地球科学与环境工程学院准平验书课程名差理与测量平差基础》学姓

号：名：

黄黎东指导老师：日

期：2015年12月日Word完美格式

可编辑版一、任务概述利用MATLAB或者C++程间接平差程序过该程序读取观测数据文件，并算出平差结果。二、计算结果截图：图

图Word完美格式

可编辑版图三、水准网图

图Word完美格式

可编辑版四、输入的数据格式数据格式为TXT文件如图所示：TXT文件格式说明：（1）一行格式第一行分别表示观测个数点数知点3个，已知点1个所有数据用英文逗号隔开（2）知点数据格式第二行开始是已知点点号和高程，一行列一个已知点点号和高程，由于该水准网只有一个已知点，所有只能列出一行。图中表示已知点点号为，高程为237.483m（3）站起始点号格式Word完美格式

可编辑版该部分表示测站的起始点点号（4）站终点点号格式该部分表示测站的终点点号（5）差格式该部分表示各测站的高差（6）离格式Word完美格式

可编辑版该部分表示各测站的距离五、流程图开始读取数据提取、P、l矩平差输出结束六、附件代码Word完美格式

可编辑版function[FileName,PathName]=uigetfile('*.txt','打开准观测数据'打开文件f=csvread(打文件并存在矩阵f中point=f(1,2);%获取所有水准点个数n=f(1,1);%得观测个数t=f(1,3);%得必要观测个数ty=f(1,4);%获得已知点个数yXX=zeros(point,1);%初化XX阵于0方便下面把已知点高程和未知点参数估值放到XX阵初化B阵，方便下面求中系数阵Bforj=1:yXX(j,1)=f(j+1,2);%把已知高程放到XX阵从文件中获取观测数据，并放到data阵中h=data(:,3);%从data中获取观测高差，并放到中初化权阵Pforj=1:nP(j,j)=10/data(j,4);%以10km观值为单位权误差计算权阵Pfor通过循环求B阵point1=data(i,1);%取某个测站的起始点号point2=data(i,2);%取某个测站的终点点号ifpoint1>y&&point2>y%当测站起始点和终点高程都未知时，求B第i行B(i,point1-y)=-1;B(i,point2-y)=1;elseifpoint1<=y&&point2>y%起始点高程已知和终点高程未知时，求第i行B(i,point2-y)=1;求i个数估值elseifpoint1>y&&point2<=y%起始点高程未知和终点高程已知时，求第i行B(i,point1-y)=-1;XX(point1,1)=XX(point2,1)-h(i,1);%第i个参数估值l=zeros(n,1);%初始小l阵方便下面求V=Bx-l的系数阵lfor通过循环求小ll(i,1)=-(XX(point2,1)-XX(point1,1)-h(i,1));%入间接平差数学模型公式进行算：Word完美格式

可编辑版r=n-t;%多余观测数x=N\W;V=B*x-l;Qxx=inv(N);Dxx=a0*a0*inv(N);%出计算结果：