CN102508268A - 一种地基导航信号观测数据生成方法 - Google Patents
一种地基导航信号观测数据生成方法 Download PDFInfo
- Publication number
- CN102508268A CN102508268A CN2011103543985A CN201110354398A CN102508268A CN 102508268 A CN102508268 A CN 102508268A CN 2011103543985 A CN2011103543985 A CN 2011103543985A CN 201110354398 A CN201110354398 A CN 201110354398A CN 102508268 A CN102508268 A CN 102508268A
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- mover
- satellite
- msup
- 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
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000012937 correction Methods 0.000 claims abstract description 4
- 239000005433 ionosphere Substances 0.000 claims description 15
- 238000004364 calculation method Methods 0.000 claims description 9
- 230000008054 signal transmission Effects 0.000 claims description 4
- 239000004576 sand Substances 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Chemical compound O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 239000005436 troposphere Substances 0.000 abstract description 8
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 abstract description 6
- 238000012795 verification Methods 0.000 abstract description 4
- 238000004088 simulation Methods 0.000 description 13
- 238000010586 diagram Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 238000009432 framing Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 239000013078 crystal Substances 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000007774 longterm Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
Images
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种地基导航信号观测数据生成方法,步骤为(1)设定地面信标站的位置参数、时间参数和环境参数;(2)根据预设的用户卫星轨道类型、用户卫星的初轨和系统时间来设置卫星轨道,生成相应的轨道点位,之后设置每一个轨道点位对应的卫星速度信息(3)计算钟差参数;(4)计算电离层延迟量;(5)计算对流层延迟量;(6)计算地面信标站发射信号时刻的信标位置坐标(7)根据步骤(3)、步骤(4)、步骤(5)和步骤(6)中得到的计算结果计算得到观测数据:伪距、伪距率和载波相位。本发明观测数据生成方法可运用到地基信号的射频信号生成中,为星载精密定轨接收机提供可靠的验证手段,提升在轨运行能力。
Description
技术领域
本发明涉及一种地基导航信号观测数据生成方法,为地基导航接收机的研究和开发提供了必要条件。
背景技术
地基导航技术,作为高分辨率对地观测系统的核心技术,可以进行精密定轨,提升高精度卫星对地观测能力。当今的激光三维测绘卫星、干涉雷达卫星等对精密定轨提出了很高的要求。
目前国际上通过激光测轨、GPS等方式对卫星进行定轨。地基导航是通过多普勒精密测速技术对卫星轨道位置进行精确测定,并成功的把测速误差提升到毫米级,径向定轨误差提升到厘米级。
我国目前尚不具备卫星精密定轨系统,对基于地基导航的精密定轨技术处于研制开发阶段。对地基导航信号的模拟,为卫星的在轨测试提供了有效的验证手段,提升了地基导航设备的在轨运行质量,是目前比较认可的一种地基导航的研究手段。但是对于地基信号得生成技术的研究多处于实验室内的软件仿真状态,通过软件接收机进行解算,这种方法的缺点是不能真实再现实际的地基导航信号,不能对实际接收设备进行调试与测试,不能对地基导航接收机设备进行性能指标评价。
发明内容
本发明解决的技术问题:克服现有技术的不足,本发明提出了一种地基导航信号生成方法,真正实现了用于精密定轨的地基导航数据以及物理射频信号模拟。
本发明所采用的技术解决方案:
一种地基导航信号观测数据生成方法,所述观测数据包括伪距、伪距率和载波相位,步骤如下:
(1)设定地面信标站的位置参数、时间参数和环境参数,其中,环境参数包括所述地面信标站所处位置的大气压力、温度和该信标站的工作状态;时间参数是指地面信标站开始工作的时间;根据地面信标站的位置参数、时间参数和环境参数生成导航电文;
(2)根据预设的用户卫星轨道类型、用户卫星的初轨和系统时间来设置卫星轨道,生成相应的轨道点位(Xu,Yu,Zu),之后设置每一个轨道点位对应的卫星速度信息;
(4)计算电离层延迟量dion;
(5)计算对流层延迟量dtrop;
(6)根据公式 计算得到地面信标站发射信号时刻的信标位置坐标(Xs,Ys,Zs),其中,(X,Y,Z)为卫星接收信号时刻的地面标信站位置坐标,修正角度α=ωe(tR-tT);ωe为地球自转角速度;tR-tT为从地面信标站发射信号到卫星接收信号的时间差;
(7)根据步骤(3)、步骤(4)、步骤(5)和步骤(6)中得到的计算结果通过如下方式计算得到观测数据:
根据公式 计算得到伪距ρ,
根据 计算伪距率其中,
步骤(4)中所述计算电离层延迟量dion采用8参数Klobuchar电离层模型,具体方法为:
通过公式dion=C×F×IZ(t)计算电离层延迟量dion,
倾斜因子F=1.0+16.0×(0.53-E)3,E为卫星高度角,
A2为白天余弦曲线的幅度,φM为穿刺点M的磁纬,φM=φi+0.064cos(λi-1.617),φi和λi为接收机到卫星之间的连线与电离层交点处的大地纬度和大地经度;A3为余弦曲线极点的地方时,且A3=50400S,A4为余弦曲线的周期且有
步骤(5)中所述计算对流层延迟量dtrop采用Hopfield对流层延迟模型,具体方法为:
通过公式 计算对流层延迟,其中E为卫星高度角, Ps和Ts分别为地面信标站的大气压与热力学温度,hd=40136+148.72×(Ts-273.16),hs为地面信标站高度, es为地面信标站处水汽分压。
本发明与现有技术相比的有益效果在于:
(1)本发明在原有对地基导航信号研究的基础上,建立了一种地基导航信号观测数据生成方法,并可运用到地基信号的射频信号仿真生成,为星载精密定轨接收机提供可靠的验证手段,提升在轨运行能力。
(2)本发明通过仿真场景设置、误差量计算、观测量与电文生成的流程把相应的数学模型转化成信号生成系统可识别的数据,保证了系统的健壮性,该数据生成方法也可作为算法验证的有力工具
(3)本发明采用统一灵活的电文处理手段,使得多频点、非整秒电文处理以及不同电文速率的调制应用,都能得到统一的处理。能够通过对载波频率的精确控制,获得高精度载波多普勒的仿真,提供高精度载波仿真方法。
附图说明
图1为地基导航方式示意图;
图2为地基导航信号生成方法流程;
图3为地基导航电文组帧方式;
图4为地基导航观测数据组帧方式;
图5地基导航电文处理方式;
图6地基导航信号源DSP部分处理结构图;
图7地基导航信号源FPGA部分处理结构图;
具体实施方式
下面结合附图对本发明具体实施方式做进一步介绍。
图1是地基导航应用于卫星定轨的示意图,信号模拟覆盖卫星运动轨迹的信标站,星载接收机接收相应的信标信号,接收机要分别测量P、L、S频点的导航信号,其中一般的测量方法是,捕获跟踪P、L频点,并优选出一个频点的信号引导对S频点信号的捕获,并根据捕获到的三个频点的信号测量它们的多普勒频移,获得多普勒频移后,就可以根据选取P、S或L、S的频率,通过双频方式消除电离层影响。根据卫星运动模型以及测定的多普勒频移就可以获得星上接收机的轨道位置。
本地基导航信号数据生成方法的基础在于其可应用于射频信号的仿真应用,该信号的结构为BPSK方式调制,包括载波、伪码和导航电文三种信息分量,信号生成的表达式如下:
如图2所示,本发明提供了一种地基导航信号观测数据生成方法,观测数据包括伪距、伪距率和载波相位,步骤如下:
(1)设定地面信标站的位置参数、时间参数和环境参数,其中,环境参数包括所述地面信标站所处位置的大气压力、温度和该信标站的工作状态;时间参数是指地面信标站开始工作的时间;根据地面信标站的位置参数、时间参数和环境参数生成导航电文;
传输数据的组帧方式见图3,具体包括帧头、信标标志和信标状态(电压、电流,温湿压传感器、超稳晶振及3台发射机工作状态);以2010年1月1日为起点的年、日和秒数;温度、湿度、气压等信息。其中,信标状态包括电压值、电流值;温度、湿度和压传感器;超稳晶振及3台发射机工作状态。
导航电文的发送周期设计为30秒,前10秒设计按导航电文编码格式1发送,中间10设计按导航电文编码格式2发送,后10秒设计按导航电文编码格式3发送,30秒为一个周期。
(2)根据预设的用户卫星轨道类型、用户卫星的初轨和系统时间来设置卫星轨道,生成相应的轨道点位(Xu,Yu,Zu),之后设置每一个轨道点位对应的卫星速度信息;
(4)计算电离层延迟量dion;
根据系统运算量和信号生成精度分析,采用8参数Klobuchar模型计算电离层延迟,能够模拟真实自然条件下电离层对电磁波传播的影响,与真实数据具有相同的变化趋势及比例关系,长期吻合程度达到60%以上。可根据如下叙述对电离层模型进行仿真:
通过公式dion=C×F×IZ(t)计算电离层延迟量dion,
倾斜因子F=1.0+16.0×(0.53-E)3,E为卫星高度角,
垂直方向延迟t为以秒为单位的接收机到卫星之间的连线与电离层交点处的地方时;A1=5×10-9s为夜间值的垂直延迟常数,
A2为白天余弦曲线的幅度,φM为穿刺点M的磁纬,φM=φi+0.064cos(λi-1.617),φi和λi为接收机到卫星之间的连线与电离层交点处的大地纬度和大地经度;A3为余弦曲线极点的地方时,且A3=50400S,A4为余弦曲线的周期且有
(5)计算对流层延迟量dtrop;
对流层延迟量:模拟真实自然条件下对流层对电磁波传播的影响,可以实现与真实数据具有相同的变化趋势及比例关系,长期吻合程度达到60%以上,对六层延迟量生成有Hopfield对流层延迟模型、改进Hopfield对流层延迟模型、Saastamoinen对流层延迟模型、常数模型等,综合计算量和精度考虑采用Hopfield对流层延迟模型,具体仿真时根据相应模型进行,并且可以利用标准大气模型模拟大气中的温度、湿度、气压等参数。
对流层延迟量dtrop采用Hopfield对流层延迟模型,具体方法为:
通过公式 计算对流层延迟,其中E为卫星高度角, Ps和Ts为地面信标站处大气压与热力学温度,hd=40136+148.72×(Ts-273.16),hs为地面信标站高度,
下面是地球自转数据生成:
(6)根据公式 计算得到地面信标站发射信号时刻的信标位置坐标(Xs,Ys,Zs),其中,(X,Y,Z)为卫星接收信号时刻的地面标信站位置坐标,修正角度α=ωe(tR-tT);ωe为地球自转角速度;tR-tT为从地面信标站发射信号到卫星接收信号的时间差;
根据用户卫星轨道数据、地面信标站位置坐标和空间环境效应仿真数据,生成地面信标站信号的伪距、伪距率和载波相位观测数据。考虑到射频信号精度仿真需求,仿真出伪距率。观测数据组帧方式如图4,每20毫秒组成一包观测数据包,包含观测数据索引、时间和每个信标站的伪距参量、伪距率参量和载波参量。
(7)根据步骤(3)、步骤(4)、步骤(5)和步骤(6)中得到的计算结果通过如下方式计算得到观测数据:
根据公式 计算得到伪距ρ,
这样就得到了地基导航信号的观测数据,载波相位、伪距和伪距率。
地基导航电文通过DSP进行格式转换,由FPGA进行调制输出,见图5。DSP根据观测数据来维护地面信标站和时间信息,并控制地基导航信号功率,见图6。FPGA根据DSP的控制信息和观测量生成地基导航信号数字中频,见图7。地基导航数字中频仿真信号输出给DAC模块,经DAC生成模拟信号,通过上变频即可生成地基导航射频信号进行发射。该信号可用于地基导航接收设备的调试与测试。
本发明未详细说明部分属本领域技术人员公知常识。
Claims (3)
1.一种地基导航信号观测数据生成方法,所述观测数据包括伪距、伪距率和载波相位,其特征在于步骤如下:
(1)设定地面信标站的位置参数、时间参数和环境参数,其中,环境参数包括所述地面信标站所处位置的大气压力、温度和该信标站的工作状态;时间参数是指地面信标站开始工作的时间;根据地面信标站的位置参数、时间参数和环境参数生成导航电文;
(2)根据预设的用户卫星轨道类型、用户卫星的初轨和系统时间来设置卫星轨道,生成相应的轨道点位(Xu,Yu,Zu),之后设置每一个轨道点位对应的卫星速度信息;
(4)计算电离层延迟量dion;
(5)计算对流层延迟量dtrop;
(6)根据公式 计算得到地面信标站发射信号时刻的信标位置坐标(Xs,Ys,Zs),其中,(X,Y,Z)为卫星接收信号时刻的地面标信站位置坐标,修正角度α=ωe(tR-tT);ωe为地球自转角速度;tR-tT为从地面信标站发射信号到卫星接收信号的时间差;
(7)根据步骤(3)、步骤(4)、步骤(5)和步骤(6)中得到的计算结果通过如下方式计算得到观测数据:
根据公式 计算得到伪距ρ,
2.根据权利要求1所述的一种地基导航信号观测数据生成方法,其特征在于:步骤(4)中所述计算电离层延迟量dion采用8参数Klobuchar电离层模型,具体方法为:
通过公式dion=C×F×IZ(t)计算电离层延迟量dion,
倾斜因子F=1.0+16.0×(0.53-E)3,E为卫星高度角,
A2为白天余弦曲线的幅度,φM为穿刺点M的磁纬,φM=φi+0.064cos(λi-1.617),φi和λi为接收机到卫星之间的连线与电离层交点处的大地纬度和大地经度;A3为余弦曲线极点的地方时,且A3=50400S,A4为余弦曲线的周期且有
3.根据权利要求1所述的一种地基导航信号观测数据生成方法,其特征在于:步骤(5)中所述计算对流层延迟量dtrop采用Hopfield对流层延迟模型,具体方法为:
通过公式 计算对流层延迟,其中E为卫星高度角, Ps和Ts分别为地面信标站的大气压与热力学温度,hd=40136+148.72×(Ts-273.16),hs为地面信标站高度, es为地面信标站处水汽分压。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110354398.5A CN102508268B (zh) | 2011-11-10 | 2011-11-10 | 一种地基导航信号观测数据生成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110354398.5A CN102508268B (zh) | 2011-11-10 | 2011-11-10 | 一种地基导航信号观测数据生成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102508268A true CN102508268A (zh) | 2012-06-20 |
CN102508268B CN102508268B (zh) | 2014-10-08 |
Family
ID=46220374
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110354398.5A Active CN102508268B (zh) | 2011-11-10 | 2011-11-10 | 一种地基导航信号观测数据生成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102508268B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103199949A (zh) * | 2013-04-10 | 2013-07-10 | 中国电子科技集团公司第五十四研究所 | 伪卫星定位系统中的零基线自闭环星地时间同步方法 |
CN105158777A (zh) * | 2015-07-31 | 2015-12-16 | 上海卫星工程研究所 | 用于无源测向定位的数据源生成方法 |
CN105182367A (zh) * | 2015-09-22 | 2015-12-23 | 桂林电子科技大学 | 一种提取bds卫星电离层穿刺点电子浓度的新方法 |
CN105607094A (zh) * | 2014-11-21 | 2016-05-25 | 航天恒星科技有限公司 | Gnss接收机模拟器及重力场反演系统 |
CN111323791A (zh) * | 2020-04-10 | 2020-06-23 | 航天恒星科技有限公司 | 导航增强产品的质量评估方法和装置 |
CN113204042A (zh) * | 2021-05-21 | 2021-08-03 | 北京交通大学 | 一种基于精密单点定位的多星座联合列车定位方法 |
CN113419223A (zh) * | 2021-05-24 | 2021-09-21 | 西安空间无线电技术研究所 | 一种模拟电离层对同源双频电磁波信号影响的方法和系统 |
CN115567133A (zh) * | 2022-10-11 | 2023-01-03 | 上海精密计量测试研究所 | 一种usb天基测控通信体制多普勒频率参数校准方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030154025A1 (en) * | 1999-04-23 | 2003-08-14 | Global Locate, Inc. | Method and apparatus for performing timing synchronization |
CN101206258A (zh) * | 2007-12-20 | 2008-06-25 | 上海伽利略导航有限公司 | 一种高灵敏度辅助定位系统及其数据处理的方法 |
CN101334463A (zh) * | 2008-07-29 | 2008-12-31 | 北京航空航天大学 | Gnss导航信号模拟器中时间参数的模拟方法 |
US20100231443A1 (en) * | 2009-03-11 | 2010-09-16 | Whitehead Michael L | Removing biases in dual frequency gnss receivers using sbas |
-
2011
- 2011-11-10 CN CN201110354398.5A patent/CN102508268B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030154025A1 (en) * | 1999-04-23 | 2003-08-14 | Global Locate, Inc. | Method and apparatus for performing timing synchronization |
CN101206258A (zh) * | 2007-12-20 | 2008-06-25 | 上海伽利略导航有限公司 | 一种高灵敏度辅助定位系统及其数据处理的方法 |
CN101334463A (zh) * | 2008-07-29 | 2008-12-31 | 北京航空航天大学 | Gnss导航信号模拟器中时间参数的模拟方法 |
US20100231443A1 (en) * | 2009-03-11 | 2010-09-16 | Whitehead Michael L | Removing biases in dual frequency gnss receivers using sbas |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103199949A (zh) * | 2013-04-10 | 2013-07-10 | 中国电子科技集团公司第五十四研究所 | 伪卫星定位系统中的零基线自闭环星地时间同步方法 |
CN103199949B (zh) * | 2013-04-10 | 2015-12-23 | 中国电子科技集团公司第五十四研究所 | 伪卫星定位系统中的零基线自闭环星地时间同步方法 |
CN105607094A (zh) * | 2014-11-21 | 2016-05-25 | 航天恒星科技有限公司 | Gnss接收机模拟器及重力场反演系统 |
CN105158777A (zh) * | 2015-07-31 | 2015-12-16 | 上海卫星工程研究所 | 用于无源测向定位的数据源生成方法 |
CN105158777B (zh) * | 2015-07-31 | 2017-08-29 | 上海卫星工程研究所 | 用于无源测向定位的数据源生成方法 |
CN105182367A (zh) * | 2015-09-22 | 2015-12-23 | 桂林电子科技大学 | 一种提取bds卫星电离层穿刺点电子浓度的新方法 |
CN111323791A (zh) * | 2020-04-10 | 2020-06-23 | 航天恒星科技有限公司 | 导航增强产品的质量评估方法和装置 |
CN113204042A (zh) * | 2021-05-21 | 2021-08-03 | 北京交通大学 | 一种基于精密单点定位的多星座联合列车定位方法 |
CN113419223A (zh) * | 2021-05-24 | 2021-09-21 | 西安空间无线电技术研究所 | 一种模拟电离层对同源双频电磁波信号影响的方法和系统 |
CN113419223B (zh) * | 2021-05-24 | 2023-11-10 | 西安空间无线电技术研究所 | 一种模拟电离层对同源双频电磁波信号影响的方法和系统 |
CN115567133A (zh) * | 2022-10-11 | 2023-01-03 | 上海精密计量测试研究所 | 一种usb天基测控通信体制多普勒频率参数校准方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102508268B (zh) | 2014-10-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102508268B (zh) | 一种地基导航信号观测数据生成方法 | |
CN102736520B (zh) | 一种卫星导航系统原理仿真方法和卫星信号模拟器 | |
Morales et al. | DGPS, RTK-GPS and StarFire DGPS performance under tree shading environments | |
CN101334463B (zh) | Gnss导航信号模拟器中时间参数的模拟方法 | |
CN105372685A (zh) | 一种区域高精度位置增强系统及方法 | |
KR20090036863A (ko) | 소프트웨어 기반 gnss 검증 신호 생성 장치 및 방법 | |
CN101975955B (zh) | 一种gnss模拟器中的通用三维载体运动轨迹生成方法 | |
CN110007326A (zh) | 一种用于星基增强系统的双频测距误差参数生成方法 | |
CN112596077B (zh) | 一种针对终端载体为低轨卫星的卫星导航信号仿真方法 | |
CN103278826A (zh) | 一种北斗b1频点中频信号仿真方法 | |
Ogaja | Introduction to GNSS Geodesy | |
CN111123345B (zh) | 一种基于gnss测量的经验电离层模型数据驱动方法 | |
Du et al. | Experimental study on GPS non-linear least squares positioning algorithm | |
CN102914781A (zh) | 一种格洛纳斯卫星信号星历电文生成方法及装置 | |
JP2004309307A (ja) | 人工衛星シミュレーション装置 | |
CN101419274A (zh) | 电离层延迟误差的获取方法及系统 | |
Tolman et al. | The GPS toolkit | |
Sakai et al. | Augmentation performance of QZSS L1-SAIF signal | |
Al-Franek | Ionospheric imaging for Canadian polar regions | |
Beran et al. | High-accuracy point positioning with low-cost GPS receivers: How good can It get? | |
Tao | Near real-time GPS PPP-inferred water vapor system development and evaluation | |
CN102830410A (zh) | 卫星导航中结合多普勒测速的定位方法 | |
Reckeweg | Integer ambiguity resolution for multi-GNSS and multi-signal raw phase observations | |
Zhang et al. | SiRF InstantFix II Technology | |
CN103389502A (zh) | 基于多个地面基站高精度确定载体加速度的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |