后方交会MATLAB程序实习报告

inv(A'*A);m0=sqrt(V'*V/2);mi=m0*sqrt(Qii);m=diag(mi);m=[m(1),m(2),m(3),m(4)/pi*180*3600,m(5)/pi*180*3600,m(6)/pi*180*3600];mXs=m(1);mYs=m(2);mZs=m(3);mfai=m(4);momig=m(5);mka=m(6);fp=fopen('后方交会计算结果.txt','wt');fprintf(fp,'迭代次数:%d\n',v);fprintf(fp,'\n旋转矩阵R:\n');[m,n]=size(R);fori=1:1:mforj=1:1:nifj==nfprintf(fp,'%g\n',R(i,j));elsefprintf(fp,'%g\t',R(i,j));endendendfprintf(fp,'\n外方位元素解:\n');fprintf(fp,'Xs=%g\t',Xs);fprintf(fp,'Ys=%g\t',Ys);fprintf(fp,'Zs=%g\n',Zs);fprintf(fp,'φ=%g\t',fai);fprintf(fp,'ω=%g\t',omig);fprintf(fp,'κ=%g\n',ka);fprintf(fp,'\n单位权中误差的值:mo=%g\n',m0);fprintf(fp,'\n外方位元素中误差为:\n');fprintf(fp,'mXs=%g米\t',mXs);fprintf(fp,'mYs=%g米\t',mYs);fprintf(fp,'mZs=%g米\n',mZs);fprintf(fp,'mφ=%g秒\t',mfai);fprintf(fp,'mω=%g秒\t',momig);fprintf(fp,'mκ=%g秒\n',mka);实习结果1、起算数据2、运行结果3、几次迭代的数据