CN110109156A - 一种卫星导航多路复用调制下boc信号数字畸变分析方法 - Google Patents
一种卫星导航多路复用调制下boc信号数字畸变分析方法 Download PDFInfo
- Publication number
- CN110109156A CN110109156A CN201910288217.XA CN201910288217A CN110109156A CN 110109156 A CN110109156 A CN 110109156A CN 201910288217 A CN201910288217 A CN 201910288217A CN 110109156 A CN110109156 A CN 110109156A
- Authority
- CN
- China
- Prior art keywords
- signal
- satellite
- distortion
- spectral density
- power spectral
- 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
- 238000004458 analytical method Methods 0.000 title claims abstract description 15
- 230000003595 spectral effect Effects 0.000 claims abstract description 23
- 238000000034 method Methods 0.000 claims abstract description 12
- 238000005311 autocorrelation function Methods 0.000 claims description 8
- 230000000630 rising effect Effects 0.000 claims description 8
- 230000009466 transformation Effects 0.000 claims description 7
- 238000007476 Maximum Likelihood Methods 0.000 claims description 4
- 230000002159 abnormal effect Effects 0.000 claims description 2
- 238000001228 spectrum Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 5
- 238000013461 design Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 108010076504 Protein Sorting Signals Proteins 0.000 description 1
- 125000002015 acyclic group Chemical group 0.000 description 1
- 238000013475 authorization Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 229910002056 binary alloy Inorganic materials 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 210000003739 neck Anatomy 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 239000000546 pharmaceutical excipient Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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/23—Testing, monitoring, correcting or calibrating of receiver elements
-
- 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/24—Acquisition or tracking or demodulation of signals transmitted by the system
- G01S19/30—Acquisition or tracking or demodulation of signals transmitted by the system code related
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Monitoring And Testing Of Transmission In General (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种卫星导航多路复用调制下BOC信号数字畸变分析方法，包括以下步骤：S100：获取卫星信号；S200：根据卫星信号求解出卫星信号的功率谱密度函数；S300：根据卫星信号的功率谱密度函数表达式，绘制出卫星信号的功率谱密度图；S400：根据卫星信号的功率谱密度图上是否出现尖峰，判断卫星信号是否发生数字畸变；S500：若卫星信号发生畸变，采用过零点计算方法精确计算出卫星信号的数字畸变量。本发明有效快捷的发现信号中发生畸变情况，并且及时精确的计算出畸变量。
Description
技术领域
本发明涉及卫星导航信号的数字畸变分析方法，尤其涉及一种卫 星导航多路复用调制下BOC信号数字畸变分析方法。
背景技术
卫星导航从基础研究领域(天文学、力学、物理学、地球动力学 等)到工程技术领域(信息传递、深空探测、空间飞行器、时间传递、 测速授时等)，以及关系国民经济建设和国家安全的诸多重要部门和 领域(海事、交通、救援、精准农业、地震监测、电子通讯等)等各个 方面，均已占据重要地位，发挥着重要作用。可以说，卫星导航技 术是各种基础理论、最新科技成果、国家科技水平和综合实力的集中 体现。全球卫星导航***(GlobalNavigation Satellite System,GNSS)以其高精度、实时性、不受气候和地域等条件限制的特 点，已发展成为在陆地交通、航海和航空中的通用导航工具，成为 应用最广泛的导航技术。世界各发达国家长期以来一直致力于卫星 导航领域的研究工作，以确保向用户提供更加精确的实时位置、速度 和时间等信息。
卫星导航信号是卫星导航***与接收机的唯一接口，完善的信号 体制设计是终端设备开发和芯片研制的首要输入条件。可以说，信号 的潜在性能决定了整个卫星导航***的性能极限。如果卫星发射的信 号本身存在缺陷，即使***其它环节的设计再优，整个***在定位、 测速、授时性能、抗干扰能力、与其它***的兼容与互操作性等关键 性能上仍然存在难以弥补的不足。
然而在现有技术中，新型GNSS导航信号虽然有着众多优点，但 同时也增加了信号设计及实现过程中的复杂性。信号的各类微小异 常，都有可能对GNSS的PVT性能产生致命影响。导航信号在星上产 生到信号发射整个过程中，任何故障或异常都会可能使得待发射信号 波形产生畸变。最终会导致接收机跟踪过程中相关峰曲线相比理想情 况下出现异常，从而对不同跟踪环路和不同鉴相器的接收机产生不同 程度的跟踪误差和测距误差。而现有技术中，对于卫星导航多路复用 调制下BOC信号的分析中，大多不能有效快捷的发现信号中发生畸变 情况，并且不能及时精确的计算出畸变量。
因此，有必要提供一种卫星导航多路复用调制下BOC信号数字畸 变分析方法，以解决上述问题。
发明内容
为解决上述技术问题，本发明提供一种卫星导航多路复用调制下 BOC信号畸变分析方法。
本发明提供一种卫星导航多路复用调制下BOC信号数字畸变分 析方法，包括以下步骤：
S100：获取卫星信号；
S200：根据卫星信号求解出卫星信号的功率谱密度函数；
S300：根据卫星信号的功率谱密度函数表达式，绘制出卫星信号 的功率谱密度图；
S400：根据卫星信号的功率谱密度图上是否出现尖峰，判断卫星 信号是否发生数字畸变；
S500：若卫星信号发生畸变，采用过零点计算方法或极大似然发 精确计算出卫星信号的数字畸变量。
进一步的，在S200中，所述根据卫星信号求解出卫星信号的功 率谱密度函数具体步骤如下：
S201：求得畸变信号的时域表达式；
S202：求得信号自相关函数；
S203：对信号自相关函数进行傅里叶变换得到卫星信号的功率谱 密度函数。
进一步的，在S500中，所述采用过零点计算方法精确计算出卫 星信号的数字畸变量的步骤如下：
S501：求得以采样点为单位的上升沿过零点的表达式；
S502：将码片的下降沿过零点减去上升沿过零点可得信号单个正 码片数字失真数值。
与相关技术相比，本发明具有如下有益效果：
在现有技术中，对于卫星导航多路复用调制下BOC信号的分析 中，大多不能有效快捷的发现信号中发生畸变情况，并且不能及时精 确的计算出畸变量。而本发明可以有效快捷的发现信号中发生畸变情 况，并且及时可以精确的计算出畸变量。
具体实施方式
下面将结合实施方式对本发明作进一步说明。
本发明提供一种卫星导航多路复用调制下BOC信号畸变分析方 法，包括以下步骤：
S100：获取多路复用调制下BOC信号；
S200：根据多路复用调制下BOC信号求解出多路复用调制下BOC 信号的功率谱密度函数；
BOC信号发生数字畸变时，存在三种不同的情形，由于现阶段 GNSS***采用BOC调制的信号的典型代表为GPS L1 M信号，因 此本文仍以BOC(10,5)信号为分析对象，详细分析在不同的数字畸变 原理下BOC信号的功率谱密度解析式。BOC信号的一般表达式为：
式中，ak为伪码信号序列，为伪码码片赋型信号，在BOC 调制信号中为矩形脉冲信号，Tc为伪码码宽，χ(t)为副载波信号，它 通常是个周期信号。
S202：求得信号自相关函数；
对sx(t)信号求解自相关函数Rx(τ)：
上式中Ra(k)为伪码序列ak的自相关函数。
S203：对信号自相关函数进行傅里叶变换得到多路复用调制下 BOC信号的功率谱密度函数：
对上式进行傅里叶变换，得到sx(t)的功率谱密度Sx(f)如下：
在此，我们将扩频符号的傅里叶变换记为副载波符号 χ(l)的傅里叶变换记为X(f)，伪码序列ak的傅里叶变换记为Sa(f)， 将上述三个傅里叶表达式代入上式可得到BOC信号的功率谱密度解 析式为：
假定二进制伪码序列ak呈独立同分布，那么这个无限长的理想伪 码序列具有如下的自相关特性：
那么其功率谱密度Sa(f)为
因此BOC信号功率谱表达式可进一步表达为：
S300：根据多路复用调制下BOC信号的功率谱密度函数表达式， 绘制出多路复用调制下BOC信号的功率谱密度图；
S400：根据多路复用调制下BOC信号的功率谱密度图上是否出现 尖峰，判断多路复用调制下BOC信号是否发生信号畸变，发生畸变后 的信号会出现尖峰谱现象，在零频率、主瓣和旁瓣以及旁瓣之间出现 了尖峰，而主瓣中心频率处没有出现尖峰。
S500：若多路复用调制下BOC信号发生畸变，采用过零点计算方 法或极大似然发精确计算出多路复用调制下BOC信号的数字畸变量。
在S500中，所述采用过零点计算方法精确计算出多路复用调制 下BOC信号的数字畸变量的步骤如下：
S501：求得以采样点为单位的上升沿过零点的表达式：
假设信号正码片发生了数字失真，且失真量Δ为正，信号的波形 中的正码片下降沿过零点发生偏移，正码片宽度增大，相邻的负码片 宽度减小了Δ。因而测量信号波形中码片的过零点之差可以获得正码 片和负码片宽度，码片的过零点可以利用边沿上符号发生改变的两个 点(x0,y0)和(x1,y1)线性拟合得到，以采样点为单位的上升沿过零点xrising求解表达式如下：
S502：将码片的下降沿过零点减去上升沿过零点可得信号单个正 码片数字失真数值：
信号的下降沿过零点xfailing也可通过线性拟合求解，将码片的下 降沿过零点减去上升沿过零点可得信号单个正码片数字失真数值：
其中，fs为信号的采样率，round为四舍五入取整运算。当数字 失真量的绝对值大于Tc/2时，信号功率谱将会产生十分严重的伪码时 钟泄露现象，信号质量严重下降，在目前GNSS***中未发现该种现 象，因而本文不作考虑。需要注意的是，基带信号已调制了导航电文， 必须根据导航电文符号分开计算数字失真。同时，地面接收信号受噪 声影响，信号实测的过零点存在抖动现象，数字失真量近似呈正态分 布，本文通过平均后基带信号信噪比远高于原始信号，通过测量累加 平均后的基带信号减少噪声对数字失真量测量的影响，对测量结果求 平均，可获得最终精确的数字失真量。
同时，求解多路复用调制下具有非周期性的授权信号的BOC信号 的数字畸变量可以采用极大似然法。由于目前只有GPS***存在数字 畸变，卫星载荷采用了差分传输技术，单个支路信号发生了数字畸变， 因而数字畸变参数仅为2个，即：
上式中ηI和ηQ为I/Q支路信号的数字畸变量，y为接收信号，X为 理想输入信号，为通道传输函数估计值。在该算法中，数字畸变量 的分辨率决定信号的采样率，可设置数字畸变量的分辨率为0.1ns， 因而数据采样率应为10GHz，上式中通过简单的二维搜索进行求解。
以上所述仅为本发明的实施例，并非因此限制本发明的专利范 围，凡是利用本发明说明书内容所作的等效结构或等效流程变换，或 直接或间接运用在其它相关的技术领域，均同理包括在本发明的专利 保护范围内。
Claims (3)
1.一种卫星导航多路复用调制下BOC信号数字畸变分析方法，其特征在于包括以下步骤：
S100：通过高增益天线获取高信噪比卫星下行导航信号；
S200：根据卫星信号求解出卫星信号的功率谱密度函数；
S300：根据卫星信号的功率谱密度函数表达式，绘制出卫星信号的功率谱密度图；
S400：根据卫星信号的功率谱密度图上是否出现尖峰，判断卫星信号是否发生数字畸变；
S500：若卫星信号发生畸变，采用过零点计算方法或极大似然法精确计算出卫星信号的数字畸变量。
2.如权利要求1所述的一种卫星导航多路复用调制下BOC信号畸变分析方法，其特征在于：在S200中，所述根据卫星信号求解出卫星信号的功率谱密度函数具体步骤如下：
S201：求得畸变信号的时域表达式；
S202：求得信号自相关函数；
S203：对信号自相关函数进行傅里叶变换得到卫星信号的功率谱密度函数。
3.如权利要求1所述的一种卫星导航多路复用调制下BOC信号畸变分析方法，其特征在于：在S500中，所述采用过零点计算方法精确计算出卫星信号的数字畸变量的步骤如下：
S501：求得以采样点为单位的上升沿过零点的表达式；
S502：将码片的下降沿过零点减去上升沿过零点可得信号单个正码片数字失真数值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910288217.XA CN110109156B (zh) | 2019-04-11 | 2019-04-11 | 一种卫星导航多路复用调制下boc信号数字畸变分析方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910288217.XA CN110109156B (zh) | 2019-04-11 | 2019-04-11 | 一种卫星导航多路复用调制下boc信号数字畸变分析方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110109156A true CN110109156A (zh) | 2019-08-09 |
| CN110109156B CN110109156B (zh) | 2023-11-21 |
Family
ID=67485292
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910288217.XA Active CN110109156B (zh) | 2019-04-11 | 2019-04-11 | 一种卫星导航多路复用调制下boc信号数字畸变分析方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110109156B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110824511A (zh) * | 2019-10-30 | 2020-02-21 | 北京航空航天大学 | 一种小型化的卫星导航信号码片畸变监测***及监测方法 |
| CN114337870A (zh) * | 2021-12-02 | 2022-04-12 | 杭州北斗时空研究院 | 一种基于Bayes准则的BOC信号质量监测方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102901973A (zh) * | 2012-09-25 | 2013-01-30 | 郑州威科姆科技股份有限公司 | 基于北斗卫星的信号实时并行快速捕获方法 |
| EP2708917A1 (fr) * | 2012-09-18 | 2014-03-19 | Centre National d'Etudes Spatiales | Authentification de signaux GNSS |
| CN109100758A (zh) * | 2018-09-26 | 2018-12-28 | 中国科学院国家授时中心 | 新一代卫星导航信号的波形畸变特性评估方法 |
-
2019
- 2019-04-11 CN CN201910288217.XA patent/CN110109156B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2708917A1 (fr) * | 2012-09-18 | 2014-03-19 | Centre National d'Etudes Spatiales | Authentification de signaux GNSS |
| CN102901973A (zh) * | 2012-09-25 | 2013-01-30 | 郑州威科姆科技股份有限公司 | 基于北斗卫星的信号实时并行快速捕获方法 |
| CN109100758A (zh) * | 2018-09-26 | 2018-12-28 | 中国科学院国家授时中心 | 新一代卫星导航信号的波形畸变特性评估方法 |
Non-Patent Citations (5)
| Title |
|---|
| 刘建成等: "卫星导航信号畸变导致的测距偏差的估计方法", 《宇航学报》 * |
| 刘瑞华等: "北斗卫星导航***空间信号数字畸变的研究", 《航天控制》 * |
| 康立等: "北斗***导航信号标称失真研究", 《电子学报》 * |
| 杨德进等: "Galileo E1导航信号质量分析与评估", 《时间频率学报》 * |
| 熊定喜等: "Galileo卫星星钟故障前后信号质量分析", 《飞行器测控学报》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110824511A (zh) * | 2019-10-30 | 2020-02-21 | 北京航空航天大学 | 一种小型化的卫星导航信号码片畸变监测***及监测方法 |
| CN114337870A (zh) * | 2021-12-02 | 2022-04-12 | 杭州北斗时空研究院 | 一种基于Bayes准则的BOC信号质量监测方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110109156B (zh) | 2023-11-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Motella et al. | Method for assessing the interference impact on GNSS receivers | |
| CN103033824B (zh) | 高性能导航卫星空间信号质量评估方法 | |
| JPH02105746A (ja) | デイジタル復調装置及び差動位相偏移キーイング復調装置ならびに低信号対雑音比入力信号復調方法 | |
| Hein et al. | A candidate for the Galileo L1 OS optimized signal | |
| US9581700B2 (en) | Method and apparatus tracking global navigation satellite system (GNSS) | |
| CN110109156A (zh) | 一种卫星导航多路复用调制下boc信号数字畸变分析方法 | |
| US10018731B2 (en) | Processor for a radio receiver | |
| Sombrin | On the formal identity of EVM and NPR measurement methods: Conditions for identity of error vector magnitude and noise power ratio | |
| CN102508263B (zh) | 一种采用多通道并行相关峰采样监测信号质量的方法 | |
| CN102571652A (zh) | 一种gnss基带信号的评估方法 | |
| CN112213742A (zh) | 一种卫星导航***信号质量监测方法 | |
| Elango et al. | Design of complete software GPS signal simulator with low complexity and precise multipath channel model | |
| Thoelert et al. | Characterization of nominal signal distortions and impact on receiver performance for GPS (IIF) L1/L5 and Galileo (IOV) E1/E5a signals | |
| CN105044734B (zh) | 一种高精度载噪比估算方法 | |
| Fan et al. | Discrete chirp‐Fourier transform‐based acquisition algorithm for weak global positioning system L5 signals in high dynamic environments | |
| Rouabah et al. | Unambiguous multipath mitigation technique for BOC (n, n) and MBOC‐modulated GNSS signals | |
| Fantino et al. | Performance analysis of MBOC, AltBOC and BOC modulations in terms of multipath effects on the carrier tracking loop within GNSS receivers | |
| CN109100758A (zh) | 新一代卫星导航信号的波形畸变特性评估方法 | |
| CN104199064B (zh) | 基于相关波动的boc信号参数盲估计方法 | |
| Pagot et al. | Threat models design for new GNSS signals | |
| CN106291612A (zh) | 一种导航卫星星间链路无线信号高性能捕获判决方法 | |
| CN108957492A (zh) | 一种gps的l1c/a和l1c联合捕获方法及装置 | |
| CN104007425A (zh) | 一种星间时差测量方法及*** | |
| Lo et al. | Loran data modulation: A primer [AESS tutorial IV] | |
| Motella et al. | An interference impact assessment model for GNSS signals |
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 | ||
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information |
Inventor after: Rao Yongnan Inventor after: Wang Xue Inventor after: Lu Xiaochun Inventor after: Kang Li Inventor before: Kang Li Inventor before: Rao Yongnan Inventor before: Wang Xue Inventor before: Lu Xiaochun |