一种抗单音干扰的统一载波测控信号捕获方法技术

本发明专利技术提供一种抗单音干扰的统一载波测控信号捕获方法，包括以下步骤：计算基带信号能量中心坐标k

【技术实现步骤摘要】
一种抗单音干扰的统一载波测控信号捕获方法
本专利技术涉及一种抗单音干扰的统一载波测控信号捕获方法，特别适用于有抗单音干扰要求的捕获，属于航天测控通信

技术介绍
统一载波测控体制由于具有在一个主载波上通过调相调制多个测控通信副载波的特点，实现了跟踪测轨、遥测、遥控、测速和测距等通信功能，完成了多种测控使命，在航天测控通信中被广泛应用。信号的同步是实现信号解调的基础，而信号的捕获则是信号的粗同步过程，是实现信号精确同步的前提。根据信号调制特点，统一载波测控系统信号捕获的关键是对主载波进行捕获。早期比较成熟的统一载波测控系统如美国的Apollo登月系统和我国的统一S波段(USB)测控系统采用的是频率扫描捕获方法，该捕获方法对多普勒进行搜索，搜索步进受限于环路滤波器带宽，扫描速度慢，捕获时间长。随着数字信号处理技术发展，信号的频谱特性得到更多关注，用来提高信号捕获性能。统一载波测控信号的频谱特性有两方面：一是调制度决定了信号能量的分配，调制度越大主载波能量越小，副载波能量越大；二是副载波信号都关于主载波呈对称分布。调制度小于1.4时调制信号的谱峰位于主载波处，可利用判最大值法对主载波进行捕获，但在大调制度下时不再适用。基于判频谱对称性的伪对称点最小误差判定法可以不受调制度的影响，但是算法需要对多普勒范围内的频点逐一判定对称性，计算量巨大，不利于工程实现。基于频谱能量中心的捕获方法计算量适中，但是只有在载波能量中心和频谱对称中心一致的前提下适用，当存在干扰时该前提不成立。目前统一载波测控体制处于军、商共用状态，技术体制公开，大部分卫星不具备抗干扰能力；且采用统一载波测控体制航天器的测控天线采用宽波束天线，无法实现空间过滤，极易受到干扰，造成信号捕获的误捕。其中人为的单音干扰是一种易于产生、干扰性强的压制式干扰。一旦单音干扰信号功率大于系统的抗干扰容限，系统捕获性能会急剧恶化，捕获时间延长、虚警概率增加，无法正确捕获到信号。所以研究具有单音干扰抑制功能的捕获算法很具实际意义。
技术实现思路
本专利技术提供一种抗单音干扰的统一载波测控系统捕获方法，该方法利用了能量中心和对称中心的位置关系确认单音干扰是否存在，在单音干扰存在的情况下，充分考虑测控信号调制度、幅值特性，依次对疑似单音干扰存在的频点进行判决和剔除，在判定单音干扰同时得到残余载波频点位置，进而实现主载波的捕获。实现本专利技术的技术方案如下：一种抗单音干扰的统一载波测控信号捕获方法，包括以下步骤：步骤一：对输入的基带信号做N点快速傅里叶变换，再进行模平方得到信号频谱幅值的平方E(k)，其中k为频点的坐标；然后计算信号能量中心坐标kc；步骤二：对基带信号的正弦副载波的调相指数mf进行判断，当1.4＜mf＜2时，直接选取出E(k)最大的2根谱线的坐标k1和k2，计算对称中心kc'＝(k1+k2)/2；当mf＜1.4时，则选出最大的1根谱线的坐标作为kc'；步骤三：若kc'＝kc，kc就是主载波的位置，进入步骤六；若kc'≠kc，进行步骤四；步骤四：当1.4≤mf＜2时，选取出E(k)最大的6根谱线，记录下坐标为p1～p6；当mf＜1.4时，则选出最大的4根谱线，记录下坐标为p1～p4；设i的初始值为1；步骤五：假定坐标pi处为疑似单音干扰，从频谱上剔除，判决剔除后的谱线是否满足对称性：若满足，输出干扰频点的频谱，进入步骤六；若不满足，令i＝i+1重复步骤五，直到满足对称性为止；步骤六：频谱对称中心的频率值即为载波频率，实现信号捕获。进一步地，本专利技术所述对称性的判决条件如下：若mf＜1.4，针对剔除疑似单音干扰后的3个频点按频率进行排序，分别记为pmin、pmedian、pmax：第一，用最大值pmax和最小值pmin计算出中值(pmin+pmax)/2，比较其与中间值pmedian是否一致；第二，比较E(pmin)和E(pmax)值是否近似相等；同时满足以上两个条件，则认为满足对称性；若1.4≤mf＜2，针对剔除疑似单音干扰后的5个频点按频率进行排序，分别记为pmin、psub_min、pmedian、psub_max、pmax：第一，用最大值pmax和最小值pmin计算出中值pestimate_1＝(pmin+pmax)/2，再用次大值psub_min和次小值psub_max计算出中值pestimate_2＝(psub_min+psub_max)/2，比较pestimate_1、pestimate_2和pmedian三者是否一致；第二，比较E(pmin)和E(pmax)、E(psub_min)和E(psub_max)的值是否分别近似相等；同时满足以上两个条件，则认为满足对称性。有益效果首先，本专利技术是基于能量中心的一种算法，所以它继承了该类算法对调制度的鲁棒性强的优势，相比最大值检测算法，本专利技术在大调制度下同样适用，相比伪对称点最小误差判定法，计算量小，实时性高。其次，统一载波测控信号频谱具有对称性，单音干扰会破坏该对称性，本专利技术利用这一特点对单音干扰进行检测和消除，能够实现单音干扰下目标信号主载波的成功捕获，相比传统能量中心算法，提高了捕获算法抗单音干扰性能。具体实施方式下面结合具体实例对本专利技术进行详细说明。下面对本专利技术的设计原理进行阐述：以上行链路信号为例，副载波为一个测距码和一个正弦波。在单音干扰和高斯白噪声影响下，输入的基带信号可表示为：s(t)＝Acos[ωct+bc(t)+mfsinωmt]+j(t)+v(t)(1)其中，A为载波幅度；ωc为载波角频率；b为测距码对载波的调相指数；c(t)为码波形，取值为±1；mf为正弦副载波的调相指数，取值范围为0～2rad；ωm为正弦副载波的角频率；j(t)＝Ajcos(ωjt)为单音干扰信号，此处只考虑一个频点的单音干扰，Aj为单音干扰幅值，ωj为单音干扰角频率；v(t)为高斯白噪声。因为测距码具有扩频特性，频谱湮没在噪声中，对判断对称性没有影响，所以此处令b＝0，看作没有测距码调制，以简化信号(1)形式。s(t)＝Acos[ωct+mfsinωmt]+j(t)+v(t)(2)将式(2)展开得到式中Jn(mf)为变量mf的第n阶贝塞尔函数；对式(3)做傅里叶变换，得到信号的频域表达式为：由式(4)可知，信号频谱由载波分量ωc和无数边频分量ωc±nωm组成。其中n＝0时只有载波分量ωc，幅值为AJ0(mf)；当n≠0时，边频分量ωc±nωm对称分布在载波两侧，幅度为AJn(mf)，相邻边频之间的间隔为ωm，所以USB信号频谱存在对称性。贝塞尔函数的性质，当n＞mf时，Jn(mf)的值迅速下降，所以当调制度mf≤1时，只考虑n＝1的边频分量；当mf＞1时，取边频数基于上述USB信号的频谱特性，本专利技术提出一种抗单音干扰的统一载波测控信号捕获方法，包括以下步骤：步骤一：对输入的基带信号做N点快速傅里叶变换(FFT)，再进行模平方得到信号频谱幅值的平方E(k)，其中k为频点的坐标。根据公式(5)算出信号能量中心坐标kc。步骤二：对基带信号的正弦副载波的调相指数mf进行判断，当1.4＜mf＜2时，峰值为第一边频的一对副载波，所以直接选取出幅值E(k)最大的2根谱线的坐标k1和k2，计算对称中心kc'＝(k1+k2)本文档来自技高网...

