CN111257911A - 一种基于数字波束成形的码片级脉冲跳时导航信号生成与播发实现方法 - Google Patents

一种基于数字波束成形的码片级脉冲跳时导航信号生成与播发实现方法 Download PDF

Info

Publication number
CN111257911A
CN111257911A CN202010148020.9A CN202010148020A CN111257911A CN 111257911 A CN111257911 A CN 111257911A CN 202010148020 A CN202010148020 A CN 202010148020A CN 111257911 A CN111257911 A CN 111257911A
Authority
CN
China
Prior art keywords
chip
navigation
digital
signal
hopping
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
CN202010148020.9A
Other languages
English (en)
Other versions
CN111257911B (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.)
Xian Institute of Space Radio Technology
Original Assignee
Xian Institute of Space Radio Technology
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 Xian Institute of Space Radio Technology filed Critical Xian Institute of Space Radio Technology
Priority to CN202010148020.9A priority Critical patent/CN111257911B/zh
Publication of CN111257911A publication Critical patent/CN111257911A/zh
Application granted granted Critical
Publication of CN111257911B publication Critical patent/CN111257911B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/02—Details of the space or ground control segments
    • 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/21—Interference related issues ; Issues related to cross-correlation, spoofing or other methods of denial of service
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00—Reducing energy consumption in communication networks
    • Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks

Abstract

本发明一种基于数字波束成形的码片级脉冲跳时导航信号生成与播发实现方法，步骤如下：(1)生成信号伪码和导航电文，进行基带调制，得到基带调制信号；对基带调制信号进行数字低通滤波，生成数字滤波后的基带导航信号
Figure DDA0002401446230000011
(2)根据阵列天线形成的波束数为，生成M个跳时伪码序列；(3)读取数字波束成形系数表，获得N个加权系数；将N个加权系数分别与步骤(1)数字滤波后的基带导航信号
Figure DDA0002401446230000012
相乘，得到加权后的N路数字导航信号，分别为，完成数字波束成形；(4)将步骤(3)加权后的N路数字导航信号，依次经过DAC、上变频、滤波、功率放大后，经N个天线阵元播发，降低实现复杂度，提升导航信号的抗干扰能力。

Description

