平差计算过程

《平差计算过程》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、一秩亏自由网平差N=B'*BN=4-1-1-1-10-13-1-100-1-1300-1-1-10200-10001000-1001>>NN=N*NNN=20-5-6-5-51-512-5-411-6-51221-4-5-42610-51112011-4002>>NM=pinv(NN)NM=0.10190.0386-0.07210.08240.1471-0.29780.03860.1523-0.00530.1284-0.1246-0.1893-0.0721-0.00530.1731-0.0995-0.29780.30160.08240.1284-0.09950.336

2、3-0.0600-0.38770.1471-0.1246-0.2978-0.06001.0255-0.6902-0.2978-0.18930.3016-0.3877-0.69021.2634>>NNM=N*NMNNM=0.21180.0035-0.05900.02430.0451-0.22570.00350.2951-0.01740.0660-0.1632-0.1840-0.0590-0.01740.2951-0.1215-0.22570.12850.02430.0660-0.12150.4618-0.1424-0.28820.0451-0.1632-0.2257-0

3、.14240.8785-0.3924-0.2257-0.18400.1285-0.2882-0.39240.9618>>X=NNM*B'*L（即：x=N-m*Bt*l）X=0.35420.10422.2292-0.2708-2.64580.2292>>V=B*X-0.001*LV=0.62402.99701.8750-1.9980-2.1210-0.25200.3750二拟稳平差>>B2=[B(1:7,4:6)]B2=100100-10100-101-1-110010]B1=[B(1:7,1:3)B1=-1000-10000001000000-100>>N=B'*BN

4、=4-1-1-1-10-13-1-100-1-1300-1-1-10200-10001000-1001>>N22=[N(4:6,4:6)]N22=4-1-1-13-1-1-13>>N21=[N(4:6,1:3)]N21=-1-10-10000-1>>N11=[N(1:3,1:3)]N11=200010001>>N12=[N(1:3,4:6)]N12=-1-10-10000-1>>N111=pinv(N11)（N11的逆矩阵）>>M1=N21*N11*N12（为））>>M=N22-M1（为））>>K1=B2'-N21*N111*B1'（K1即为αT）>>MM=M*M>>

5、MMM=pinv(MM)>>MMMM=M*MMM(MMMM即为Mm->>K2=MMMM*K1（K2即为α-）>>KP=B1'-N12*K2（kp即为）>>K3=N111*KP（K3即为β）>>X1=K3*L结果为：V=-0.6260-0.37800.50002.0020-2.8710-2.3770-3.0000三转换H1=I-ST'*ST*PX/3H1=1.000000-0.3333-0.3333-0.333301.00000-0.3333-0.3333-0.3333001.0000-0.3333-0.3333-0.33330000.6667-0.3333-0.333

6、3000-0.33330.6667-0.3333000-0.3333-0.33330.6667>>XN=H1*XXN=1.25001.00003.12500.6250-1.75001.1250>>V=B*XN-0.001*LV=0.62402.99701.8750-1.9980-2.1210-0.25200.3750残差v的结果与拟稳平差的结果一致四假设一点a的高程，利用经典平差，进行平差计算。>>N=A*Q*A'N=3223>>K=inv(N)*WK=2.00000>>V=Q*A'*KV=002.000002.0000-2.00000