水准网平差实例.pdf

该【水准网平差实例 】是由【青山代下】上传分享，文档一共【4】页，该文档可以免费在线阅读，需要了解更多关于【水准网平差实例 】的内容，可以使用淘豆网的站内搜索功能，选择自己适合的文档，以下文字是截取该文章内的部分文字，如需要获得完整电子版，请下载此文档到您的设备，方便您编辑和打印。:..tiaojian%A=[1-100100;001-1100;0010011;010-1000];W=[786-3]';%协因数阵的组成Q=diag([]);P=inv(Q);%法方程组成Naa=A*Q*A';%求解联系数K和改正数K=-inv(Naa)*W;V=Q*A'*K;%计算观测值的平差值h=[-]';adjh=h+V/1000;%精度评定VTPV=V'*P*V;r=4;var0=sqrt(VTPV/r);%计算单位权方差%评定观测值平差值函数的精度fs=[0000100]';%函数的系数Qss=fs'*Q*fs-(A*Q*fs)'*inv(Naa)*A*Q*fs;%观测值平差值函数的协因数varS=Qss*var0;%观测值平差值函数的方差%写输出结果文件fid=fopen(陈岑桌面 ','w+');fprintf(fid,'%s','######水准网平差算例######');fprintf(fid,'%s','观测值改正数');fprintf(fid,'');fori=1:7fprintf(fid,'%s%','V',i);fprintf(fid,'%',V(i));endfprintf(fid,'');fprintf(fid,'%s','观测值平差值');fori=1:7fprintf(fid,'%s%','h',i);fprintf(fid,'%',h(i));fprintf(fid,'');endfprintf(fid,'');:..fprintf(fid,'%s','');fprintf(fid,'');fprintf(fid,'%s%','m0=',var0);fprintf(fid,'');fclose(fid);######水准网平差算例######观测值改正数V1-------=--------------------------------------------------------------------------------------间接平差(水准网)function[B,l,P,h]=Finput%误差方程的组成B=[10;01;10;01;-11;-10];l=[004372]';%协因数阵Q的组成Q=diag([]);P=inv(Q);%观测值:..functionjianjie1[B,l,P,h]=Finput;%fafangchengdezuchengNbb=B'*P*B;%adjxVadjx=inv(Nbb)*B'*P*l;V=B*adjx-l;%adjh=h+V/1000;%精度评定VTPV=V'*P*V;r=4;var0=sqrt(VTPV/r);%计算单位权方差%写输出结果文件fid=fopen(陈岑桌面,'w+');fprintf(fid,'%s','######水准网间接平差算例######');fprintf(fid,'%s','观测值改正数');fprintf(fid,'');fori=1:6fprintf(fid,'%s%','V',i);fprintf(fid,'%',V(i));endfprintf(fid,'');fprintf(fid,'%s','观测值平差值');fori=1:6fprintf(fid,'%s%','adjh',i);fprintf(fid,'%',adjh(i));fprintf(fid,'');endfprintf(fid,'');fprintf(fid,'%s','单位权中误差');fprintf(fid,'');fprintf(fid,'%s%','m0=',var0);fprintf(fid,'');fclose(fid);######水准网间接平差算例######观测值改正数:..-----=