测量平差MATLAB方法-Toy模板网

这篇具有很好参考价值的文章主要介绍了测量平差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)

（本方法不提供观测方程的线性结算）系数矩阵，常数矩阵需要自行解决。文章来源地址https://www.toymoban.com/news/detail-509499.html

到了这里，关于测量平差MATLAB方法的文章就介绍完了。如果您还想了解更多内容，请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章，希望大家以后多多支持TOY模板网！