武汉大学测量平差基础-第五章、条件平差

function Excise5_3_38 %本程序用于计算误差理论与测量平差基础习题集，武汉大学测绘学院与测量平差教学组编著的教材 %S5-1条件平差原理中练习5_2_27的计算P31 %本程序版权所有:王胜利 %联系方式： %QQ： 12782600 %MSN： [email protected] %E-mail: [email protected] L(1)=dms2rad(63,14,25.02); L(2)=dms2rad(23,28,50.06); L(3)=dms2rad(23,31,29.31); L(4)=dms2rad(69,45,14.74); L(5)=dms2rad(61,40,57.38); L(6)=dms2rad(25,2,19.23); L(7)=dms2rad(27,24,8.77); L(8)=dms2rad(65,52,35.08); n=8;t=4;r=4; p0=180*3600/pi; A=zeros(4,8); w(1)=L(1)+L(2)+L(3)+L(4)-pi; w(2)=L(5)+L(6)+L(7)+L(8)-pi; w(3)=L(3)+L(4)-L(7)-L(8); A(1,1)=1;A(1,2)=1;A(1,3)=1;A(1,4)=1; A(2,5)=1;A(2,6)=1;A(2,7)=1;A(2,8)=1; A(3,3)=1;A(3,4)=1;A(3,7)=-1;A(3,8)=-1; w(4)=1-(sin(L(1)+L(8))*sin(L(3))*sin(L(4))*sin(L(6)+L(7)))/(sin(L(7))*sin(L(2)+L(3))*sin(L(4)+L(5))*sin(L(8))); w=p0*w; A(4,1)=cos(L(1)+L(8))/sin(L(1)+L(8)); A(4,2)=-cos(L(2)+L(3))/sin(L(2)+L(3)); A(4,3)=cos(L(3))/sin(L(3))-cos(L(2)+L(3))/sin(L(2)+L(3)); A(4,4)=cos(L(4))/sin(L(4))-cos(L(4)+L(5))/sin(L(4)+L(5)); A(4,5)=-cos(L(4)+L(5))/sin(L(4)+L(5)); A(4,6)=cos(L(6)+L(7))/sin(L(6)+L(7)); A(4,7)=cos(L(6)+L(7))/sin(L(6)+L(7))-cos(L(7))/sin(L(7)); A(4,8)=cos(L(1)+L(8))/sin(L(1)+L(8))-1/atan(L(8)); P=eye(n); Naa=A*inv(P)*A'; K=-inv(Naa)*w'; V=inv(P)*A'*K; V=V/p0; L=L+V'; Q=inv(P); PL=P-P*A'*inv(Naa)*A*P; f=zeros(1,8); f(1)=cos(L(1)+L(8))/sin(L(1)+L(8)); f(2)=-cos(L(2)+L(3))/sin(L(2)+L(3)); f(3)=-cos(L(2)+L(3))/sin(L(2)+L(3)); f(4)=cos(L(4))/sin(L(4)); f(7)=-cos(L(7))/sin(L(7)); f(8)=cos(L(1)+L(8))/sin(L(1)+L(8)); Pf=f*PL*f'