哪位大侠能提供一个秩亏自由平差改正数和精度评定的源程序呢？

flywoud

谢了！自己编了一个，但觉得结果不正确，哪位大下给修改一下？
load shujua1.txt;load shujuH0.txt;load shujuh.txt;load shujup.txt;
A=shujua1;H0=shujuH0;h=shujuh;p=shujup;
p=diag(p);
l=A*H0-h;
l=1000*l
A1=A(1:7,1:4);A2=A(1:7,5);
ATPL=A1'*p*l;
N11=A1'*p*A1
N12=A1'*p*A2
N21=N12'
N11*N11+N12*N21
invNNNN=inv(N11*N11+N12*N21)
K=-invNNNN*ATPL
N1112=[N11;N21]
N2111=[N11 N12]
X=N1112*K
V=A*X+l
m0=sqrt(V'*p*V/(24-15))
Q=N1112*invNNNN*N11*invNNNN*N2111
Qxx=diag(Q)
m=m0*sqrt(Qxx)