实验三利用matlab程序设计语言完成某工程导线网平差计算Word文件下载.docx

1、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

2、）;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/360

3、0）-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 a15X20=Xb+s10*cosd（a1）;X30=X20+s20*cosd（a2）;X40=X30+s30*cosd（a3）;X50a=X40+s40*cosd（a4

4、）;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

5、 X50eY20=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=Y15

6、0+s140*sind（a15）; %个点从坐标近似值Y0=Y20 Y30 Y40 Y60 Y70 Y80 Y90 Y100 Y130 Y140 Y150 Y50a Y50c Y50eP=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（

7、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=

8、-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-a

9、b14+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-s14e1=（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）

10、/s60;e7=（abs（X80-X70）/s70;e8=（abs（X90-X80）/s80;e9=（abs（X100-X90）/s90;e10=（abs（X50-X100）/s1;e11=（abs（X130-Xf）/s110;e12=（abs（X140-X130）/s120;e13=（abs（X150-X140）/s130;e14=（abs（X50-X150）/s14;e=e1 e2 e3 e4 e5 e6 e7 e8 e9 e10 e11 e12 e13 e14m1=（abs（Y20-Yb）/s10;m2=（abs（Y30-Y20）/s20;m3=（abs（Y40-Y30）/s30;m4=

11、（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（Y50-Y100）/s1;m11=（abs（Y130-Yf）/s110;m12=（abs（Y140-Y130）/s120;m13=（abs（Y150-Y140）/s130;m14=（abs（Y50-Y150）/s14;m=m1 m2 m3 m4 m5 m6 m7 m8 m9 m10 m11 m12 m13 m14 %以上为求得误差方程系

12、数B=-0.6851 0.5271 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 00.3872 0.0289 -1.0723 -0.556 0 0 0 0 0 0 0 0 0 0 0 0 0 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 0 0 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 0 0 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0.8

13、256 -0.1536 0 0 0 0 0 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 -0.6191 -0.7809 -0.2065 0.9345 0 0 0 0 0 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.109 0 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 -

14、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.1339 1.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 -

15、0.7588 0.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 0 0 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

16、0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0-0.4603 -0.8878 0.4603 0.8878 0 0 0 0 0 0 0 0 0 0 0 0 0 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 0 0 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 0 0 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 0 0 0 0 0 0 00 0 0 0

17、 0 0 0 0 -0.9764 -0.2158 0.9764 0.2158 0 0 0 0 0 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 0 0 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.9442 0 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.1461 0.9893 0 0 0 0 0 00 0 0 0 0 0 0.9574 0.2888 0 0 0

18、 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.2631 0.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 0 0 0 0 0 -0.0594 -0.9982

19、 %系数矩阵BP=blkdiag（1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,100/s10,100/s20,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.41 4.779 4.9396669193517

20、7 4.24372266973269 4.72376669705670 4.21449589983501 4.606 4.51 4.277 3.66643202793325 4.42847073890118 3.71293519499601 4.468 3.98387916563715 4.38851642122011 4.21197375648003 4.35564105689259 4.43455120690799 4.925 5.905 4.57985606835358 5.518 4.59452770624570 4.734点名 坐 标（m） X Y 2 5018.25824092.1

21、405 3 4939.66694243.7227 4 4723.76674214.4959 5 4606.29214510.1872 6 4218.04233666.4320 7 4428.47073712.9352 8 4468.07143983.8792 9 4388.51644211.9738 10 4355.64114434.5512 13 4925.18735046.9904 14 4579.85615019.2538 15 4595.52774707.5059Qx1=1.488 Qy1=1.5596 Qx2=2.2058 Qy2=5.1342 Qx15=1.1836 Qy15=4.9268