CN109471138A - 一种多模卫星导航信号的合成方法 - Google Patents

一种多模卫星导航信号的合成方法 Download PDF

Info

Publication number
CN109471138A
CN109471138A CN201910028327.2A CN201910028327A CN109471138A CN 109471138 A CN109471138 A CN 109471138A CN 201910028327 A CN201910028327 A CN 201910028327A CN 109471138 A CN109471138 A CN 109471138A
Authority
CN
China
Prior art keywords
signal
satellite
digital
phase
satellite navigation
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
Application number
CN201910028327.2A
Other languages
English (en)
Other versions
CN109471138B (zh
Inventor
段章山
王徐华
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xi'an Reqing Electronic Technology Co Ltd
Original Assignee
Xi'an Reqing Electronic Technology Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xi'an Reqing Electronic Technology Co Ltd filed Critical Xi'an Reqing Electronic Technology Co Ltd
Priority to CN201910028327.2A priority Critical patent/CN109471138B/zh
Publication of CN109471138A publication Critical patent/CN109471138A/zh
Application granted granted Critical
Publication of CN109471138B publication Critical patent/CN109471138B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/35Constructional details or hardware or software details of the signal processing chain
    • G01S19/37Hardware or software details of the signal processing chain

Abstract

本发明公开了一种多模卫星导航信号的合成方法,包括如下步骤:导航信息时空统一、多载波仿真的信号合成、仿真信号调频、调相和调幅、微信导航信号合成;本发明的有益效果在于,首先,本发明解决多卫星导航系统间的时空差异问题;通过将同一卫星导航系统下的时钟差按照对应模型实时计算,并在卫星伪距计算中予以修正,解决同系统多导航卫星间时间差异问题;其次,本发明采用直接数字频率合成、多载波信号合成技术,实现了对多模卫星信号的合成,实现生成GPS卫星系统,GLONASS卫星系统和BD(北斗)卫星系统及其混合卫星轨道、导航电文即观测数据等信息,并由此生成任意指定位置和速度的卫星导航模拟信号。

Description

一种多模卫星导航信号的合成方法
技术领域
本发明涉及电子信息技术领域,具体是一种多模卫星导航信号的合成方法。
背景技术
由于无人机成本低廉、操作方便、功能强大、容易获取,易形成“黑飞”。通常这些无人机由于飞行速度较慢,体积小、飞行高度低,导致雷达反射截面积小,不易被发现和拦截。某些无资质、未经审批的个人和组织利用无人机进行飞行活动也易对国家社会、空中航道、重点目标、大型活动等造成严重安全威胁。更为严重的是,一些不法分子可利用无人机携带炸药、相机、毒品等,对诸如监狱、机场、核设施、军用设施基地、边境等重要区域造成极大安全隐患,严重危及国家和社会公共安全。因此,反无人机系统应运而生。目前,反无人机系统主要分为物理打击、强压制干扰和欺骗式诱导三种。
物理空中打击费用较高、准确度低,容易形成碎片,对地面人群可能造成伤害。使用强电磁脉冲对无人机进行压制扰乱,使无人机失控,自动返航,或者坠落,这种强电磁信号极易干扰市民正常生活,另外失控的无人机也易造成二次伤害。最有效的方式是采用全球导航卫星系统(GNSS)欺骗的功能,将无人机欺骗到指定位置,再进行各种针对性处理,不但增强处理的效果,同时有效降低处置带来的民众影响。
但是,由于不同卫星导航系统由不同国家或组织所开发,因而其所使用的空间基准和时间基准也有所不同,同时,由于各系统中导航卫星个体的特性,其卫星时间基准也不统一,而是相对各系统时间基准有一个偏差,因此,在启用GNSS欺骗功能之前需要将多模卫星导航信号进行合成。
发明内容
1.需要解决的技术问题
本发明的目的在于提供一种多模卫星导航信号的合成方法,在短时间内实现对GPS、GLONASS、BD(北斗)三模导航信号的合成。
2.技术方案
为实现上述目的,本发明提供如下技术方案:一种多模卫星导航信号的合成方法,包括如下步骤:
步骤1,导航信号时空统一:确定基准空间坐标系、基准时间轴与各卫星系统的空间坐标系、时间轴的转换关系,将仿真时刻与位置参量按照仿真计算周期实时转换到对应时空基准下;
步骤2,多载波仿真的信号合成:采用直接数字频率合成技术实现多载波仿真信号的合成;
步骤3,仿真信号调频、调相和调幅:将步骤2的仿真信号在仿真程序中完成连续频率调制,在卫星信号合成通道内完成相位调制和幅值调制;
步骤4,卫星导航信号合成:通过公式(1)对卫星导航信号进行合成:
式中,k为离散的采样点(采样周期为fs),Ts为采样周期,Tc为仿真计算周期;A(k)为该路信号幅值,实际上是仿真载体坐标系下载体到卫星极矢量的函数,同时由于A(k)是一个慢变过程,因此每个计算周期更新一次,并将其在一个计算周期内视为常值;分别为k采样时刻扩频码、载波相位;分别为k采样时刻的扩频码和导航信息,取值空间为-1或1;FCWPRN(n)、FCWCarrier(n)分别为k采样时刻所在计算周期计算所得的DDS频率控制控制字;分别为扩频码初相位和载波初相位。
优选的是,步骤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)进行波形变换:
其中,bweight为累加器字位权重,AI,AQ分别为两路正交波形I、Q的幅值,设ROM表输入相位为 为波形变换后的相位,ROM输入相位量化位数WP,ROM输出位数为WQ
优选的是,所述累加器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为仿真计算周期;A(k)为该路信号幅值,实际上是仿真载体坐标系下载体到卫星极矢量的函数,同时由于A(k)是一个慢变过程,因此每个计算周期更新一次,并将其在一个计算周期内视为常值;分别为k采样时刻扩频码、载波相位;分别为k采样时刻的扩频码和导航信息,取值空间为-1或1;FCWPRN(n)、FCWCarrier(n)分别为k采样时刻所在计算周期计算所得的DDS频率控制控制字;分别为扩频码初相位和载波初相位。
优选的是,步骤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)进行波形变换:
其中,bweight为累加器字位权重,AI,AQ分别为两路正交波形I、Q的幅值,设ROM表输入相位为 为波形变换后的相位,ROM输入相位量化位数WP,ROM输出位数为WQ
优选的是,所述累加器ACC和数字波形表ROM之间引入数字加法器,实现数字相位调制。
优选的是,所述数字相位表ROM和数模转换器DAC之间引入数字乘法器,实现数字幅值调制。
以上所述仅为本发明的优选实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书内容所作的等效结构变换,或直接或间接运用附属在其他相关产品的技术领域,均同理包括在本发明的专利保护范围内。

