平差实验报告及完整matlab程序..doc

秩亏自由网平差综合计算 一、实验名称：秩亏自由网平差综合计算 二、实验目的：掌握应用MATLAB程序设计语言进行秩亏自由网平差计算：具体包括建立数学模型，点位稳定性分析、平差基准选取等。 三、实验设备：电脑 四、实验数据 序号 hi(m) Si(km) XO 1 86.809 204.1 -0.004 2 25.714 188.7 86.806 3 31.225 344.8 14.846 4 -71.952 149.2 25.706 5 -61.084 142.9 31.216 6 -44.178 250.0 42.626 7 10.847 128.2 8 16.350 98.0 9 11.409 196.1 五、计算过程 1、经典自由网平差（以6号点为例） 由已知可列出误差方程 （1） 其中 根据最小二乘准则 （2） 组成法方程 （3） 解得参数的最佳估值为 参数的协因数阵为 2、秩亏自由网平差 由已知可列出误差方程 （2-1） 其中 根据重心基准的约束条件 (2-2) 得出 (2-3) 求得参数的最佳估值为 参数的协因数阵为 3、拟稳平差（选1、2、3、4号点为稳定点） 由已知可列出误差方程 （3-1） 其中 拟稳平差的约束条件为 （3-2） 其中 (3-3) (3-4) 得出 （3-5） 求得参数的最佳估值为 参数的协因数阵为 六、综合分析 从结果中我们可以看出，三种平差方法所得高差的改正数是相同的，所得参数的最佳估值是不同的。 经典自由网平差适用于起算数据充足的情况，而秩亏自由网平差适用于起算数据不足，需要加一个基准条件才能进行平差计算。拟稳平差则是选择一些相对稳定的点进行平差，是那些不稳定的点不参与平差。根据实际情况我们选择适宜的平差方法。 七、MATLAB程序代码 1、经典自由网平差 [f,p1]=uigetfile(*.txt,选择高差文件) fnal=strcat(p1,f);fpath=p1; fm=fopen(fnal,r); while (~feof(fm)) l1=fscanf(fm,%f,inf); end [f,p1]=uigetfile(*.txt,选择距离文件) fnal=strcat(p1,f);fpath=p1; fm=fopen(fnal,r); while (~feof(fm)) s=fscanf(fm,%f,inf); end %------------------从文件获取数据-------------------------------- b=[-1 1 0 0 0 0;-1 0 0 1 0 0;-1 0 0 0 1 0;0 -1 1 0 0 0;0 -1 0 1 0 0;0 -1 0 0 0 1;0 0 -1 1 0 0;0 0 -1 0 1 0;0 0 0 0 -1 1]; x0=[-0.004;86.806;14.846;25.706;31.216;42.626]; %------------------从文本框输入数据------------------------------- %m=get(findobj(tag,edit1),string); %l1=str2num(m); %l1=l1; %-----------------从文本框输入数据------------------------------- l2=l1-b*x0; dh_index=get(findobj(tag,dh),value); %----------选择点号--------------------------------------------- switch dh_index case 1 b=b(:,2:6); x0=x0(2:6,:); case 2