CN109471138B - 一种多模卫星导航信号的合成方法 - Google Patents
一种多模卫星导航信号的合成方法 Download PDFInfo
- Publication number
- CN109471138B CN109471138B CN201910028327.2A CN201910028327A CN109471138B CN 109471138 B CN109471138 B CN 109471138B CN 201910028327 A CN201910028327 A CN 201910028327A CN 109471138 B CN109471138 B CN 109471138B
- Authority
- CN
- China
- Prior art keywords
- satellite
- carrier
- signal
- digital
- synthesizing
- 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
Images
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/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/35—Constructional details or hardware or software details of the signal processing chain
- G01S19/37—Hardware or software details of the signal processing chain
Landscapes
- Engineering & Computer Science (AREA)
- Signal Processing (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
本发明公开了一种多模卫星导航信号的合成方法,包括如下步骤:导航信息时空统一、多载波仿真的信号合成、仿真信号调频、调相和调幅、微信导航信号合成;本发明的有益效果在于,首先,本发明解决多卫星导航系统间的时空差异问题;通过将同一卫星导航系统下的时钟差按照对应模型实时计算,并在卫星伪距计算中予以修正,解决同系统多导航卫星间时间差异问题;其次,本发明采用直接数字频率合成、多载波信号合成技术,实现了对多模卫星信号的合成,实现生成GPS卫星系统,GLONASS卫星系统和BD(北斗)卫星系统及其混合卫星轨道、导航电文即观测数据等信息,并由此生成任意指定位置和速度的卫星导航模拟信号。
Description
技术领域
本发明涉及电子信息技术领域,具体是一种多模卫星导航信号的合成方法。
背景技术
由于无人机成本低廉、操作方便、功能强大、容易获取,易形成“黑飞”。通常这些无人机由于飞行速度较慢,体积小、飞行高度低,导致雷达反射截面积小,不易被发现和拦截。某些无资质、未经审批的个人和组织利用无人机进行飞行活动也易对国家社会、空中航道、重点目标、大型活动等造成严重安全威胁。更为严重的是,一些不法分子可利用无人机携带炸药、相机、毒品等,对诸如监狱、机场、核设施、军用设施基地、边境等重要区域造成极大安全隐患,严重危及国家和社会公共安全。因此,反无人机系统应运而生。目前,反无人机系统主要分为物理打击、强压制干扰和欺骗式诱导三种。
物理空中打击费用较高、准确度低,容易形成碎片,对地面人群可能造成伤害。使用强电磁脉冲对无人机进行压制扰乱,使无人机失控,自动返航,或者坠落,这种强电磁信号极易干扰市民正常生活,另外失控的无人机也易造成二次伤害。最有效的方式是采用全球导航卫星系统(GNSS)欺骗的功能,将无人机欺骗到指定位置,再进行各种针对性处理,不但增强处理的效果,同时有效降低处置带来的民众影响。
但是,由于不同卫星导航系统由不同国家或组织所开发,因而其所使用的空间基准和时间基准也有所不同,同时,由于各系统中导航卫星个体的特性,其卫星时间基准也不统一,而是相对各系统时间基准有一个偏差,因此,在启用GNSS欺骗功能之前需要将多模卫星导航信号进行合成。
发明内容
1.需要解决的技术问题
本发明的目的在于提供一种多模卫星导航信号的合成方法,在短时间内实现对GPS、GLONASS、BD(北斗)三模导航信号的合成。
2.技术方案
为实现上述目的,本发明提供如下技术方案:一种多模卫星导航信号的合成方法,包括如下步骤:
步骤1,导航信号时空统一:确定基准空间坐标系、基准时间轴与各卫星系统的空间坐标系、时间轴的转换关系,将仿真时刻与位置参量按照仿真计算周期实时转换到对应时空基准下;
步骤2,多载波仿真的信号合成:采用直接数字频率合成技术实现多载波仿真信号的合成;
步骤3,仿真信号调频、调相和调幅:将步骤2的仿真信号在仿真程序中完成连续频率调制,在卫星信号合成通道内完成相位调制和幅值调制;
步骤4,卫星导航信号合成:通过公式(1)对卫星导航信号进行合成:
式中,k为离散的采样点,其采样周期为fs,Ts为采样周期,Tc为仿真计算周期;S(k)为第k个采样时刻的信号幅值;A(k)为第k个采样时刻单载波信号幅值,同时由于A(k)是一个慢变过程,因此每个计算周期更新一次,并将其在一个计算周期内视为常值;分别为k采样时刻扩频码、载波相位;分别为k采样时刻的扩频码和导航信息,取值空间为-1或1;分别为第个采样时刻所在计算周期计算所得的码频率和载波频率控制字,表示对计算所得值进行向下取整;分别为扩频码初相位和载波初相位;n(k)为第k个采样时刻的噪声;为第k个采样时刻载体到卫星之间的距离,为第个采样周期时载体到卫星之间的距离;表示在载体到卫星之间的距离为时单载波信号所有采样点幅值的均值,表示在载体到卫星之间的距离为时单载波信号所有采样点幅值的均值进行量化后的数值。
优选的是,步骤1中卫星系统包括GPS卫星系统、GLONSASS卫星系统和BD卫星系统中的一种或多种。
优选的是,步骤2包括如下步骤:
步骤2.1,通过累加器ACC完成多载波频率在时间上积分,得到积分波;
步骤2.2,通过数字波形表ROM完成步骤2.1的积分波由相位到波形的变换,得到数字波形;
步骤2.3,通过数模转换器DAC完成步骤2.2变换后的数字波形到模拟波形的转换,得到模拟波形;
步骤2.4,通过滤波器Filter完成步骤2.3转换后的模拟波形的平滑和滤波,得到多载波的仿真信号。
根据权利要求3所述的一种多模卫星导航信号的合成方法,其特征在于:步骤2采用二进制直接数字频率合成技术实现多载波仿真信号的合成,其中,信号合成的工作时钟fs为100MHz,频率分辨率小于0.05Hz,合成信号寄生功率小于-60dBC。
优选的是,步骤2.2中采用公式(2)进行波形变换:
优选的是,所述累加器ACC和数字波形表ROM之间引入数字加法器,实现数字相位调制。
优选的是,所述数字相位表ROM和数模转换器DAC之间引入数字乘法器,实现数字幅值调制。
3.有益效果
综上所述,本发明的有益效果在于:
(1)本发明解决多卫星导航系统间的时空差异问题;通过将同一卫星导航系统下的时钟差按照对应模型实时计算,并在卫星伪距计算中予以修正,解决同系统多导航卫星间时间差异问题;
(2)本发明采用直接数字频率合成、多载波信号合成技术,实现了对多模卫星信号的合成,实现生成GPS卫星系统,GLONASS卫星系统和BD卫星系统及其混合卫星轨道、导航电文即观测数据等信息,并由此生成任意指定位置和速度的卫星导航模拟信号。
附图说明
图1为一种多模卫星导航信号的合成方法的原理图;
图2为一种多模卫星导航信号的合成方法步骤2的工作原理图。
具体实施方式
下面结合具体实施方式对本专利的技术方案作进一步详细地说明。
请参阅附图1及附图2所示,一种多模卫星导航信号的合成方法,包括如下步骤:
步骤1,导航信号时空统一:确定基准空间坐标系、基准时间轴与各卫星系统的空间坐标系、时间轴的转换关系,将仿真时刻与位置参量按照仿真计算周期实时转换到对应时空基准下;
步骤2,多载波仿真的信号合成:采用直接数字频率合成技术实现多载波仿真信号的合成;
步骤3,仿真信号调频、调相和调幅:将步骤2的仿真信号在仿真程序中完成连续频率调制,在卫星信号合成通道内完成相位调制和幅值调制;
步骤4,卫星导航信号合成:通过公式(1)对卫星导航信号进行合成:
式中,k为离散的采样点,其采样周期为fs,Ts为采样周期,Tc为仿真计算周期;S(k)为第k个采样时刻的信号幅值;A(k)为第k个采样时刻单载波信号幅值,同时由于A(k)是一个慢变过程,因此每个计算周期更新一次,并将其在一个计算周期内视为常值;分别为k采样时刻扩频码、载波相位;分别为k采样时刻的扩频码和导航信息,取值空间为-1或1;分别为第个采样时刻所在计算周期计算所得的码频率和载波频率控制字,表示对计算所得值进行向下取整;分别为扩频码初相位和载波初相位;n(k)为第k个采样时刻的噪声;为第k个采样时刻载体到卫星之间的距离,为第个采样周期时载体到卫星之间的距离;表示在载体到卫星之间的距离为时单载波信号所有采样点幅值的均值,表示在载体到卫星之间的距离为时单载波信号所有采样点幅值的均值进行量化后的数值。
优选的是,步骤1中卫星系统包括GPS卫星系统、GLONSASS卫星系统和BD卫星系统中的一种或多种。
优选的是,步骤2包括如下步骤:
步骤2.1,通过累加器ACC完成多载波频率在时间上积分,得到积分波;
步骤2.2,通过数字波形表ROM完成步骤2.1的积分波由相位到波形的变换,得到数字波形;
步骤2.3,通过数模转换器DAC完成步骤2.2变换后的数字波形到模拟波形的转换,得到模拟波形;
步骤2.4,通过滤波器Filter完成步骤2.3转换后的模拟波形的平滑和滤波,得到多载波的仿真信号。
根据权利要求3所述的一种多模卫星导航信号的合成方法,其特征在于:步骤2采用二进制直接数字频率合成技术实现多载波仿真信号的合成,其中,信号合成的工作时钟fs为100MHz,频率分辨率小于0.05Hz,合成信号寄生功率小于-60dBC。
优选的是,步骤2.2中采用公式(2)进行波形变换:
优选的是,所述累加器ACC和数字波形表ROM之间引入数字加法器,实现数字相位调制。
优选的是,所述数字相位表ROM和数模转换器DAC之间引入数字乘法器,实现数字幅值调制。
以上所述仅为本发明的优选实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书内容所作的等效结构变换,或直接或间接运用附属在其他相关产品的技术领域,均同理包括在本发明的专利保护范围内。
Claims (7)
1.一种多模卫星导航信号的合成方法,其特征在于,包括如下步骤:
步骤1,导航信号时空统一:确定基准空间坐标系、基准时间轴与各卫星系统的空间坐标系、时间轴的转换关系,将仿真时刻与位置参量按照仿真计算周期实时转换到对应时空基准下;
步骤2,多载波仿真的信号合成:采用直接数字频率合成技术实现多载波仿真信号的合成;
步骤3,仿真信号调频、调相和调幅:将步骤2的仿真信号在仿真程序中完成连续频率调制,在卫星信号合成通道内完成相位调制和幅值调制;
步骤4,卫星导航信号合成:通过公式(1)对卫星导航信号进行合成:
式中,k为离散的采样点,其采样周期为fs,Ts为采样周期,Tc为仿真计算周期;S(k)为第k个采样时刻的信号幅值;A(k)为第k个采样时刻单载波信号幅值,同时由于A(k)是一个慢变过程,因此每个计算周期更新一次,并将其在一个计算周期内视为常值;分别为k采样时刻扩频码、载波相位;分别为k采样时刻的扩频码和导航信息,取值空间为-1或1; 分别为第个采样时刻所在计算周期计算所得的码频率和载波频率控制字,表示对计算所得值进行向下取整;分别为扩频码初相位和载波初相位;n(k)为第k个采样时刻的噪声;为第k个采样时刻载体到卫星之间的距离,为第个采样周期时载体到卫星之间的距离;表示在载体到卫星之间的距离为时单载波信号所有采样点幅值的均值,表示在载体到卫星之间的距离为时单载波信号所有采样点幅值的均值进行量化后的数值。
2.根据权利要求1所述的一种多模卫星导航信号的合成方法,其特征在于:步骤1中卫星系统包括GPS卫星系统、GLONSASS卫星系统和BD卫星系统中的一种或多种。
3.根据权利要求1所述的一种多模卫星导航信号的合成方法,其特征在于:步骤2包括如下步骤:
步骤2.1,通过累加器ACC完成多载波频率在时间上积分,得到积分波;
步骤2.2,通过数字波形表ROM完成步骤2.1的积分波由相位到波形的变换,得到数字波形;
步骤2.3,通过数模转换器DAC完成步骤2.2变换后的数字波形到模拟波形的转换,得到模拟波形;
步骤2.4,通过滤波器Filter完成步骤2.3转换后的模拟波形的平滑和滤波,得到多载波的仿真信号。
4.根据权利要求3所述的一种多模卫星导航信号的合成方法,其特征在于:步骤2采用二进制直接数字频率合成技术实现多载波仿真信号的合成,其中,信号合成的工作时钟fs为100MHz,频率分辨率小于0.05Hz,合成信号寄生功率小于-60dBC。
6.根据权利要求4所述的一种多模卫星导航信号的合成方法,其特征在于:所述累加器ACC和数字波形表ROM之间引入数字加法器,实现数字相位调制。
7.根据权利要求3所述的一种多模卫星导航信号的合成方法,其特征在于:所述数字波形表ROM和数模转换器DAC之间引入数字乘法器,实现数字幅值调制。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910028327.2A CN109471138B (zh) | 2019-01-11 | 2019-01-11 | 一种多模卫星导航信号的合成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910028327.2A CN109471138B (zh) | 2019-01-11 | 2019-01-11 | 一种多模卫星导航信号的合成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109471138A CN109471138A (zh) | 2019-03-15 |
CN109471138B true CN109471138B (zh) | 2023-01-31 |
Family
ID=65678527
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910028327.2A Active CN109471138B (zh) | 2019-01-11 | 2019-01-11 | 一种多模卫星导航信号的合成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109471138B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114255626A (zh) * | 2021-12-16 | 2022-03-29 | 商飞软件有限公司 | 一种航空模拟通信仿真系统及其仿真方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6823170B1 (en) * | 2000-07-26 | 2004-11-23 | Ericsson Inc. | Satellite communications system using multiple earth stations |
CN101969690A (zh) * | 2010-10-09 | 2011-02-09 | 东莞市泰斗微电子科技有限公司 | 一种利用导航卫星进行授时的方法及装置 |
US8067738B1 (en) * | 2009-08-25 | 2011-11-29 | Ball Aerospace & Technologies Corp. | Space based calibration transfer spectroradiometer |
CN102841362A (zh) * | 2012-09-24 | 2012-12-26 | 桂林电子科技大学 | 三模卫星信号模拟方法及模拟器 |
US8930706B2 (en) * | 2010-06-30 | 2015-01-06 | Antonio Pujante Cuadrupani | Method, device and network for authenticating the position of a navigation receiver |
CN108073455A (zh) * | 2016-11-11 | 2018-05-25 | 南京航空航天大学 | 一种基于gpu的卫星导航信号并行模拟方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2592077C1 (ru) * | 2015-07-17 | 2016-07-20 | Открытое акционерное общество "Российская корпорация ракетно-космического приборостроения и информационных систем" (ОАО "Российские космические системы") | Способ измерения взаимного положения объектов по сигналам глобальных навигационных спутниковых систем |
-
2019
- 2019-01-11 CN CN201910028327.2A patent/CN109471138B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6823170B1 (en) * | 2000-07-26 | 2004-11-23 | Ericsson Inc. | Satellite communications system using multiple earth stations |
US8067738B1 (en) * | 2009-08-25 | 2011-11-29 | Ball Aerospace & Technologies Corp. | Space based calibration transfer spectroradiometer |
US8930706B2 (en) * | 2010-06-30 | 2015-01-06 | Antonio Pujante Cuadrupani | Method, device and network for authenticating the position of a navigation receiver |
CN101969690A (zh) * | 2010-10-09 | 2011-02-09 | 东莞市泰斗微电子科技有限公司 | 一种利用导航卫星进行授时的方法及装置 |
CN102841362A (zh) * | 2012-09-24 | 2012-12-26 | 桂林电子科技大学 | 三模卫星信号模拟方法及模拟器 |
CN108073455A (zh) * | 2016-11-11 | 2018-05-25 | 南京航空航天大学 | 一种基于gpu的卫星导航信号并行模拟方法 |
Non-Patent Citations (2)
Title |
---|
An experiment of GPS+GLONASS common-view time transfer using new multi-system receivers;J. Nawrocki et.al;《Proceedings of the 20th European Frequency and Time Forum》;20120705;第562-566页 * |
基于强跟踪滤波的多模卫星容错导航算法;赵欣 等;《现代雷达》;20090930;第47-51页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109471138A (zh) | 2019-03-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Gonzalez et al. | Maximum kick from nonspinning black-hole binary inspiral | |
Barr et al. | ELF and VLF radio waves | |
Sahr et al. | The Manastash Ridge radar: A passive bistatic radar for upper atmospheric radio science | |
Hunsucker | Atmospheric gravity waves generated in the high‐latitude ionosphere: A review | |
CN105261265B (zh) | 一种电子对抗技术实验教学系统 | |
Oberheide et al. | The geospace response to variable inputs from the lower atmosphere: a review of the progress made by Task Group 4 of CAWSES-II | |
CN104898107B (zh) | 一种多发多收合成孔径激光雷达信号处理方法 | |
CN101033974A (zh) | 基于软件的导航卫星信号生成系统 | |
CN102141627A (zh) | 一种猝发式的导航信号体制与接收方法 | |
KR101534167B1 (ko) | 실시간 위성항법 재밍 효과도 분석 장치 및 방법 | |
Dizgah et al. | Neutrino properties with ground-based millimeter-wavelength line intensity mapping | |
Chen et al. | Fast and refined processing of radar maneuvering target based on hierarchical detection via sparse fractional representation | |
CN113671537B (zh) | 一种三频信标信号电离层信道仿真方法 | |
CN109471138B (zh) | 一种多模卫星导航信号的合成方法 | |
Dzib et al. | Radio emission from the Be/black hole binary MWC 656 | |
Zhao et al. | The Second Plateau in X-Ray Afterglow Providing Additional Evidence for Rapidly Spinning Magnetars as the GRB Central Engine | |
Pevtsov et al. | Long-term studies of photospheric magnetic fields on the Sun | |
Yang et al. | Characteristics of thunderstorm structure and lightning activity causing negative and positive sprites | |
Yuan et al. | Coordinated investigation of summer time mid‐latitude descending E layer (Es) perturbations using Na lidar, ionosonde, and meteor wind radar observations over Logan, Utah (41.7° N, 111.8° W) | |
CN209446772U (zh) | 一种北斗三代卫星信号模拟系统 | |
CN113985461B (zh) | 一种无人机探测、诱骗系统及方法 | |
CN103954978B (zh) | 基于pxi架构动态导航信号源的实现方法 | |
CN109597100A (zh) | 一种多模卫星导航诱骗反无人机系统 | |
RU151663U1 (ru) | Имитатор радиолокационной обстановки с синтезатором сигналов радиотехнических средств | |
Koval et al. | Morphology of Solar Type II Bursts Caused by Shock Propagation through Turbulent and Inhomogeneous Coronal Plasma |
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 | ||
PE01 | Entry into force of the registration of the contract for pledge of patent right |
Denomination of invention: A Method for Synthesizing Multimode Satellite Navigation Signals Effective date of registration: 20230831 Granted publication date: 20230131 Pledgee: Xi'an Caijin Financing Guarantee Co.,Ltd. Pledgor: XI'AN THUNDER DRIVEN ELECTRONIC TECHNOLOGY CO.,LTD. Registration number: Y2023980054869 |
|
PE01 | Entry into force of the registration of the contract for pledge of patent right |