实验三利用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)-180a2=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(a6);Y80=Y70+s70*sind(a7);Y90=Y80+s80*sind(a8);Y100=Y90+s90*sind(a9);Y50c=Y100+s100*sind(a10);Y130=Yf+s110*sind(a12);Y140=Y130+s120*sind(a13);Y150=Y140+s130*sind(a14);Y50e=Y150+s140*sind(a15); %个点从坐标近似值Y0=[Y20 Y30 Y40 Y60 Y70 Y80 Y90 Y100 Y130 Y140 Y150 Y50a Y50c Y50e]'P=[X0 Y0];X50=(X50a+X50c+X50e)/3Y50=(Y50a+Y50c+Y50e)/3s4=sqrt((Y40-Y50)^2+(X40-X50)^2);s1=sqrt((Y100-Y50)^2+(X100-X50)^2);s14=sqrt((Y150-Y50)^2+(X150-X50)^2);A1=[cosd(a1) cosd(a2) cosd(a3) cosd(a4) cos(a5) cosd(a6) cosd(a7) cosd(a8) cosd(a9) cosd(a10) cosd(a12) cosd(a13) cosd(a14) cosd(a15)]';B11=[sind(a1) sind(a2) sind(a3) sind(a4) sin(a5) sind(a6) sind(a7) sind(a8) sind(a9) sind(a10) sind(a12) sind(a13) sind(a14) sind(a15)]';s=blkdiag(s10,s20,s30,s4,s50,s60,s70,s80,s90,s10',s110, s120,s130,s14);a=206.2648026*inv(s)*B11b=-206.2648026*inv(s)*A1ab4=atand((Y50-Y40)/(X50-X40))+180;ab10=atand((Y50-Y100)/(X50-X100));ab14=atand((Y50-Y150)/(X50-X150))+360;m4=ab4-a3+180;m10=ab10-a9+180;m11=ab4-ab10;m15=ab14-a14+180;m16=ab10-ab14+360;m04=103+57/60+34/3600;m010=98+22/60+4/3600;m011=94+53/60+50/3600;m015=180+41/60+18/3600;m016=ab10-ab14+360;l=[0 0 0 m4-103-57/60-34/3600 0 0 0 0 0 m10-98-22/60-4/3600 m11-94-53/60-50/3600 0 0 0 m15-180-41/60-18/3600 m16-103-23/60-8/3600 0 0 0 s40-s4 0 0 0 0 0 s100-s1 0 0 0 s140-s14]';e1=(abs(X20-Xb))/s10;e2=(abs(X30-X20))/s20;e3=(abs(X40-X30))/s30;e4=(abs(X50-X40))/s4;e5=(abs(X60-Xd))/s50;e6= (abs(X70-X60))/s60;e7=(abs(X80-X70))/s70;e8=(abs(X90-X80))/s80;e9=(abs(X100-X90))/s90;e10=(abs(X 50-X100))/s1;e11=(abs(X130-Xf))/s110;e12=(abs(X140-X130 ))/s120;e13=(abs(X150-X140))/s130;e14=(abs(X50-X150))/s 14;e=[e1 e2 e3 e4 e5 e6 e7 e8 e9 e10 e11 e12 e13 e14]'m1=(abs(Y20-Yb))/s10;m2=(abs(Y30-Y20))/s20;m3=(abs(Y40-Y30))/s30;m4=(abs(Y50-Y40))/s4;m5=(abs(Y60-Yd))/s50;m6= (abs(Y70-Y60))/s60;m7=(abs(Y80-Y70))/s70;m8=(abs(Y90-Y80))/s80;m9=(abs(Y100-Y90))/s90;m10=(abs(Y 50-Y100))/s1;m11=(abs(Y130-Yf))/s110;m12=(abs(Y140-Y130 ))/s120;m13=(abs(Y150-Y140))/s130;m14=(abs(Y50-Y150))/s 14;m=[m1 m2 m3 m4 m5 m6 m7 m8 m9 m10 m11 m12 m13 m14]' %以上为求得误差方程系数B=[-0.6851 0.5271 0 0 0 0 0 0 0 0 0 0 0 0 0 0 00 0 0 0 0 0 00.3872 0.0289 -1.0723 -0.556 0 0 0 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 00 0 0.9453 1.4942 0.127 -0.9382 -0.9382 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 00 0 0.127 -0.9382 0.4756 1.1776 -0.6026 -0.2394 0 0 0 0 0 00 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0.8256 -0.1536 0 0 0 0 0 0 0 00 0 0 0 0 00 0 0 0 0 0 0 0 -0.6191 -0.7809 -0.2065 0.9345 0 0 0 00 0 0 0 0 0 0 00 0 0 0 0 0 0 0 -0.2065 0.9345 0.9518 -1.0435 -0.7453 0.1090 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 0 -0.7453 0.109 1.5516 0.1722 -0.8063-0.2812 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 0 0 0 -0.8063 -0.2812 1.7132 0.4151-0.9069 -0.1339 0 0 0 0 0 00 0 0 0 0 0 -0.2491 0.8277 0 0 0 0 0 0 -0.9069 0.13391.156 -0.6938 0 0 0 0 0 00 0 0 0 0.6026 0.2394 0.8517 -0.5863 0 0 0 0 0 0 0 0-0.2491 0.8277 0 0 0 0 0 00 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.7111-0.194 0 0 0 00 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.75880.7875 0.0477 -0.5935 0 00 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.0477-0.5935 -0.7078 0.6246 0.6601 -0.03110 0 0 0 0 0 1.0417 -0.0616 0 0 0 0 0 0 0 0 0 00 0 0.6601 -0.0311 -1.7018 0.09270 0 0 0 0 0 -1.2809 0.7661 0 0 0 0 0 0 0 0 0.2491-0.8277 0 0 0 0 1.0417 -0.0666160.6097 0.7926 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 00 0 0 0 0 0-0.4603 -0.8878 0.4603 0.8878 0 0 0 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 00 0 -0.991 -0.1342 0.991 0.1342 0 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 00 0 0 0 -0.3692 -0.9293 0.3692 0.9293 0 0 0 0 0 0 0 00 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0.3996 0.9167 0 0 0 0 0 0 0 00 0 0 0 0 00 0 0 0 0 0 0 0 -0.9764 -0.2158 0.9764 0.2158 0 0 0 00 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 0 -0.1447 -0.9895 0.1447 0.9895 0 00 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 0 0 0 -0.3293 -0.9442 0.3293 0.94420 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.1461 -0.9893 0.14610.9893 0 0 0 0 0 00 0 0 0 0 0 0.9574 0.2888 0 0 0 0 0 0 0 0 -0.9574-0.2888 0 0 0 0 0 00 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.26310.9648 0 0 0 00 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.9968-0.0801 0.9968 0.0801 0 00 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0-0.047 -0.9989 0.047 0.99890 0 0 0 0 0 0.0594 0.9982 0 0 0 0 0 0 0 0 0 00 0 0 0 -0.0594 -0.9982] %系数矩阵BP=blkdiag(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,100/s10,100/s 20,100/s30,100/s40,100/s50,100/s60,100/s70,100/s80,100/ s90,100/s100,100/s110,100/s120,100/s130,100/s140); %定义权矩阵Nbb=B'*P*BW=B'*P*l;x=inv(Nbb)*WV=B*x-l;inv(Nbb);Y=V'*P*V;O=sqrt(Y/6)*3600 %精度评定计算结果：平差值坐标X：1.0e+003 *5.414.7794.939666919351774.243722669732694.723766697056704.214495899835014.6064.514.2773.666432027933254.428470738901183.712935194996014.4683.983879165637154.388516421220114.211973756480034.35564105689259 4.434551206907994.9255.9054.579856068353585.5184.59452770624570 4.734Qx1=1.488 Qy1=1.5596 Qx2=2.2058 Qy2=5.1342 ……Qx15=1.1836 Qy15=4.926811 / 11。