一种基于数字波束成形的码片级脉冲跳时导航信号生成与播 发实现方法
技术领域
本发明涉及一种基于数字波束成形的码片级脉冲跳时导航信号生成与播发实现方法，属于卫星导航技术领域。
背景技术
近年来，全球卫星导航系统(GNSS)信号的脆弱性问题引起来广泛的关注。一方面，卫星导航信号的落地功率低(约-160dBW)，容易受到有意或者无意的干扰，另一方面，GNSS开放服务信号结构公开，容易受到欺骗干扰的攻击。
为了缓解GNSS信号的脆弱性问题，出现了多种方法。对于GNSS开放服务信号，Galileo和GPS研究引入了信号认证功能，采用导航电文认证(NMA)和扩频码认证(SCA)，提升开放服务信号的抗欺骗能力。对于GNSS授权服务信号，一种方式是采用点波束功率增强。专利“一种基于Ka频段多波束天线的低轨卫星导航信号功率增强方法”(专利号201811068149.8)给出了一种导航信号功率增强方法，将低轨卫星多波束天线形成的所有波束指向导航信号功率增强目标区域，提高信号落地功率。点波束功率增强方法，缩小了覆盖范围。
为了提升导航信号的抗干扰能力，另一方法是采用抗干扰信号体制。专利“一种混合扩频信号的生成方法、生成装置及发送装置”(专利号201810141670.3)给出了一种混合直扩跳频信号生成方法。专利“一种DS-FH-TH混合扩频系统”(专利号201620428971.0)给出了一种直扩跳频跳时信号生成系统。然而，直扩跳频信号的载波相位的高精度测量存在一定问题，传统的直扩跳时信号在时域上划分多个时隙，按照时隙进行播发，信号测量不连续。
另一种抗干扰信号体制是码片级脉冲跳时信号体制，在码片层面进行跳时，将码序列以每Nc个码片分组，每Nc个码片只播发一个码片，播发码片的位置是伪随机跳动的，Nc个码片的功率集中于一个码片上，实现将低平均功率的连续信号变为高瞬时功率的随机脉冲位置准连续信号，在同等平均功率前提下抗干扰能力提升(10lgNc)dB。然而，对于这种码片级脉冲跳时信号，需要采用脉冲固放实现信号放大，存在两个问题，一是难以实现高瞬时功率，二是脉冲切换时间难以保证，工程实现复杂度高。
发明内容
本发明解决的技术问题为：克服上述现有技术的不足，提供一种基于数字波束成形的码片级脉冲跳时导航信号生成与播发实现方法，降低实现复杂度，提升导航信号的抗干扰能力。
本发明解决的技术方案为：一种基于数字波束成形的码片级脉冲跳时导航信号生成与播发实现方法，步骤如下：
(1)生成信号伪码和导航电文，进行基带调制，得到基带调制信号sbase(t)；对基带调制信号sbase(t)进行数字低通滤波，生成数字滤波后的基带导航信号
Figure BDA0002401446210000021
(2)根据阵列天线形成的波束数为Nbeam，生成M个跳时伪码序列，
Figure BDA0002401446210000022
(3)读取数字波束成形系数表，获得N个加权系数；将N个加权系数分别与步骤(1)数字滤波后的基带导航信号
Figure BDA0002401446210000023
相乘，得到加权后的N路数字导航信号，分别为Sn(t),n＝1,2,…,N，完成数字波束成形；
(4)将步骤(3)加权后的N路数字导航信号，依次经过DAC、滤波、上变频、功率放大后，经N个天线阵元播发，实现码片级脉冲跳时导航信号生成与播发。
优选的，步骤(1)中所述的数字滤波后的基带导航信号
Figure BDA0002401446210000024
生成，通过以下方法得到：
(1)信号伪码发生器生成信号伪码序列记为{cs}，cs∈{0,1}，s＝0,1,2,…，信号伪码序列的码速率为Rchip，信号伪码序列的码片宽度为Tchip＝1/Rchip；对码片进行赋形，得到码波形，实现扩频码(即伪码)波形生成；
码波形为：
Figure BDA0002401446210000031
式中，t表示时间，s表示码片的序号，(1-2cs)将码片值{0,1}映射为{1,-1}，
Figure BDA0002401446210000032
为信号伪码序列的码片波形，可以采用当前GNSS信号的任意一种波形；对于矩形码片波形，有：
Figure BDA0002401446210000033
(2)将导航电文经过组帧、信号编码后，得到符号速率为Rd的电文数据{dl}，dl∈{1,-1}，l＝0,1,2,…，导航电文符号宽度Td＝1/Rd，且Td是Tchip的整数倍，将电文数据写为连续波形形式D(t)，即导航电文生成，则有：
Figure BDA0002401446210000034
式中，pd(t)为宽度Td的矩形函数：
Figure BDA0002401446210000035
(3)将连续波形形式的电文数据调制到伪码波形上，得到基带调制信号：
sbase(t)＝D(t)·Cs(t)；
(4)将sbase(t)进行FIR数字低通滤波，得到低通滤波后的信号
Figure BDA0002401446210000036
完成数字滤波后的基带导航信号生成。
优选的，步骤(2)的跳时伪码序列生成，满足以下要求：
(1)跳时伪码发生器生成M个跳时伪码序列，分别为
Figure BDA0002401446210000041
Figure BDA0002401446210000042
Figure BDA0002401446210000043
m＝1,2,……,M，s＝0,1,2,…；
Figure BDA0002401446210000044
Nbeam为阵列天线形成的波束数；
(2)M个跳时伪码序列的码速率与信号伪码序列的码速率一样，都是Rchip；跳时伪码序列与信号伪码序列在时域上是对齐的；
(3)M个跳时伪码序列的第s个码片值的组合，决定了导航信号的第s个码片值的播发波束(即决定了导航信号的第s个码片值在哪个波束播发)。
优选的，步骤(3)所述的数字波束成形，通过以下方法得到：
(1)数字波束成形系数表存储了加权系数；数字波束成形系数表是一个Nbeam×N二维表，记为{wb,k}，b＝1,2,……,Nbeam，k＝1,2,……,N；数字波束成形系数表通过数值优化生成，与阵列天线的布局以及波束增益方向图设计有关；使用该表中第nbeam行的系数进行数字波束成形，能够让信号在第nbeam个波束播发，nbeam＝1,2，……，Nbeam；
(2)使用M个跳时伪码序列，控制数字波束成形系数表输出加权系数，具体为：
在第s个码片时刻，M个跳时伪码序列的值分别为：
Figure BDA0002401446210000045
转换为十进制数，即
Figure BDA0002401446210000046
将NTH,s对Nbeam进行求模加1，得到属于{1,2,…,Nbeam}的数Ns，实现数字波束成形系数选择，即
Ns＝(NTH,s mod Nbeam)+1
在第s个码片时刻，信号在第Ns个波束播发，数字波束成形系数表输出第Ns行加权系数
Figure BDA0002401446210000047
k＝1,2,……,N，
Figure BDA0002401446210000048
是一个复系数，表示为实部和虚部的形式，有
Figure BDA0002401446210000049
Figure BDA00024014462100000410
是
Figure BDA00024014462100000411
的实部，
Figure BDA00024014462100000412
是
Figure BDA00024014462100000413
的虚部；
(3)将生成的基带导航信号
Figure BDA0002401446210000051
与N个加权系数分别相乘，实现数字波束成形，得到第s个码片的持续时间内N个加权后的数字导航信号，分别为：
Figure BDA0002401446210000052
实现数字波束成形。
优选的，生成信号伪码和导航电文，进行基带调制，具体为：
生成信号伪码和导航电文，将导航电文调制到伪码上。
优选的，数字波束成形时，阵列天线的阵元数为N，形成的波束数为Nbeam，数字波束成形系数表存储了加权系数。
优选的，数字波束成形系数表是一个Nbeam×N二维表将生成基带导航信号。
优选的，M个跳时伪码序列，控制不同的码片伪随机的在不同的波束里面播发。
优选的，将N个加权系数分别与步骤(1)数字滤波后的基带导航信号
Figure BDA0002401446210000053
相乘，得到加权后的N路数字导航信号，具体如下：
在跳时伪码序列控制下，每一个信号伪码期间，选择N个加权系数，分别与基带导航信号相乘，得到加权后的N路数字导航信号，分别为Sn(t),n＝1,2,…,N。
优选的，(4)将步骤(3)加权后的N路数字导航信号，依次经过DAC、滤波、上变频、功率放大后，经N个天线阵元播发，具体为
第n路数字导航信号Sn(t),n＝1,2,…,N经第n个DAC、第n个滤波器、第n个上变频器、第n个功率放大器、第n个天线阵元播发，实现码片级脉冲跳时导航信号生成与播发。
本发明与现有技术相比的优点在于：
(1)本发明公开了一种数字波束成形的码片级脉冲跳时导航信号的实现方法，实现将连续导航信号的不同码片在不同的波束中播发，每个波束下收到的都是等效的码片级脉冲跳时导航信号，提升了导航信号的抗干扰能力。
(2)当前的码片级脉冲跳时导航信号生成方案，需要采用脉冲固放，将低平均功率的连续信号变为高瞬时功率的随机脉冲位置准连续信号，但是存在两个问题，一是难以实现高瞬时功率，二是脉冲切换时间难以保证，工程实现复杂度高。而本专利提出的基于数字波束成形的码片级脉冲跳时导航信号实现方法，不需要使用脉冲固放，采用连续工作的固放，每个固放的功率要求低，工程实现性好，导航信号的质量有保障。
(3)传统的混合扩频抗干扰系统，采用直扩+跳频+扩频体制实现抗干扰能力的提升，信号体制与当前卫星导航信号体制不同，接收复杂度高，而本专利提供的基于数字波束成形的码片级脉冲跳时导航信号，即使在不使用跳时伪码序列的情况下，可以采用当前的接收方法进行接收，兼容性好。
(4)本发明提出的基于数字波束成形的码片级脉冲跳时导航信号实现方法，从时域上上看，导航信号的不同码片在不同波束播发，每个波束下的用户都能接收到导航信号，服务区域不缩小，不缩小服务区域。在平均功率一定的前提下，对于有Nbeam个波束的系统，抗干扰能力能够提升(10lgNbeam)dB。
附图说明
图1为本发明公开的基于数字波束成形的码片级脉冲跳时导航信号实现方法的技术方案流程图；
图2为跳时伪码序列与信号伪码序列的时域关系示意图；
图3为不同波束下接收到的码片级脉冲跳时信号的示意图；
图4为当前卫星导航信号与本发明导航信号的干扰容限对比图。
具体实施方式
下面结合附图和具体实施例对本发明做进一步详细描述。
本发明一种基于数字波束成形的码片级脉冲跳时导航信号生成与播发实现方法，步骤如下：(1)生成信号伪码和导航电文，进行基带调制，得到基带调制信号；对基带调制信号进行数字低通滤波，生成数字滤波后的基带导航信号
Figure BDA0002401446210000071
(2)根据阵列天线形成的波束数为，生成M个跳时伪码序列；(3)读取数字波束成形系数表，获得N个加权系数；将N个加权系数分别与步骤(1)数字滤波后的基带导航信号
Figure BDA0002401446210000072
相乘，得到加权后的N路数字导航信号，分别为，完成数字波束成形；(4)将步骤(3)加权后的N路数字导航信号，依次经过DAC、上变频、滤波、功率放大后，经N个天线阵元播发，降低实现复杂度，提升导航信号的抗干扰能力。
导航卫星是一个能量受限的系统，导致卫星导航信号的落地功率低，通常只有-160dBW左右，抗干扰能力弱。因此，在复杂电磁环境下，存在干扰(如匹配谱干扰)时，卫星导航服务可能不可用。为提高卫星导航信号抗干扰能力，提升服务可用性，本发明提出了一种基于数字波束成形的码片级脉冲跳时导航信号的生成与播发方法，由基带信号生成、跳时伪码发生器、数字波束成形、N个DAC、N个上变频器、N个滤波器、N个天线阵元组成，N个天线阵元形成多波束阵列天线，每一个阵元前面都有个功率放大器。通过数字波束成形，控制码片在不同的波束播发，在每个波束覆盖区域内，收到的都是码片级脉冲跳时导航信号。由于功率放大器的输入和输出是连续的，能够有效保障导航信号质量，每个固放的功率要求不高，降低工程实现难度。采用本法明的方法，能够在导航卫星能量受限，平均功率一定的情况下，提升导航信号的抗干扰能力，同时，不减小服务区域。能够在有干扰的环境下，提供导航定位与授时服务，提高卫星导航系统服务的可用性。
本发明的一种码片级脉冲跳时导航信号的生成与播发方法，优选步骤包括基带导航信号生成、跳时伪码序列、数字波束成形、射频滤波放大与播发，如图1所示，具体如下：
(1)生成信号伪码和导航电文，信号伪码发生器生成信号伪码序列记为{cs}，cs∈{0,1}，s＝0,1,2,…，信号伪码序列的码速率为Rchip，信号伪码序列的码片宽度为Tchip＝1/Rchip。对码片进行赋形，得到码波形，实现扩频码(即伪码)波形生成；
码波形为：
Figure BDA0002401446210000081
式中，t表示时间，s表示码片的序号，(1-2cs)将码片值{0,1}映射为{1,-1}，
Figure BDA0002401446210000082
为信号伪码序列的码片波形，可以采用当前GNSS信号的任意一种波形。对于矩形码片波形，有：
Figure BDA0002401446210000083
导航电文由地面站上注到卫星上，主要包含了导航卫星的轨道参数(即导航卫星星历)和钟差参数，也可以包含电离层改正参数等。将导航电文经过组帧、信号编码后，得到符号速率为Rd的电文数据{dl}，dl∈{1,-1}，l＝0,1,2,…，导航电文符号宽度Td＝1/Rd，且Td是Tchip的整数倍，将电文数据写为连续波形形式D(t)，即导航电文生成，则有：
Figure BDA0002401446210000084
式中，pd(t)为宽度Td的矩形函数：
Figure BDA0002401446210000085
将连续波形形式的电文数据调制到伪码波形上，进行基带调制，得到基带调制信号：
sbase(t)＝D(t)·Cs(t)；
将sbase(t)进行FIR数字低通滤波，得到低通滤波后的信号
Figure BDA0002401446210000086
完成数字滤波后的基带导航信号生成。
导航电文宽度Td是信号伪码序列的码片宽度Tchip的整数倍，因此，基带导航信号的最小时间间隔为码片级，即码片宽度Tchip。
(2)根据阵列天线形成的波束数为Nbeam，生成M个跳时伪码序列，
Figure BDA0002401446210000091
跳时伪码发生器生成M个跳时伪码序列，分别为
Figure BDA0002401446210000092
Figure BDA0002401446210000093
Figure BDA0002401446210000094
m＝1,2,……,M，s＝0,1,2,…。
Figure BDA0002401446210000095
Nbeam为阵列天线形成的波束数。
M个跳时伪码序列的码速率与信号伪码序列的码速率一样，都是Rchip。跳时伪码序列与信号伪码序列在时域上是对齐的。
M个跳时伪码序列的第s个码片值的组合，决定了导航信号的第s个码片值的播发波束(即决定了导航信号的第s个码片值在哪个波束播发)。
由于M个跳时伪码序列的码速率与信号伪码序列的码速率一样，都是Rchip，因此，跳时伪码序列的最小时间间隔也是码片级，即码片宽度Tchip。跳时伪码序列与信号伪码序列的时域关系示意图如图2所示。
(3)读取数字波束成形系数表，获得N个加权系数；将N个加权系数分别与步骤(1)数字滤波后的基带导航信号
Figure BDA0002401446210000096
相乘，得到加权后的N路数字导航信号，分别为Sn(t),n＝1,2,…,N，完成数字波束成形；
数字波束成形系数表存储了加权系数。数字波束成形系数表是一个Nbeam×N二维表，记为{wb,k}，b＝1,2,……,Nbeam，k＝1,2,……,N。数字波束成形系数表通过数值优化生成，与阵列天线的布局以及波束增益方向图设计有关。
N个天线阵元的相对位置坐标记为dn＝[xn,yn,zn]T,n＝1,2,L,N，对于阵列天线下方的一个点，天顶角记为
Figure BDA0002401446210000097
方位角记为θ，则其单位方向向量
Figure BDA0002401446210000098
Figure BDA0002401446210000099
N个阵元在该方向下的阵列矢量为：
Figure BDA0002401446210000101
式中，fc为载波频率。将
Figure BDA0002401446210000102
与加权系数矢量点积，得到天顶角
Figure BDA0002401446210000103
方位角θ时的增益，即
Figure BDA0002401446210000104
对于数字波束成形系数表的第b行的加权系数，设置目标函数：
Figure BDA0002401446210000105
其中，
Figure BDA0002401446210000106
表示第b个波束的天顶角范围，θb表示第b个波束的方位角范围，
Figure BDA0002401446210000107
表示除去第b个波束外的天顶角范围，
Figure BDA0002401446210000108
表示除去第b个波束外的方位角范围，GTH为设置的门限值，表示波束间的增益抑制差。
通过已有的数值优化算法，如拟牛顿法，遗传算法或者粒子群算法，优化计算出每一行的加权系数。
使用数字波束成形系数表中第nbeam行的系数进行数字波束成形，能够让信号在第nbeam个波束播发，nbeam＝1,2，……，Nbeam。
使用M个跳时伪码序列，控制数字波束成形系数表输出加权系数，具体为：在第s个码片时刻，M个跳时伪码序列的值分别为：
Figure BDA0002401446210000109
转换为十进制数，即
Figure BDA00024014462100001010
将NTH,s对Nbeam进行求模加1，得到属于{1,2,…,Nbeam}的数Ns，实现数字波束成形系数选择，即
Ns＝(NTH,s mod Nbeam)+1
在第s个码片时刻，信号在第Ns个波束播发，数字波束成形系数表输出第Ns行加权系数
Figure BDA0002401446210000111
k＝1,2,……,N，
Figure BDA0002401446210000112
是一个复系数，表示为实部和虚部的形式，有
Figure BDA0002401446210000113
Figure BDA0002401446210000114
是
Figure BDA0002401446210000115
的实部，
Figure BDA0002401446210000116
是
Figure BDA0002401446210000117
的虚部。
将生成的基带导航信号
Figure BDA0002401446210000118
与N个加权系数分别相乘，实现数字波束成形，得到第s个码片的持续时间内N个加权后的数字导航信号，分别为：
Figure BDA0002401446210000119
实现数字波束成形。
由于M个跳时伪码序列码速率与信号伪码序列码速率一样，因此，M个跳时伪码序列的值以码片级变化，N个波束成形系数
Figure BDA00024014462100001110
也是以码片级变化，从而实现让导航信号以码片级在不同的波束播发，即生成了码片级脉冲跳时导航信号。
(4)将步骤(3)加权后的N路数字导航信号，依次经过DAC、滤波、上变频、功率放大后，经N个天线阵元播发；第n路数字导航信号Sn(t),n＝1,2,…,N经第n个DAC、第n个上变频器、第n个滤波器、第n个功率放大器、n个天线阵元播发，实现码片级脉冲跳时导航信号生成与播发。在每个波束中，播发的导航信号是脉冲的，脉冲的位置是伪随机跳动，脉冲的最小宽度是码片宽度Tchip。因此，是码片级脉冲跳时导航信号。
本发明的设计思路如下：导航卫星是能量受限系统，平均功率一定，为了提升导航信号的抗干扰能力，可以将信号集中一个小区域内，提升信号的落地功率，但是，这样会降低服务区域。为此，本发明采用基于数字波束成形的码片级脉冲跳时导航信号方案，将服务区域划分为Nbeam个波束，每一个码片持续时间内，只在一个波束播发，不同码片由跳时伪码序列控制，在不同的波束播发。在每一个波束内，都是码片级脉冲跳时导航信号，瞬时功率提升Nbeam，平均功率保持不变。干扰信号由于不知道跳时伪码序列，只能发连续干扰信号。码片级脉冲跳时导航信号的接收终端具有跳时伪码序列，在时域上只接收对应的码片，因此，干扰信号功率等效降低Nbeam倍，抗干扰能力增加10lg(Nbeam)dB。
1)本发明实现信号生成时，信号伪码可采用周期性的伪码，码长为Lc，为保证信号伪码的自相关特性与互相关特性，需要满足Lc/Nbeam≥100，Nbeam为波束个数。
2)本发明实现信号生成时，为了保证足够的抗干扰能力提升，波束Nbeam≥10，且log2(Nbeam)为整数，保证码片在Nbeam个波束间均匀的跳动。
3)本发明在数字波束成形系数表生成时，为了保证足够的波束间隔离度，数值优化时目标函数的门限值满足GTH≥15dB。
(1)基带导航信号生成
基带导航信号生成包括信号伪码生成，导航电文生成、电文调制、以及数字低通滤波。信号伪码的码序列为周期伪码序列，码序列长度为20460，记为{cs}，cs∈{0,1}，s＝0,1,2,…，20460,信号伪码的码速率为Rchip＝10.23Mcps，码片宽度为Tchip＝1/Rchip。对码片进行赋形，得到一个周期的码波形：
Figure BDA0002401446210000121
其中，(1-2cs)将码片值{0,1}映射为{1,-1}，
Figure BDA0002401446210000122
为信号伪码码片波形，采用矩形码片波形，有：
Figure BDA0002401446210000123
导航电文符号速率为Rd＝100sps的电文数据流{dl}，dl∈{1,-1}，l＝0,1,2,…，电文符号宽度Td＝1/Rd。将电文数据写为连续波形形式，有：
Figure BDA0002401446210000131
式中，pd(t)为宽度Td的矩形函数：
Figure BDA0002401446210000132
将电文调制到伪码上，得到：
sbase(t)＝D(t)·Cs(t)
将sbase(t)进行低通数字滤波得到低通滤波后的信号
Figure BDA0002401446210000133
完成基带导航信号生成。低通数字滤波采用FIR滤波器，单边3dB带宽为30.69MHz。
(2)跳时伪码生成
阵列天线形成的波束数为Nbeam＝128，为了控制不同的码片伪随机的在不同的波束里面播发，跳时伪码序列的个数为
Figure BDA0002401446210000134
跳时伪码发生器生成7个跳时伪码序列，分别为
Figure BDA0002401446210000135
s＝0,1,2,…。7个跳时伪码序列的码速率以信号伪码序列一样，都是Rchip＝10.23Mcps。7个跳时伪码序列的第s个码片值的组合，决定了导航信号的第s个码片值在哪个波束播发。
(3)数字波束成形
阵列天线的阵元数为N＝256，形成的波束数为Nbeam＝128，数字波束成形系数表存储了加权系数。数字波束成形系数表是一个100×200的二维表，记为{wb,k}，b＝1,2,……,128，k＝1,2,……,256。数字波束成形系数表通过优化生成，与阵列天线的布局以及波束增益方向图设计有关。
200个天线阵元的相对位置坐标记为dn＝[xn,yn,zn]T,n＝1,2,L,256，对于阵列天线下方的一个点，天顶角记为
Figure BDA0002401446210000136
方位角记为θ，则其单位方向向量
Figure BDA0002401446210000137
Figure BDA0002401446210000138
N个阵元在该方向下的阵列矢量为：
Figure BDA0002401446210000141
式中，fc为载波频率,如1575.42MHz。将
Figure BDA0002401446210000142
与加权系数矢量点积，得到天顶角
Figure BDA0002401446210000143
方位角θ时的增益，即
Figure BDA0002401446210000144
对于数字波束成形系数表的第b行的加权系数，设置目标函数：
Figure BDA0002401446210000145
其中，
Figure BDA0002401446210000146
表示第b个波束的天顶角范围，θb表示第b个波束的方位角范围，
Figure BDA0002401446210000147
表示除去第b个波束外的天顶角范围，
Figure BDA0002401446210000148
表示除去第b个波束外的方位角范围，GTH为设置的门限值，表示波束间的增益抑制差，为保证足够高的波束隔离度，GTH＝20dB。
通过已有的数值优化算法，如拟牛顿法，遗传算法或者粒子群算法，优化计算出每一行的加权系数。
使用该表中第nbeam行的系数进行数字波束成形，能够让信号在第nbeam个波束播发，nbeam＝1,2，……，128。
使用7个跳时伪码序列，控制数字波束成形系数表输出系数。在第s个码片时刻，7个跳时伪码序列的值分别为：
Figure BDA0002401446210000149
转换为十进制数，即
Figure BDA00024014462100001410
将NTH,s对Nbeam＝128进行求模加1，得到属于{1,2,…,128}的数Ns，即
Ns＝(NTH,s mod Nbeam)+1
因此，在第s个码片时刻，信号在第Ns个波束播发，数字波束成形系数表输出第Ns行系数
Figure BDA0002401446210000151
k＝1,2,……,N，
Figure BDA0002401446210000152
是一个复系数，表示为实部和虚部的形式，有
Figure BDA0002401446210000153
Figure BDA0002401446210000154
是
Figure BDA0002401446210000155
的实部，
Figure BDA0002401446210000156
是
Figure BDA0002401446210000157
的虚部。
将生成的基带导航信号
Figure BDA0002401446210000158
与256个加权系数分别相乘，得到第s个码片期间，200个数字导航信号，分别为：
Figure BDA0002401446210000159
(4)射频滤波放大与播发
将256个数字导航信号，分别经DAC、滤波、放大后，经256个天线阵元播发。第n个数字导航信号Sn(t),n＝1,2,…,256经第n个DAC、第n个滤波器、第n个放大器、第n个天线阵元播发。
单个波束下接收到的码片级脉冲跳时信号，与基带导航信号的示意图如图3所示。可以看到，单个波束下受到的信号是码片级脉冲跳时信号，平均功率与传统导航信号一样，但是码片脉冲峰值功率高了128倍。
本发明通过理论与仿真分析，采用干扰容限指标，对比本发明与当前卫星导航信号的抗干扰能力，干扰信号采为匹配谱干扰，噪声功率谱密度为-204dBW/Hz，载波跟踪门限为28dB-Hz。当平均发射功率相同时，当前卫星导航信号与本发明的干扰容限如图4所示。可以看到，同样的平均功率下，本发明抗干扰能力提升了20dB。

