MATLAB处理测量平差

B=input('请输入B的值:');
B
P=input('请输入P的值:');
P
l=input('请输入l的值:');
l
C=B';
N=C*P*B;
M=inv(N);
%待定参数x的平差值：
x=M*C*P*l;
x
%观测值平差值
V=B*x-l;
V
A=V';
E=A*P*V;
%单位权中误差
r=input('请输入多余观测r的值:');
H=E/r;
O=sqrt(H);
O
F=input('请输入F的值:');
F
G=F';
Q=F*M*G;
%函数的协方差阵
D=H*Q;
D
%函数的中误差阵