测量平差MATLAB方法

条件平差

P=input('请输入权阵 ');
A=input('请输入矩阵A ');
w=input('请输入矩阵w ');
n=A*inv(P)*A.';
k=-inv(n)*w;
v=inv(P)*A.'*k
disp('单位权中误差的估值为 ');
sqrt(v.' *P*v/size(w,1))
disp('观测值平差值的协因数阵为')
Qh=inv(P)-inv(P)*A.'*inv(n)*A*inv(P)
f=input('请输入平差值函数的权函数 ')
disp('平差值的权倒数为');
q=f*Qh*f.'

有参数的条件平差

P=input('请输入权阵 ');
A=input('请输入矩阵A');
B=input('请输入矩阵B');
w=input('请输入矩阵w ');
[c,u]=size(B);
Q=inv(P);
Na=A*Q*A.';
Nb=B.'*inv(Na)*B;
disp('参数解算')
x=-inv(Nb)*B.'*inv(Na)*w
K=-inv(Na)*(B*x+w);
disp('改正数v ');
v=Q*A.'*K
disp('单位权中误差的估值o为 ');
o=sqrt(v.' *P*v/(c-u))
disp('观测值平差值的协因数阵为');
Qll=Q-Q*A.'*(inv(Na)-inv(Na)*B*inv(Nb)*B.'*inv(Na))*A*Q
disp('观测值平差值的中误差oLL为');
oLL=o*sqrt(diag(Qll))

间接平差

P=input('请输入权阵 ');
B=input('请输入矩阵B');
l=input('请输入矩阵误差方程常数项向量值l ');
[n,t]=size(B);
Nb=B.'*P*B;
W=B.'*P*l;
disp('参数解算')
x=inv(Nb)*W
disp('改正数v ');
v=B*x-l
disp('单位权中误差的估值o为 ');
o=sqrt(v.' *P*v/(n-t))
disp('观测值平差值的协因数阵为');
QLL=B*inv(Nb)*B.'

有限制条件的间接平差

P=input('请输入权阵 ');
B=input('请输入矩阵B');
c=input('请输入矩阵c');
w=input('请输入矩阵w');
l=input('请输入矩阵误差方程常数项向量值l ');
[n,u]=size(B);
s=size(c,1);
Nb=B.'*P*B;
W=B.'*P*l;
Nc=c*inv(Nb)*c.';
K=inv(Nc)*(c*inv(Nb)*W+w);
disp('参数解算')
x=inv(Nb)*(c.'*K-W)
disp('改正数v ');
v=B*x-l
disp('单位权中误差的估值o为 ');
o=sqrt(v.' *P*v/(n-u+s))
disp('参数平差值的协因数阵为');
Qx=inv(Nb)-inv(Nb)*c.'*inv(Nc)*c*inv(Nb)

（本方法不提供观测方程的线性结算）系数矩阵，常数矩阵需要自行解决。