本发明涉及利用卫星导航系统对载体进行姿态解算的方法,具体涉及到一种双天线gps测姿方法。
背景技术:
为了满足不同应用平台对姿态信息的需求,目前已有大量的姿态测量设备问世。例如,用于空间载体姿态测量的星敏感器、地平跟踪器、太阳敏感器和地磁仪等;用于为陆上或水下载体提供航向的磁罗盘、电磁罗盘;用于为各种飞行器、陆上运输设备、舰船、潜器、空间载体提供姿态信息的惯性器件等等。然而,在对姿态角及姿态变化率进行测量的过程中,上述各种姿态测量设备均存在一定的问题。光学和红外线仪器主要用于低动态和静态平台监测,例如星敏感器可以同时得到多颗观测星的视赤经、赤纬以及对应在像平面的位置信息,解算出载体在惯性系中的三轴姿态,并且可靠性高,隐蔽性好,精度不随时间变化。但星敏感器易受天气、地形或者其它客观因素的影响,不能实时地给出姿态或方位值。目前,高精度的实时姿态测量系统大都采用由陀螺和加速度计组成的惯性系统来实现,这种系统是自主导航系统,不受外界干扰,具有良好的隐蔽性和较强的高频动态测量能力,但其结构复杂、价格昂贵,当工作时间较长时惯性器件误差会随时间积累而引起测姿精度的降低。
gps在全球范围内的三维精确定位、测时、测速和测定姿态的能力,使其成为先进的导航系统,受到了各国军方的重视,在卫星导航应用技术上投入了大量研究资源,以期提高武器装备的信息含量。gps测姿系统只需要利用低成本的接收机即可以提供较高精度的姿态信息,以此来取代造价昂贵的传统测姿设备,同时还可以完成运载体的定位和授时,并且受环境的影响小,能够长时间进行高精度测姿任务,这使得gps姿态测量成为近几年研究的热点之一。利用gps载波相位观测量进行定向和定姿的研究为gps开辟了又一个新的应用领域。
基于gps双天线测姿系统具有传统的姿态系统所没有的优点,不仅能提供载体的位置、航向、速度等基本信息,而且能提供载体的偏航角和俯仰角。但是,在现有技术gps双天线测姿方法中,没有对gps定位数据进行修正,也没有对整周模糊度进行限制,以此得到的测姿数据精度相对较低,并且,整周模糊度的搜索范围较大,需要花费较长的搜索时间。显然,现有技术gps双天线测姿方法存在着测姿数据精度相对较低和花费较长的整周模糊度搜索时间等问题。
技术实现要素:
为解决现有技术gps双天线测姿方法存在的测姿数据精度相对较低和花费较长的整周模糊度搜索时间等问题,本发明提出一种双天线gps测姿方法。
本发明双天线gps测姿方法,在测姿载体上呈几何分布设置两个gps信号接收天线且对应设置两台接收机,采用gps系统3r定位法确定天线ⅰ的位置坐标,将天线ⅰ的位置坐标定为概略坐标并对其进行修正,确定整周模糊度的值,然后,解算天线ⅱ的位置坐标,最后,根据定位出的天线ⅰ的位置坐标和天线ⅱ的位置坐标确定测姿载体的偏航角和俯仰角,即进行测姿载体的姿态角解算。
进一步的,本发明双天线gps测姿方法,包括以下步骤:
s1、在测姿载体上呈几何分布设置两个gps信号接收天线且对应设置两台接收机,即天线ⅰ和天线ⅱ,接收机ⅰ和接收机ⅱ;
s2、利用gps系统3r定位法确定天线ⅰ的位置坐标;
s3、将天线ⅰ的位置坐标设定为概略坐标,根据伪距线性误差方程对天线ⅰ的位置坐标进行修正;
s4、根据天线ⅰ与天线ⅱ之间的基线长度与载波波长的关系确定整周模糊度搜索范围,再进一步确定整周模糊度的值;
s5、将修正后的天线ⅰ的位置坐标作为天线ⅱ位置坐标的概略坐标,根据双差,最小二乘法定位出天线ⅱ的位置坐标;
s6、根据定位出的天线ⅰ的位置坐标和天线ⅱ位置坐标确定测姿载体的偏航角和俯仰角,即进行测姿载体的姿态角解算。
进一步的,所述步骤s2中利用gps系统3r定位法确定天线ⅰ的位置坐标,包括,通过空间中三颗卫星为原点以各自到接收机ⅰ的距离为半径得到三个球面的相交点,确定天线ⅰ位置坐标在空间中的位置;设天线ⅰ利用到达时间分别测定gps系统中的三颗卫星s1、s2和s3到天线ⅰ位置的分别距离为p1、p2和p3;天线ⅰ的接收机ⅰ利用导航电文计算得到s1、s2和s3三颗卫星的位置分别为(x1,y1,z1),(x2,y2,z2)和(x3,y3,z3);根据距离交汇法求解天线ⅰ三维坐标的观测方程为:
其中,(x,y,z)为天线ⅰ的坐标位置;根据上式解算出的天线ⅰ位置坐标为(x1,y1,z1)。
进一步的,所述步骤s3中,将天线ⅰ的位置坐标设定为概略坐标,根据伪距线性误差方程对天线ⅰ的位置坐标进行修正,包括,
将步骤s2解算出的天线ⅰ的位置坐标(x1,y1,z1)设定为概略坐标(x'1,y'1,z'1),设概略位置与真实位置之间的偏移记为(δx,δy,δz),可得伪距定位线性误差方程为:
其中,
解算该方程组得到修正后的天线ⅰ的位置坐标(x'1,y'1,z'1)。
进一步的,所述步骤s4中,根据天线ⅰ与天线ⅱ之间的基线长度与载波波长的关系确定整周模糊度搜索范围,再进一步确定整周模糊度的值,包括,设在时刻t,接收机ⅰ和接收机ⅱ同步观测到卫星j,得到两个载波观测方程;
式中,
再让接收机ⅰ和接收机ⅱ同时观测卫星k,将再次得到的单差方程与上面的公式进行作差,得到的双差方程如下公式:
式中,
由此得到:
由于天线ⅰ和天线ⅱ之间的基线长度小于载波波长的1/2,
对
进一步的,所述步骤s5中,将修正后的天线ⅰ的位置坐标作为天线ⅱ位置坐标的概略坐标,根据双差,最小二乘法定位出天线ⅱ的位置坐标,包括,设接收机ⅰ为基准站,天线ⅰ为主天线,接收机ⅱ为流动站,天线ⅱ为从天线,将修正后的天线ⅰ的位置坐标(x'1,y'1,z'1)作为天线ⅱ的概略坐标,设天线ⅱ的概略坐标为(x'2,y'2,z'2)其中(x'2,y'2,z'2)=(x'1,y'1,z'1),真实坐标为(x2,y2,z2),概略坐标与真实坐标的偏移记为(δx',δy',δz');将天线ⅱ与卫星j和k之间的距离
因为基线长度远远小于天线与卫星的距离,因此,选用天线ⅰ的坐标作为天线ⅱ的概略坐标,式中(x'2,y'2,z'2)为天线ⅱ的概略坐标,即为步骤s3中求出的天线ⅰ的修正后位置坐标(x'1,y'1,z'1),
令
整理得到:
所以得出
其中,
进一步的,所述步骤s6中,根据定位出的天线ⅰ的位置坐标和天线ⅱ位置坐标确定测姿载体的偏航角和俯仰角,即进行姿态解算,包括,设天线ⅰ为点a,天线ⅱ为点b,天线ⅱ与天线ⅱ组成向量
δx=x2-x1
δy=y2-y1
δz=z2-z1
由上述公式可以解算出偏航角
本发明双天线gps测姿方法的有益技术效果是先对天线ⅰ的位置坐标进行了修正,在一定程度上提高了定位精度;其次,在整周模糊度值的确定上,采用波长对整周模糊度的搜索区间进行限制,降低了整周模糊度的搜索范围和时间复杂度;在此基础上完成测姿载体的姿态角解算,从而提高了测姿的精确度和及时性。
附图说明
图1为测姿载体的三维姿态示意图;
图2为测姿载体的姿态角示意图。
下面结合附图对本发明双天线gps测姿方法作进一步的说明。
具体实施方式
显而易见地,下面描述中的附图仅仅是本发明的一些实施例,应当理解,此处所描述的优选实施例仅用于说明和解释本发明,并不用于限定本发明。
附图1为测姿载体的三维姿态示意图,附图2为测姿载体的姿态角示意图。由图可知,本发明双天线gps测姿方法,在测姿载体上呈几何分布设置两个gps信号接收天线且对应设置两台接收机,采用gps系统3r定位法确定天线ⅰ的位置坐标,将天线ⅰ的位置坐标定为概略坐标并对其进行修正,确定整周模糊度的值,然后,解算天线ⅱ的位置坐标,最后,根据定位出的天线ⅰ的位置坐标和天线ⅱ的位置坐标确定测姿载体的偏航角和俯仰角,即进行测姿载体的姿态角解算。
本发明双天线gps测姿方法,包括以下步骤:
s1、在测姿载体上呈几何分布设置两个gps信号接收天线且对应设置两台接收机,即天线ⅰ和天线ⅱ,接收机ⅰ和接收机ⅱ;
s2、利用gps系统3r定位法确定天线ⅰ的位置坐标;
s3、将天线ⅰ的位置坐标设定为概略坐标,根据伪距线性误差方程对天线ⅰ的位置坐标进行修正;
s4、根据天线ⅰ与天线ⅱ之间的基线长度与载波波长的关系确定整周模糊度搜索范围,再进一步确定整周模糊度的值;
s5、将修正后的天线ⅰ的位置坐标作为天线ⅱ位置坐标的概略坐标,根据双差,最小二乘法定位出天线ⅱ的位置坐标;
s6、根据定位出的天线ⅰ的位置坐标和天线ⅱ位置坐标确定测姿载体的偏航角和俯仰角,即进行测姿载体的姿态角解算。
作为优选,所述步骤s2中利用gps系统3r定位法确定天线ⅰ的位置坐标,包括,通过空间中三颗卫星为原点以各自到接收机ⅰ的距离为半径得到三个球面的相交点,确定天线ⅰ位置坐标在空间中的位置;设天线ⅰ利用到达时间分别测定gps系统中的三颗卫星s1、s2和s3到天线ⅰ位置的分别距离为p1、p2和p3;天线ⅰ的接收机ⅰ利用导航电文计算得到s1、s2和s3三颗卫星的位置分别为(x1,y1,z1),(x2,y2,z2)和(x3,y3,z3);根据距离交汇法求解天线ⅰ三维坐标的观测方程为:
其中,(x,y,z)为天线ⅰ的坐标位置;根据上式解算出的天线ⅰ位置坐标为(x1,y1,z1)。
作为优选,所述步骤s3中,将天线ⅰ的位置坐标设定为概略坐标,根据伪距线性误差方程对天线ⅰ的位置坐标进行修正,包括,
将步骤s2解算出的天线ⅰ的位置坐标(x1,y1,z1)设定为概略坐标(x'1,y'1,z'1),设概略位置与真实位置之间的偏移记为(δx,δy,δz),可得伪距定位线性误差方程为:
其中,
解算该方程组得到修正后的天线ⅰ的位置坐标(x'1,y'1,z'1)。
作为优选,所述步骤s4中,根据天线ⅰ与天线ⅱ之间的基线长度与载波波长的关系确定整周模糊度搜索范围,再进一步确定整周模糊度的值,包括,设在时刻t,接收机ⅰ和接收机ⅱ同步观测到卫星j,得到两个载波观测方程;
式中,
再让接收机ⅰ和接收机ⅱ同时观测卫星k,将再次得到的单差方程与上面的公式进行作差,得到的双差方程如下公式:
式中,
由此得到:
由于天线ⅰ和天线ⅱ之间的基线长度小于载波波长的1/2,
对
作为优选,所述步骤s5中,将修正后的天线ⅰ的位置坐标作为天线ⅱ位置坐标的概略坐标,根据双差,最小二乘法定位出天线ⅱ的位置坐标,包括,设接收机ⅰ为基准站,天线ⅰ为主天线,接收机ⅱ为流动站,天线ⅱ为从天线,将修正后的天线ⅰ的位置坐标(x'1,y'1,z'1)作为天线ⅱ的概略坐标,设天线ⅱ的概略坐标为(x'2,y'2,z'2),其中(x'2,y'2,z'2)=(x'1,y'1,z'1),真实坐标为(x2,y2,z2),概略坐标与真实坐标的偏移记为(δx',δy',δz');将天线ⅱ与卫星j和k之间的距离
因为基线长度远远小于天线与卫星的距离,因此,选用天线ⅰ的坐标作为天线ⅱ的概略坐标,式中(x'2,y'2,z'2)为天线ⅱ的概略坐标,即为步骤s3中求出的天线ⅰ的修正后位置坐标(x'1,y'1,z'1),
令
整理得到:
所以得出
其中,
作为优选,所述步骤s6中,根据定位出的天线ⅰ的位置坐标和天线ⅱ位置坐标确定测姿载体的偏航角和俯仰角,即进行姿态解算,包括,设天线ⅰ为点a,天线ⅱ为点b,天线ⅱ与天线ⅱ组成向量
δx=x2-x1
δy=y2-y1
δz=z2-z1
由上述公式可以解算出偏航角
显然,本发明双天线gps测姿方法的有益技术效果是先对天线ⅰ的位置坐标进行了修正,在一定程度上提高了定位精度;其次,在整周模糊度值的确定上,采用波长对整周模糊度的搜索区间进行限制,降低了整周模糊度的搜索范围和时间复杂度;在此基础上完成测姿载体的姿态角解算,从而提高了测姿的精确度和及时性。