實(shí)驗(yàn)三-利用matlab程序設(shè)計(jì)語言完成某工程導(dǎo)線網(wǎng)平差計(jì)算

1、精選優(yōu)質(zhì)文檔-傾情為你奉上實(shí)驗(yàn)三 利用matlab程序設(shè)計(jì)語言完成某工程導(dǎo)線網(wǎng)平差計(jì)算實(shí)驗(yàn)數(shù)據(jù)；某工程項(xiàng)目按城市測量規(guī)范（CJJ8-99）不設(shè)一個(gè)二級導(dǎo)線網(wǎng)作為首級平面控制網(wǎng)，主要技術(shù)要求為：平均邊長200cm，測角中誤差8,導(dǎo)線全長相對閉合差1/10000,最弱點(diǎn)的點(diǎn)位中誤差不得大于5cm，經(jīng)過測量得到觀測數(shù)據(jù)，設(shè)角度為等精度觀測值、測角中誤差為m=8秒，鞭長光電測距、測距中誤差為m=0.8smm，根據(jù)所學(xué)的誤差理論與測量平差基礎(chǔ)提出一個(gè)最佳的平差方案，利用matlab完成該網(wǎng)的嚴(yán)密平差級精度評定計(jì)算;平差程序設(shè)計(jì)思路:1采用間接平差方法,12個(gè)點(diǎn)的坐標(biāo)的平差值作為參數(shù).利用matlab進(jìn)行

2、坐標(biāo)反算，求出已知坐標(biāo)方位角；根據(jù)已知圖形各觀測方向方位角；2計(jì)算各待定點(diǎn)的近似坐標(biāo), 然后反算出近似方位角，近似邊.計(jì)算各邊坐標(biāo)方位角改正數(shù)系數(shù)；3確定角和邊的權(quán),角度權(quán)Pj=1；邊長權(quán)Ps=100/S；4計(jì)算角度和邊長的誤差方程系數(shù)和常數(shù)項(xiàng)，列出誤差方程系數(shù)矩陣B，算出Nbb=BPB,W=BPl,參數(shù)改正數(shù)x=inv(Nbb)*W;角度和邊長改正數(shù)V=Bx-l；6建立法方程和解算x,計(jì)算坐標(biāo)平差值,精度計(jì)算；程序代碼以及說明:s10=238.619;s20=170.759;s30=217.869;s40=318.173;s50=245.635;s60=215.514;s70=273.829

3、;s80=241.560;s90=224.996;s100=261.826;s110=279.840;s120=346.443; s130=312.109;s140=197.637; %已知點(diǎn)間距離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; %已知點(diǎn)坐標(biāo)值a0=atand(Yb-Ya)/(Xb-Xa)+180;d0=atand(Yd-Yc)/(Xd-X

4、c);f0=atand(Yf-Ye)/(Xf-Xe)+360; %坐標(biāo)反算方位角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/36

5、00)-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; %推算個(gè)點(diǎn)方位角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(a

6、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); %各點(diǎn)橫坐標(biāo)近似值X0=X20 X30 X40 X60 X70 X80 X90 X100 X130 X140 X150 X50a X50

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

8、50+s140*sind(a15); %個(gè)點(diǎn)從坐標(biāo)近似值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) cos

9、d(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.*inv(s)*B11b=-20

10、6.*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=

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

12、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=(abs(Y50-

13、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 %以上為求得誤差方程系數(shù)B=-0.685

14、10.527100000000000000000000000.38720.0289-1.0723-0.55600000000000000000000000.94531.49420.127-0.9382-0.938200000000000000000000.127-0.93820.47561.1776-0.6026-0.23940000000000000000000000000.8256-0.15360000000000000000000000-0.6191-0.7809-0.20650.934500000000000000000000-0.20650.93450.9518-1.0435-0.7

15、4530.10900000000000000000000-0.74530.1091.55160.1722-0.8063-0.281200000000000000000000-0.8063-0.28121.71320.4151-0.9069-0.1339000000000000-0.24910.8277000000-0.90690.13391.156-0.693800000000000.60260.23940.8517-0.586300000000-0.24910.82770000000000000000000000000.7111-0.1940000000000000000000000-0.7

16、5880.78750.0477-0.5935000000000000000000000.0477-0.5935-0.70780.62460.6601-0.03110000001.0417-0.06160000000000000.6601-0.0311-1.70180.0927000000-1.28090.7661000000000.2491-0.827700001.0417-0.0.60970.79260000000000000000000000-0.4603-0.88780.46030.88780000000000000000000000-0.991-0.13420.9910.1342000

17、0000000000000000000-0.3692-0.92930.36920.92930000000000000000000000000.39960.91670000000000000000000000-0.9764-0.21580.97640.21580000000000000000000000-0.1447-0.98950.14470.98950000000000000000000000-0.3293-0.94420.32930.94420000000000000000000000-0.1461-0.98930.14610.98930000000000000.95740.2888000

18、00000-0.9574-0.28880000000000000000000000000.26310.96480000000000000000000000-0.9968-0.08010.99680.08010000000000000000000000-0.047-0.99890.0470.99890000000.05940.998200000000000000-0.0594-0.9982 %系數(shù)矩陣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