一种基于码周期扩频码认证的导航信号安全性增强方法与流程

1.本发明涉及一种基于码周期扩频码认证的导航信号安全性增强方法，属于卫星导航领域。
背景技术：
：：2.全球卫星导航系统(gnss)正越来来广泛的参与到社会经济生产中，在民用领域中，广泛依赖gnss提供的开放服务信号进行定位导航与授时。一方面，以电力、金融、通信为代表的基础设施的授时主要依赖于gnss开放服务信号，据2019年美国商务部估计，gps系统30天丢失对美国经济的影响将为10亿美元/天，在特定的时间发生，则可能会高出50％。另一方面，以共享单车、在线点餐为代表的新兴领域，依赖于准确的定位实现经济利益的最大化。3.而伴随着gnss技术的发展与应用，针对gnss的干扰技术也越加成熟，其中，欺骗式干扰的威胁逐渐增加。欺骗攻击生成虚假的gnss信号，诱导被攻击者获得错误的定位授时解，甚至导致灾难性后果。gnss开放服务信号面临的欺骗攻击更加普遍，因为开放服务信号的信号体制是公开的，任何欺骗者都可以按照协议生成合规的gnss信号。4.因此，增强开放服务的安全性，是gnss发展的新需求，在gnss开放服务信号中引入安全认证特征正成为一种共识，通过进行gnss信号认证，接受者能够验证所接收的gnss信号是否来自真正的gnss卫星，而不是欺骗者。目前，研究者已针对galileo、gps、qzss信号认证进行了研究，信号认证可在电文和扩频码进行，分别称为导航电文认证(nma)和扩频码认证(sca)。5.导航电文认证是信息层面的认证，通常采用数字签名认证方案，选取导航电文中随时间变化的比特，作为参考认证导航电文数据，经密钥生成签名数据，与导航电文同时播发。接收机接收播发的导航电文、数字签名与密钥，进行导航电文认证。扩频码认证通常是在开放扩频码中，插入不可预测的扩频码序列，实现测量层面的认证，即对开放扩频码序列的部分码片进行了加密，为了保证不对非认证用户造成影响，需要提升信号的发射功率。6.相比于导航电文认证，扩频码认证具有更好的抗再生式欺骗攻击的能力。galileo将在e1b提供osnma信号认证服务，在e6c同时提供扩频码认证服务，并正在开始测试；gps计划在l1c上引入chimera(chips‑messagerobustauthentication)认证方案，同时实现扩频码认证与导航电文认证。7.然而，现有的扩频码认证方法，需要在原来开放信号的扩频码上引入认证特征，每个码周期都需要引入认证特征，为避免对非认证用户造成影响，引入认证特征的码周期功率要提升，导致功率增加。同时，现有的扩频码认证方法，是通过对开放扩频码的部分码片(5％或者10％)进行加密实现的，需要同时生成加密码片的位置和对应的码片值，信号生成相对复杂。此外，现有方法扩频码认证方法，由于是对部分码片进行加密，因此，在接收机端，需要存储每一个采样值，直到获取密钥进行认证，接收机实现扩频码认证所需的存储资源多。针对上述问题，本发明提出一种基于码周期扩频码认证的导航信号安全性增强方法。技术实现要素：8.本发明解决的技术问题是：克服现有技术的不足，针对卫星导航开放服务信号的安全认证需要，在不影响非认证用户使用的约束下，提出一种基于码周期扩频码认证的导航信号安全性增强方法，实现扩频码认证功能，使得欺骗攻击者难以通过接收该导航认证信号恢复出正确的认证信号，实现安全性增强。9.本发明解决技术的方案是：一种基于码周期扩频码认证的导航信号安全性增强方法，该方法中，导航卫星执行如下步骤：10.s1、生成码周期认证序列；所述码周期认证序列的长度等于历元的码周期数，第n个历元的码周期认证序列中第m个元素ccpa(n)(m)为0或者1；11.s2、当码周期认证序列中的元素为1时，生成扩频码认证序列csca(n)，并针对第n个历元的第m个码周期，根据扩频码认证序列生成加入扩频码认证特征的基带波形；否则，针对第n个历元的第m个码周期，生成常规的导航基带波形；12.所述扩频码认证序列的长度等于扩频码长度nperiod·ls，nperiod为历元内的码周期数，ls为一个扩频码周期内的码片数，扩频码认证序列中第l个元素为0或者1，当扩频码认证序列中的元素为1时，第n个历元第m个码周期的第l个码片的幅度为第一幅度a1，否则，第n个历元第m个码周期的第l个码片的幅度为第二幅度a2，n∈[1,n],m∈[1,nperiod],l∈[1,ls]，n为导航电文中的历元数，即：[0013][0014]其中，a1为第一幅度和a2为第二信号幅度；[0015]s3、将导航电文调制到基带波形上，得到基带信号；[0016]s4、将基带信号然后调制到数字中频上，得到数字中频调制信号；[0017]s5、将数字中频调制信号经数模转换，变为模拟中频信号，再经上变频、滤波、放大后得到导航增强射频信号，将导航增强射频信号经天线播发。[0018]所述第一幅度a1为和第二信号幅度a2的计算公式为：[0019][0020]其中，γchip为加入扩频码认证特征的码片数占码周期内码片数的比例，β为正常码周期的平均功率与认证码周期的平均功率之比，0.5≤γchip≤β。[0021]所述码周期认证序列通过如下方法得到；[0022]s1.1a、生成码周期认证序列查找表tablecpa，所述扩频码认证序列的查找表tablecpa是一个二维查找表，大小为rowcpa行，colcpa列，表中的每一行都是由{‘0’‘1’}中的元素组成的序列，且每一行中‘1’所占的比例都是γcode，每一行的序列不相等；[0023]s1.2a、随机生成码周期认证序列查找表tablecpa的行序号，根据行序号提取码周期认证序列查找表tablecpa中的一行；[0024]s1.3a、重复步骤1.2a，直到得到长度为nperiod的码周期认证序列。[0025]步骤1.2a采用如下方法生成码周期认证序列查找表tablecpa的行序号：[0026]s1.2.1a、采用标准哈希算法，根据导航电文数据生成第n个历元的哈希序列hash(n)∈{0,1}，长度为nhash比特；[0027]s1.2.2a、第n个历元的哈希序列的nhash比特按照比特划分为组，每比特对应的数值即为码周期认证序列查找表tablecpa的行序号，用于对应提取码周期认证序列查找表tablecpa中的一行，[0028]所述扩频码认证序列根据如下方法确定：[0029]s1.1b、生成扩频码认证序列查找表tablesca，所述扩频码认证序列的查找表tablesca是一个二维查找表，大小为rowsca行，colsca列，表中的每一行都是由{‘0’‘1’}中的元素组成的序列，且每一行中‘1’所占的比例都是γchip，每一行的序列不相等；[0030]s1.2b、随机生成扩频码认证序列查找表tablesca的行序号，根据行序号提取扩频码认证序列查找表tablesca中的一行；[0031]s1.3b、重复步骤1.2b，直到得到长度为ls的扩频码认证序列。[0032]所述步骤1.2b采用如下方法生成扩频码认证序列查找表tablesca的行序号：[0033]s1.2.1b、采用标准哈希算法，根据导航电文数据生成第n个历元的哈希序列hash(n)∈{0,1}，长度为nhash比特；[0034]s1.2.2b、第n个历元的哈希序列的nhash比特按照比特划分为组，每比特对应的数值即为扩频码认证序列查找表tablesca的行序号，用于对应提取扩频码认证序列查找表tablesca中的一行，[0035]用户终端执行如下步骤：[0036]r1‑1、导航增强信号经接收天线接收，经低噪放放大后进行下变频、滤波，变为模拟中频信号；[0037]r1‑2、模拟中频信号经采样和模数转换后变为数字中频信号；[0038]r1‑3、对数字信号接收处理，完成捕获、跟踪、解调与数据存储，得到导航电文与导航增强信号；[0039]r1‑4、对导航电文与导航增强信号进行扩频码认证，如果扩频码认证通过，则认为收到的导航增强信号是真实的，如果扩频码认证未通过则认为导航增强信号可能为欺骗信号。[0040]所述步骤r4扩频码认证的方法为：[0041]r4.1、跟踪稳定后，剥离载波，得到同相支路的基带数据基带数据按照一个码片宽度ts进行累加，得到第n个历元的第m个码周期的第w个码片累积值的均值[0042]r4.2、存储第n个历元的所有nperiod个码周期的所有码片累加值的均值1≤m≤nperiod，1≤w≤ls，表示第m个码周期的第w个码片样点数量；[0043]r4.3、将认证码周期的码片累积值的均值，分为两组：[0044]第一组：[0045]第二组：[0046]r4.4、对第一组码片累加值均值求平均，得到第一累加平均值ip1(n)；将第二组码片累加值均值求平均，得到第二累加平均值[0047]r4.5、计算第一累加平均值与第二累加平均值之差如果th为检测门限值，则扩频码认证成功，反之则认证失败。[0048]所述检测门限值th按照如下方式确定：[0049]其中，pfa为给定的最大虚警概率。[0050]所述导航电文为导航认证电文，通过如下步骤生成：[0051]s2.1、生成第n个历元的导航电文数据导航电文数据包括同步头、周内秒、周记数、轨道参数、钟差参数、crc等，长度为nnav比特，[0052]s2.2、采用标准哈希算法，根据导航电文数据生成第n个历元的哈希序列hash(n)∈{0,1}，长度为nhash比特；[0053]s2.3、根据哈希序列hash(n)，采用非对称加密算法生成第n个历元的私钥及其对应的公钥为第n个历元私钥对应的公钥长度为nkey,public比特；[0054]s2.4、根据哈希序列hash(n)、第n个历元的私钥、第n个历元的私钥长度为nkey,private比特，生成第n个历元的数字签名ds(n)∈{0,1}，长度为nds比特；[0055]s2.4、将第n个历元的导航电文数据第n个历元的数字签名ds(n)，以及第n‑1个历元的组成第n个历元的导航认证电文长度为比特，电文符号速率为rnav。[0056]用户终端执行如下步骤：[0057]r2‑1、导航增强信号经接收天线接收，经低噪放放大后进行下变频、滤波，变为模拟中频信号；[0058]r2‑2、模拟中频信号经采样和模数转换后变为数字中频信号；[0059]r2‑3、对数字信号接收处理，完成捕获、跟踪、解调与数据存储；[0060]r2‑4、当接收到第n个历元的公钥后，进行导航电文认证，导航电文认证通过后，再进行扩频码认证；[0061]r2‑5、如果导航电文认证和扩频码认证都通过，则认证结果为第n个历元的导航增强信号是真实的，如果扩频码认证未通过或者导航电文认证未通过则认为第n个历元的导航增强信号可能为欺骗信号。[0062]所述导航电文认证步骤如下：[0063]r2‑4.1、使用第n个历元的公钥解密第n个历元的数字签名，得到[0064]r2‑4.2、采用步骤s2.2相同的哈希算法，对进行hash计算，得到[0065]r2‑4.3、如果则导航电文认证成功，否则，导航电文认证失败。[0066]本发明与现有技术相比的有益效果是：[0067](1)、本发明由于采用了码周期认证序列，实现了只在部分码周期引入认证特征，相比于现有的chimera方法，在不影响信号认证能力的前提下，降低了平均发生功率。[0068](2)、本发明的采用对开放扩频码的幅度进行调制的方法，相比于现有的osnma和chimera方法，在新的纬度实现了扩频码认证，因此，能够与现有扩频码认证方法结合使用，进一步提升安全性。[0069](3)、本发明由于采用了查找表结合哈希序列的方式，实现了不可预测的码周期认证序列和扩频码认证序列的生成，相比于现有的chimera方法，在不影响信号认证能力的前提下，简化了生成方法。[0070](4)、本发明由于采用对码片幅度进行调制的扩频码认证方法，不对码片值进行加密，实现了在接收机端进行认证只需要存储码片累加值均值，即一个码片存储一个值，相比于现有的chimera方法，在保持认证性能不变的前提下，大幅降低所需要的存储资源。附图说明[0071]图1为本发明实施例码周期扩频码认证的导航信号安全性增强方法；[0072]图2为本发明导航认证电文结构示意图；[0073]图3为本发明码周期扩频码认证导航信号基带波形示意图；[0074]图4为本发明扩频码认证中的概率密度函数曲线示意图；[0075]图5为本发明扩频码认证在不同载噪比下的检测概率。具体实施方式[0076]下面结合实施例对本发明作进一步阐述。[0077]为实现上述目的，本发明公开了一种基于码周期扩频码认证的导航信号安全性增强方法。[0078]一种基于码周期扩频码认证的导航信号安全性增强方法包括：[0079]码周期扩频码认证导航信号生成与播发。导航认证电文生成，生成第n个历元的导航认证电文；基带波形调制，加入码周期认证和扩频码认证特征，生成第n个历元的基带波形；电文调制与数字中频调制，将第n个历元的导航认证电文调制到第n个历元的基带波形上，然后调制到数字中频上；数字中频调制后的信号经dac、变为模拟中频信号，再经上变频、滤波、放大后变为模拟射频信号，再经天线播发。[0080]码周期扩频码认证导航信号接收与认证。天线接收放大、滤波、下变频、adc，发射的信号经接收天线接收，经低噪放放大后进行下变频、滤波，变为模拟中频信号，再经adc采样变为数字中频信号，接收到的第n个历元的数字中频信号记为tsample为采样速率，l表示第l个采样点；对第n个历元的采样信号进行数字信号接收处理，完成捕获、跟踪、解调与数据存储；当接收到第n个历元的公钥后，进行导航电文认证，导航电文认证通过后，再进行扩频码认证；最后输出认证结果，如果导航电文认证和扩频码认证都通过，则认证结果为第n个历元的信号是真实的，如果扩频码认证未通过或者导航电文认证未通过则表明可能存在欺骗信号。[0081]上述方法也可以采用普通的导航信号，导航卫星只做码周期和扩频码认证：[0082]生成码周期认证序列；所述码周期认证序列的长度等于历元的码周期数，第n个历元的码周期认证序列中第m个元素ccpa(n)(m)为0或者1；[0083]当码周期认证序列中的元素为1时，生成扩频码认证序列csca(n)，并针对第n个历元的第m个码周期，根据扩频码认证序列生成加入扩频码认证特征的基带波形；否则，针对第n个历元的第m个码周期，生成常规的导航基带波形；[0084]用户终端执行如下步骤：[0085]r1‑1、导航增强信号经接收天线接收，经低噪放放大后进行下变频、滤波，变为模拟中频信号；[0086]r1‑2、模拟中频信号经采样和模数转换后变为数字中频信号；[0087]r1‑3、对数字信号接收处理，完成捕获、跟踪、解调与数据存储，得到导航电文与导航增强信号；[0088]r1‑4、对导航电文与导航增强信号进行扩频码认证，如果扩频码认证通过，则认为收到的导航增强信号是真实的，如果扩频码认证未通过则认为导航增强信号可能为欺骗信号。[0089]步骤1中(1)所述导航认证电文生成，包括如下步骤：[0090]生成第n个历元的导航电文数据导航电文数据包括同步头、周内秒、周记数、轨道参数、钟差参数、crc等，长度为nnav比特，[0091]采用标准哈希算法，可以是国内标准sm3密码杂凑算法，根据导航电文数据生成第n个历元的哈希序列hash(n)∈{0,1}，长度为nhash比特；[0092]根据哈希序列hash(n)，采用非对称加密算法，可以是rsa或者ecdsa算法，根据哈希序列hash(n)，采用非对称加密算法生成第n个历元的私钥及其对应的公钥为第n个历元私钥对应的公钥为对应的公钥为长度为nkey,public比特。[0093]使用第n个历元的私钥使用第n个历元的私钥长度为nkey,private比特，生成第n个历元的数字签名ds(n)∈{0,1}，长度为nds比特。[0094]将第n个历元的导航电文数据第n个历元的数字签名ds(n)，以及第n‑1个历元的组成第n个历元的导航认证电文长度为比特，电文符号速率为rnav。[0095]步骤1中(1)所述的基带波形生成，包括如下步骤：[0096]信号扩频码生成。生成信号的扩频码序列cs∈{1,0}，扩频码长度为ls，扩频码速率为rs，码片宽度为ts＝1/rs，扩频码周期为tcode＝ls/rs。则第n个历元包括的码周期数为。则第n个历元包括的码周期数为是一个整数，表示一个电文符号内的码周期数。[0097]码周期认证序列生成。根据hash(n)，生成第n个历元的码周期认证序列ccpa(n)∈{1,0}，ccpa(n)长度为nperiod。当码周期认证序列为0时，表示正常码周期，当码周期认证序列为1时，表示认证码周期，认证码周期所占比例，等于码周期认证序列中为1的比例，记为γcode，0≤γcode≤1。特殊情况，γcode＝0表示无认证码周期，退化为传统的卫星导航信号，γcode＝1表示全部都是认证码周期。[0098]扩频码认证序列生成。根据hash(n)，生成第n个历元的扩频码认证序列csca(n)∈{1,0}，csca(n)长度为nperiod·ls。在认证码周期中，当扩频码认证序列为1时，信号幅度为a1，当扩频码认证序列为0时，信号幅度为a2。扩频码认证序列中为1的比例记为γchip。可选择a1≥1，0≤a2≤1，则1≥γchip≥0.5。[0099]基带波形调制。对于第n个历元的第m个码周期，如果是正常码周期，如ccpa(n)(m)＝0，则基带波形调制为：[0100][0101]其中，pchip(t)是码片波形，t表示数字域时间。[0102]如果是认证码周期，如ccpa(n)(m)＝1，则基带波形调制为：[0103][0104]其中，第n个历元第m个码周期的第l个码片的幅度为：[0105][0106]也就是说，所述扩频码认证序列的长度等于扩频码长度nperiod·ls，nperiod为历元内的码周期数，ls为一个扩频码周期内的码片数，扩频码认证序列中第l个元素为0或者1，当扩频码认证序列中的元素为1时，第n个历元第m个码周期的第l个码片的幅度为第一幅度a1，否则，第n个历元第m个码周期的第l个码片的幅度为第二幅度a2，n∈[1,n],m∈[1,nperiod],l∈[1,ls]，n为导航电文中的历元数，即：[0107][0108]其中，a1为第一幅度和a2为第二信号幅度；[0109]式中a1和a2的值由波形调制参数确定方法确定。[0110]步骤1中(1)所述的电文调制与数字中频调制，包括如下步骤：[0111]将基带波形调制后的信号调制电文，第n个历元，电文调制后得到的信号为：[0112][0113]将电文调制后得到的信号调制到数字中频fif，得到第n个历元的数字中频信号：[0114][0115]步骤2中2)所述的码周期认证序列生成，按照如下方式生成：[0116]s1.1a、生成码周期认证序列查找表tablecpa，所述扩频码认证序列的查找表tablecpa是一个二维查找表，大小为rowcpa行，colcpa列，表中的每一行都是由{‘0’‘1’}中的元素组成的序列，且每一行中‘1’所占的比例都是γcode，每一行的序列不相等；[0117]s1.2a、随机生成码周期认证序列查找表tablecpa的行序号，根据行序号提取码周期认证序列查找表tablecpa中的一行；[0118]s1.3a、重复步骤1.2a，直到得到长度为nperiod的码周期认证序列。[0119]步骤1.2a采用第n个历元得到的哈希序列hash(n)生成码周期认证序列查找表tablecpa的行序号：[0120]s1.2.1a、采用标准哈希算法，根据导航电文数据生成第n个历元的哈希序列hash(n)∈{0,1}，长度为nhash比特；[0121]s1.2.2a、第n个历元的哈希序列的nhash比特按照比特划分为组，每比特对应的数值即为码周期认证序列查找表tablecpa的行序号，用于对应提取码周期认证序列查找表tablecpa中的一行，[0122]第n个历元得到的哈希序列hash(n)共有nhash比特，将nhash比特按照比特划分为每比特作为地址查表生成码周期认证序列。即哈希序列hash(n)的每比特，对应查找表中的一行colcpa比特，依次得到colcpa比特码周期认证序列。在第n个历元，得到的码周期认证序列长度为[0123]步骤2中3)所述的扩频码认证序列生成，按照如下方式生成：[0124]s1.1b、生成扩频码认证序列查找表tablesca，所述扩频码认证序列的查找表tablesca是一个二维查找表，大小为rowsca行，colsca列，表中的每一行都是由{‘0’‘1’}中的元素组成的序列，且每一行中‘1’所占的比例都是γchip，每一行的序列不相等；[0125]s1.2b、随机生成扩频码认证序列查找表tablesca的行序号，根据行序号提取扩频码认证序列查找表tablesca中的一行；[0126]s1.3b、重复步骤1.2b，直到得到长度为ls的扩频码认证序列。[0127]所述步骤1.2b采用如下方法生成扩频码认证序列查找表tablesca的行序号：[0128]s1.2.1b、采用标准哈希算法，根据导航电文数据生成第n个历元的哈希序列hash(n)∈{0,1}，长度为nhash比特；[0129]s1.2.2b、第n个历元的哈希序列的nhash比特按照比特划分为组，每比特对应的数值即为扩频码认证序列查找表tablesca的行序号，用于对应提取扩频码认证序列查找表tablesca中的一行，[0130]利用第n个历元得到的哈希序列hash(n)生成扩频码认证序列。在第n个历元，得到的哈希序列hash(n)是随机的，共有nhash比特。将nhash比特按照比特划分为组，每比特作为地址查表生成扩频码认证序列。即哈希序列hash(n)的每比特，对应查找表中的一行colsca比特，依次得到colsca比特码周期认证序列。[0131]在第n个历元，共有nperiod个码周期，每个码周期，扩频码认证序列长度为ls。第n个历元得到的扩频码认证序列长度为历元得到的扩频码认证序列长度为[0132]步骤2中4)所述的波形调制参数确定方法，按照如下方式进行：[0133]考虑到不影响非授权用户，要求非授权用户使用正常码周期的信号扩频码，与接收信号的相关值保持不变，则有[0134]a1γchip a2(1‑γchip)＝1[0135]认证码周期相比正常码周期平均发射功率增加，正常码周期的平均功率与认证码周期的平均功率之比β，0<β≤1，则[0136][0137]设置β和γchip的值，且β和γchip的值满足：[0138][0139]求解得到a1和a2的值，即[0140][0141]其中，γchip为加入扩频码认证特征的码片数占码周期内码片数的比例，β为正常码周期的平均功率与认证码周期的平均功率之比，0.5≤γchip≤β。[0142]步骤1中(2)所述的数字信号接收处理，包括如下步骤：[0143]对数字中频信号，接收机复现正常码周期信号，按照常规的捕获、跟踪、解调方法，对接收的信号进行捕获、载波跟踪与码跟踪，并解调出导航认证电文。[0144]存储第n个历元解调出的导航认证电文为共有nauthentication比特。包括第n个历元的导航电文数据长度为nnav，第n个历元的数字签名共有nds比特，以及第n‑1个历元的公钥长度为nkey,public。[0145]跟踪稳定后，剥离载波，得到同相支路的基带数据基带数据按照一个码片宽度ts进行累加。第n个历元的第m个码周期的第w个码片累积值的均值为：[0146][0147](m‑1)tcode (w‑1)ts≤l·tsample<(m‑1)tcode wts[0148]存储第n个历元的所有nperiod个码周期的所有码片累加值的均值1≤m≤nperiod，1≤w≤ls，表示第m个码周期的第w个码片样点数量。[0149]步骤1中(2)所述的导航电文认证，包括如下步骤：[0150]在n 1个历元，解调得到第n个历元的公钥[0151]使用第n个历元的公钥解密第n个历元的数字签名，得到[0152]采用相同的哈希算法，对进行hash计算，得到[0153]如果则导航电文认证成功，否则，导航电文认证失败。[0154]步骤1中(2)所述的扩频码认证，包括如下步骤：[0155]根据和码周期认证序列的查找表tablecpa，生成第n个历元的码周期认证序列ccpa(n)∈{1,0}，ccpa(n)长度为nperiod。[0156]根据和扩频码认证码的查找表tablesca，生成第n个历元的扩频码认证序列csca(n)∈{1,0}，csca(n)长度为nperiod·ls。[0157]将认证码周期的码片累积值的均值，分为两组：[0158]第一组：[0159]第二组：[0160]进行扩频码认证。对第一组码片累加值均值求平均，得到ip1(n)；将第二组码片累加值均值求平均，得到计算的值，如果th为检测门限值，则扩频码认证成功，反之则认证失败。[0161]步骤10中(4)所述的监测门限值，按照如下方式确定：[0162]对于真实的信号，在高斯白噪声信道下，ip1(n)满足高斯分布均值tobs是选择的认证码周期的信号长度。满足高斯分布均值均值表示正常码周期的载噪比。此时，满足高斯分布满足高斯分布满足高斯分布概率密度为[0163]对于可能的欺骗信号，此时，ip1(n)满足高斯分布均值均值均值满足高斯分布均值此时，满足高斯分布概率密度为[0164]对于给定的最大虚警概率为pfa，根据可求得门限值th。[0165]本发明从码周期层面和扩频码层面分别引入认证特性，实现高安全的信号认证服务。将开放服务信号的码周期划分为非认证码周期与认证码周期，由码周期认证序列决定是否为认证码周期，只要在认证码周期才添加扩频码认证特性。扩频码认证特性不是对码片值进行加密，而是对扩频码幅度进行二进制幅度调制。相对于非认证码周期的码片幅度为1，认证码周期的一部分码片的幅度增加到a1，另一部分码片的幅度减少到a2，码片的幅度值由扩频码认证序列决定。a1和a2的取值保证非认证用户接收到的信号功率保持不变，不对非认证用户造成影响。由于没有对码片值进行加密，因此可以和现有的导航电文认证(nma)和扩频码认证(sca)方法结合使用，进一步提升安全性。[0166]实施例：[0167]为使本发明的目的、技术方案及优点更加清楚明白，以下参考附图，对本发明某一具体实施例进一步详细说明。[0168]本发明公开的基于码周期扩频码认证的导航信号安全性增强方法如图1所示，包括：[0169](a)码周期扩频码认证导航信号生成与播发。[0170](1)导航认证电文生成。[0171]1)生成第n个历元的导航电文数据导航电文数据包括同步头、周内秒、周记数、轨道参数、钟差参数、crc等，长度为nnav比特，[0172]2)利用标准哈希算法，可以是国内标准sm3密码杂凑算法，根据导航电文数据，生成第n个历元的哈希序列hash(n)∈{0,1}，长度为nhash比特；[0173]3)根据哈希序列hash(n)，采用非对称加密算法，如rsa或者ecdsa算法，使用第n个历元的私钥历元的私钥长度为nkey,private比特，生成第n个历元的数字签名ds(n)∈{0,1}，长度为nds比特。第n个历元私钥对应的公钥为公钥为长度为nkey,public比特。[0174]4)将第n个历元的导航电文数据第n个历元的数字签名ds(n)，以及第n‑1个历元的组成第n个历元的导航认证电文长度为比特，电文符号速率为rnav。图2给出了导航电文数据的结构示意图。[0175](2)基带波形生成。[0176]1)信号扩频码生成。生成信号的扩频码序列cs∈{1,0}，扩频码长度为ls，扩频码速率为rs，码片宽度为ts＝1/rs，扩频码周期为tcode＝ls/rs。则第n个历元包括的码周期数为是一个整数，表示一个电文符号内的码周期数。[0177]2)码周期认证序列生成。根据hash(n)，生成第n个历元的码周期认证序列ccpa(n)∈{1,0}，ccpa(n)长度为nperiod。[0178]当码周期认证序列为0时，表示正常码周期，当码周期认证序列为1时，表示认证码周期，认证码周期所占比例，等于码周期认证序列中为1的比例，记为γcode，0≤γcode≤1。特殊情况，γcode＝0表示无认证码周期，退化为传统的卫星导航信号，γcode＝1表示全部都是认证码周期。码周期认证序列具有不可预测性，可采用基于查找表的方式，根据哈希序列hash(n)生成，0,1所占的比例可重置。[0179]码周期认证序列的查找表tablecpa是一个二维查找表，大小为rowcpa行，colcpa列。表中的每一行都是由{‘0’‘1’}组成的序列，且每一行中‘1’所占的比例都是γcode，每一行的序列不相等。例如，表为[0180]表1码周期认证序列的查找表示意[0181][0182]在第n个历元，得到的哈希序列hash(n)是随机的，共有nhash比特。将nhash比特按照比特划分为每比特作为地址查表生成码周期认证序列。即哈希序列hash(n)的每比特，对应查找表中的一行colcpa比特，依次得到colcpa比特码周期认证序列。在第n个历元，得到的码周期认证序列长度为[0183]3)扩频码认证序列生成。根据hash(n)，生成第n个历元的扩频码认证序列csca(n)∈{1,0}，csca(n)长度为nperiod·ls。[0184]在认证码周期中，当扩频码认证序列为1时，信号幅度为a1，当扩频码认证序列为0时，信号幅度为a2。扩频码认证序列中为1的比例记为γchip。可选择a1≥1，0≤a2≤1，则1≥γchip≥0.5。扩频码认证序列具有不可预测性，也可以采用基于查找表的方式，根据哈希序列hash(n)生成，0,1所占的比例可重置。[0185]扩频码认证序列的查找表tablesca是一个二维查找表，大小为rowsca行，colsca列。表中的每一行都是由{‘0’‘1’}组成的序列，且每一行中‘1’所占的比例都是γchip，每一行的序列不相等。例如，表的前4行的示意如下表所示。[0186]表3扩频码认证序列的查找表示意[0187][0188]在第n个历元，得到的哈希序列hash(n)是随机的，共有nhash比特。将nhash比特按照比特划分为组，每比特作为地址查表生成扩频码认证序列。即哈希序列hash(n)的每比特，对应查找表中的一行colsca比特，依次得到colsca比特码周期认证序列。[0189]在第n个历元，共有nperiod个码周期，每个码周期，扩频码认证序列长度为ls。第n个历元得到的扩频码认证序列长度为因此，要求[0190]4)基带波形调制[0191]对于第n个历元的第m个码周期，如果是正常码周期，如ccpa(n)(m)＝0，则基带波形调制为：[0192][0193]其中，pchip(t)是码片波形，t表示数字域时间。[0194]如果是认证码周期，如ccpa(n)(m)＝1，则基带波形调制为：[0195][0196]其中，第n个历元第m个码周期的第l个码片的幅度为：[0197][0198]图3给出了基带波形的示意图。[0199]波形调制参数确定方法如下：考虑到不影响非授权用户，要求非授权用户使用正常码周期的信号扩频码，与接收信号的相关值保持不变，则有[0200]a1γchip a2(1‑γchip)＝1[0201]认证码周期相比正常码周期平均发射功率增加，正常码周期的平均功率与认证码周期的平均功率之比β，0<β≤1，则[0202][0203]设置β和γchip的值，则可以求解得到a1和a2的值，即[0204][0205]且β和γchip的值满足：[0206][0207]5)电文调制与数字中频调制[0208]将基带波形调制后的信号调制电文，第n个历元，电文调制后得到的信号为：[0209][0210]将电文调制后得到的信号调制到数字中频fif，得到第n个历元的数字中频信号：[0211][0212]6)dac、上变频、滤波、放大、播发[0213]数字中频信号经dac变为模拟中频信号，再经上变频、滤波、放大后变为模拟射频信号，再经天线播发。[0214](b)码周期扩频码认证导航信号生成与播发。[0215](1)发射的信号经接收天线接收，经低噪放放大后进行下变频、滤波，变为模拟中频信号，再经adc采样变为数字中频信号，接收到的第n个历元的数字中频信号记为tsample为采样速率，l表示第l个采样点。[0216](2)数字信号接收处理。[0217]对数字中频信号，接收机复现正常码周期信号，按照常规的捕获、跟踪、解调方法，对接收的信号进行捕获、载波跟踪与码跟踪，并解调出导航认证电文。[0218]存储第n个历元解调出的导航认证电文为共有nauthentication比特。包括第n个历元的导航电文数据长度为nnav，第n个历元的数字签名共有nds比特，以及第n‑1个历元的公钥长度为nkey,public。[0219]跟踪稳定后，剥离载波，得到同相支路的基带数据基带数据按照一个码片宽度ts进行累加。第n个历元的第m个码周期的第w个码片累积值的均值为：[0220][0221]存储第n个历元的所有nperiod个码周期的所有码片累加值的均值1≤m≤nperiod，1≤w≤ls，表示第m个码周期的第w个码片样点数。[0222](3)导航电文认证。[0223]1)在n 1个历元，解调得到第n个历元的公钥[0224]2)使用第n个历元的公钥解密第n个历元的数字签名，得到[0225]3)采用相同的哈希算法，对进行hash计算，得到[0226]4)如果则导航电文认证成功，否则，导航电文认证失败。[0227](4)扩频码认证。导航电文认证成功后，则可以进行扩频码认证。[0228]1)根据和码周期认证序列的查找表tablecpa，生成第n个历元的码周期认证序列ccpa(n)∈{1,0}，ccpa(n)长度为nperiod。[0229]2)根据和扩频码认证序列的查找表tablesca，生成第n个历元的扩频码认证序列csca(n)∈{1,0}，csca(n)长度为nperiod·ls[0230]3)将认证码周期的码片累积值的均值，分为两组[0231]第一组：[0232]第二组：[0233]4)进行扩频码认证。对第一组码片累加值均值求平均，得到ip1(n)；将第二组码片累加值均值求平均，得到[0234]5)计算的值，如果th为检测门限值，则扩频码认证成功，反之则认证失败。[0235](5)门限值th确定方法。[0236]1)对于真实的信号，在高斯白噪声信道下，ip1(n)满足高斯分布均值tobs是选择的认证码周期的信号长度。满足高斯分布均值均值表示正常码周期的载噪比。此时，满足高斯分布满足高斯分布满足高斯分布概率密度为[0237]2)对于可能的欺骗信号，若采用再生式欺骗，电文是真实的，但都是正常码周期，无认证码周期。此时，ip1(n)满足高斯分布均值均值均值满足高斯分布均值此时，满足高斯分布概率密度为tobs＝100ms，c/n0＝40db‑hz时的概率密度函数曲线示意图如图4所示。[0238]3)对于给定的虚警概率要求为pfa，根据可求得门限值th，并计算出对应的检测概率为：设定选择则不同载噪比下的检测概率如图5所示。[0239]本发明虽然已以较佳实施例公开如上，但其并不是用来限定本发明，任何本领域技术人员在不脱离本发明的精神和范围内，都可以利用上述揭示的方法和技术内容对本发明技术方案做出可能的变动和修改，因此，凡是未脱离本发明技术方案的内容，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化及修饰，均属于本发明技术方案的保护范围。当前第1页12当前第1页12