CN102937717B - 一种优选导航卫星扩频码的方法 - Google Patents
一种优选导航卫星扩频码的方法 Download PDFInfo
- Publication number
- CN102937717B CN102937717B CN201210453838.7A CN201210453838A CN102937717B CN 102937717 B CN102937717 B CN 102937717B CN 201210453838 A CN201210453838 A CN 201210453838A CN 102937717 B CN102937717 B CN 102937717B
- Authority
- CN
- China
- Prior art keywords
- spreading code
- value
- calculated value
- maximum
- code
- 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
Landscapes
- Mobile Radio Communication Systems (AREA)
Abstract
本发明针对卫星导航系统提出了一种优选导航卫星扩频码的方法。本方法计算获得扩频码组中各个扩频码的码平衡性、最大自相关旁峰值、最大互相关峰值、功率谱最大值,并依据这四个指标对扩频码进行排序和加权计算值的计算,根据加权计算值选取扩频码组中抗干扰性能较好的扩频码。采用本方法可在不改变系统软硬件设计方案和不增加系统建设成本的条件下,通过扩频码的优选有效提高了卫星导航系统的抗干扰性能性能。
Description
技术领域
本发明涉及卫星导航技术领域,特别是一种优选导航卫星扩频码的方法。
背景技术
在卫星导航系统中,导航卫星发射调制有扩频码的射频信号。地面接收机接收此射频信号,经过变频、解调获得卫星扩频码,并使用本地扩频码与卫星扩频码进行码相关计算,依据相关计算的结果获得测距结果。扩频码的优劣直接决定导航卫星的测距精度。另外,扩频码对导航接收机的抗干扰能力也有重要影响。因此,扩频码的选取是导航卫星系统设计的重要工作之一。
目前扩频码的选取主要为对扩频码码组的选取。一个码组是指由相同生成多项式产生,但相位选择器不同的移位寄存器产生的扩频码所组成的集合。由于属于同一码组的各个扩频码在码平衡性、最大自相关旁峰、最大互相关峰值和功率谱最大值方面的差别,造成扩频码在测距性能、抗干扰能力存在一定的差别。现有方法一般采取在同一码组中随机选取扩频码的方式,忽略了同一码组中扩频码的差异,不能优选出性能好的扩频码,导致卫星导航系统性能一定程度的降低。
发明内容
本发明的技术解决问题是:克服现有技术的不足,提供了一种优选导航卫星扩频码的方法,该方法在不改变卫星导航系统软硬件设计方案和不增加系统建设成本的条件下,能够选取出扩频码组中抗干扰性能较好的扩频码。
本发明的技术解决方案是:一种优选导航卫星扩频码的方法,其步骤如下:
(3)将步骤(2)中得到的码平衡性Di的值从小到大排序,数值相同的Di划入同一集合形成M1个集合,其中1≤j≤M1;数值最小的Di对应的集合为数值最大的Di对应的集合为依照排序结果得出码平衡性Di的扩频码计算值,对应的扩频码计算值为1,对应的扩频码计算值为2,依次类推,对应的扩频码计算值为M1;所有扩频码计算值构成集合其中为第i个扩频码的计算值;
(6)将步骤(5)中得出的最大自相关函数旁峰值从小到大排序,数值相同的划入同一集合形成M2个集合,其中1≤j≤M2,数值最小的对应的集合为数值最大的对应的集合为依照排序结果得出最大自相关函数旁峰值的扩频码计算值,对应的扩频码计算值为1,对应的扩频码计算值为2,依次类推,对应的扩频码计算值为M2;所有扩频码的计算值构成集合其中为第i个扩频码的计算值;
(9)将步骤(8)中获得的最大互相关峰值从小到大排序,数值相同的划入同一集合形成M3个集合,其中1≤j≤M3,数值最小的对应的集合为数值最大的对应的集合为依照排序结果得出最大互相关峰值的扩频码计算值,对应的扩频码计算值为1,对应的扩频码计算值为2,依次类推,对应的扩频码计算值为M3。所有扩频码的计算值构成集合其中为第i个扩频码的计算值;
(12)将步骤(11)中获得的功率谱最大值从小到大排序,数值相同的划入同一集合形成M4个集合,其中1≤j≤M4,数值最小的对应的集合为数值最大的对应的集合为依照排序结果得出功率谱最大值的扩频码计算值,对应的扩频码计算值为1,对应的扩频码计算值为2,依次类推,对应的扩频码计算值为M4;所有扩频码的计算值构成集合 为第i个扩频码的计算值;
(13)分别给码平衡性Di、最大自相关旁峰值最大互相关峰值功率谱最大值的扩频码计算值集合S1、S2、S3和S4赋予权值w1,w2,w3,w4,权值w1,w2,w3,w4应满足0≤w1≤1,0≤w2≤1,0≤w3≤1和0≤w4≤1,且w1+w2+w3+w4=1;对于带宽、发射功率受限的系统取w1=0.4,w2=w3=w4=0.2;对于要求接收机假锁概率的系统取w2=0.4,w1=w3=w4=0.2;对于要求抗多址干扰能力的系统取w3=0.4,w1=w2=w4=0.2;对于要求抗窄带干扰能力的系统w4=0.4,w1=w2=w3=0.2;其余系统取w1=w2=w3=w4=0.25;
(14)根据权值和扩频码计算值集合S1、S2、S3和S4,计算各个扩频码的加权计算值S,即S=w1·S1+w2·S2+w3·S3+w4·S4;加权计算值S值最小的扩频码为最优。
本发明与现有技术相比的有益效果是:
给出了一种优选导航卫星扩频码的方法,该方法通过分别计算扩频码组中各个扩频码的码平衡性、自相关旁峰值、互相关峰值和功率谱最大值,并依据这四个指标对扩频码进行排序和计算加权计算值,根据加权计算值选取扩频码组中抗干扰性能较好的扩频码,在不改变卫星导航系统软硬件设计方案和不增加系统建设成本的条件下,提高了卫星导航系统的抗干扰性能性能。
附图说明
图1为本发明方法的流程图;
图2为循环移位示意图。
具体实施方式
码平衡性好的扩频码在传输时占用带宽少,所需发射功率低。使用最大自相关旁峰小的扩频码有利于降低接收机假锁的概率。最大互相关峰值小的扩频码抗多址干扰的能力较强。功率谱最大值较小的扩频码抗窄带干扰的能力较强。因此,根据码平衡性好、最大自相关旁峰、最大互相关峰值和功率谱最大值这四个指标有减少降低系统占用的带宽和功率,降低接收机假锁的概率,提高抗干扰能力,综合提升卫星导航系统的性能。本发明方法主要包括以下具体步骤:
(1)对于包含N个扩频码的码组,计算第i个扩频码中数值等于1的码片个数1≤i≤N。由于码片的值只能等于0或者等于1,因此将所有L个码片的值相加即获得扩频码中数值等于1的码片个数计算第i个扩频码中数值等于0的个数使用L减去即为
(3)将码平衡性Di从小到大排序,数值相同的Di划入同一集合形成M1个集合,其中1≤j≤M1,最小值对应的集合为最大值对应的集合为依照排序结果得出码平衡性Di的扩频码计算值,对应的扩频码计算值为1,对应的扩频码计算值2,依次类推,对应的扩频码为M1。所有扩频码的计算值构成集合S1, 为第i个扩频码的计算值。
(4)依次计算码组中第i个扩频码的循环移位自相关值1≤l≤L,Ci(l)为第i个扩频码的第l个码片,k为码片循环移位的偏移量,0≤k≤L-1,Ci(l+k)为第i个扩频码的第l+k个码片;循环移位的具体实现方法如图1所示:两个长度为L的扩频码在进行相关运算时需进行循环移位操作。左图表示偏移量k等于0时两个扩频码码片的位置关系,其中外圈表示第1个扩频码,内圈表示第2个扩频码,每格代表1个码片,可见两个扩频码的第1到L个码片分别对齐。右图表示偏移量k等于3时两个扩频码码片的位置关系。相对于左图,右图外圈沿箭头所示方向旋转3格。对于任意1≤k≤L的值,外圈沿图中箭头所示方向旋转k格即实现偏移量为k的循环移位。若第1个扩频码和第2个扩频码相同,则通过循环移位获得扩频码的自相关值。若第1个扩频码和第2个扩频码不相同,则通过循环移位获得扩频码的互相关值。
(6)将最大自相关函数旁峰值按照从小到大排序,将最大自相关旁峰值从小到大排序,数值相同的划入同一集合形成M2个集合,其中1≤j≤M2,最小值对应的集合为最大值对应的集合为依照排序得出最大自相关函数旁峰值按的扩频码计算值,对应的扩频码计算值为1,对应的扩频码计算值为2,依次类推,对应的扩频码为M2。所有扩频码的计算值构成集合S2, 为第i个扩频码的计算值。
(9)将最大互相关峰值按照从小到大排序,数值相同的划入同一集合形成M3个集合,其中1≤j≤M3,最小值对应的集合为最大值对应的集合为依照排序结果得出最大互相关峰值的扩频码计算值,对应的扩频码计算值为1,对应的扩频码计算值为2,依次类推,对应的扩频码计算值为M3。所有扩频码的计算值构成集合S3, 为第i个扩频码的计算值。
(10)依次根据第i个扩频码的循环移位自相关值Ri(k),计算码组中第i个扩频码的功率谱值Θi,即其中Ri为Ri(k)组成的序列,1≤k≤L,为傅里叶变换,具体为:Θi为Θi(m)组成的序列,0≤m≤L-1。
(12)将功率谱最大值按照从小到大排序,数值相同的划入同一集合形成M4个集合,其中1≤j≤M4,最小值对应的集合为最大值对应的集合为依照排序结果得出功率谱最大值的扩频码计算值,对应的扩频码计算值为1,对应的扩频码计算值为2,依次类推,对应的扩频码为M4。所有扩频码的计算值构成集合S4, 为第i个扩频码的计算值。
(13)分别给码平衡性、最大自相关旁峰值,最大互相关峰值,功率谱最大值的扩频码的计算值集合S1、S2、S3和S4赋予权值w1,w2,w3,w4;w1,w2,w3,w4应满足0≤w1≤1,0≤w2≤1,0≤w3≤1和0≤w4≤1,并且w1+w2+w3+w4=1。对于带宽、发射功率受限的系统取w1=0.4,w2=w3=w4=0.2。对于要求接收机假锁概率的系统取w2=0.4,w1=w3=w4=0.2。对于要求抗多址干扰能力的系统取w3=0.4,w1=w2=w4=0.2。对于要求抗窄带干扰能力的系统w4=0.4,w1=w2=w3=0.2。其余取系统w1=w2=w3=w4=0.25。
(14)根据权值和扩频码计算值集合S1、S2、S3和S4,计算各个扩频码的加权计算值S,即S=w1·S1+w2·S2+w3·S3+w4·S4;加权计算值S值最小的扩频码为最优。
本发明说明书中未作详细描述的内容属于本领域专业技术人员的公知技术。
Claims (1)
1.一种优选导航卫星扩频码的方法,其特征在于步骤如下:
(3)将步骤(2)中得到的码平衡性Di的值从小到大排序,数值相同的Di划入同一集合形成M1个集合,其中1≤j≤M1;数值最小的Di对应的集合为数值最大的Di对应的集合为依照排序结果得出码平衡性Di的扩频码计算值,对应的扩频码计算值为1,对应的扩频码计算值为2,依次类推,对应的扩频码计算值为M1;所有扩频码计算值构成集合其中为第i个扩频码的计算值;
(4)依次计算码组中第i个扩频码的循环移位自相关值1≤l≤L,L为扩频码的长度;Ci(l)为第i个扩频码的第l个码片,k为码片循环移位的偏移量,0≤k≤L-1,Ci(l+k)为第i个扩频码的第l+k个码片;
(6)将步骤(5)中得出的最大自相关函数旁峰值从小到大排序,数值相同的划入同一集合形成M2个集合,其中1≤j≤M2,数值最小的对应的集合为数值最大的对应的集合为依照排序结果得出最大自相关函数旁峰值的扩频码计算值,对应的扩频码计算值为1,对应的扩频码计算值为2,依次类推,对应的扩频码计算值为M2;所有扩频码的计算值构成集合其中为第i个扩频码的计算值;
(9)将步骤(8)中获得的最大互相关峰值从小到大排序,数值相同的划入同一集合形成M3个集合,其中1≤j≤M3,数值最小的对应的集合为数值最大的对应的集合为依照排序结果得出最大互相关峰值的扩频码计算值,对应的扩频码计算值为1,对应的扩频码计算值为2,依次类推,对应的扩频码计算值为M3;所有扩频码的计算值构成集合其中为第i个扩频码的计算值;
(10)根据步骤(4)中获得的第i个扩频码的循环移位自相关值Ri(k),依次计算码组中第i个扩频码的功率谱值Θi=F(Ri),其中F为傅里叶变换,Ri为由Ri(k)组成的序列;
(11)根据步骤(10)中获得的第i个扩频码的功率谱Θi,计算得到功率谱最大值
(12)将步骤(11)中获得的功率谱最大值从小到大排序,数值相同的划入同一集合形成M4个集合,其中1≤j≤M4,数值最小的对应的集合为数值最大的对应的集合为依照排序结果得出功率谱最大值的扩频码计算值,对应的扩频码计算值为1,对应的扩频码计算值为2,依次类推,对应的扩频码计算值为M4;所有扩频码的计算值构成集合 为第i个扩频码的计算值;
(13)分别给码平衡性Di、最大自相关旁峰值最大互相关峰值功率谱最大值的扩频码计算值集合S1、S2、S3和S4赋予杈值w1,w2,w3,w4,权值w1,w2,w3,w4应满足0≤w1≤1,0≤w2≤1,0≤w3≤1和0≤w4≤1,且w1+w2+w3+w4=1;对于带宽、发射功率受限的系统取w1=0.4,w2=w3=w4=0.2;对于要求接收机假锁概率的系统取w2=0.4,w1=w3=w4=0.2;对于要求抗多址干扰能力的系统取w3=0.4,w1=w2=w4=0.2;对于要求抗窄带干扰能力的系统w4=0.4,w1=w2=w3=0.2;其余系统取w1=w2=w3=w4=0.25;
(14)根据权值和扩频码计算值集合S1、S2、S3和S4,计算各个扩频码的加权计算值S,即S=w1·S1+w2·S2+w3·S3+w4·S4;加权计算值S值最小的扩频码为最优。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210453838.7A CN102937717B (zh) | 2012-11-13 | 2012-11-13 | 一种优选导航卫星扩频码的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210453838.7A CN102937717B (zh) | 2012-11-13 | 2012-11-13 | 一种优选导航卫星扩频码的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102937717A CN102937717A (zh) | 2013-02-20 |
CN102937717B true CN102937717B (zh) | 2014-07-02 |
Family
ID=47696627
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210453838.7A Active CN102937717B (zh) | 2012-11-13 | 2012-11-13 | 一种优选导航卫星扩频码的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102937717B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103580717B (zh) * | 2013-10-17 | 2015-04-15 | 华中科技大学 | 一种导航信号扩频码的优选方法 |
CN104143997B (zh) * | 2014-07-23 | 2016-06-01 | 西安空间无线电技术研究所 | 一种基于扩频体制的多址干扰抗错锁判决方法 |
CN108880732B (zh) * | 2018-04-18 | 2019-10-22 | 中国电子科技集团公司第五十四研究所 | 一种基于分治算法的卫星导航系统扩频码互相关优化方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102638319A (zh) * | 2012-04-23 | 2012-08-15 | 北京空间飞行器总体设计部 | 一种导航卫星二进制偏移载波信号的调制性能测试方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007116578A (ja) * | 2005-10-24 | 2007-05-10 | Sony Corp | 衛星信号受信装置および衛星信号受信方法 |
JP5375774B2 (ja) * | 2010-08-30 | 2013-12-25 | セイコーエプソン株式会社 | コード位相誤差算出方法及びコード位相算出方法 |
-
2012
- 2012-11-13 CN CN201210453838.7A patent/CN102937717B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102638319A (zh) * | 2012-04-23 | 2012-08-15 | 北京空间飞行器总体设计部 | 一种导航卫星二进制偏移载波信号的调制性能测试方法 |
Non-Patent Citations (5)
Title |
---|
JP特开2011-7808A 2011.01.13 |
一种新的混沌扩频序列优选算法;罗冬梅等;《计算机工程》;20101231;第36卷(第24期);99-101 * |
数字化移动通信系统中扩频码优选;邬国扬;《无线电工程》;19921231;第22卷(第4期);1-8 * |
罗冬梅等.一种新的混沌扩频序列优选算法.《计算机工程》.2010,第36卷(第24期),99-101. |
邬国扬.数字化移动通信系统中扩频码优选.《无线电工程》.1992,第22卷(第4期),1-8. |
Also Published As
Publication number | Publication date |
---|---|
CN102937717A (zh) | 2013-02-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105282021B (zh) | 信号集中器设备 | |
CN102540187B (zh) | 编队卫星sar的一种正交波形设计方法 | |
CN103424754B (zh) | 基于合成相关函数的mboc调制信号无模糊多径抑制方法 | |
CN102937717B (zh) | 一种优选导航卫星扩频码的方法 | |
CN105577229B (zh) | Cpu辅助gpu实现扩频信号快速捕获方法 | |
CN103018721B (zh) | 一种mimo-oth雷达波形的生成方法 | |
CN101072044A (zh) | 一种长周期扩频码的频率折叠时频并行搜索方法 | |
CN100553169C (zh) | 基于频率折叠的短周期扩频码时频并行搜索方法 | |
CN108401489B (zh) | 一种量子混沌波包数字信号生成方法 | |
CN100562762C (zh) | 用于定位全球定位系统相关峰值信号的方法和系统 | |
CN105842713A (zh) | 一种基于sft的ins辅助bds信号快速捕获方法 | |
CN103259755A (zh) | 一种全域覆盖多波束卫星lte的主同步序列设计方法 | |
CN104468451A (zh) | 一种基于变换域通信的多级循环移位键控调制解调方法 | |
CN105137454A (zh) | 一种基于协方差矩阵特征分解的抗干扰算法的fpga实现方法及实现装置 | |
CN101478331A (zh) | 一种双时频非重复非完全跳频序列的构造方法 | |
CN111398946A (zh) | 一种基于lfm信号调频率调制的共享信号设计方法 | |
CN105372678B (zh) | 一种正弦boc调制信号的无模糊跟踪方法 | |
CN105553635A (zh) | 基于循环自相关的boc信号参数盲估计方法 | |
Qi et al. | Fast acquisition method of navigation receiver based on folded PMF-FFT | |
CN102801459A (zh) | 多模卫星信号接收装置及其方法 | |
CN103763233A (zh) | 基于数字调制解调的gsm外辐射雷达信号处理方法 | |
CN105608057A (zh) | 一种分时复用硬件资源的信号子空间分解的fpga实现模块及其fpga实现方法 | |
CN103152070A (zh) | 一种基于可变位标序列的扩频通信方法 | |
CN103701539A (zh) | 一种cdma用户数量测定方法 | |
CN103698751A (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 |