CN105403903B - 一种载体姿态测量方法及设备 - Google Patents
一种载体姿态测量方法及设备 Download PDFInfo
- Publication number
- CN105403903B CN105403903B CN201510732259.XA CN201510732259A CN105403903B CN 105403903 B CN105403903 B CN 105403903B CN 201510732259 A CN201510732259 A CN 201510732259A CN 105403903 B CN105403903 B CN 105403903B
- Authority
- CN
- China
- Prior art keywords
- fuzziness
- observation
- carrier
- satellite
- error
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/53—Determining attitude
- G01S19/54—Determining attitude using carrier phase measurements; using long or short baseline interferometry
- G01S19/55—Carrier phase ambiguity resolution; Floating ambiguity; LAMBDA [Least-squares AMBiguity Decorrelation Adjustment] method
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种载体姿态测量方法及设备,方法包括下述步骤:获得当前历元的卫星观测值;利用伪距差分技术对观测值进行粗差检查,如果没通过检查,退出当前历元,将下一历元设为当前历元,返回上一步,如果通过检查,建立基线约束定向的观测模型和站星双差随机模型;利用最小二乘法求解基线约束定向的观测模型和站星双差随机模型的模糊度浮点解及其协因数阵;使用LAMBDA方法、二次约束和累积中误差算法,获得正确的模糊度;然后解算基线向量,获得载体姿态。本发明的载体姿态测量方法,采用LAMBDA、二次约束、累积中误差算法,提高模糊度解算成功率,提高获得载体姿态的成功率;本发明的载体姿态测量设备实现实时动态高精度定向。
Description
技术领域
本发明属于姿态测量技术领域,具体地说,是涉及一种载体姿态测量方法及载体姿态测量设备。
背景技术
姿态测量是航空、航天、航海以及陆地导航系统的关键技术之一。北斗卫星导航系统(BeiDou Navigation Satellite System,简称为BDS,早期也称作COMPASS)是我国自行研制发展、独立运行的全球卫星定位与通信系统,是继美国(GPS)、俄国(GLONASS)之后第三个已经投入运行的成熟的卫星导航系统。
随着北斗二代导航卫星系统的快速发展,利用北斗卫星组成姿态测量系统来测量载体的姿态迅速成为一个重要的研究方向。
目前使用的姿态测量方法主要是组建双差观测方程,使用LAMBDA方法解算模糊度,从而求解基线向量,进而获得载体姿态,但是使用LAMBDA方法解算模糊度,模糊度解算成功率较低,从而影响了载体姿态的获得,不利于姿态测量技术的发展。
发明内容
本发明提供了一种载体姿态测量方法,提高了模糊度解算成功率。
为解决上述技术问题,本发明采用下述技术方案予以实现:
一种载体姿态测量方法,包括下述步骤:
(1)获得当前历元下多个观测卫星的观测值;
(2)利用伪距差分技术对观测值进行粗差检查;
如果没通过检查,退出当前历元,将下一历元设为当前历元,返回(1);
如果通过检查,进入(3);
(3)建立基线约束定向的观测模型和站星双差随机模型;
(4)利用最小二乘法求解基线约束定向的观测模型和站星双差随机模型的模糊度浮点解及其协因数阵;
(5)使用LAMBDA方法进行模糊度搜索,获得多组模糊度备选值;
(6)对每一组模糊度备选值进行二次约束;
(7)对所有通过二次约束的模糊度备选值使用累积中误差算法,获得正确的模糊度;
(8)使用正确的模糊度解算基线向量;
(9)通过基线向量获得载体姿态。
进一步的,所述步骤(2)具体包括下述步骤:
(21)建立站间伪距单差模型;
(22)利用最小二乘法求解中误差,并判断中误差是否大于阈值;
如果中误差大于阈值,说明存在粗差,进入(23);
如果中误差不大于阈值,说明不存在粗差,进入(24);
(23)判断当前历元的有效观测卫星数量是否大于1;
如果大于1,求解观测卫星的残差,找出残差最大的观测卫星,将残差最大的观测卫星设为当前历元的无效观测卫星,并丢弃该观测卫星的观测值,返回(22);
如果不大于1,进入(25);
(24)判断当前历元的有效观测卫星的数量是否小于预设卫星数;
如果小于,进入(25);
如果不小于,进入(26);
(25)失败,说明没通过检查,退出;
(26)成功,说明通过检查,退出。
又进一步的,所述建立基线约束定向的观测模型包括下述步骤:
(31)建立载波相位站星双差观测方程;
(32)建立伪距双差观测方程;
(33)建立基线约束虚拟观测方程;
(34)综合上述三个观测方程,获得基线约束定向的观测模型。
更进一步的,在所述步骤(4)之前还运行模糊度桥接步骤,如果桥接成功,则直接进入步骤(8),如果桥接失败,则进入步骤(4)。
优选的,所述模糊度桥接具体包括下述步骤:
(a1)判断是否成功获取历史历元的正确模糊度;
如果没有成功获取,桥接失败,退出;
如果成功获取,进入(a2);
(a2)判断参考卫星是否发生变化;
如果发生变化,进入(a3);
如果没发生变化,进入(a4);
(a3)将获取的模糊度按照当前历元的参考卫星进行参数变换;
(a4)获取共视卫星;
(a5)判断共视卫星的数量是否小于预设卫星数;
如果小于,桥接失败,退出;
如果不小于,进入(a6);
(a6)判断当前历元与获取的模糊度对应的历元是否有相同的共视卫星;
如果没有,桥接失败,退出;
如果有,进入(a7);
(a7)判断当前历元与获取的模糊度对应的历元是否有不相同的共视卫星;
如果有,计算出新的共视卫星的模糊度,然后输出当前历元的正确模糊度,桥接成功,退出;
如果没有,输出当前历元的正确模糊度,桥接成功,退出。
进一步的,在使用LAMBDA方法进行模糊度搜索中,通过扩大搜索体积和备选值数量的方法,获得多组模糊度备选值。
又进一步的,所述对每一组模糊度备选值进行二次约束是指对每一组模糊度备选值进行基线长度约束、姿态约束和中误差约束。
再进一步的,对所有通过二次约束的模糊度备选值使用累积中误差算法获得正确的模糊度是指对所有通过二次约束的模糊度备选值按照累积中误差进行排序,找出累积中误差最小的模糊度备选值,获得该模糊度备选值的累积成功次数,如果该累积成功次数大于设定次数,则该模糊度备选值为正确的模糊度备选值。
优选的,所述预设卫星数为5。
本发明还提出了一种载体姿态测量设备,包括数字信号处理器、两个北斗B1天线、两个北斗B1单频接收机,第一北斗B1天线安装在载体一端,第二北斗B1天线安装在载体另一端;第一北斗B1单频接收机接收所述第一北斗B1天线发送的数据,并传输至所述数字信号处理器;第二北斗B1单频接收机接收所述第二北斗B1天线发送的数据,并传输至所述数字信号处理器;所述数字信号处理器运行权利要求1至8中任一项所述的载体姿态测量方法,输出载体姿态。
与现有技术相比,本发明的优点和积极效果是:本发明的载体姿态测量方法,通过采用LAMBDA、二次约束以及累积中误差算法,提高了模糊度解算成功率,进而提高了获得载体姿态的成功率,促进了姿态测量技术的发展;本发明的载体姿态测量设备实现了实时动态高精度定向,且成本低、集成度高。
结合附图阅读本发明的具体实施方式后,本发明的其他特点和优点将变得更加清楚。
附图说明
图1是本发明提出的载体姿态测量方法的一个实施例的流程图;
图2是本发明提出的载体姿态测量方法的又一个实施例的流程图;
图3是图2中利用伪距单差技术对观测值进行粗差检查的一个实施例的流程图;
图4是图2中模糊度桥接的一个实施例的流程图;
图5是本发明提出的载体姿态测量方法的再一个实施例的流程图;
图6是本发明提出的载体姿态测量设备的一个实施例的结构框图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下将结合附图和实施例,对本发明作进一步详细说明。
在被测载体的两端安装有第一天线和第二天线,第一天线接收的数据传输给第一接收机,第二天线接收的数据传输给第二接收机,第一接收机和第二接收机将接收到的数据传输给数据处理器,数据处理器进行数据处理,计算出载体的姿态。第一天线和第一接收机组成第一测站(测站1),第二天线和第二接收机组成第二测站(测站2)。
实施例一、参见图1所示,本实施例的载体姿态测量方法包括下述步骤:
步骤101:流程开始。
步骤102:获得当前历元下多个观测卫星的观测值。
第一测站和第二测站分别获得当前历元下多个观测卫星的载波相位、伪距、多普勒、星历数据等观测值。
步骤103:利用伪距单差技术对观测值进行粗差检查。
步骤104:判断是否通过粗差检查。
如果没有通过检查,进入步骤105;
如果通过检查,进入步骤106。
步骤105:退出当前历元,将下一历元设为当前历元,返回步骤102;
步骤106:建立基线约束定向的观测模型和站星双差随机模型。
步骤107:求解模糊度浮点解及其协因数阵。
使用最小二乘法求解基线约束定向的观测模型和站星双差随机模型的模糊度浮点解及其协因数阵。
步骤108:使用LAMBDA方法进行模糊度搜索,获得多组模糊度备选值。
在使用LAMBDA方法进行模糊度搜索中,通过扩大搜索体积和备选值数量的方法,获得多组模糊度备选值。在本实施例中,搜索体积为5,备选值数量为100,可以得到100组模糊度备选值。
步骤109:对每一组模糊度备选值进行二次约束。
步骤110:对所有通过二次约束的模糊度备选值使用累积中误差算法,获得正确的模糊度。
累积中误差算法的核心是用时间换取成功率,随着时间的增加,卫星位置发生变化,错误的模糊度在新的方程中误差会越来越大,甚至无法通过二次约束,只有唯一正确的模糊度才能保证在所有的时间内都能出现,且其基线中误差相对较小,其累积的中误差最小。
对所有通过二次约束的模糊度备选值按照累积中误差进行排序,找出累积中误差最小的模糊度备选值,获得该模糊度备选值的累积成功次数,如果该累积成功次数大于设定次数,则该模糊度备选值为正确的模糊度备选值。
通过采用LAMBDA、二次约束以及累积中误差算法,提高了模糊度解算成功率,进而提高了获得载体姿态的成功率,促进了姿态测量技术的发展。
步骤111:使用正确的模糊度解算基线向量。
将正确的模糊度代入观测模型,求得基线向量。
步骤112:通过基线向量获得载体姿态。
获得基线向量后,可求解两个天线之间的航向和纵向角度等,即可获得当前历元下载体的姿态。
步骤113:流程结束。
实施例二、参见图2所示,本实施例的载体姿态测量方法包括下述步骤:
步骤201:流程开始。
步骤202:获得当前历元下多个观测卫星的观测值。
第一测站和第二测站分别获得当前历元下多个观测卫星的载波相位、伪距、多普勒、星历数据等观测值。
步骤203:利用伪距单差技术对观测值进行粗差检查。
参见图3所示,进行粗差检查的具体步骤是:
步骤203-1:建立站间伪距单差模型。
ρ12=l2dΔX12+m2dΔY12+n2dΔZ12+M12+ε12
式中,ρ12为测站1和测站2之间的伪距观测值之差,l2,m2,n2分别为伪距单差方程线性化后的系数,dΔX12,dΔY12,dΔZ12为测站坐标差形式的改正参数,M12为伪距单差方程线性化后的常数项,ε12为随机误差。
步骤203-2:利用最小二乘法求解中误差。
步骤203-3:判断中误差是否大于阈值。
如果中误差大于阈值,说明存在粗差,进入步骤203-4;
如果中误差不大于阈值,说明不存在粗差,进入步骤203-6。
在本实施例中,阈值选择为5。
步骤203-4:判断当前历元的有效观测卫星数量是否大于1;
如果大于1,进入步骤203-5;
如果不大于1,进入步骤203-7;
步骤203-5:求解观测卫星的残差,找出残差最大的观测卫星,将残差最大的观测卫星设为当前历元的无效观测卫星,并丢弃该观测卫星的观测值,返回步骤203-2;
步骤203-6:判断当前历元的有效观测卫星的数量是否小于预设卫星数;
如果小于,进入步骤203-7;
如果不小于,进入步骤203-8。
在本实施例中,预设卫星数为5。
步骤203-7:失败,说明没通过检查,退出。
步骤203-8:成功,说明通过检查,退出。
步骤204:判断是否通过粗差检查。
如果没有通过检查,进入步骤205;
如果通过检查,进入步骤206。
步骤205:退出当前历元,将下一历元设为当前历元,返回步骤202;
步骤206:建立基线约束定向的观测模型和站星双差随机模型。
(1)线性化后的北斗载波相位模型为:
式中,为接收机i对卫星j的载波观测值,λ为波长,分别为线性化后的系数,dXi,dYi,dZi为测站坐标未知参数,()0为线性化后的常数项,为整周模糊度,c为光速,δti为接收机钟差,δtj为卫星钟差,分别为对流层误差和电离层误差修正值,为随机误差。i=1,2,接收机1即为第一接收机,接收机2即为第二接收机。j=1,2,3,...N,N为观测卫星的数量。
第一测站和第二测站同时观测卫星j和k,j和k是多个观测卫星中的两个,当然,其中一个是参考卫星。
建立载波相位站星双差观测方程如下:
λ为波长,为载波相位观测值的差分,dΔX12,dΔY12,dΔZ12同上,为测站坐标差形式的改正参数,为分别为卫星j和k的线性化方程系数之差,为双差模糊度,为双差方程移位后剩下的常数项,ε1为随机误差。
(2)建立伪距双差观测方程
式中,指的是卫星j和卫星k的伪距单差观测值之差,dΔX12,dΔY12,dΔZ12和同上,为伪距双差观测方程线性化后的常数项,ε2为随机误差。
(3)建立基线约束虚拟观测方程
l2=(ΔX12+dΔX12)2+(ΔY12+dΔY12)2+(ΔZ12+dΔZ12)2+ε3
式中,l为基线长度,ΔX12,ΔY12,ΔZ12分别为三个方向的基线向量的近似值,dΔX12,dΔY12,dΔZ12为基线改正数向量,也是双差定位中的未知数向量,ε3为随机误差。
令该数值由初始基线精度确定。展开后,可以得到:
l0基线长度的初始估计值。
(4)综合上述三个观测方程,即可获得基线约束定向的观测模型。
其中,B=[2ΔX12 2ΔY12 2ΔZ12],N为观测卫星的数量,dΔX=[dΔX12 dΔY12 dΔZ12]T,v=[ε1 ε2ε3]T。
(5)建立站星双差随机模型。
在随机模型中,一般认为,北斗伪距观测与相位观测值不相关,不同北斗卫星之间的观测是相互独立的,且基线长度虚拟观测值与北斗观测信息不相关,于是,其随机模型为:
其中,
Ql=m
k=0.0001,l=4,m=0.0001。
步骤207:进行模糊度桥接。
参见图4所示,模糊度桥接具体包括下述步骤:
步骤207-1:尝试获取历史历元的正确模糊度。
在本实施例中,历史历元是指上个历元。
步骤207-2:判断是否成功获取历史历元的正确模糊度。
如果没有成功获取,进入步骤207-11;
如果成功获取,进入步骤207-3。
步骤207-3:判断参考卫星是否发生变化。
如果发生变化,进入步骤207-4;
如果没有发生变化,进入步骤207-5。
步骤207-4:参数变换。
将获取的模糊度按照当前历元的参考卫星进行参数变换,获得变换后的模糊度。
步骤207-5:获取共视卫星。
步骤207-6:判断共视卫星的数量是否小于预设卫星数。
如果小于,进入步骤207-11;
如果不小于,进入步骤207-7。
在本实施例中,预设卫星数为5。
步骤207-7:判断当前历元与获取的模糊度对应的历元是否有相同的共视卫星。
在本实施例中,获取的模糊度对应的历元是指上一历元。
如果没有,进入步骤207-11;
如果有,进入步骤207-8。
步骤207-8:判断当前历元与获取的模糊度对应的历元是否有不相同的共视卫星。
在本实施例中,获取的模糊度对应的历元是指上一历元。
如果有,进入步骤207-9;
如果没有,进入步骤207-10。
步骤207-9:计算出新的共视卫星的模糊度。
新的共视卫星,即步骤207-8中判断出的不相同的共视卫星。
步骤207-10:输出当前历元的正确模糊度,桥接成功,退出。
步骤207-11:桥接失败,退出。
步骤208:判断模糊度桥接是否成功。
如果桥接成功,进入步骤213。
如果桥接失败,进入步骤209。
步骤209:求解模糊度浮点解及其协因数阵。
使用最小二乘法求解基线约束定向的观测模型和站星双差随机模型的模糊度浮点解及其协因数阵。
步骤210:使用LAMBDA方法进行模糊度搜索,获得多组模糊度备选值。
在使用LAMBDA方法进行模糊度搜索中,通过扩大搜索体积和备选值数量的方法,获得多组模糊度备选值。在本实施例中,搜索体积为5,备选值数量为100,可以得到100组模糊度备选值。
步骤211:对每一组模糊度备选值进行二次约束。
在本实施例中,二次约束包括基线长度约束、姿态约束和中误差约束。
也就是说,对每一组模糊度备选值进行基线长度约束、姿态约束和中误差约束。
步骤212:对所有通过二次约束的模糊度备选值使用累积中误差算法,获得正确的模糊度。
累积中误差算法的核心是用时间换取成功率,随着时间的增加,卫星位置发生变化,错误的模糊度在新的方程中误差会越来越大,甚至无法通过二次约束,只有唯一正确的模糊度才能保证在所有的时间内都能出现,且其基线中误差相对较小,其累积的中误差最小。
对所有通过二次约束的模糊度备选值按照累积中误差进行排序,找出累积中误差最小的模糊度备选值,获得该模糊度备选值的累积成功次数,如果该累积成功次数大于设定次数,则该模糊度备选值为正确的模糊度备选值。
通过采用LAMBDA、二次约束以及累积中误差算法,提高了模糊度解算成功率,进而提高了获得载体姿态的成功率,促进了姿态测量技术的发展。
步骤213:使用正确的模糊度解算基线向量。
将正确的模糊度代入观测模型,求得基线向量。
步骤214:通过基线向量获得载体姿态。
获得基线向量后,可求解两个天线之间的航向和纵向角度等,即可获得当前历元下载体的姿态。
步骤215:流程结束。
在本实施例中,进行模糊度解算之前,先进行模糊度桥接,如果模糊度桥接成功,可以直接得到当前历元的正确模糊度,而不需要进行模糊度的解算,从而提高了获得基线向量的效率,进而提高了获得载体姿态的效率,促进了姿态测量技术的发展;而且由于不需要进行模糊度的解算,极大地节省了处理器资源。
实施例三、参见图5所示,本实施例的载体姿态测量方法包括下述步骤:
步骤301:流程开始。
步骤302:获得当前历元下多个观测卫星的观测值。
第一测站和第二测站分别接收当前历元下多个观测卫星的载波相位、伪距、多普勒、星历数据等观测值。
步骤303:利用伪距单差技术对观测值进行粗差检查。
进行粗差检查的具体步骤可参照实施例二,此处不再赘述。
步骤304:判断是否通过粗差检查。
如果没有通过检查,进入步骤305;
如果通过检查,进入步骤306。
步骤305:退出当前历元,将下一历元设为当前历元,返回步骤302;
步骤306:建立基线约束定向的观测模型和站星双差随机模型。
具体步骤可参照实施例二,此处不再赘述。
步骤307:进行模糊度桥接。
具体步骤可参照实施例二,此处不再赘述。
步骤308:判断模糊度桥接是否成功。
如果桥接成功,进入步骤314。
如果桥接失败,进入步骤309。
步骤309:求解模糊度浮点解及其协因数阵。
使用最小二乘法求解基线约束定向的观测模型和站星双差随机模型的模糊度浮点解及其协因数阵。
步骤310:使用LAMBDA方法进行模糊度搜索,获得多组模糊度备选值。
在使用LAMBDA方法进行模糊度搜索中,通过扩大搜索体积和备选值数量的方法,获得多组模糊度备选值。在本实施例中,在使用LAMBDA方法进行模糊度搜索时,搜索体积为5,备选值数量为100,可以得到100组模糊度备选值。
上个历元的第二组合中的模糊度备选值以及当前历元的上述100组模糊度备选值构成当前历元的第一组合。将第一组合中的模糊度备选值代入观测模型求解基线向量。
下面,依次对第一组合中的每一组模糊度备选值进行二次约束。
在本实施例中,二次约束包括基线长度约束、姿态约束和中误差约束。
步骤311-1:判断当前组模糊度备选值是否满足基线长度约束条件。
根据事先测定的基线长度精度确定约束条件。在本实施例中,基线长度约束条件为基线长度<=0.01m。
如果不满足,进入步骤311-2;
如果满足,进入步骤311-5。
步骤311-2:在第一组合中剔除该组模糊度备选值,进入311-3。
步骤311-3:判断该组模糊度备选值是否是第一组合的最后一组。
如果不是,进入步骤311-4。
如果是,进入步骤312。
步骤311-4:将第一组合的下一组模糊度备选值设为当前组模糊度备选值,返回311-1。
步骤311-5:判断当前组模糊度备选值是否满足姿态约束条件。
使用上个历元的姿态约束条件或已知的姿态约束条件。在本实施例中,姿态约束条件为载体的纵向角小于20°。
如果不满足,进入步骤311-2;
如果满足,进入步骤311-6。
步骤311-6:计算基线中误差。
步骤311-7:判断当前组模糊度备选值是否满足中误差约束条件。
在本实施例中,中误差约束条件为中误差<=0.05。
如果不满足,进入步骤311-2;
如果满足,进入步骤311-8。
步骤311-8:将该组模糊度备选值加入第二组合,并计算该组模糊度备选值的累积中误差以及累积成功次数,进入311-3。
步骤312:对第二组合中的模糊度备选值按照累积中误差进行排序,找出累积中误差最小的模糊度备选值,获得该模糊度备选值的累积成功次数。
步骤313:判断该累积成功次数是否大于设定次数。
如果不大于,则失败,没有找到正确的模糊度,返回步骤305;
如果大于,则成功,该组模糊度备选值为正确的模糊度,进入步骤314;
在本实施例中,设定次数为60。
步骤314:使用正确的模糊度解算基线向量。
将正确的模糊度代入观测模型,求得基线向量。
步骤315:通过基线向量获得载体姿态。
获得基线向量后,可求解两个天线之间的航向和纵向角度等,即可获得当前历元下载体的姿态。
步骤316:流程结束。
在本实施例中,通过采用LAMBDA、二次约束以及累积中误差算法,提高了模糊度解算成功率,进而提高了获得载体姿态的成功率,促进了姿态测量技术的发展。
在本实施例中,进行模糊度解算之前,先进行模糊度桥接,如果模糊度桥接成功,可以直接得到当前历元的正确模糊度,而不需要进行模糊度的解算,提高了模糊度的搜索效率,从而提高了获得基线向量的效率,进而提高了获得载体姿态的效率,促进了姿态测量技术的发展;而且由于不需要进行模糊度的解算,极大地节省了处理器资源。
实施例四、本实施例提出了一种姿态测量设备,主要包括数字信号处理器DSP、两个北斗B1天线、两个北斗B1单频接收机,第一北斗B1天线安装在载体一端,第二北斗B1天线安装在载体另一端;第一北斗B1单频接收机接收第一北斗B1天线发送的数据,并传输至数字信号处理器DSP;第二北斗B1单频接收机接收第二北斗B1天线发送的数据,并传输至数字信号处理器DSP;数字信号处理器DSP运行实施例一至实施例三提出的载体姿态测量方法,输出载体姿态,参见图6所示。
在本实施例中,数字信号处理器DSP采用低功耗、低尺寸、低成本的DSP嵌入式处理器OMAPL137,该处理器支持高速浮点运算,可以满足姿态测量中模糊度实时动态搜索等算法中的大运算量要求。
所述姿态测量设备还包括外部SDRAM,DSP在运算过程中将大量数据存储在SDRAM,减小了内部RAM的存储压力。
本实施例的姿态测量设备实现了实时动态高精度定向,且具有成本低、集成度高、体积小、安装灵活等优点,可广泛应用于航天测量、车船姿态测量、精密着陆等方面。
以上实施例仅用以说明本发明的技术方案,而非对其进行限制;尽管参照前述实施例对本发明进行了详细的说明,对于本领域的普通技术人员来说,依然可以对前述实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或替换,并不使相应技术方案的本质脱离本发明所要求保护的技术方案的精神和范围。
Claims (8)
1.一种载体姿态测量方法,其特征在于:包括下述步骤:
(1)获得当前历元下多个观测卫星的观测值;
(2)利用伪距差分技术对观测值进行粗差检查;具体包括下述步骤:
(21)建立站间伪距单差模型;
(22)利用最小二乘法求解中误差,并判断中误差是否大于阈值;
如果中误差大于阈值,说明存在粗差,进入(23);
如果中误差不大于阈值,说明不存在粗差,进入(24);
(23)判断当前历元的有效观测卫星数量是否大于1;
如果大于1,求解观测卫星的残差,找出残差最大的观测卫星,将残差最大的观测卫星设为当前历元的无效观测卫星,并丢弃该观测卫星的观测值,返回(22);
如果不大于1,进入(25);
(24)判断当前历元的有效观测卫星的数量是否小于预设卫星数;
如果小于,进入(25);
如果不小于,进入(26);
(25)失败,说明没通过检查,退出;
(26)成功,说明通过检查,退出;
判断是否通过粗差检查;
如果没通过检查,退出当前历元,将下一历元设为当前历元,返回(1);
如果通过检查,进入(3);
(3)建立基线约束定向的观测模型和站星双差随机模型;
(4)利用最小二乘法求解基线约束定向的观测模型和站星双差随机模型的模糊度浮点解及其协因数阵;
(5)使用LAMBDA方法进行模糊度搜索,获得多组模糊度备选值;
(6)对每一组模糊度备选值进行二次约束;
(7)对所有通过二次约束的模糊度备选值使用累积中误差算法,获得正确的模糊度:
对所有通过二次约束的模糊度备选值按照累积中误差进行排序,找出累积中误差最小的模糊度备选值,获得该模糊度备选值的累积成功次数,如果该累积成功次数大于设定次数,则该模糊度备选值为正确的模糊度;
(8)使用正确的模糊度解算基线向量;
(9)通过基线向量获得载体姿态。
2.根据权利要求1所述的载体姿态测量方法,其特征在于:所述建立基线约束定向的观测模型包括下述步骤:
(31)建立载波相位站星双差观测方程;
(32)建立伪距双差观测方程;
(33)建立基线约束虚拟观测方程;
(34)综合上述三个观测方程,获得基线约束定向的观测模型。
3.根据权利要求1所述的载体姿态测量方法,其特征在于:在所述步骤(4)之前还运行模糊度桥接步骤,如果桥接成功,则直接进入步骤(8),如果桥接失败,则进入步骤(4)。
4.根据权利要求3所述的载体姿态测量方法,其特征在于:所述模糊度桥接具体包括下述步骤:
(a1)判断是否成功获取历史历元的正确模糊度;
如果没有成功获取,桥接失败,退出;
如果成功获取,进入(a2);
(a2)判断参考卫星是否发生变化;
如果发生变化,进入(a3);
如果没发生变化,进入(a4);
(a3)将获取的模糊度按照当前历元的参考卫星进行参数变换;
(a4)获取共视卫星;
(a5)判断共视卫星的数量是否小于预设卫星数;
如果小于,桥接失败,退出;
如果不小于,进入(a6);
(a6)判断当前历元与获取的模糊度对应的历元是否有相同的共视卫星;
如果没有,桥接失败,退出;
如果有,进入(a7);
(a7)判断当前历元与获取的模糊度对应的历元是否有不相同的共视卫星;
如果有,计算出新的共视卫星的模糊度,然后输出当前历元的正确模糊度,桥接成功,退出;
如果没有,输出当前历元的正确模糊度,桥接成功,退出。
5.根据权利要求1所述的载体姿态测量方法,其特征在于:在使用LAMBDA方法进行模糊度搜索中,通过扩大搜索体积和备选值数量的方法,获得多组模糊度备选值。
6.根据权利要求5所述的载体姿态测量方法,其特征在于:所述对每一组模糊度备选值进行二次约束是指对每一组模糊度备选值进行基线长度约束、姿态约束和中误差约束。
7.根据权利要求1所述的载体姿态测量方法,其特征在于:所述预设卫星数为5。
8.一种载体姿态测量设备,其特征在于:包括数字信号处理器、两个北斗B1天线、两个北斗B1单频接收机,第一北斗B1天线安装在载体一端,第二北斗B1天线安装在载体另一端;第一北斗B1单频接收机接收所述第一北斗B1天线发送的数据,并传输至所述数字信号处理器;第二北斗B1单频接收机接收所述第二北斗B1天线发送的数据,并传输至所述数字信号处理器;所述数字信号处理器运行权利要求1至6中任一项所述的载体姿态测量方法,输出载体姿态。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510732259.XA CN105403903B (zh) | 2015-11-02 | 2015-11-02 | 一种载体姿态测量方法及设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510732259.XA CN105403903B (zh) | 2015-11-02 | 2015-11-02 | 一种载体姿态测量方法及设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105403903A CN105403903A (zh) | 2016-03-16 |
CN105403903B true CN105403903B (zh) | 2018-01-23 |
Family
ID=55469501
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510732259.XA Active CN105403903B (zh) | 2015-11-02 | 2015-11-02 | 一种载体姿态测量方法及设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105403903B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107632316A (zh) * | 2017-09-11 | 2018-01-26 | 北京数码汇博科技有限公司 | 一种利用卫星定位测量车身姿态的方法及结构 |
CN109613585A (zh) * | 2018-12-14 | 2019-04-12 | 中国科学院国家授时中心 | 一种对基站天线超短基线gnss双天线实时测向的方法 |
CN109613565A (zh) * | 2019-01-14 | 2019-04-12 | 中国人民解放军战略支援部队信息工程大学 | 基于多星座gnss的电离层层析方法及系统 |
CN113253322B (zh) * | 2021-04-12 | 2023-09-01 | 中国人民解放军军事科学院国防科技创新研究院 | 基于双天线的移动载体相对位置解算方法 |
CN113551673B (zh) * | 2021-07-23 | 2024-02-13 | 山东泉清通信有限责任公司 | 一种惯导系统及其设计方法、惯导方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101825717B (zh) * | 2010-04-16 | 2012-10-24 | 北京航空航天大学 | 一种基于载波平滑码伪距技术的动态定姿方法 |
CN103163542B (zh) * | 2013-03-26 | 2014-12-10 | 东南大学 | 一种gnss基线解算中基于观测常量的粗差探测方法 |
CN104483688A (zh) * | 2014-11-03 | 2015-04-01 | 中国人民解放军63961部队 | 基于北斗卫星导航系统的高精度基线解算方法 |
-
2015
- 2015-11-02 CN CN201510732259.XA patent/CN105403903B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN105403903A (zh) | 2016-03-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106168672B (zh) | 一种gnss多模单频rtk周跳探测方法及装置 | |
CN105403903B (zh) | 一种载体姿态测量方法及设备 | |
CN107102346B (zh) | 一种基于北斗系统的多天线测姿方法 | |
CN114518586B (zh) | 一种基于球谐展开的gnss精密单点定位方法 | |
CN105158783B (zh) | 一种实时动态差分定位方法及其设备 | |
CN108981559B (zh) | 基于北斗地基增强系统的实时变形监测方法及系统 | |
CN107272039B (zh) | 一种基于双天线gps的定位测姿方法 | |
CN107193029B (zh) | 北斗三频信号的网络rtk基准站间模糊度快速确定方法 | |
CN102565834B (zh) | 一种单频gps测向系统及其测向定位方法 | |
CN110531392A (zh) | 一种基于ppp算法的高精度定位方法和系统 | |
CN111559372B (zh) | 一种泊车系统的车位检测性能的评价方法和电子设备 | |
CN106526629B (zh) | 卫星导航设备及其定向方法及装置 | |
CN105842719B (zh) | 一种顾及对流层影响的cors基准站网基线模糊度解算方法 | |
CN105158778B (zh) | 多系统联合实施载波相位差分故障卫星剔除方法及其系统 | |
CN108196281A (zh) | 一种基于位置域曲线约束的单频动态周跳探测与修复方法 | |
CN104215977B (zh) | 一种基于卫星导航系统的精度评估方法及系统 | |
CN1967282B (zh) | 全站仪与gps单频实时动态组合测量方法及其系统 | |
JP2009264844A (ja) | 相対位置検知装置、及び相対位置検知システム | |
CN109738902B (zh) | 一种基于同步信标模式的水下高速目标高精度自主声学导航方法 | |
CN106547002A (zh) | 一种提高rtk定位精度的选星方法 | |
CN107462905A (zh) | 一种北斗gps双模定位方法、定位终端及定位系统 | |
MX2011005844A (es) | Un sistema y metodo para determinacion de posicion de referencia. | |
CN107607032A (zh) | 一种gnss形变监测系统 | |
CN116819580B (zh) | 一种惯性辅助的双天线gnss海上船舶姿态确定方法 | |
CN106772483A (zh) | 一种基于cors系统的数据后处理方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |