《地面坐标系中某一点的地形矢在北东坐标系中的投影》由会员分享，可在线阅读，更多相关《地面坐标系中某一点的地形矢在北东坐标系中的投影（4页珍藏版）》请在金锄头文库上搜索。
1、一、计算地面坐标系中某一点m的地形矢 在北东坐标系中的投影？1、地面坐标系：原点O在厦门，其经纬度从网上查询为北纬24度27 分、东经118度 6分，天文瞄准方位角为45度。所以 , , ，a=24.5TB.8T45TA6378140, b=6356755北东坐标系：其原点m在地面坐标系中的坐标为 =(10km,10km,0)OMgVT；2、求在地面坐标系下的 gsOVO点的地心纬度为： 1.5sin(2)TTB所以根据 22|*(ta)/tas ssbb可求得 在地心坐标系中的长度，从而求得在地心坐标系中的坐标：s(|co,|coin,|i)sS S SOOsTOsTOsVVV再在进行坐标变
2、换，把 通过坐标变换矩阵Cgs得到 ：s gsOg=Cs*ssOOS其中， 213()()2sTTABC3、求地心大地直角坐标系下的 sSOMV有 ，在通过坐标变换矩阵Cgs得到 ：SOMggOsgVV sSOMVs=*SSMC4、求在东北坐标系中的 sm根据在地心大地直角坐标系下所求的 的向量求得如下：sSOMV22ss=arcos(qt(1)+()/SSsmOMOV 2sss/1()SSSs Ort5、最后求出在北天东坐标系中的 =n*SSMMsVC其中， 213()(snsmsC仿真的结果为： = SOMnV0m, 6384547.7m, 0.1072432780m基于如上解题过程，编写
3、Matlab程序如下：a=6378140; %地球长轴b=6356755; %地球短轴lamda_T=118.1/180*pi; %发射点O天文经度B_T=24.25/180*pi; %发射点O地理纬度A_T=pi/4; %天文瞄准方位角phi_s=B_T-11.5/60/180*pi*sin(2*B_T); %发射点O地心纬度Vom_g=10000,10000,0; %向量Om在地面坐标系中的坐标OsO=a*b*sqrt(1+tan(phi_s)2)/(a*tan(phi_s)2+b2);%地心坐标系下OsO的长度VOsO_s=OsO*cos(phi_s)*cos(lamda_T);OsO*
4、cos(phi_s)*sin(lamda_T);OsO*sin(phi_s); %向量OsO 在地心坐标系中的坐标 VOsO_g=Cgs(lamda_T,B_T,A_T)*VOsO_s; %向量OsO 在地面坐标系中的坐标VOsM_g=VOsO_g+Vom_g; %向量OsM在地面坐标系中的坐标VOsM_s=Cgs(lamda_T,B_T,A_T)*VOsM_g; %向量OsM在地心大地直角坐标系中的坐标OsM=sqrt(VOsM_s(1)2+VOsM_s(2)2+VOsM_s(3)2);%地心坐标系下OsM的长度phi_sm=acos(sqrt(VOsM_s(1)2+VOsM_s(2)2)/
5、OsM);%东北坐标系中的phi_smlamda_sm=acos(VOsM_s(1)/sqrt(VOsM_s(1)2+VOsM_s(2)2); %东北坐标系中的 lamda_smVOsM_n=Csn(lamda_sm,phi_sm)*VOsM_s; %向量OsM在北天东坐标系中的坐标 VOsM_n=vpa(VOsM_n,8) %保留8位有效数字编写Cgs函数如下：（即地心大地直角坐标系转换成地面坐标系的转换矩阵）function result=Cgs(lamda_T,B_T,A_T)a=B_T;b=-pi/2-A_T;r=lamda_T-pi/2;C1=1 0 0;0 cos(a) sin(a
6、);0 -sin(a) cos(a);C2=cos(b) 0 -sin(b);0 1 0;sin(b) 0 cos(b);C3=cos(r) sin(r) 0;-sin(r) cos(r) 0;0 0 1;result=C2*C1*C3;end编写Csn 函数如下： (即地心大地直角坐标系转换成地面坐标系的转换矩阵)function result=Csn(lamda_sm,phi_sm)a=pi/2-lamda_sm;b=pi/2;r=phi_sm;C1=1 0 0;0 cos(a) sin(a);0 -sin(a) cos(a);C2=cos(b) 0 -sin(b);0 1 0;sin(b) 0 cos(b);C3=cos(r) sin(r) 0;-sin(r) cos(r) 0;0 0 1;result=C2*C1*C3;end仿真结果如下：VOsM_n =0.6384547.70.10724327e-80
