开发工具:
文件大小: 11mb
下载次数: 0
上传时间: 2014-05-06
详细说明： unction [Err,Ai,Bj,perO,perN,Enew]=balancing(A,E,Flag,J,bJ) %计算测网焦点差值的平差值。 %输出：Ai-主测网平差值，应用（主测线-Ai）；Bj-联络测网平差值，应用（联络测网+Bj） % perO-平差前测网精度,perN-平差后测网精度 % Err-平差是否正确完成，‘0’-正确，其他-出错信息 %输入：A-测网焦点构造的矩阵，E-测网焦点的差异值（主-辅） % Flag-平差算法选择，0-规则测网平差（要求交点都存在），1-不规则测网最小二乘平差 % （J，bJ）当Flag==0时无用，当Flage==1时，如果J=0表示最小方差平差，bJ无用 % 当0(M+N) Err='测线号指定错误，超出测线范围。'; fprintf(Err); return; end b=calE(E); j=0; b1=zeros(M+N-1,1); for i=1:M+N if i~=J j=j+1; b1(j)=b(i)-A(i,J)*bJ; end end A(J,:)=[]; A(:,J)=[]; AB=A\b1; if J<=M Ai=zeros(1,M); j=0; Ai(J)=bJ; for i=1:M-1 j=j+1; if i==J j=j+1; end Ai(j)=AB(i); end Bj=AB(M:M+N-1)'; else Ai=AB(1:M)'; Bj=zeros(1,N); j=0; Bj(J-M)=bJ; for i=M+1:M+N-1 j=j+1; if i==J j=j+1; end Bj(j)=AB(i); end end end for i=1:M for j=1:N if A(i,M+j)==1 E(i,j)=E(i,j)-Ai(i)-Bj(j); end end end W=sum(sum(E.^2)); perN=sqrt(W/Z/2);%平差后测网精度计算 end Enew=E; end function b=calE(E) %计算AM=b的b值,M-主测线条数；N-联络测线条数 [M,N]=size(E); b=zeros(M+N,1); for i=1:M b(i)=sum(E(i,:)); end for i=1:N b(i+M)=sum(E(:,i)); end end ...展开收缩
(系统自动生成,下载前可以参看下载内容)
下载文件列表
相关说明
- 本站资源为会员上传分享交流与学习,如有侵犯您的权益,请联系我们删除.
- 本站是交换下载平台，提供交流渠道，下载内容来自于网络，除下载问题外，其它问题请自行百度。
- 本站已设置防盗链，请勿用迅雷、QQ旋风等多线程下载软件下载资源，下载后用WinRAR最新版进行解压.
- 如果您发现内容无法下载，请稍后再次尝试；或者到消费记录里找到下载记录反馈给我们.
- 下载后发现下载的内容跟说明不相乎，请到消费记录里找到下载记录反馈给我们，经确认后退回积分.
- 如下载前有疑问，可以通过点击"提供者"的名字，查看对方的联系方式，联系对方咨询.