Claims (10)

1.一种基于数字波束成形的码片级脉冲跳时导航信号生成与播发实现方法，其特征在于步骤如下：
(1)生成信号伪码和导航电文，进行基带调制，得到基带调制信号sbase(t)；对基带调制信号sbase(t)进行数字低通滤波，生成数字滤波后的基带导航信号
Figure FDA0002401446200000011
(2)根据阵列天线形成的波束数为Nbeam，生成M个跳时伪码序列，
Figure FDA0002401446200000012
(3)读取数字波束成形系数表，获得N个加权系数；将N个加权系数分别与步骤(1)数字滤波后的基带导航信号
Figure FDA0002401446200000013
相乘，得到加权后的N路数字导航信号，分别为Sn(t),n＝1,2,…,N，完成数字波束成形；
(4)将步骤(3)加权后的N路数字导航信号，依次经过DAC、滤波、上变频、功率放大后，经N个天线阵元播发，实现码片级脉冲跳时导航信号生成与播发。
2.根据权利要求1所述的一种基于数字波束成形的码片级脉冲跳时导航信号生成与播发实现方法，其特征在于：步骤(1)中所述的数字滤波后的基带导航信号
Figure FDA0002401446200000014
生成，通过以下方法得到：
(1)信号伪码发生器生成信号伪码序列记为{cs}，cs∈{0,1}，s＝0,1,2,…，信号伪码序列的码速率为Rchip，信号伪码序列的码片宽度为Tchip＝1/Rchip；对码片进行赋形，得到码波形，实现扩频码(即伪码)波形生成；
码波形为：
Figure FDA0002401446200000015
式中，t表示时间，s表示码片的序号，(1-2cs)将码片值{0,1}映射为{1,-1}，
Figure FDA0002401446200000021
为信号伪码序列的码片波形，可以采用当前GNSS信号的任意一种波形；对于矩形码片波形，有：
Figure FDA0002401446200000022
(2)将导航电文经过组帧、信号编码后，得到符号速率为Rd的电文数据{dl}，dl∈{1,-1}，l＝0,1,2,…，导航电文符号宽度Td＝1/Rd，且Td是Tchip的整数倍，将电文数据写为连续波形形式D(t)，即导航电文生成，则有：
Figure FDA0002401446200000023
式中，pd(t)为宽度Td的矩形函数：
Figure FDA0002401446200000024
(3)将连续波形形式的电文数据调制到伪码波形上，得到基带调制信号：
sbase(t)＝D(t)·Cs(t)；
(4)将sbase(t)进行FIR数字低通滤波，得到低通滤波后的信号
Figure FDA0002401446200000025
完成数字滤波后的基带导航信号生成。
3.根据权利要求1所述的一种基于数字波束成形的码片级脉冲跳时导航信号生成与播发实现方法，其特征在于：步骤(2)的跳时伪码序列生成，满足以下要求：
(1)跳时伪码发生器生成M个跳时伪码序列，分别为
Figure FDA0002401446200000026
Figure FDA0002401446200000027
Nbeam为阵列天线形成的波束数；
(2)M个跳时伪码序列的码速率与信号伪码序列的码速率一样，都是Rchip；跳时伪码序列与信号伪码序列在时域上是对齐的；
(3)M个跳时伪码序列的第s个码片值的组合，决定了导航信号的第s个码片值的播发波束(即决定了导航信号的第s个码片值在哪个波束播发)。
4.根据权利要求1所述的一种基于数字波束成形的码片级脉冲跳时导航信号生成与播发实现方法，其特征在于：步骤(3)所述的数字波束成形，通过以下方法得到：
(1)数字波束成形系数表存储了加权系数；数字波束成形系数表是一个Nbeam×N二维表，记为{wb,k}，b＝1,2,……,Nbeam，k＝1,2,……,N；数字波束成形系数表通过数值优化生成，与阵列天线的布局以及波束增益方向图设计有关；使用该表中第nbeam行的系数进行数字波束成形，能够让信号在第nbeam个波束播发，nbeam＝1,2，……，Nbeam；
(2)使用M个跳时伪码序列，控制数字波束成形系数表输出加权系数，具体为：
在第s个码片时刻，M个跳时伪码序列的值分别为：
Figure FDA0002401446200000031
转换为十进制数，即
Figure FDA0002401446200000032
将NTH,s对Nbeam进行求模加1，得到属于{1,2,…,Nbeam}的数Ns，实现数字波束成形系数选择，即
Ns＝(NTH,s mod Nbeam)+1
在第s个码片时刻，信号在第Ns个波束播发，数字波束成形系数表输出第Ns行加权系数
Figure FDA0002401446200000033
Figure FDA0002401446200000034
是一个复系数，表示为实部和虚部的形式，有
Figure FDA0002401446200000035
Figure FDA0002401446200000036
是
Figure FDA0002401446200000037
的实部，
Figure FDA0002401446200000038
是
Figure FDA0002401446200000039
的虚部；
(3)将生成的基带导航信号
Figure FDA00024014462000000310
与N个加权系数分别相乘，实现数字波束成形，得到第s个码片的持续时间内N个加权后的数字导航信号，分别为：
Figure FDA0002401446200000041
实现数字波束成形。
5.根据权利要求1所述的一种基于数字波束成形的码片级脉冲跳时导航信号生成与播发实现方法，其特征在于：生成信号伪码和导航电文，进行基带调制，具体为：
生成信号伪码和导航电文，将导航电文调制到伪码上。
6.根据权利要求1所述的一种基于数字波束成形的码片级脉冲跳时导航信号生成与播发实现方法，其特征在于：数字波束成形时，阵列天线的阵元数为N，形成的波束数为Nbeam，数字波束成形系数表存储了加权系数。
7.根据权利要求1所述的一种基于数字波束成形的码片级脉冲跳时导航信号生成与播发实现方法，其特征在于：数字波束成形系数表是一个Nbeam×N二维表将生成基带导航信号。
8.根据权利要求1所述的一种基于数字波束成形的码片级脉冲跳时导航信号生成与播发实现方法，其特征在于：M个跳时伪码序列，控制不同的码片伪随机的在不同的波束里面播发。
9.根据权利要求1所述的一种基于数字波束成形的码片级脉冲跳时导航信号生成与播发实现方法，其特征在于：将N个加权系数分别与步骤(1)数字滤波后的基带导航信号
Figure FDA0002401446200000042
相乘，得到加权后的N路数字导航信号，具体如下：
在跳时伪码序列控制下，每一个信号伪码期间，选择N个加权系数，分别与基带导航信号相乘，得到加权后的N路数字导航信号，分别为Sn(t),n＝1,2,…,N。
10.根据权利要求1所述的一种基于数字波束成形的码片级脉冲跳时导航信号生成与播发实现方法，其特征在于：(4)将步骤(3)加权后的N路数字导航信号，依次经过DAC、滤波、上变频、功率放大后，经N个天线阵元播发，具体为
第n路数字导航信号Sn(t),n＝1,2,…,N经第n个DAC、第n个滤波器、第n个上变频器、第n个功率放大器、第n个天线阵元播发，实现码片级脉冲跳时导航信号生成与播发。
CN202010148020.9A 2020-03-05 2020-03-05 一种基于数字波束成形的码片级脉冲跳时导航信号生成与播发实现方法 Active CN111257911B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010148020.9A CN111257911B (zh) 2020-03-05 2020-03-05 一种基于数字波束成形的码片级脉冲跳时导航信号生成与播发实现方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010148020.9A CN111257911B (zh) 2020-03-05 2020-03-05 一种基于数字波束成形的码片级脉冲跳时导航信号生成与播发实现方法