Claims (7)

1.一种多模卫星导航信号的合成方法,其特征在于,包括如下步骤:
步骤1,导航信号时空统一:确定基准空间坐标系、基准时间轴与各卫星系统的空间坐标系、时间轴的转换关系,将仿真时刻与位置参量按照仿真计算周期实时转换到对应时空基准下;
步骤2,多载波仿真的信号合成:采用直接数字频率合成技术实现多载波仿真信号的合成;
步骤3,仿真信号调频、调相和调幅:将步骤2的仿真信号在仿真程序中完成连续频率调制,在卫星信号合成通道内完成相位调制和幅值调制;
步骤4,卫星导航信号合成:通过公式(1)对卫星导航信号进行合成:
式中,k为离散的采样点(采样周期为fs),Ts为采样周期,Tc为仿真计算周期;A(k)为该路信号幅值,实际上是仿真载体坐标系下载体到卫星极矢量的函数,同时由于A(k)是一个慢变过程,因此每个计算周期更新一次,并将其在一个计算周期内视为常值;分别为k采样时刻扩频码、载波相位;分别为k采样时刻的扩频码和导航信息,取值空间为-1或1;FCWPRN(n)、FCWCarrier(n)分别为k采样时刻所在计算周期计算所得的DDS频率控制控制字;分别为扩频码初相位和载波初相位。
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。
5.根据权利要求3所述的一种多模卫星导航信号的合成方法,其特征在于:步骤2.2中采用公式(2)进行波形变换:
其中,bweight为累加器字位权重,AI,AQ分别为两路正交波形I、Q的幅值,设ROM表输入相位为为波形变换后的相位,ROM输入相位量化位数WP,ROM输出位数为WQ
6.根据权利要求4所述的一种多模卫星导航信号的合成方法,其特征在于:所述累加器ACC和数字波形表ROM之间引入数字加法器,实现数字相位调制。
7.根据权利要求2所述的一种多模卫星导航信号的合成方法,其特征在于:所述数字相位表ROM和数模转换器DAC之间引入数字乘法器,实现数字幅值调制。
CN201910028327.2A 2019-01-11 2019-01-11 一种多模卫星导航信号的合成方法 Active CN109471138B (zh)

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 true CN109471138A (zh) 2019-03-15
CN109471138B 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)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114255626A (zh) * 2021-12-16 2022-03-29 商飞软件有限公司 一种航空模拟通信仿真系统及其仿真方法

