空间后方交会报告.doc

《空间后方交会报告.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、任务已知f=153.24mm，m=10000，限差0.1’各点坐标点号像点坐标地面坐标x（mm）y（mm）X（m）Y（m）Z（m）1-86.15-68.9936589.4125273.322195.172-53.4082.2137631.0831324.51728.693-14.78-76.6339100.9724934.982386.50410.4664.4340426.5430319.81757.31求近似垂直摄影情况下后方交会解设计任务1、确定未知数的初始值：Φ0=ω0=К0=0,内方位元素x0=0.0，y0=0.0，f=153.24mm。Zs0=H+14

2、i=14Zi=mf+14i=14Zi=3049.3m；Xs0=14i=14Xi=38437m；Ys0=14i=14Yi=27963.16m2、计算旋转矩阵R利用角元素的近似值计算方向余弦值，组成R阵根据《摄影测量学》P32中的公式（3-9），初步计算R阵R[0][0]=cos(Φ)*cos(K)-sin(Φ)*sin(W)*sin(K);R[0][1]=-cos(Φ)*sin(K)-sin(Φ)*sin(W)*cos(K);R[0][2]=-sin(Φ)*cos(W);R[1][0]=cos(W)*sin(K);R[1][1]=cos(W)*cos(K);R[1

3、][2]=-sin(W);R[2][0]=sin(Φ)*cos(K)+cos(Φ)*sin(W)*sin(K);R[2][1]=-sin(Φ)*sin(K)+cos(Φ)*sin(W)*cos(K);R[2][2]=cos(Φ)*cos(W);得初始R阵1000100013、逐点计算近似值（x），（y）：带入《摄影测量学》P61的公式（5-1）；得(x1)1(y1)1(x2)1(y2)1(x3)1(y3)1(x4)1(y4)1-0.33147-0.48258-0.053220.221960.15351-0.700980.133020.157564、组成误差方程式

4、：按（5-8）；（5-9b）、（5-4）式逐点计算误差方程式的系数和常数项根据Lx=x-(x)；Ly=y-(y)得Lx1Ly1Lx2Ly2Lx3Ly3Lx4Ly40.245320.41359-0.00018-0.13975-0.168290.62347-0.12256-0.09313解得A阵为-0.000100.000056219-0.201678-0.038785-0.068990-0.00010.0000450209-0.00000387855-0.18430.08615-0.000100.0000348473-0.171848-0.28648-0.0822

5、10-0.0001-0.0000536479-0.0000028648-0.1973440.0534534-0.000100.000009645-0.154666-0.00739096-0.076630-0.00010.0000500065-0.000000739096-0.191560.01478-0.00010-0.00000682589-0.153954-0.004397920.064430-0.0001-0.0000420452-0.000000439792-0.18033-0.01046根据公式X=（ATA）-1ATL解除X阵，检验Φ1、ω1、К1并不

6、在限差0.1’之内;利用解得的Φ1、ω1、К1，带入2中循环；知道解出合格结果！1、结果，通过C++程序解得结果2、源程序代码#include#include#include#include#includeusingnamespacestd;constintn=6;voidinverse(doublec[n][n]);templatevoidtranspose(T1*mat1,T2*mat2,inta,intb);template<

7、typenameT1,typenameT2>voidmulti(T1*mat1,T2*mat2,T2*result,inta,intb,intc);templatevoidinput(T*mat,inta,intb);templatevoidoutput(T*mat,char*s,inta,intb);intmain(){ofstreamoutFile;cout.precision(5);doublex0=0.0,y0=0.0;doublefk=0.15324;//内方位元素doublem=10000;//估算比例尺d

8、oubleB[4][5]