Publications (2)

Publication Number Publication Date
CN111257911A true CN111257911A (zh) 2020-06-09
CN111257911B CN111257911B (zh) 2022-08-12

Family

ID=70947601

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010148020.9A Active CN111257911B (zh) 2020-03-05 2020-03-05 一种基于数字波束成形的码片级脉冲跳时导航信号生成与播发实现方法

Country Status (1)

Country Link
CN (1) CN111257911B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111781621A (zh) * 2020-06-19 2020-10-16 西安空间无线电技术研究所 一种基于低轨宽带互联网星座的导航方法及系统
CN112821925A (zh) * 2021-01-21 2021-05-18 重庆邮电大学 一种基于大规模mimo天线阵列的混合波束赋形算法

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1991011732A1 (en) * 1990-01-30 1991-08-08 Nauchno-Issledovatelsky Institut Kosmicheskogo Priborostroenia Method and device for radio-navigational determinations using artificial earth-satellites
JP2002064331A (ja) * 2000-08-21 2002-02-28 Nec Corp マルチビーム空中線装置
US20150236751A1 (en) * 2014-02-19 2015-08-20 Airbus Ds Gmbh Receiver for Acquiring and Tracking Spread Spectrum Navigation Signals with Changing Subcarriers
CN106556849A (zh) * 2016-11-11 2017-04-05 深圳市天弓导航科技有限公司 一种卫星导航接收设备的数字波束形成方法
WO2017153985A1 (en) * 2016-03-07 2017-09-14 Satixfy Uk Limited Digital beam forming system and method
CN107431508A (zh) * 2015-03-03 2017-12-01 平流层平台有限公司 用于与用户天线进行通信的过程和装置
CN107666333A (zh) * 2017-01-20 2018-02-06 西南电子技术研究所（中国电子科技集团公司第十研究所） 扩/跳频信号体制阵列天线自适应波束形成系统
BR112017009539A2 (pt) * 2014-11-17 2018-07-03 Mediatek Inc arquitetura de transceptor para sistemas de múltiplas antenas
US10031234B1 (en) * 2015-04-03 2018-07-24 Interstate Electronics Corporation Global navigation satellite system beam based attitude determination
CN109116306A (zh) * 2018-07-26 2019-01-01 河海大学 多载波宽带信号的数字波束形成方法
CN110208822A (zh) * 2019-05-28 2019-09-06 西安空间无线电技术研究所 一种基于低轨移动通信卫星的通信方法
CN110703279A (zh) * 2019-09-16 2020-01-17 西安空间无线电技术研究所 一种基于码片级脉冲跳时的卫星导航信号生成方法

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1991011732A1 (en) * 1990-01-30 1991-08-08 Nauchno-Issledovatelsky Institut Kosmicheskogo Priborostroenia Method and device for radio-navigational determinations using artificial earth-satellites
JP2002064331A (ja) * 2000-08-21 2002-02-28 Nec Corp マルチビーム空中線装置
US20150236751A1 (en) * 2014-02-19 2015-08-20 Airbus Ds Gmbh Receiver for Acquiring and Tracking Spread Spectrum Navigation Signals with Changing Subcarriers
BR112017009539A2 (pt) * 2014-11-17 2018-07-03 Mediatek Inc arquitetura de transceptor para sistemas de múltiplas antenas
CN107431508A (zh) * 2015-03-03 2017-12-01 平流层平台有限公司 用于与用户天线进行通信的过程和装置
US10031234B1 (en) * 2015-04-03 2018-07-24 Interstate Electronics Corporation Global navigation satellite system beam based attitude determination
WO2017153985A1 (en) * 2016-03-07 2017-09-14 Satixfy Uk Limited Digital beam forming system and method
CN106556849A (zh) * 2016-11-11 2017-04-05 深圳市天弓导航科技有限公司 一种卫星导航接收设备的数字波束形成方法
CN107666333A (zh) * 2017-01-20 2018-02-06 西南电子技术研究所（中国电子科技集团公司第十研究所） 扩/跳频信号体制阵列天线自适应波束形成系统
CN109116306A (zh) * 2018-07-26 2019-01-01 河海大学 多载波宽带信号的数字波束形成方法
CN110208822A (zh) * 2019-05-28 2019-09-06 西安空间无线电技术研究所 一种基于低轨移动通信卫星的通信方法
CN110703279A (zh) * 2019-09-16 2020-01-17 西安空间无线电技术研究所 一种基于码片级脉冲跳时的卫星导航信号生成方法

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
U. GIROLA, A. PICCIRIELLO, D. VINCENZONI: "Smart Antenna Receiver Based on a Single Chip Solution for GSM/DCS Baseband Processing", 《 PROCEEDINGS DESIGN, AUTOMATION AND TEST IN EUROPE CONFERENCE AND EXHIBITION 2000 (CAT. NO. PR00537)》 *
ZHANGJIAN LU, YANSONG MENG: "A Short Range Two-way Frequency Transfer Experiment in Free Space", 《2019 IEEE 11TH INTERNATIONAL CONFERENCE ON COMMUNICATION SOFTWARE AND NETWORKS》 *
戴卫恒: "利用敌方军用导航信号实现定位的新方法", 《 第一届中国卫星导航学术年会论文集》 *
沈大海，蒙艳松: "基于低轨通信星座的全球导航增强系统", 《太赫兹科学与电子信息学报》 *
王华等: "一种基于现场可编程门阵列的卫星导航信号生成方法", 《导航定位学报》 *
王鹏毅等: "宽带扩频阵列信号处理技术研究", 《飞行器测控学报》 *
谢完成: "智能天线的波束形成算法与实现", 《通信技术》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111781621A (zh) * 2020-06-19 2020-10-16 西安空间无线电技术研究所 一种基于低轨宽带互联网星座的导航方法及系统
CN111781621B (zh) * 2020-06-19 2023-06-06 西安空间无线电技术研究所 一种基于低轨宽带互联网星座的导航方法及系统
CN112821925A (zh) * 2021-01-21 2021-05-18 重庆邮电大学 一种基于大规模mimo天线阵列的混合波束赋形算法

