trel.m

利用ＭＡＴＬＡＢ实现了对测量中的平面网

%函数 [maa,ms]=trel(uw,rr,tt)是根据平差后的方差-协方差阵计算任意两点rr、
%tt之间平差计算后的方位和边长误差maa、ms。其中uw为单位权中误差，c为法方程求逆后的协因素阵。
function [maa,ms]=trel(uw,rr,tt)   
global ed dd x0 y0 pn c
n=2*(dd+ed);
dx=x0(tt)-x0(rr);
dy=y0(tt)-y0(rr);
ss=sqrt(dx*dx+dy*dy);
a0=alfa(dx,dy);
aij=2062.648*sin(a0)/ss;
bij=-2062.648*cos(a0)/ss;
a=2*rr-1;
b=a+1;
f=2*tt-1;
g=f+1;
da=(a-1)*(n-a/2.0);
db=(b-1)*(n-b/2.0);
df=(f-1)*(n-f/2.0);
dg=(g-1)*(n-g/2.0);
q1=c(da+a)+c(df+f)-2.0*c(da+f);
q2=c(db+b)+c(dg+g)-2.0*c(db+g);
q3=c(da+b)+c(df+g)-c(da+g)-c(db+f);
x=q1-q2;
y=2.0*q3;
z=q1+q2;
qs=q1*cos(a0)^2+q2*sin(a0)^2+q3*sin(2.0*a0);
qa=q1*aij*aij+q2*bij*bij+2.0*q3*aij*bij;
ms=sqrt(qs)*uw;
maa=sqrt(qa)*uw;
return