一种单星伪距的测姿标校方法技术

本发明专利技术具体公开了一种单星伪距的测姿标校方法，所述方法包括：S1、基于卫星导航接收机获取伪距观测值，并求出所述伪距观测值的二阶导数；S2、计算出所述二阶导数中的直流分量；S3、基于所述直流分量计算出所述二阶导数中的交流分量；S4、对所述交流分量进行频域分析以找出最大频谱值对应的信号点，进而获取旋转载体的旋转频率，从而估算出旋转载体的姿态。本发明专利技术首先对伪距观测值进行二阶求导，然后求出二阶导数中的直流分量和交流分量，并基于交流分量做频域分析找出最大频谱值对应的信号点，最后获取旋转载体的旋转频率，进而估算出旋转载体的姿态，既能够有效降低成本，又保证了转速测量的准确性。速测量的准确性。速测量的准确性。

【技术实现步骤摘要】
一种单星伪距的测姿标校方法


[0001]本专利技术涉及卫星导航
，尤其涉及一种单星伪距的测姿标校方法。

技术介绍

[0002]目前，获取飞行器姿态信息对飞行器的飞行轨迹估计以及飞行过程控制有着重大意义，若是飞行器姿态信息不准确，将严重影响到控制系统的控制精度，甚至导致飞行器不稳定。而部分飞行器为保证其飞行稳定性，会在其径向运动时保持自身旋转，而旋转状态给飞行器载体姿态测量带来了一定程度的影响。
[0003]现有获取飞行器姿态信息的手段包括惯性导航系统、地磁传感器、陀螺仪、卫星导航接收机等，均能有效为飞行器提供后续控制机构所需的信息。但各类传感器在载体高速自旋过程中，其适应性有所区别。地磁传感器因其耐高过载、功耗低等优点在姿态测量中被广泛使用，地磁传感器具有转速波形相位关系明确，转速适应范围宽，检测得到的转速变化精细、稳定等特点，但是存在使用成本高，易受电机及外界磁场变化的影响以及弹体、发射器等铁磁物质的磁屏蔽干扰，同时由于地磁传感器成本造价较高、测姿算法相对复杂，如果与卫星导航接收机组合工作，又增加了器件机械结构的复杂度。
[0004]鉴于此，研究一种成本低、标定便捷且测量准确度高的单星伪距测姿标校方法是本
人员需要解决的技术问题。

技术实现思路

