资源描述：
《2016测量平差实习报告.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、误差理论与测量平差基础（MATLAB）实习报告学号：姓名：班级：Z专业：测绘工程课程名称：误差理论与测量平差基础任课老师：2016年5月一.水准网间接平差1.列出误差方程设P1和P2点高程平差值为和，相应的近似值取为按已知数据及观测数据表列出观测方程后，将有关观测数据带入即得误差方程V1=V2=V3=V4=V5=V6=式中常数项以m为单位。2.列出权函数式P1至P2间高差平差值的权函数式为3.组成法方程以1km水准测量的观测高差为单位权观测值，各观测值互相独立，定权式为Pi=1/Si，得权阵为p=0.90910000000.58820000000.
2、43480000000.37040000000.41670000000.2500由此组成法方程为NBB-1=0.53070.16080.16080.77581.计算VV=2.精度评定单位权中误差P1，P2点高程中误差P1至P2点高差平差值中误差3.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;%改正数d
3、isp(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);Nbbn=inv(Nbb);%求逆矩阵disp('高程平差值中误差')cx=c0*sqrt(diag(Nbbn));%高程平差值中误差disp(cx);returnloadB.txt;loadd.txt;loads.txt;loadx0.txt;loadh.txt;loadn.txt;loadt.txtl
4、oadf.txt[v,ch,cx]=szw(s,h,B,x0,d,n,t,f);计算结果展示二.导线网间接平差本题n=7，既有7个误差方程，其中有4个角度误差方程，3个边长误差方程。必要观测数t=4。现选取待定点坐标平差值为参数，即1.计算待定点近似坐标根据坐标正算公式，算得E,F俩点的近似坐标XE0=.366YE0=-59253.0948XF0=.813YF0=-59451.60042列出误差方程3.确定边和角的权已知测角中误差，则角度观测值的权为各导线边的权为P=0.487900000000.499700000000.289700000001.
5、000000000001.000000000001.000000000001.00004.计算边长和角度误差方程系数和常数项，V列为边长和角度改正数，在解出坐标改正数后给出5.法方程的组成和解算NBB*x-BT*P*L=0x=NBB-1*BT*P*L3.平差值计算（1）坐标平差值将坐标改正数加上近似值，即得平差值（2）观测值的平差值将改正数与观测值相加，即得观测量的平差值7.精度计算（1）单位权中误差（2）待定点点位中误差由NBB-1中可得未知数的权倒数（3）各点位中误差为8.MATLAB解算过程loadn.txt%观测总数loadt.txt%必要
6、观测数loaddws.txt%导入点位数loadcsgs.txt%导入参数个数loadbs.txt%导入待测边数loaddcdh1.txt%%第一个待测点号loadx1y1.txt%导入已知点位一的坐标loadxydws.txt%导入末点位的坐标loaddeg.txt%导入由观测角计算而得的方位角loads.txt%导入边长观测值loadgcj.txt%各点位观测角digits(12);digits(15);x1=x1y1(1,1);y1=x1y1(1,2);x0=zeros(dws,1);y0=zeros(dws,1);a=zeros(1,3);
7、x0(1)=x1;y0(1)=y1;fori=2:(dws-1)a=deg(i-1,:);[x2,y2]=zbzs(x1,y1,s(i-1),a);x1=x2;y1=y2;x0(i)=x1;y0(i)=y1;endx1=xydws(1,1);y1=xydws(1,2);x0(dws)=x1;%各点位x的近似值y0(dws)=y1;%各点位y的近似值s0=zeros(dws-1,1);%各各点位之间距离的近似值fori=1:(dws-1)s0(i)=sqrt((x0(i)-x0(i+1)).^2+(y0(i)-y0(i+1)).^2);end[B1,
8、L1]=jsB1L1(x0,y0,s0,s,bs,dcdh1,csgs);[B2,L2,ajs]=jsB2L2(x0,y0
