卫星导航GPS典型例题编程报告讲述.doc

GPS Matlab编程报告 一、通过星历计算卫星坐标（修正） 计算归化时间 计算计算观测历元的平近点角 计算偏近点角 利用不动点迭代法求解此方程： 得出 不动点迭代法的迭代公式为： 计算卫星的地心矢径 计算真近点角 计算升交点角距 计算摄动改正项： 计算经过摄动改正的升交点角距，卫星矢径，和轨道面倾角 计算观测历元的升交点经度 计算卫星在轨道平面坐标系中的位置 计算卫星在地固坐标系下的坐标 Matlab源程序clear; %---------t时刻，卫星星历数据--------- t=1.728128741751984e+005; %当前时间（接收到卫星信号时间） Omega0= -3.1122016819; %按参考时间计算的升交点经度 I0=0.3*pi; %非GEO卫星 %I0= 0.0879969068; %GEO卫星 SqrtA= 6493 %长半轴平方根 Ecc= 0 %偏心率 Small_Omega= -1.7409228484; %近地点角距 Mu0= 0.0005884201; %参考时间的平近点角 Delta_n= -2.4123e-009; %卫星平均运动速度与计算值之差 I_Dot= -3.4333e-010; %轨道倾角变化率 Omega_dot= 3.47676e-009; %升交点经度变化率 C_uc= 4.692e-006; %纬度幅角的余弦调和改正项的振幅 C_us= -1.10036e-005; %纬度幅角的正弦调和改正项的振幅 C_ic= 9.01e-008; %轨道倾角的余弦调和改正项的振幅 C_is= 2.01e-008; %轨道倾角的正弦调和改正项的振幅 C_rc= 325.73368; %轨道半径的余弦调和改正项的振幅 C_rs= 146.91074; %轨道半径的正弦调和改正项的振幅 Toe= 172800; %星历参考时间 IODC= 10; %钟差数据龄期 URAI= 0; %用户距离精度标志 IODE= 10; %星历数据龄期 Toc = 172800; %本时段钟差参数参考时间 a0 = -2.72893e-005; %卫星钟差改正0阶多项式系数 a1 = -6.7531e-014; %卫星钟差改正1阶多项式系数 a2 = 5.787e-018; %卫星钟差改正2阶多项式系数 %---------定义常量--------- c=28; %光速 omegae=7.2921151467e-5; %地球自转角速度 mu=3.986004418e14; %地球引力常数GM %---------1、计算GPS卫星运行的平均速度n--------- a=SqrtA*SqrtA; n=sqrt(mu/(a^3))+Delta_n; %---------2、计算归化时间Delta_t--------- Delta_t=t-Toe; %---------3、计算观测历元t的平近点角M--------- M=Mu0+n*Delta_t; %---------4、计算偏近点角E--------- eps=1e-20; E=M; tol=1; while (toleps) %不动点迭代法 E0=E; E=M+(Ecc)*sin(E0); tol=abs(E-E0); end %---------5、计算卫星的地心矢径r0--------- r0=a*(1-Ecc*cos(E)); %---------6、计算真近点角f--------- %f=2*atan(sqrt((1+Ecc)/(1-Ecc))*tan(E/2)); f=atan((sqrt(1-Ecc^2)*sin(E))/(cos(E)-Ecc)); %---------7、计算升交点角距Phi0--------- Phi0=Small_Omega+f; %---------8、计算摄动改正项：Delta_u,Delta_r,Delta_i-----