Also Published As

Publication number Publication date
CN111257911B (zh) 2022-08-12

Similar Documents

Publication Publication Date Title
Hassanien et al. Dual-function radar communication systems: A solution to the spectrum congestion problem
CN111257911B (zh) 一种基于数字波束成形的码片级脉冲跳时导航信号生成与播发实现方法
JP2763099B2 (ja) 衛星または遠隔中継器を使用するスペクトル拡散多元接続通信システム
Chen et al. Monocycle shapes for ultra wideband system
CA2224572C (en) Radio with peak power and bandwidth efficient modulation
TW392394B (en) A system for generating and using global radio frequency maps
US10312984B2 (en) Distributed airborne beamforming system
JP2005510169A (ja) ウルトラワイドバンド・アンテナアレイ
CN107666333A (zh) 扩/跳频信号体制阵列天线自适应波束形成系统
US9300354B2 (en) Receiver for acquiring and tracking spread spectrum navigation signals with changing subcarriers
US6486827B2 (en) Sparse frequency waveform radar system and method
CN110208822B (zh) 一种基于低轨移动通信卫星的通信方法
CN106680852B (zh) 基于卫星导航信号欺骗式干扰的定位方法及系统
CN110703279B (zh) 一种基于码片级脉冲跳时的卫星导航信号生成方法
WO2023168734A1 (zh) 数字全宽带直接转发式雷达应答器系统及其工作方法
CN113556205A (zh) 基于单边带时间调制阵列的物理层保密通信系统
US20190181563A1 (en) Per-element power control for array based communications
CA3072809C (en) Long range navigation system having direct sequence spread spectrum (dsss) rf signal
Lucente et al. Analysis and design of a point-to-point radio-link at W band for future satellite telecommunication experiments
Rabideau et al. Mitigation of digital array nonlinearities
Shi et al. A novel ionospheric oblique-incidence sounding network consisting of the ionospheric oblique backscatter sounder and the parasitic oblique-incidence sounder
Kamio et al. Adaptive sidelobe control for clutter rejection of atmospheric radars
CN113726699B (zh) 抗频偏的通信卫星天线实时校准方法及装置
Garmatyuk et al. Radar and GPS system inter-operability with UWB-OFDM signals
Floch et al. S-band signal design considering interoperability and spectral separation

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