【技术保护点】
一种抗单音干扰的统一载波测控信号捕获方法，其特征在于，包括以下步骤：步骤一：对输入的基带信号做N点快速傅里叶变换，再进行模平方得到信号频谱幅值的平方E(k)，其中k为频点的坐标；然后计算信号能量中心坐标k

【技术特征摘要】
1.一种抗单音干扰的统一载波测控信号捕获方法，其特征在于，包括以下步骤：步骤一：对输入的基带信号做N点快速傅里叶变换，再进行模平方得到信号频谱幅值的平方E(k)，其中k为频点的坐标；然后计算信号能量中心坐标kc；步骤二：对基带信号的正弦副载波的调相指数mf进行判断，当1.4＜mf＜2时，直接选取出E(k)最大的2根谱线的坐标k1和k2，计算对称中心kc'＝(k1+k2)/2；当mf＜1.4时，则选出最大的1根谱线的坐标作为kc'；步骤三：若kc'＝kc，进入步骤六；若kc'≠kc，进入步骤四；步骤四：当1.4≤mf＜2时，选取出E(k)最大的6根谱线，记录下坐标为p1～p6；当mf＜1.4时，则选出最大的4根谱线，记录下坐标为p1～p4；设i的初始值为1；步骤五：假定坐标pi处为疑似单音干扰，从频谱上剔除，判决剔除后的谱线是否满足对称性：若满足，输出干扰频点的频谱，进入步骤六；若不满足，令i＝i+1重复步骤五，直到满足对称性为止；步骤六：将频谱对称中心的频率值确定为载波频率，实现信号捕获。2.根据权利要求1所述抗单音干扰的统一载波测控信号捕获方法，其特征在于，所述对称性的判决条件如下：若mf＜1.4，针对剔除疑似单音干扰后的3个频点按频率进行...

【专利技术属性】
技术研发人员：申宇瑶，吴芸芸，王永庆，吴嗣亮，于秀丽，马淑芬，刘东磊，
申请(专利权)人：北京理工大学，
类型：发明
国别省市：北京,11