测量平差第一章.docVIP

测量平差第一章 3 .题目：独立观测值、的中误差分别为1、2，求下列函数的中误差： 解答： 因、是独立观测值，则 ∴ 对此函数式进行全微分，得： 对此函数式进行全微分，得： 所以 4 .解答： Matlab 实现如下： function [q]=No4(p) [mp,np]= size(p); if mp~=np warning( P error); return end q=inv(p); end 5 . 题目：已知观测量及其协方差阵，组成一下函数 为常数阵。 求协方差阵、和。 解答： 6 .解答： 由题意可知: 7 .题目：设P1点及P2点的坐标为： ， 向量的协方差阵为： 试求坐标差函数与的方差—协方差阵； 求两点间边长与坐标方位角的方差—协方差阵。 解答： 向量的协方差阵得 同理 ∴ 对函数式进行全微分： 由协方差传播方程 = = 所以方差—协方差阵： MATLAB程序编码样例 （以该题第一问第一个公式为例，其他与之类似）： save saf k1=[1 -1];k2=k1; D=[3 2 2 0;2 4 -1 3;2 -1 6 -2;0 3 -2 8] D = 3 2 2 0 2 4 -1 3 2 -1 6 -2 0 3 -2 8 d2x1= D(1,1);d2x2=D(3,3);dx1x2=dx2x1=D(1,2); d=k1*[d2x2 dx1x2;dx1x2 d2x1]*k2 d = 5 clc 8 .解答： Matlab实现如下： function [Q1]=No8(a,q,b) nl=length(a); n2=length(b); [mp,np]= size(q); if mp~=nl warning(q error); return end if mp~=n2 warning(q error); return end Q1=a*q*b; end 10 .解答： 11 .题目：设半侧回观测一次方向中误差为，试问一测回的角度中误差等于多少。 解答：设上半侧回测得角度，下半测回，A为方向角观测值 同理 所以一测回的角度中误差 12 ．解答： 解 设路线总长S公里，按照测量学上的附合路线计算步骤，则路线闭合差 由于是路线中点，故 则线路中点高程 设每公里高差观测中误差为，则按误差传播定律 Matlab实现如下： function Ox=No12(a,b) nl=length(a); [mp,np]= size(b); if mp~=nl warning(b error); return end Ox=a*b*a; end 13. 题目： 有一角度测２０测回，得中误差±０．４２＇＇，问再增加多少测回，其中误差为±０.２８＇＇。 解答：　 设每个测回的中误差为σ0，需要再增加n个测回，则： , , 由上式可解出n,即： 14. 解答： 15 .题目：已知独立观测值的方差阵及单位权方差，试求权阵PLL及权P1P2。 解答： ∴ 16 解答：由题意可求出即： 同理可得： Matlab实现如下: function Qff=No16(Pa,Pb,L,F) %///Pa A的权阵 Pb B的权阵 L 观测值 F函数系数 nl=length(Pa); L1=diag(L); [ma,na]=(L1); if ma~= n1 warning (L error); return end if na~= n1 warning (L error); return end Qxx=Pa*L1*Pa; Qyy=Pb*L1*Pb; Qxy=Pa*L1*Pb; Qtt=[Qxx,Qxy;Qxy,Qyy]; Qff=F*Qtt*F; 17 .题目： 已知观测值向量，其协因数阵为单位阵。有方程： 式中：为已知的系数阵， 为可逆矩阵。 求协因数阵、; 证明与和 均互不相关。 解答： （1）的协因数阵为单位阵。 (2) 因为 所以 因为 所以 18 .解： ⑴每