实验三利用matlab程序设计语言完成某工程导线网平差计算.docVIP

实验三 利用matlab程序设计语言完成某工程导线网平差计算 实验数据； 某工程项目按城市测量规范（CJJ8-99）不设一个二级导线网作为首级平面控制网，主要技术要求为：平均边长200cm，测角中误差±8,导线全长相对闭合差1/10000,最弱点的点位中误差不得大于5cm，经过测量得到观测数据，设角度为等精度观测值、测角中误差为m=±8秒，鞭长光电测距、测距中误差为m=±0.8√smm，根据所学的‘误差理论与测量平差基础’提出一个最佳的平差方案，利用matlab完成该网的严密平差级精度评定计算; 平差程序设计思路: 1采用间接平差方法,12个点的坐标的平差值作为参数.利用matlab进行坐标反算，求出已知坐标方位角；根据已知图形各观测方向方位角； 2计算各待定点的近似坐标, 然后反算出近似方位角，近似边.计算各边坐标方位角改正数系数； 3确定角和边的权,角度权Pj=1；边长权Ps=100/S； 4计算角度和边长的误差方程系数和常数项，列出误差方程系数矩阵B，算出Nbb=B’PB,W=B’Pl,参数改正数x=inv(Nbb)*W;角度和边长改正数V=Bx-l； 6 建立法方程和解算x,计算坐标平差值, 精度计算； 程序代码以及说明: s10=238.619;s20=170.759; s30=217.869;s40=318.173; s50=245.635;s60=215.514; s70=273.829;s80=241.560; s90=224.996;s100=261.826; s110=279.840;s120=346.443; s130=312.109;s140=197.637; %已知点间距离 Xa=5256.953;Ya=4520.068; Xb=5163.752;Yb=4281.277; Xc=3659.371;Yc=3621.210; Xd=4119.879;Yd=3891.607; Xe=4581.150;Ye=5345.292; Xf=4851.554;Yf=5316.953; %已知点坐标值 a0=atand((Yb-Ya)/(Xb-Xa))+180; d0=atand((Yd-Yc)/(Xd-Xc)); f0=atand((Yf-Ye)/(Xf-Xe))+360; %坐标反算方位角 a1=a0+(163+45/60+4/3600)-180 a2=a1+(64+58/60+37/3600)-180; a3=a2+(250+18/60+11/3600)-180; a4=a3+(103+57/60+34/3600)-180; a5=d0+(83+8/60+5/3600)+180; a6=a5+(258+54/60+18/3600)-180-360; a7=a6+(249+13/60+17/3600)-180; a8=a7+(207+32/60+34/3600)-180; a9=a8+(169+10/60+30/3600)-180; a10=a9+(98+22/60+4/3600)-180; a12=f0+(111+14/60+23/3600)-180; a13=a12+(79+20/60+18/3600)-180; a14=a13+(268+6/60+4/3600)-180; a15=a14+(180+41/60+18/3600)-180; %推算个点方位角 aa=[a1 a2 a3 a4 a5 a6 a7 a8 a9 a10 a12 a13 a14 a15] X20=Xb+s10*cosd(a1); X30=X20+s20*cosd(a2); X40=X30+s30*cosd(a3); X50a=X40+s40*cosd(a4); X60=Xd+s50*cosd(a5); X70=X60+s60*cosd(a6); X80=X70+s70*cosd(a7); X90=X80+s80*cosd(a8); X100=X90+s90*cosd(a9); X50c=X100+s100*cosd(a10); X130=Xf+s110*cosd(a12); X140=X130+s120*cosd(a13); X150=X140+s130*cosd(a14); X50e=X150+s140*cosd(a15); %各点横坐标近似值 X0=[X20 X30 X40 X60 X70 X80 X90 X100 X130 X140 X150 X50a X50c X50e] Y20=Yb+s10*sind(a1); Y30=Y20+s20*sind(a2); Y40=Y30+s30*sind(a3); Y50a=Y40+s40*sind(a4); Y60=Yd+s50*sind(a5); Y70=Y60+s60*sind