资源描述：
《2016测量平差实习报告》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、EASTCHINAINSTITUTEOFTECHNOLOGY误差理论与测量平差基础(MATLAB)实习报告学号：姓名:班级：1420501Z专业：测绘工程课程名称.•误差理论与测量平差基础任课老师：2016年5月一.水准网间接平差1.列出误差方程设P1和P2点高程平差值为和，相应的近似值取按已知数据及观测数据表列出观测方程后，将有关观测:据带入即得误差方程V2=V3=V5:V6=式中常数项以m为单位。2.列出权函数式P1至P2间高差平差值的权函数式为3.组成法方程p=0.909100000.58820000以lkm水准测量的观测高差为
2、单位权观测值，各观测宣互相独立，定权式为Pi=l/Si，得权阵为000000.4348000000.37040000.4167000.2500由此组成法方程为Nbb'1=0.53070.16080.16080.77581.计算VV=2.精度评定单位权中误差Pl,P2点高程中误差P1至P2点高差平差值中误差6.MATLAB解算过程function[v,ch,ex]=szw(s,h,B,xO,d,n,t,f)%改正数，商差屮误差，商程屮误差p=diag(1./s);%定义权阵l=h-B*xO-d；W=B,*p*l；Nbb=B'*p*B；x
3、=(inv(Nbb))*W；disp(•改正数1)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);%高左〒•方ffi中识方disp(ch)；Nbbn=inv(Nbb);%求逆):1(昨disp(•高程平差值中误差i)cx=cO*sqrt(diag(Nbbn));%髙权平差值中误差disp(ex)；returnloadB.txt；loadd.txt
4、；loads.txt；loadxO•txt;loadh.txt；loadn.txt；loadt.txtloadf.txt[v,ch,ex]=szw(s,h,B,xO,d,n,t,f);计算结果展示ap当前文件夹ProgramFiles►MATLAB►R2015b►bin►zfCB网阅2.0®r命令行窗口>>shuizhunwang改正教-0.00020.0029-0.0042-0.0001-0.0039-0.0018S差平差值中误差0.0022高程平差值中误差0.00160.0020IA»二.导线网间接平差本题n=7,既有7个误差方程
5、，其中有4个角度误差方程，3个边长误差方程。必要观测数t=4。现选取待定点坐标平差值为参1.计算待定点近似坐标根据坐标正算公式算得E，F俩点的近似坐标XE0=203046.366XF0=203071.813YEO=-59253.0948YF0=-59451.60042列出误差方程3.确定边和角的权已知测角中误差各导线边的权为,则角度观测值的权为000-289700000001.000000000001.000000000001.00000000.49970000000000000000L0000P=0.4879TJ4.计算边长和角度误
6、差方程系数和常数项，V列为边长和j度改正数，在解出坐标改正数后给出5.法方程的组成和解算Nbb*x-Bt*P*L=Ox=Nbb4*Bt*P*L6.平差值计算(1)坐标平差值将坐标改正数加上近似值，即得平差值(2)观测值的平差值将改正数与观测值相加，即得观测量的平差值7.精度计算(X)单位权中误差(2)待定点点位中误差由Nbb4中可得未知数的权倒(3)各点位中误差为8.MATLAB解算过程load11.七父七％观测兑数loadt.txt%必要观测数loaddws.txt%¥入，,:位数loadcsgs.txt%导入参数个数loadbs.
7、txt%导入彳、Y测边数loaddcdhl.txt%%第、个待测点号loadxlyl.txt%导入已知点位-、的坐标loadxydws.txt%4•入末/乂位的坐标loaddeg.txt%导入山观测角计17而得的77位角loads.txt%@入边长观测估loadgcj.txt%各位观测ffjdigits(12)；digits(15);xl=xlyl(1,1);yl=xlyl(1,2)；xO=zeros(dws,1);yO=zeros(dws,1)；a=zeros(1,3);xO(1)=xl；yO(1)=yi；fori=2:(dws-1
8、)a=deg(i-lz:)；[x2,y2]=zbzs(xl,yl,s(i-1),a)；xl=x2；yl=y2;xO(i)=xl;yO(i)=y1;endxl=xydws(1,1);yl=xydws(1,2);xO(dws
