CN107656300A - 基于北斗/gps双模软件接收机的卫星/惯性超紧组合系统及方法 - Google Patents
基于北斗/gps双模软件接收机的卫星/惯性超紧组合系统及方法 Download PDFInfo
- Publication number
- CN107656300A CN107656300A CN201710695052.9A CN201710695052A CN107656300A CN 107656300 A CN107656300 A CN 107656300A CN 201710695052 A CN201710695052 A CN 201710695052A CN 107656300 A CN107656300 A CN 107656300A
- Authority
- CN
- China
- Prior art keywords
- mtd
- mrow
- msub
- mover
- mfrac
- 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.)
- Granted
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/42—Determining position
- G01S19/45—Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement
- G01S19/47—Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement the supplementary measurement being an inertial measurement, e.g. tightly coupled inertial
-
- 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/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/33—Multimode operation in different systems which transmit time stamped messages, e.g. GPS/GLONASS
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种基于北斗/GPS双模软件接收机的卫星/惯性超紧组合系统及方法,系统包括有源卫星天线模块、卫星射频接收模块、卫星基带处理模块、惯性导航模块和组合导航模块。方法包括如下步骤:北斗/GPS双模软件接收机完成首次定位后对惯性导航系统进行初始化;进入超紧组合模式,估计出I、Q信息IINS、QINS;北斗/GPS双模软件接收机得到的I、Q信息与IINS、QINS的差值作为EKF的观测量;对惯性导航系统进行校正;估计出卫星信号的载波多普勒频率和伪码码相位,分别反馈到载波NCO和码NCO,形成闭环回路。本发明可以改善卫星/惯性组合导航系统的定位精度、提高系统的鲁棒性,改善卫星接收机在高动态、微弱信号、强干扰条件下的跟踪性能。
Description
技术领域
本发明涉及北斗技术应用、多传感器数据融合及应用领域,尤其是一种基于北斗/GPS双模软件接收机的卫星/惯性超紧组合系统及方法。
背景技术
卫星导航系统和惯性导航系统具有高度的互补特性,这一特性已经被广泛应用于卫星/惯性组合导航系统。传统的卫星/惯性组合导航系统通常采用松组合和紧组合两种信息融合方式。在松组合工作方式下,卫星导航系统和惯性导航系统独立工作,分别输出位置、速度、姿态等导航信息,然后再经由导航滤波器进行数据融合。而紧组合则利用了更加原始的导航信息,将惯性导航系统结合星历估计出的伪距、伪距率与卫星接收机得到的伪距、伪距率之间的差值作为观测量。虽然这两种组合方式都可以得到比单个定位系统更加优良的定位精度,系统也更加鲁棒,但是信息融合之后得到的误差校正量都是针对惯性导航系统的,对于卫星接收机的性能改善并没有帮助。这种退耦带来的结果是接收机跟踪环路对噪声和载体运动高度敏感,极易在微弱信号及载体高动态环境下造成卫星失锁,使得接收机跟踪环成为整个系统当中最为薄弱的环节。
超紧组合就是针对上述问题而提出的新型信息融合方式。其主要思想就是要打破传统接收机的跟踪环路,而将组合滤波器估计出的跟踪环路误差反馈到接收机内部的载波和伪码发生器,形成新的闭合环路。通过惯性导航系统来辅助卫星接收机的跟踪环路可以极大改善接收机在高动态、微弱信号、强干扰等环境下的性能,进一步提高了系统的鲁棒性,真正意义上实现卫星和惯性导航系统的双向辅助。
北斗二代导航系统是我国拥有自主知识产权的卫星导航系统,对于争取我国空间信息领域的话语权、保证国内位置服务的安全性具有重要意义。除了军用领域不再受制于人,卫星导航系统在民用领域同样有着广阔的市场前景和巨大的经济价值,这也成为GPS、北斗、GLONASS、GALILEO四大卫星导航系统最为激烈的争夺方向。因而加速北斗相关产品的开发,拓宽北斗导航系统的应用领域成为一项迫切的需求。
多模接收机相比于单模接收机来说,可以同时利用多个卫星星座的导航信息,并从中选择出空间分布最好的几颗卫星参与定位,从而提高接收机的定位精度。同时,多个载波频率有利于载波相位整周模糊度的快速解算和减弱电离层效应的影响。而软件接收机相比于以往基于ASIC平台的接收机来说,其突出的可编程特性为修改接收机跟踪环路提供了可能。
以往的超紧组合研究存在以下几个问题:
一、研究的开展往往在松组合和紧组合的基础之上,通过增加惯性导航结果到卫星接收机的反馈通道来实现,没有进一步利用更加原始的导航信息。
二、由于国内外对GPS的研究相对成熟,因而对超紧组合的研究也通常是针对GPS卫星导航系统的超紧组合研究,对于多模接收机,尤其是涉及北斗接收机的超紧组合研究相对欠缺。
三、国内对于超紧组合的研究往往停留在理论和方法层面的论证和仿真,缺乏在工程实现方面的探究。
发明内容
本发明所要解决的技术问题在于,提供一种基于北斗/GPS双模软件接收机的卫星/惯性超紧组合系统及方法,能够改善卫星/惯性组合导航系统的定位精度、提高系统的鲁棒性,改善卫星接收机在高动态、微弱信号、强干扰条件下的跟踪性能。
为解决上述技术问题,本发明提供一种基于北斗/GPS双模软件接收机的卫星/惯性超紧组合系统,基于FPGA+DSP架构,包括:有源卫星天线模块、卫星射频接收模块、卫星基带处理模块、惯性导航模块和组合导航模块;有源卫星天线模块同时接收GPS L1和北斗B1信号,经功分器分为两路,分别接入卫星射频接收模块的两个通道;卫星射频接收模块具有GPS L1和北斗B1两个通道,信号经过低噪声放大、滤波和混频后输出近似零中频的I、Q两路基带信号,其中I表示中频信号同相输出,Q表示中频信号正交输出。中频信号进行AD采样转为数字量后送入卫星基带处理模块;卫星基带处理模块包括FPGA和DSP两个处理单元;组合导航模块由卫星基带处理模块中的DSP处理单元实现;惯性导航模块将位置、速度、加速度、姿态信息送入DSP处理单元,同时接收DSP处理单元发回的位置、速度、加速度、姿态的误差校正量。
优选的,FPGA处理单元分为多个独立的通道,分别对北斗和GPS卫星进行捕获跟踪。
优选的,DSP处理单元负责完成卫星接收机的算法,包括:通过相关累加结果判断是否成功捕获卫星,实现捕获决策,载波环和码环的鉴相,环路滤波及整个环路控制过程,位同步与帧同步,获取导航电文,完成PVT解算。
优选的,惯性导航模块包括3个单轴的加速度计和3个单轴的陀螺仪。
相应的,一种基于北斗/GPS双模软件接收机的卫星/惯性超紧组合方法,包括如下步骤:
(1)北斗/GPS双模软件接收机同时对北斗和GPS卫星进行捕获、跟踪、位同步和帧同步、双模定位,完成首次定位后对惯性导航系统进行初始化;
(2)完成步骤(1)后进入超紧组合模式,根据惯性导航系统得到的导航信息,结合星历估计出I、Q信息IINS、QINS,其中INS是惯性导航系统(Inertial Navigation System)的简写;
(3)北斗/GPS双模软件接收机得到的I、Q信息IGNSS、QGNSS与IINS、QINS的差值作为EKF的观测量,表示为:
zk=[dI1+ηI1,dI2+ηI2,…,dIm+ηIm,dQ1+ηQ1,dQ2+ηQ2,…,dQm+ηQm]T
其中I表示中频信号同相输出,Q表示中频信号正交输出,η代表噪声,下标数字代表卫星通道号,EKF是扩展卡尔曼滤波器(Extended Kalman Filter)的简写,GNSS是全球导航卫星系统(Global Navigation Satellite System)的简写;
(4)通过步骤(3)得到的误差信息,对惯性导航系统进行校正;
(5)通过步骤(4)校正后的位置、速度、加速度信息结合卫星星历估计出卫星信号的载波多普勒频率和伪码码相位,分别反馈到载波NCO和码NCO,形成闭环回路。其中NCO是数字控制振荡器(Numerically Controlled Oscillator)的简写。
优选的,步骤(1)中的双模定位过程将GPS时GPST和北斗时BDT转换到协调时间时UTC,并将转换误差作为未知数进行解算。
优选的,步骤(1)中的双模定位过程将计算各组卫星的加权GDOP值,并从中选取GDOP值最小的一组卫星进行定位。
优选的,步骤(3)中EKF的状态向量表示为:
其中,这17个状态分别为三个方向上的位置误差δR=[δx,δy,δz]T、速度误差加速度计漂移δa=[δax,δay,δaz]T、姿态误差陀螺仪随机漂移ε=[εx,εy,εz]T以及卫星钟差δtu和钟漂δtru,其中,前15维状态量为INS的状态误差,第16、17维状态量为卫星接收机的状态误差。
优选的,步骤(3)中EKF的观测矩阵表示为:
其中观测矩阵各因子可以表示为下式,i表示通道号,取值为1~m,m为接收机的最大通道数;
其余各因子都可以按照上式列出,其中载波角频率误差ωe和载波相位误差θe可以表示为:
E[I],E[Q]是I、Q信号的期望,ω表示载波角频率,c表示光速。Ru和分别为接收机到地心距离矢量的量测值和估计值;vu和分别为接收机到地心速度矢量的量测值和估计值;au和分别为接收机到地心加速度矢量的量测值和估计值。Re、ve和ae分别代表位置误差、速度误差和加速度误差的大小。
优选的,观测矩阵耦合了加速度信息与多普勒频率的变化率。
本发明的有益效果为:本发明使用的超紧组合方案利用了更加原始的导航信息I、Q值,相较于单纯的速度与位置信息,本方法在计算观测矩阵时耦合了加速度和多普勒频率变化率信息,可以改善卫星/惯性组合导航系统的定位精度、提高系统的鲁棒性,改善卫星接收机在高动态、微弱信号、强干扰条件下的跟踪性能;北斗/GPS双模定位可以同时利用北斗和GPS星座的卫星进行定位,通过卫星优选选择空间分布更好的卫星来改善几何精度因子,进一步提高系统的定位精度和可靠性。
附图说明
图1为本发明的超紧组合系统结构示意图。
图2为本发明的超紧组合方法流程示意图。
具体实施方式
如图1所示,一种基于北斗/GPS双模软件接收机的卫星/惯性超紧组合系统,在硬件上基于FPGA+DSP架构,包括:有源卫星天线模块、卫星射频接收模块、卫星基带处理模块、惯性导航模块和组合导航模块。
有源卫星天线模块可以同时接收GPS L1频点和北斗B1频点的信号,经功分器将信号分为两路,分别接入卫星射频接收模块的两个通道进行带通滤波和低噪声放大。GPS L1通道和北斗B1通道分别产生1575MHz和1561MHz的本振信号,与输入信号进行混频,从而将GPS L1信号下变频到0.42MHz近似零中频的I、Q两路基带信号,北斗B1信号下变频到0.098MHz近似零中频的I、Q两路基带信号,随后进行AD采样转为数字量并送入卫星基带处理模块。
卫星基带处理模块与卫星射频接收模块相连,包括FPGA和DSP两个处理单元。
其中FPGA处理单元接收到AD采样后的数字I、Q信号,随后完成载波剥离和伪码剥离,在完成相关累加后将累加结果送入DSP处理单元。
DSP处理单元与FPGA处理单元相连,负责完成卫星接收机的算法,包括:通过相关累加结果判断是否成功捕获卫星、实现捕获决策、载波环和码环的鉴相、环路滤波及整个环路控制过程、位同步与帧同步、获取导航电文、完成PVT解算。
组合导航模块也是由卫星基带处理模块中的DSP处理单元实现的。DSP处理单元在读取惯性导航模块送入的位置、速度、加速度与姿态信息后,结合导航电文中星历信息所提供的卫星位置和速度,得到一个I、Q信息的估计IINS和QINS。同时,DSP处理单元将经过FPGA处理过后的I、Q信息IGNSS、QGNSS和经过惯性导航模块估计出来的I、Q信息IINS、QINS分别作差,将差值作为导航滤波器的观测值。
惯性导航模块包括3个单轴的加速度计和3个单轴的陀螺仪。它与DSP处理单元相连,将位置、速度、加速度与姿态信息送入DSP处理单元,同时接收DSP处理单元发回的位置、速度、加速度与姿态的误差校正量。
DSP处理单元通过校正后的惯性导航数据,结合卫星星历估计出载波多普勒频率和伪码的码相位,反馈到FPGA处理单元的载波NCO和码NCO模块,从而调节本地载波的载波频率和本地伪码的码相位。
FPGA处理单元分为32个独立的通道,每个通道都可以独立地对一颗卫星进行捕获和跟踪。其中20个通道用于北斗卫星,12个通道用于GPS卫星。每个卫星通道所包含的载波NCO模块可以通过改变频率控制字来产生不同的载波频率,且上述卫星通道包含不同的码NCO模块,分别产生北斗CB1I码和GPS C/A码,分别对北斗卫星和GPS卫星进行捕获跟踪。
如图2所示,一种基于北斗/GPS双模软件接收机的卫星/惯性超紧组合方法,包括如下步骤:
(1)北斗/GPS双模软件接收机的32个通道同时对北斗和GPS卫星进行捕获、跟踪、位同步和子帧同步、获取双模观测量,完成PVT解算后获得首次定位信息,并用该定位信息完成惯性导航系统的初始化对准。在双模定位过程当中,需要将GPS时(GPST)和北斗时(BDT)全部转换到协调时间时(UTC),并将转换误差作为未知数进行解算,从而需要五颗卫星参与定位。而GPS系统所采用的WGS-84坐标系和北斗系统所采用的CGCS 2000坐标系之间的差异则可以忽略不计。
步骤(1)中的双模定位过程将计算各组卫星的加权GDOP值,并从中选取GDOP值最小的一组卫星进行定位。具体的选星算法可以分为如下几个步骤:
a.在H矩阵中选取行列式最小的k组可见卫星,其中H是伪距残差的量测矩阵。该可见卫星当中包括了北斗卫星和GPS卫星;
b.分别计算k组卫星的加权GDOP值;
c.在k组卫星之间选取加权GDOP值最小的一组可见卫星进行定位解算。
其中k值可以自由选取,一般k越大则计算量越大,但是得到的最小加权GDOP值越接近于实际的最小加权GDOP值。
(2)完成步骤(1)后进入超紧组合模式,卫星接收机在完成导航电文的译码之后不再进行PVT计算。惯性导航系统得到的位置、速度、加速度结合导航电文中的卫星星历所提供的卫星位置、速度信息,估计出I、Q信息IINS、QINS;
进一步的,由于GPS L1信号和北斗B1信号得到I、Q信息在表达式上没有区别,在后面的阐述中将不加区分地将I、Q信号的期望表示为:
这里利用I和Q的期望值E[I]和E[Q]代替I和Q本身进行建模。其中A为卫星的信号强度,k为测量计数,T为测量间隔时间,c为光速。载波角频率误差ωe和相位误差θe满足如下关系式:
其中,Ru和分别为接收机到地心距离矢量的量测值和估计值;vu和分别为接收机到地心速度矢量的量测值和估计值;au和分别为接收机到地心加速度矢量的量测值和估计值。Re、ve和ae分别代表位置误差、速度误差和加速度误差的大小,满足如下关系,其中下标e代表误差,(x,y,z)是三个方向上的位置测量值,是三个方向上的位置估计值,上标一点代表一阶导数,上标两点代表二阶导数:
(3)北斗/GPS双模软件接收机得到的I、Q信息IGNSS和QGNSS与惯性导航系统结合星历估计出来的I、Q信息IINS和QINS分别作差得到dI=IGNSS-IINS和dQ=QGNSS-QINS,将差值作为EKF的观测量,滤波后得到17个状态量;
进一步的,步骤(3)中的观测量可以表示为:
zk=[dI1+ηI1,dI2+ηI2,…,dI32+ηI32,dQ1+ηQ1,dQ2+ηQ2,…,dQ32+ηQ32]T
其中η代表噪声,下标数字代表1到32个卫星跟踪通道。
进一步的,步骤(3)中的EKF状态向量可以表示为:
其中,这17个状态分别为三个方向上的位置误差δR=[δx,δy,δz]T、速度误差加速度计漂移δa=[δax,δay,δaz]T、姿态误差和陀螺仪随机漂移ε=[εx,εy,εz]T以及卫星钟差δtu和钟漂δtru。其中,前15维状态量为INS的状态误差,第16、17维状态量为卫星接收机的状态误差。
进一步的,步骤(3)中的EKF,其观测矩阵可以表示为:
进一步的,上述观测矩阵各因子可以表示为下式,其中i表示通道号,取值为1~32。
进一步的,上述各个因子可通过下式得到:
(4)利用步骤(3)得到的误差信息,对惯性导航系统进行校正;
(5)利用步骤(4)校正后的位置、速度、加速度信息结合卫星星历所提供的卫星的位置、速度,估计出接收到的卫星信号的载波多普勒频率和伪码码相位,分别反馈到载波NCO和码NCO,形成闭环回路。
尽管本发明就优选实施方式进行了示意和描述,但本领域的技术人员应当理解,只要不超出本发明的权利要求所限定的范围,可以对本发明进行各种变化和修改。
Claims (10)
1.一种基于北斗/GPS双模软件接收机的卫星/惯性超紧组合系统,其特征在于,基于FPGA+DSP架构,包括:有源卫星天线模块、卫星射频接收模块、卫星基带处理模块、惯性导航模块和组合导航模块;有源卫星天线模块同时接收GPS L1和北斗B1信号,经功分器分为两路,分别接入卫星射频接收模块的两个通道;卫星射频接收模块具有GPS L1和北斗B1两个通道,信号经过低噪声放大、滤波和混频后输出近似零中频的I、Q两路基带信号,进行AD采样转为数字量后送入卫星基带处理模块;卫星基带处理模块包括FPGA和DSP两个处理单元;组合导航模块由卫星基带处理模块中的DSP处理单元实现;惯性导航模块将位置、速度、加速度、姿态信息送入DSP处理单元,同时接收DSP处理单元发回的位置、速度、加速度、姿态的误差校正量。
2.如权利要求1所述的基于北斗/GPS双模软件接收机的卫星/惯性超紧组合系统,其特征在于,FPGA处理单元分为多个独立的通道,分别对北斗和GPS卫星进行捕获跟踪。
3.如权利要求1所述的基于北斗/GPS双模软件接收机的卫星/惯性超紧组合系统,其特征在于,DSP处理单元负责完成卫星接收机的算法,包括:通过相关累加结果判断是否成功捕获卫星,实现捕获决策,载波环和码环的鉴相,环路滤波及整个环路控制过程,位同步与帧同步,获取导航电文,完成PVT解算。
4.如权利要求1所述的基于北斗/GPS双模软件接收机的卫星/惯性超紧组合系统,其特征在于,惯性导航模块包括3个单轴的加速度计和3个单轴的陀螺仪。
5.一种基于北斗/GPS双模软件接收机的卫星/惯性超紧组合方法,其特征在于,包括如下步骤:
(1)北斗/GPS双模软件接收机同时对北斗和GPS卫星进行捕获、跟踪、位同步和帧同步、双模定位,完成首次定位后对惯性导航系统进行初始化;
(2)完成步骤(1)后进入超紧组合模式,根据惯性导航系统得到的导航信息,结合星历估计出I、Q信息IINS、QINS;
(3)北斗/GPS双模软件接收机得到的I、Q信息IGNSS、QGNSS与IINS、QINS的差值作为EKF的观测量,表示为:
zk=[dI1+ηI1,dI2+ηI2,…,dIm+ηIm,dQ1+ηQ1,dQ2+ηQ2,…,dQm+ηQm]T
其中I表示中频信号同相输出,Q表示中频信号正交输出,η代表噪声,下标数字代表卫星通道号;
(4)通过步骤(3)得到的误差信息,对惯性导航系统进行校正;
(5)通过步骤(4)校正后的位置、速度、加速度信息结合卫星星历估计出卫星信号的载波多普勒频率和伪码码相位,分别反馈到载波NCO和码NCO,形成闭环回路。
6.如权利要求5所述的基于北斗/GPS双模软件接收机的卫星/惯性超紧组合方法,其特征在于,步骤(1)中的双模定位过程将GPS时GPST和北斗时BDT转换到协调时间时UTC,并将转换误差作为未知数进行解算。
7.如权利要求5所述的基于北斗/GPS双模软件接收机的卫星/惯性超紧组合方法,其特征在于,步骤(1)中的双模定位过程将计算各组卫星的加权GDOP值,并从中选取GDOP值最小的一组卫星进行定位。
8.如权利要求5所述的基于北斗/GPS双模软件接收机的卫星/惯性超紧组合方法,其特征在于,步骤(3)中EKF的状态向量表示为:
其中,这17个状态分别为三个方向上的位置误差δR=[δx,δy,δz]T、速度误差加速度计漂移δa=[δax,δay,δaz]T、姿态误差陀螺仪随机漂移ε=[εx,εy,εz]T以及卫星钟差δtu和钟漂δtru,其中,前15维状态量为INS的状态误差,第16、17维状态量为卫星接收机的状态误差。
9.如权利要求5所述的基于北斗/GPS双模软件接收机的卫星/惯性超紧组合方法,其特征在于,步骤(3)中EKF的观测矩阵表示为:
<mrow>
<msub>
<mi>H</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<msub>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>h</mi>
<mrow>
<mi>I</mi>
<mi>x</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>h</mi>
<mrow>
<mi>I</mi>
<mi>y</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>h</mi>
<mrow>
<mi>I</mi>
<mi>z</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>h</mi>
<mrow>
<mi>I</mi>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>h</mi>
<mrow>
<mi>I</mi>
<mover>
<mi>y</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>h</mi>
<mrow>
<mi>I</mi>
<mover>
<mi>z</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>h</mi>
<mrow>
<mi>I</mi>
<mover>
<mi>x</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>h</mi>
<mrow>
<mi>I</mi>
<mover>
<mi>y</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>h</mi>
<mrow>
<mi>I</mi>
<mover>
<mi>z</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>h</mi>
<mrow>
<mi>I</mi>
<mi>x</mi>
<mi>m</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>h</mi>
<mrow>
<mi>I</mi>
<mi>y</mi>
<mi>m</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>h</mi>
<mrow>
<mi>I</mi>
<mi>z</mi>
<mi>m</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>h</mi>
<mrow>
<mi>I</mi>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>m</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>h</mi>
<mrow>
<mi>I</mi>
<mover>
<mi>y</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>m</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>h</mi>
<mrow>
<mi>I</mi>
<mover>
<mi>z</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>m</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>h</mi>
<mrow>
<mi>I</mi>
<mover>
<mi>x</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mi>m</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>h</mi>
<mrow>
<mi>I</mi>
<mover>
<mi>y</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mi>m</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>h</mi>
<mrow>
<mi>I</mi>
<mover>
<mi>z</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mi>m</mi>
</mrow>
</msub>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>h</mi>
<mrow>
<mi>Q</mi>
<mi>x</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>h</mi>
<mrow>
<mi>Q</mi>
<mi>y</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>h</mi>
<mrow>
<mi>Q</mi>
<mi>z</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>h</mi>
<mrow>
<mi>Q</mi>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>h</mi>
<mrow>
<mi>Q</mi>
<mover>
<mi>y</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>h</mi>
<mrow>
<mi>Q</mi>
<mover>
<mi>z</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>h</mi>
<mrow>
<mi>Q</mi>
<mover>
<mi>x</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>h</mi>
<mrow>
<mi>Q</mi>
<mover>
<mi>y</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>h</mi>
<mrow>
<mi>Q</mi>
<mover>
<mi>z</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>h</mi>
<mrow>
<mi>I</mi>
<mi>x</mi>
<mi>m</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>h</mi>
<mrow>
<mi>I</mi>
<mi>y</mi>
<mi>m</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>h</mi>
<mrow>
<mi>I</mi>
<mi>z</mi>
<mi>m</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>h</mi>
<mrow>
<mi>I</mi>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>m</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>h</mi>
<mrow>
<mi>I</mi>
<mover>
<mi>y</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>m</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>h</mi>
<mrow>
<mi>I</mi>
<mover>
<mi>z</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>m</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>h</mi>
<mrow>
<mi>I</mi>
<mover>
<mi>x</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mi>m</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>h</mi>
<mrow>
<mi>I</mi>
<mover>
<mi>y</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mi>m</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>h</mi>
<mrow>
<mi>I</mi>
<mover>
<mi>z</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mi>m</mi>
</mrow>
</msub>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<mn>2</mn>
<mi>m</mi>
<mo>&times;</mo>
<mn>17</mn>
</mrow>
</msub>
</mrow>
其中观测矩阵各因子可以表示为下式,i表示通道号,取值为1~m,m为接收机的最大通道数;
<mrow>
<msub>
<mi>h</mi>
<mrow>
<mi>I</mi>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>6</mn>
</mfrac>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>E</mi>
<mo>&lsqb;</mo>
<mi>I</mi>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&theta;</mi>
<mi>e</mi>
</msub>
</mrow>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&theta;</mi>
<mi>e</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>E</mi>
<mo>&lsqb;</mo>
<mi>I</mi>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&omega;</mi>
<mi>e</mi>
</msub>
</mrow>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&omega;</mi>
<mi>e</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>h</mi>
<mrow>
<mi>I</mi>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>6</mn>
</mfrac>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>E</mi>
<mo>&lsqb;</mo>
<mi>I</mi>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&theta;</mi>
<mi>e</mi>
</msub>
</mrow>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&theta;</mi>
<mi>e</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>E</mi>
<mo>&lsqb;</mo>
<mi>I</mi>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&omega;</mi>
<mi>e</mi>
</msub>
</mrow>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&omega;</mi>
<mi>e</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>h</mi>
<mrow>
<mi>I</mi>
<mover>
<mi>x</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>6</mn>
</mfrac>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>E</mi>
<mo>&lsqb;</mo>
<mi>I</mi>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&theta;</mi>
<mi>e</mi>
</msub>
</mrow>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&theta;</mi>
<mi>e</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mover>
<mi>x</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>E</mi>
<mo>&lsqb;</mo>
<mi>I</mi>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&omega;</mi>
<mi>e</mi>
</msub>
</mrow>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&omega;</mi>
<mi>e</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mover>
<mi>x</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>h</mi>
<mrow>
<mi>Q</mi>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>6</mn>
</mfrac>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>E</mi>
<mo>&lsqb;</mo>
<mi>Q</mi>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&theta;</mi>
<mi>e</mi>
</msub>
</mrow>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&theta;</mi>
<mi>e</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>E</mi>
<mo>&lsqb;</mo>
<mi>Q</mi>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&omega;</mi>
<mi>e</mi>
</msub>
</mrow>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&omega;</mi>
<mi>e</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>h</mi>
<mrow>
<mi>Q</mi>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>6</mn>
</mfrac>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>E</mi>
<mo>&lsqb;</mo>
<mi>Q</mi>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&theta;</mi>
<mi>e</mi>
</msub>
</mrow>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&theta;</mi>
<mi>e</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>E</mi>
<mo>&lsqb;</mo>
<mi>Q</mi>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&omega;</mi>
<mi>e</mi>
</msub>
</mrow>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&omega;</mi>
<mi>e</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>h</mi>
<mrow>
<mi>Q</mi>
<mover>
<mi>x</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>6</mn>
</mfrac>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>E</mi>
<mo>&lsqb;</mo>
<mi>Q</mi>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&theta;</mi>
<mi>e</mi>
</msub>
</mrow>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&theta;</mi>
<mi>e</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mover>
<mi>x</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>E</mi>
<mo>&lsqb;</mo>
<mi>Q</mi>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&omega;</mi>
<mi>e</mi>
</msub>
</mrow>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&omega;</mi>
<mi>e</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mover>
<mi>x</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
其余各因子都可以按照上式列出,其中载波角频率误差ωe和载波相位误差θe可以表示为:
<mrow>
<msub>
<mi>&omega;</mi>
<mi>e</mi>
</msub>
<mo>=</mo>
<mfrac>
<mi>&omega;</mi>
<mi>c</mi>
</mfrac>
<mo>&lsqb;</mo>
<mo>|</mo>
<mrow>
<msub>
<mi>v</mi>
<mi>u</mi>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>v</mi>
<mo>^</mo>
</mover>
<mi>u</mi>
</msub>
</mrow>
<mo>|</mo>
<mo>+</mo>
<mo>|</mo>
<mrow>
<msub>
<mi>a</mi>
<mi>u</mi>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>a</mi>
<mo>^</mo>
</mover>
<mi>u</mi>
</msub>
</mrow>
<mo>|</mo>
<mi>t</mi>
<mo>&rsqb;</mo>
<mo>=</mo>
<mfrac>
<mi>&omega;</mi>
<mi>c</mi>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>v</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<msub>
<mi>a</mi>
<mi>e</mi>
</msub>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&theta;</mi>
<mi>e</mi>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mi>&omega;</mi>
<mi>c</mi>
</mfrac>
<mo>&lsqb;</mo>
<mrow>
<mo>|</mo>
<mrow>
<msub>
<mi>R</mi>
<mi>u</mi>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>R</mi>
<mo>^</mo>
</mover>
<mi>u</mi>
</msub>
</mrow>
<mo>|</mo>
</mrow>
<mo>-</mo>
<mrow>
<mo>|</mo>
<mrow>
<msub>
<mi>v</mi>
<mi>u</mi>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>v</mi>
<mo>^</mo>
</mover>
<mi>u</mi>
</msub>
</mrow>
<mo>|</mo>
</mrow>
<mi>t</mi>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mrow>
<mo>|</mo>
<mrow>
<msub>
<mi>a</mi>
<mi>u</mi>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>a</mi>
<mo>^</mo>
</mover>
<mi>u</mi>
</msub>
</mrow>
<mo>|</mo>
</mrow>
<msup>
<mi>t</mi>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mi>&omega;</mi>
<mi>c</mi>
</mfrac>
<mo>&lsqb;</mo>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>-</mo>
<msub>
<mi>v</mi>
<mi>e</mi>
</msub>
<mi>t</mi>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msub>
<mi>a</mi>
<mi>e</mi>
</msub>
<msup>
<mi>t</mi>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
</mrow>
E[I],E[Q]是I、Q信号的期望,ω表示载波角频率,c表示光速;Ru和分别为接收机到地心距离矢量的量测值和估计值;vu和分别为接收机到地心速度矢量的量测值和估计值;au和分别为接收机到地心加速度矢量的量测值和估计值;Re、ve和ae分别代表位置误差、速度误差和加速度误差的大小。
10.如权利要求5所述的基于北斗/GPS双模软件接收机的卫星/惯性超紧组合方法,其特征在于,观测矩阵耦合了加速度信息与多普勒频率的变化率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710695052.9A CN107656300B (zh) | 2017-08-15 | 2017-08-15 | 基于北斗/gps双模软件接收机的卫星/惯性超紧组合方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710695052.9A CN107656300B (zh) | 2017-08-15 | 2017-08-15 | 基于北斗/gps双模软件接收机的卫星/惯性超紧组合方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107656300A true CN107656300A (zh) | 2018-02-02 |
CN107656300B CN107656300B (zh) | 2020-10-02 |
Family
ID=61128514
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710695052.9A Active CN107656300B (zh) | 2017-08-15 | 2017-08-15 | 基于北斗/gps双模软件接收机的卫星/惯性超紧组合方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107656300B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108896044A (zh) * | 2018-07-23 | 2018-11-27 | 湖南格纳微信息科技有限公司 | 一种基于惯性导航和卫星导航的定位方法及装置 |
CN109150407A (zh) * | 2018-09-25 | 2019-01-04 | 中国科学院国家授时中心 | 多路复用信号伪随机码盲估方法、误码纠错方法及装置 |
CN109212573A (zh) * | 2018-10-15 | 2019-01-15 | 东南大学 | 一种城市峡谷环境下用于测绘车辆的定位系统及方法 |
CN109307873A (zh) * | 2018-11-08 | 2019-02-05 | 北京理工大学 | 一种INS辅助的双Kalman滤波器卫星信号跟踪环路 |
CN110058288A (zh) * | 2019-04-28 | 2019-07-26 | 北京微克智飞科技有限公司 | 无人机ins/gnss组合导航系统航向误差修正方法及系统 |
CN110133699A (zh) * | 2019-05-23 | 2019-08-16 | 深圳市联和安业科技有限公司 | 一种超紧组合导航定位系统及方法 |
CN112612047A (zh) * | 2020-11-27 | 2021-04-06 | 航天恒星科技有限公司 | 一种卫星/惯性组合的微型导航系统 |
CN113701751A (zh) * | 2021-08-24 | 2021-11-26 | 西安大衡天成信息科技有限公司 | 一种基于多波束天线的导航装置 |
US11550063B2 (en) * | 2021-04-30 | 2023-01-10 | Samsung Electronics Co., Ltd | Method and apparatus for synchronizing with global navigation satellite system |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2007135115A1 (fr) * | 2006-05-19 | 2007-11-29 | Thales | Dispositif de navigation aerienne a capteurs inertiels et recepteurs de radionavigation et procede de navigation aerienne utilisant de tels elements |
CN103969672A (zh) * | 2014-05-14 | 2014-08-06 | 东南大学 | 一种多卫星系统与捷联惯性导航系统紧组合导航方法 |
CN105068102A (zh) * | 2015-08-11 | 2015-11-18 | 南京理工大学 | 一种基于dsp+fpga的超紧组合导航方法 |
CN105758401A (zh) * | 2016-05-14 | 2016-07-13 | 中卫物联成都科技有限公司 | 一种基于多源信息融合的组合导航方法和设备 |
CN106526635A (zh) * | 2016-12-29 | 2017-03-22 | 中国人民解放军国防科学技术大学 | 一种gnss信号载波跟踪和导航解算紧组合的滤波方法 |
-
2017
- 2017-08-15 CN CN201710695052.9A patent/CN107656300B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2007135115A1 (fr) * | 2006-05-19 | 2007-11-29 | Thales | Dispositif de navigation aerienne a capteurs inertiels et recepteurs de radionavigation et procede de navigation aerienne utilisant de tels elements |
CN103969672A (zh) * | 2014-05-14 | 2014-08-06 | 东南大学 | 一种多卫星系统与捷联惯性导航系统紧组合导航方法 |
CN105068102A (zh) * | 2015-08-11 | 2015-11-18 | 南京理工大学 | 一种基于dsp+fpga的超紧组合导航方法 |
CN105758401A (zh) * | 2016-05-14 | 2016-07-13 | 中卫物联成都科技有限公司 | 一种基于多源信息融合的组合导航方法和设备 |
CN106526635A (zh) * | 2016-12-29 | 2017-03-22 | 中国人民解放军国防科学技术大学 | 一种gnss信号载波跟踪和导航解算紧组合的滤波方法 |
Non-Patent Citations (2)
Title |
---|
余威 等: ""基于DSP+FPGA的超紧组合导航系统研究"", 《航天控制》 * |
谢非: ""北斗软件接收机机及惯性/卫星超紧组合导航关键技术研究"", 《中国博士学位论文全文数据库 信息科技辑》 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108896044A (zh) * | 2018-07-23 | 2018-11-27 | 湖南格纳微信息科技有限公司 | 一种基于惯性导航和卫星导航的定位方法及装置 |
CN108896044B (zh) * | 2018-07-23 | 2021-09-10 | 湖南格纳微信息科技有限公司 | 一种基于惯性导航和卫星导航的定位方法及装置 |
CN109150407A (zh) * | 2018-09-25 | 2019-01-04 | 中国科学院国家授时中心 | 多路复用信号伪随机码盲估方法、误码纠错方法及装置 |
CN109150407B (zh) * | 2018-09-25 | 2021-12-17 | 中国科学院国家授时中心 | 多路复用信号伪随机码盲估方法、误码纠错方法及装置 |
CN109212573B (zh) * | 2018-10-15 | 2020-02-07 | 东南大学 | 一种城市峡谷环境下用于测绘车辆的定位系统及方法 |
WO2020077941A1 (zh) * | 2018-10-15 | 2020-04-23 | 东南大学 | 一种城市峡谷环境下用于测绘车辆的定位系统及方法 |
CN109212573A (zh) * | 2018-10-15 | 2019-01-15 | 东南大学 | 一种城市峡谷环境下用于测绘车辆的定位系统及方法 |
CN109307873A (zh) * | 2018-11-08 | 2019-02-05 | 北京理工大学 | 一种INS辅助的双Kalman滤波器卫星信号跟踪环路 |
CN110058288A (zh) * | 2019-04-28 | 2019-07-26 | 北京微克智飞科技有限公司 | 无人机ins/gnss组合导航系统航向误差修正方法及系统 |
CN110133699A (zh) * | 2019-05-23 | 2019-08-16 | 深圳市联和安业科技有限公司 | 一种超紧组合导航定位系统及方法 |
CN112612047A (zh) * | 2020-11-27 | 2021-04-06 | 航天恒星科技有限公司 | 一种卫星/惯性组合的微型导航系统 |
CN112612047B (zh) * | 2020-11-27 | 2024-03-15 | 航天恒星科技有限公司 | 一种卫星/惯性组合的微型导航系统 |
US11550063B2 (en) * | 2021-04-30 | 2023-01-10 | Samsung Electronics Co., Ltd | Method and apparatus for synchronizing with global navigation satellite system |
CN113701751A (zh) * | 2021-08-24 | 2021-11-26 | 西安大衡天成信息科技有限公司 | 一种基于多波束天线的导航装置 |
CN113701751B (zh) * | 2021-08-24 | 2024-05-07 | 西安大衡天成信息科技有限公司 | 一种基于多波束天线的导航装置 |
Also Published As
Publication number | Publication date |
---|---|
CN107656300B (zh) | 2020-10-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107656300A (zh) | 基于北斗/gps双模软件接收机的卫星/惯性超紧组合系统及方法 | |
CN109212573B (zh) | 一种城市峡谷环境下用于测绘车辆的定位系统及方法 | |
CN103675861B (zh) | 一种基于星载gnss多天线的卫星自主定轨方法 | |
CN104297773B (zh) | 一种高精度北斗三频sins深组合导航系统 | |
CN101666650B (zh) | 一种sins/gps超紧致组合导航系统及实现方法 | |
CN105607076B (zh) | 一种北斗二代b1和b3双频接收机 | |
CN100381835C (zh) | 基于数字广播电视信号的无线电组合定位方法 | |
CN103941271B (zh) | 一种时间-空间差分的gps/sins超紧组合导航方法 | |
CN108196272A (zh) | 一种基于实时精密单点定位的卫星导航定位装置及方法 | |
CN100437143C (zh) | 区域卫星导航系统及方法 | |
CN109782313A (zh) | 地面精密单点定位数据处理方法 | |
CN106019333B (zh) | 一种基于非相干鉴别器的北斗导航信号矢量跟踪方法 | |
CN110109158A (zh) | 基于gps、glonass和bds多系统的事后超快速rtk定位算法 | |
CN104614739B (zh) | 基于抗干扰滤波器的北斗多频接收机信号联合跟踪方法 | |
CN107607971A (zh) | 基于gnss共视时间比对算法的时间频率传递方法及接收机 | |
US11409005B2 (en) | GNSS receiver adapted to produce, use, and communicate software-generated satellite signal data | |
CN102983901B (zh) | 兼容gps、bd2和glonass系统的辅助同步方法 | |
CN112748449A (zh) | 一种卫星导航接收机锁相环与锁频环结合的矢量跟踪方法 | |
CN202305821U (zh) | 精密单点定位与惯性测量紧组合导航系统 | |
Chi et al. | Enabling robust and accurate navigation for UAVs using real-time GNSS precise point positioning and IMU integration | |
US10534087B1 (en) | Differential vector phase locked loop GPS reception method | |
CN102621569A (zh) | 一种分布滤波式全球定位及捷联惯导组合导航方法 | |
Brewer et al. | Differential vector phase locked loop | |
CN105652299B (zh) | 一种基于最大相关信号能量的卫星导航定位方法 | |
TW448304B (en) | Fully-coupled positioning process and system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |