GPS从入门到放弃(八) --- GPS卫星速度解算
GPS从入门到放弃(八) — GPS卫星速度解算
在阅读这一篇强烈建议先阅读GPS卫星位置解算。
为了计算卫星速度,需要对卫星的位置求导。
计算各参变量对时间的导数如下:
M˙k=n\dot{M}_k = nM˙k=n
E˙k=M˙k1−ecosEk\dot{E}_k = \frac{\dot{M}_k}{1-e\cos E_k}E˙k=1−ecosEkM˙k
ν˙k=1−e2E˙k1−ecosEk\dot{\nu}_k = \frac{\sqrt{1-e^2}\dot{E}_k}{1-e\cos E_k}ν˙k=1−ecosEk1−e2E˙k
Φ˙k=ν˙k\dot{\Phi}_k = \dot{\nu}_kΦ˙k=ν˙k
δu˙k=2Φ˙k(Cuscos2Φk−Cucsin2Φk)\delta\dot{u}_k = 2\dot{\Phi}_k(C_{us}\cos2\Phi_k - C_{uc}\sin2\Phi_k)δu˙k=2Φ˙k(Cuscos2Φk−Cucsin2Φk)
δr˙k=2Φ˙k(Crscos2Φk−Crcsin2Φk)\delta\dot{r}_k = 2\dot{\Phi}_k(C_{rs}\cos2\Phi_k - C_{rc}\sin2\Phi_k)δr˙k=2Φ˙k(Crscos2Φk−Crcsin2Φk)
δi˙k=2Φ˙k(Ciscos2Φk−Cicsin2Φk)\delta\dot{i}_k = 2\dot{\Phi}_k(C_{is}\cos2\Phi_k - C_{ic}\sin2\Phi_k)δi˙k=2Φ˙k(Ciscos2Φk−Cicsin2Φk)
Ω˙k=Ω˙−Ω˙e\dot{\Omega}_k = \dot{\Omega} - \dot{\Omega}_eΩ˙k=Ω˙−Ω˙e
i˙k=i˙+δi˙k\dot{i}_k = \dot{i} + \delta\dot{i}_ki˙k=i˙+δi˙k
r˙k=AeE˙ksinEk+δr˙k\dot{r}_k = Ae\dot{E}_k\sin E_k + \delta\dot{r}_kr˙k=AeE˙ksinEk+δr˙k
u˙k=Φ˙k+δu˙k\dot{u}_k = \dot{\Phi}_k + \delta\dot{u}_ku˙k=Φ˙k+δu˙k
计算卫星在轨道平面内的速度和WGS-84坐标系中的速度如下:
x˙k′=r˙kcosuk−rku˙ksinuk\dot{x}'_k = \dot{r}_k\cos u_k - r_k \dot{u}_k \sin u_kx˙k′=r˙kcosuk−rku˙ksinuk
y˙k′=r˙ksinuk+rku˙kcosuk\dot{y}'_k = \dot{r}_k\sin u_k + r_k \dot{u}_k \cos u_ky˙k′=r˙ksinuk+rku˙kcosuk
x˙k=−ykΩ˙k−(y˙k′cosik−zki˙k)sinΩk+x˙k′cosΩk\dot{x}_k = -y_k\dot{\Omega}_k - (\dot{y}'_k\cos i_k - z_k\dot{i}_k)\sin\Omega_k + \dot{x}'_k\cos\Omega_kx˙k=−ykΩ˙k−(y˙k′cosik−zki˙k)sinΩk+x˙k′cosΩk
y˙k=xkΩ˙k+(y˙k′cosik−zki˙k)cosΩk+x˙k′sinΩk\dot{y}_k = x_k\dot{\Omega}_k + (\dot{y}'_k\cos i_k - z_k\dot{i}_k)\cos\Omega_k + \dot{x}'_k\sin\Omega_ky˙k=xkΩ˙k+(y˙k′cosik−zki˙k)cosΩk+x˙k′sinΩk
z˙k=y˙k′sinik+yk′i˙kcosik\dot{z}_k = \dot{y}'_k\sin i_k + y'_k\dot{i}_k\cos i_kz˙k=y˙k′sinik+yk′i˙kcosik
在之前的位置解算的过程中,我们已经求得了很多中间变量,利用以上公式,即可求得卫星速度。
代码如下:
M_k_Dot = n
E_k_Dot = M_k_Dot/(1-e*math.cos(E_k))
nu_k_Dot = math.sqrt(1-e**2)*E_k_Dot/(1-e*math.cos(E_k))
Phi_k_Dot = nu_k_Dot
delta_u_k_Dot = 2*Phi_k_Dot*(Cus*math.cos(2*Phi_k) - Cuc*math.sin(2*Phi_k))
delta_r_k_Dot = 2*Phi_k_Dot*(Crs*math.cos(2*Phi_k) - Crc*math.sin(2*Phi_k))
delta_i_k_Dot = 2*Phi_k_Dot*(Cis*math.cos(2*Phi_k) - Cic*math.sin(2*Phi_k))
Omega_k_Dot = Omega_Dot - Omega_e_Dot
i_k_Dot = i_Dot + delta_i_k_Dot
r_k_Dot = A*e*E_k_Dot*math.sin(E_k) + delta_r_k_Dot
u_k_Dot = Phi_k_Dot + delta_u_k_Dotx_p_k_Dot = r_k_Dot*math.cos(u_k) - r_k*u_k_Dot*math.sin(u_k)
y_p_k_Dot = r_k_Dot*math.sin(u_k) + r_k*u_k_Dot*math.cos(u_k)x_k_Dot = -y_k*Omega_k_Dot - (y_p_k_Dot*math.cos(i_k) - z_k*i_k_Dot)*math.sin(Omega_k) + x_p_k_Dot*math.cos(Omega_k)
y_k_Dot = x_k*Omega_k_Dot + (y_p_k_Dot*math.cos(i_k) - z_k*i_k_Dot)*math.cos(Omega_k) + x_p_k_Dot*math.sin(Omega_k)
z_k_Dot = y_p_k_Dot*math.sin(i_k) + y_p_k*i_k_Dot*math.cos(i_k)print("x_k_Dot={}".format(x_k_Dot))
print("y_k_Dot={}".format(y_k_Dot))
print("z_k_Dot={}".format(z_k_Dot))
用GPS卫星位置解算中的数据进行解算,可以得到卫星在WGS-84坐标系中的速度为:
[x˙ky˙kz˙k]=[1088.40932985282189.0005213284817−1527.1796186285617](米/秒)\left[ \begin{array}{r} \dot{x}_k \\ \dot{y}_k \\ \dot{z}_k \end{array} \right]= \left[ \begin{array}{c} 1088.4093298528 \\ 2189.0005213284817 \\ -1527.1796186285617 \end{array} \right](米/秒) ⎣⎡x˙ky˙kz˙k⎦⎤=⎣⎡1088.40932985282189.0005213284817−1527.1796186285617⎦⎤(米/秒)
GPS从入门到放弃(八) --- GPS卫星速度解算相关推荐
- GPS从入门到放弃(七) --- GPS卫星位置解算
GPS从入门到放弃(七) - GPS卫星位置解算 上一篇讲了开普勒轨道参数,根据这些参数就可以确定卫星的位置,这一篇我们来实际计算一下. WGS-84基本参数 首先给出几个WGS-84坐标系中的基本参 ...
- 【转】GPS从入门到放弃(一) --- GPS基础原理
转自:https://blog.csdn.net/tyst08/article/details/100529424 GPS从入门到放弃(一) - GPS基础原理 GPS,全称Global Positi ...
- GPS从入门到放弃(二十六) --- RTKLIB函数解析
GPS从入门到放弃(二十六) - RTKLIB函数解析 为了贴合这个系列的标题"从入门到放弃",在入门之后现在就要放弃此方向了.虽然感觉遗憾,暂时也没有办法.在此附上此系列最后一篇 ...
- GPS从入门到放弃(六) --- 开普勒轨道参数
GPS从入门到放弃(六) - 开普勒轨道参数 开普勒轨道参数是用于描述卫星轨道的,GPS卫星的无摄椭圆轨道运动就是用开普勒轨道参数来描述.参考GPS导航电文,电文中的星历和历书都是用开普勒轨道参数来描 ...
- GPS从入门到放弃(二) --- GPS时间
GPS从入门到放弃(二) - GPS时间 上一篇GPS基础原理中提到了一个钟差的概念,并没有细说.GPS涉及到的时间有不少,这一篇就详细讲一讲. 时间是一个相对的概念,最简单的例子:你的手表的时间和我 ...
- GPS从入门到放弃(十二) --- 多普勒定速
GPS从入门到放弃(十二) - 多普勒定速 多普勒效应 多普勒效应在我们日常生活中有很多,比如当一辆救护车迎面驶来的时候,听到声音比原来高:而车离去的时候声音的音高比原来低. 这个效应是为纪念奥地利物 ...
- GPS从入门到放弃(十一) --- 差分GPS
GPS从入门到放弃(十一) - 差分GPS 原理 差分GPS是一种通过消除或减小GPS测量误差从而提高GPS定位精度的方法. 差分GPS消除或减小测量误差的基础原理是利用误差的相关性.以卫星时钟误差举 ...
- GPS从入门到放弃(三) --- GPS坐标系
GPS从入门到放弃(三) - GPS坐标系 定位就需要坐标,坐标当然是相对坐标系而言的,我们描述一个物体的位置,首先就需要建立坐标系. 按大类来分,坐标系可以分为惯性坐标系和非惯性坐标系.惯性坐标系是 ...
- GPS从入门到放弃(十) --- 定位方程解算和定位精度
GPS从入门到放弃(十) - 定位方程解算和定位精度 上一篇伪距与载波相位中我们介绍了伪距的计算方法,也得到了包含 (x,y,z,δt)(x,\ y,\ z,\ \delta_t)(x, y, z, ...
最新文章
- Centos7 安装 telnet 服务
- python中的匿名函数与lambda
- iosanimationWithKeyPath
- java-容器-ArrayList
- NeurIPS 2021 助力YOLOv5涨点 Alpha-IoU:IoU Loss大一统
- linux命令行怎么注释,Bash Shell 注释多行的几种方法
- Vaadin提示:以声明方式构建UI
- ubuntu PPA
- 插入排序InsertSort
- 密码学:一文读懂常用加密技术原理及其逻辑与应用方法
- 华为IPD你学不会(转)
- 人体性神经系统分布图,女性神经分布图高清
- Word修改标题样式缩进不起作用原因
- 计算机DSU,研科d5000dsu驱动
- SpringBoot项目打包成jar后,使用ClassPathResource获取classpath(resource)下文件失败
- 餐饮管理系统哪个好用?
- 张小龙-年薪近3亿的微信之父,他是如何做到的?
- 【华为云CCE】 k8s部署nacos集群
- LeetCode题解(0657):根据操作列表判断机器人是否能回到原点(Python)
- 蓝海卓越NSWMP管理协议简介