Citations (7)

* Cited by examiner, † Cited by third party
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
RU2592077C1 (ru) * 2015-07-17 2016-07-20 Открытое акционерное общество "Российская корпорация ракетно-космического приборостроения и информационных систем" (ОАО "Российские космические системы") Способ измерения взаимного положения объектов по сигналам глобальных навигационных спутниковых систем
CN108073455A (zh) * 2016-11-11 2018-05-25 南京航空航天大学 一种基于gpu的卫星导航信号并行模拟方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
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 桂林电子科技大学 三模卫星信号模拟方法及模拟器
RU2592077C1 (ru) * 2015-07-17 2016-07-20 Открытое акционерное общество "Российская корпорация ракетно-космического приборостроения и информационных систем" (ОАО "Российские космические системы") Способ измерения взаимного положения объектов по сигналам глобальных навигационных спутниковых систем
CN108073455A (zh) * 2016-11-11 2018-05-25 南京航空航天大学 一种基于gpu的卫星导航信号并行模拟方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
J. NAWROCKI ET.AL: "An experiment of GPS+GLONASS common-view time transfer using new multi-system receivers", 《PROCEEDINGS OF THE 20TH EUROPEAN FREQUENCY AND TIME FORUM》 *
赵欣 等: "基于强跟踪滤波的多模卫星容错导航算法", 《现代雷达》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114255626A (zh) * 2021-12-16 2022-03-29 商飞软件有限公司 一种航空模拟通信仿真系统及其仿真方法

Also Published As

Publication number Publication date
CN109471138B (zh) 2023-01-31

Similar Documents

Publication Publication Date Title
EP2816371B1 (en) Method and device for steering attitude of satellite carrying synthetic aperture radar
US10803210B2 (en) Real-time electromagnetic environmental simulator
CN102937713B (zh) 一种可重构的卫星导航射频信号模拟方法及装置
CN104777464B (zh) 利用ads‑b数据实现多雷达数据半实物仿真的方法
EP2759847A1 (en) Method and apparatus for determining equivalent velocity
CN110910511B (zh) 一种基于混合现实的搜救流程仿真与交互系统及方法
CN102226841A (zh) 基于高阶多项式距离方程的同步轨道sar成像方法
CN109683147A (zh) 杂乱脉冲流信号的实时生成方法、装置及电子设备
CN102176029A (zh) 一种gps直达与多径信号模拟器及模拟方法
Chen et al. Hybrid N-inception-LSTM-based aircraft coordinate prediction method for secure air traffic
CN103105614A (zh) 基于惯导辅助的空时域联合抗干扰方法
CN109471138A (zh) 一种多模卫星导航信号的合成方法
Burczyk et al. Voice multilateration system
Wang et al. Impacts of gps spoofing on path planning of unmanned surface ships
CN102478653A (zh) 一种基于距离分割的sar回波时频混合模拟方法
CN104155652B (zh) 一种地表冻融状态的监测方法及装置
Wang et al. Approach for high‐resolution inverse synthetic aperture radar imaging of ship target with complex motion
CN105607054A (zh) 收发一体式超宽带飞行雷达定时采样探测方法
CN105891798A (zh) 变prf条件下雷达目标微动特征提取的方法
CN109597100A (zh) 一种多模卫星导航诱骗反无人机系统
Ofek Space and velocity distributions of Galactic isolated old Neutron stars
CN109471135A (zh) 一种卫星信号的真实模拟方法
CN207799083U (zh) 一种反制无人机用干扰系统
CN113985461B (zh) 一种无人机探测、诱骗系统及方法
CN109471133A (zh) 一种无缝侵入无人机导航系统的欺骗技术

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
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