[0005]为解决上述技术问题，本专利技术提供的一种单星伪距的测姿标校方法，所述方法包括以下步骤：
[0006]S1、基于卫星导航接收机获取伪距观测值，并求出所述伪距观测值的二阶导数；
[0007]S2、计算出所述二阶导数中的直流分量；
[0008]S3、基于所述直流分量计算出所述二阶导数中的交流分量；
[0009]S4、对所述交流分量进行频域分析以找出最大频谱值对应的信号点，进而获取旋转载体的旋转频率，从而估算出旋转载体的姿态。
[0010]上述的单星伪距的测姿标校方法，所述步骤S1中所述卫星伪距观测值的二阶导数计算公式为：
[0011][0012]式(1)中，表示卫星伪距观测值P(t)的一阶导数，t表示时域采样点。
[0013]上述的单星伪距的测姿标校方法，所述步骤S2中所述直流分量的计算公式为：
[0014][0015]式(2)中，A表示直流分量，K表示与旋转载体的滚转周期和信号采样率有关的历元
数，i表示时刻。
[0016]上述的单星伪距的测姿标校方法，所述步骤S3中交流分量的计算公式为：
[0017]C(t)＝D(t)
‑
A
ꢀꢀꢀ
(3)
[0018]式(3)中，C(t)表示交流分量。
[0019]上述的单星伪距的测姿标校方法，所述步骤S4中对所述交流分量进行频域分析用公式表示为：
[0020]C(jkω)＝F[C(t)]ꢀꢀꢀ
(4)
[0021]式(4)中，t表示时域采样点，k表示频域值的索引，ω表示频率，j表示复数中的虚数单位。
[0022]上述的单星伪距的测姿标校方法，所述步骤S4中找出最大频谱值对应的信号点用公式表示为：
[0023][0024]式(5)中，f
r
表示旋转载体的转速，f
s
表示采样频率，m表示最大频谱值对应的信号点，N表示FFT运算点数。
[0025]与现有技术相比，本专利技术提出了一种基于卫星导航接收机接收的单星伪距测姿标校方法，该方法对获取的卫星伪距观测值进行二阶求导，然后求出所述二阶导数中的直流分量和交流分量，进而基于交流分量做频域分析找出最大频谱值对应的信号点，根据最大频谱值对应的信号点获取旋转载体的旋转频率，进而估算出旋转载体的姿态，本专利技术中的单星伪距测姿标校方法既能够有效降低成本，又保证了转速测量的准确性。
附图说明
[0026]图1是本专利技术中一种单星伪距的测姿标校方法的流程图，
[0027]图2是本专利技术中旋转载体测向单天线旋转矢量在视线矢量上投影的示意图，
[0028]图3是本专利技术仿真实验中一种场景下旋转载体旋转速度估计结果的示意图，
[0029]图4是本专利技术仿真实验中另一种场景下旋转载体旋转速度估计结果的示意图，
[0030]图5是本专利技术仿真实验中第三种场景下旋转载体旋转速度估计结果的示意图，
[0031]图6是本专利技术仿真实验中第四种场景下旋转载体旋转速度估计结果的示意图。
具体实施方式
[0032]为了能够更清楚地理解本专利技术的上述目的、特征和优点，下面结合附图和具体实施方式对本专利技术进行进一步的详细描述。需要说明的是，在不冲突的情况下，本申请的实施例及实施例中的特征可以相互组合。
[0033]如图1所示，图1示出了一种单星伪距的测姿标校方法的具体流程。
[0034]在其中一个实施例中，一种单星伪距的测姿标校方法，所述方法包括以下步骤：
[0035]S1、基于卫星导航接收机获取伪距观测值，并求出所述伪距观测值的二阶导数；
[0036]本步骤中，所述卫星伪距观测值的二阶导数计算公式为：
[0037][0038]式(1)中，表示卫星伪距观测值P(t)的一阶导数，t表示时域采样点，如图2所示，O1表示天线安装位置，P为旋转载体质心，S为卫星位置。随着旋转载体旋转的卫星接收机所测得的伪距观测值P(t)包含两部分：一部分是旋转载体质心P到卫星位置S的位置矢量在视线矢量上的投影因旋转载体与卫星的相对运动实时发生变化；另一部分是卫星接收机天线在旋转载体横截面上旋转时位置矢量在视线矢量上的投影因卫星接收机天线的圆周运动呈正弦变化，其振荡频率与旋转载体的旋转速度一致。
[0039]对伪距观测值P(t)求导，即可得到卫星伪距观测值P(t)的一阶导数
[0040]对一阶导数P
′
(t)再次求导，得到卫星伪距观测值的二阶导数
[0041]伪距测量值P(t)中的它与旋转载体质心和卫星的相对运动速度、加速度相关，与旋转载体旋转无关，即在一定时间内，变化量为旋转载体质心和卫星的相对运动距离值，其一阶差分量为两者相对运动速度，二阶差分量为两者相对运动的加速度，当选取的观测时间段在旋转载体自由落体阶段，其加速度恒定变化时，该二阶差分量为常数，称之为直流分量；
[0042]伪距测量值P(t)中的它与旋转载体和卫星的相对运动无关，其一阶差分量为卫星接收机在旋转载体表面圆周运动的速度在视线矢量上的投影，二阶差分量为圆周运动的加速度在视线矢量上的投影，两者变化规律均呈正弦波动，称之为交流分量；
[0043]S2、计算出所述二阶导数中的直流分量；
[0044]本步骤中，所述直流分量的计算公式为：
[0045][0046]式(2)中，A表示直流分量，K表示与旋转载体的滚转周期和信号采样率有关的历元数，i表示时刻；通常情况下，历元数K应覆盖旋转载体滚转的多个周期，且覆盖滚转周期数量越多，在滚转频率估计时所获得的频谱能量越大，对滚转频率的估计就越准确，但是，由于滚转频率估本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种单星伪距的测姿标校方法，其特征在于，所述方法包括以下步骤：S1、基于卫星导航接收机获取伪距观测值，并求出所述伪距观测值的二阶导数；S2、计算出所述二阶导数中的直流分量；S3、基于所述直流分量计算出所述二阶导数中的交流分量；S4、对所述交流分量进行频域分析以找出最大频谱值对应的信号点，进而获取旋转载体的旋转频率，从而估算出旋转载体的姿态。2.如权利要求1所述的基于机器学习的单星伪距的测姿标校方法，其特征在于，所述步骤S1中所述卫星伪距观测值的二阶导数计算公式为：式(1)中，表示卫星伪距观测值P(t)的一阶导数，t表示时域采样点。3.如权利要求2所述的基于机器学习的单星伪距的测姿标校方法，其特征在于，所述步骤S2中所述直流分量的计算公式为：式(2)中，A表示直流分量，K表示与旋转载体的滚转周期和信号采样率有关的历元数，i表示时刻。4.如权利要求3所述的基于机器学习的单星伪距...

【专利技术属性】
技术研发人员：吴鹏，冯璐，黄雅静，仝海波，陈雨莺，肖学美，方新怡，颜锦辉，
申请(专利权)人：长沙学院，
类型：发明
国别省市：