CN114650108B - 一种变换域通信系统信号的检测方法及系统 - Google Patents

一种变换域通信系统信号的检测方法及系统 Download PDF

Info

Publication number
CN114650108B
CN114650108B CN202210517701.7A CN202210517701A CN114650108B CN 114650108 B CN114650108 B CN 114650108B CN 202210517701 A CN202210517701 A CN 202210517701A CN 114650108 B CN114650108 B CN 114650108B
Authority
CN
China
Prior art keywords
signal
communication system
value
autocorrelation
period
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
Application number
CN202210517701.7A
Other languages
English (en)
Other versions
CN114650108A (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.)
Naval University of Engineering PLA
Original Assignee
Naval University of Engineering PLA
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 Naval University of Engineering PLA filed Critical Naval University of Engineering PLA
Priority to CN202210517701.7A priority Critical patent/CN114650108B/zh
Publication of CN114650108A publication Critical patent/CN114650108A/zh
Application granted granted Critical
Publication of CN114650108B publication Critical patent/CN114650108B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B17/00Monitoring; Testing
    • H04B17/30Monitoring; Testing of propagation channels
    • H04B17/309Measuring or estimating channel quality parameters

Landscapes

  • Engineering & Computer Science (AREA)
  • Quality & Reliability (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Noise Elimination (AREA)
  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)

Abstract

本申请公开了一种变换域通信系统信号的检测方法及系统。该方法包括:按设定的采样率对接收的变换域通信系统信号进行采样,获取信号采样序列,求解所述信号采样序列的自相关函数,所述自相关函数的自变量为延迟时间,因变量为自相关值;将所述自相关函数的最大自相关值与比例因子相乘的结果作为谱峰判断阈值δ 1,进行谱峰判断和虚假谱峰抑制;计算每对相邻的两个谱峰对应的延迟时间的差值,将所述差值作为周期;统计每个周期出现的频次比,若频次比最大的周期占比超过预设阈值δ 2,则判断存在变换域通信系统信号,并将频次比最大的周期作为所述变换域通信系统信号的周期。本发明能够较为精准地检测变换域通信系统信号的周期。

Description

一种变换域通信系统信号的检测方法及系统
技术领域
本申请涉及通信技术领域,更具体地,涉及一种变换域通信系统信号的检测方法及系统。
背景技术
变换域通信系统(Transform Domain Communication System,TDCS)是指能够是主动感知周围的电磁频谱、规避干扰频段、利用未被占用的频段生成基础调制波形以实现通信的系统。因为其基础调制波形不含干扰频段,因而使得信号具有较强的抗干扰能力。目前对TDCS信号的检测技术还不成熟。TDCS信号本身是一种扩频信号,且产生方式与伪随机序列相关。现有常用的直接序列扩频(Direct Sequence Spread Spectrum,DSSS)信号的检测方法对于TDCS信号并不适用。例如,能量检测法虽然可以适用于TDCS信号的检测,但该方法有很大的缺陷:该方法前提是信道中除噪声外,若有则仅有TDCS信号,且该方法仅能检测TDCS信号的有无,并无法对TDCS信号的相关参数给予估计。循环谱法是针对单一载波的DSSS信号载频的估计,根据非零循环轴上谱线的有无判断DSSS信号的存在性,并根据最大值的位置估计载频。但是对TDCS信号的处理是基带信号的处理,不存在载频,因而循环谱法也不适用 TDCS信号的检测。
发明内容
针对现有技术的至少一个缺陷或改进需求,本发明提供了一种变换域通信系统信号的检测方法及系统,能够较为精准地检测变换域通信系统信号的周期。
为实现上述目的,按照本发明的第一个方面,提供了一种变换域通信系统信号的检测方法,包括:
按设定的采样率对接收信号进行采样,获取信号采样序列,求解所述信号采样序列的自相关函数,所述自相关函数的自变量为延迟时间,因变量为自相关值;
将所述自相关函数的最大自相关值与比例因子相乘的结果作为谱峰判断阈值δ 1,将大于所述谱峰判断阈值δ 1的自相关值作为谱峰;
计算每对相邻的两个谱峰对应的延迟时间的差值,将所述差值作为周期;
统计每个周期出现的频次比,若频次比最大的周期占比超过预设阈值δ 2,则判断所述接收信号包括变换域通信系统信号,并将频次比最大的周期作为所述变换域通信系统信号的周期。
进一步地,所述比例因子是与信噪比相关的函数,所述比例因子的函数采用仿真方法获得,具体包括:
设定信噪比和比例因子的变化范围,在信噪比和比例因子的变化范围内取M 1个信噪比和M 2个比例因子候选值,每个信噪比下仿真产生M 3个变换域通信系统信号;
对每个信噪比下的M 3个变换域通信系统信号,分别采用该M 2个比例因子候选值和所述的变换域通信系统信号的检测方法进行检测,判断检测的通信系统信号的周期是否正确,从该M 2个比例因子候选值中筛选周期检测正确的概率大于预设阈值δ 2的比例因子,获得M 1个信噪比下对应的M 1组比例因子;
根据M 1个信噪比下对应的M 1组比例因子拟合比例因子的函数。
进一步地,所述比例因子的函数为:
Figure 709545DEST_PATH_IMAGE001
其中,μ为比例因子,SNR为信噪比,snr(fs, np)为与信号采样率fs和采样点数np相关的函数,表示求在采样率为fs和采样点数为np的条件下,使得上述M 1×M 2×M 3次信号周期检测准确的概率大于预设阈值δ 2的最小信噪比,μ snr 为最小信噪比下的比例因子。
进一步地,当fs=1000,np≥256000时,所述比例因子的函数为:
Figure 829948DEST_PATH_IMAGE002
进一步地,所述计算每对相邻的两个谱峰对应的自相关值的差值前还包括:
执行虚假谱峰抑制步骤,所述虚假谱峰抑制步骤包括:取所有谱峰对应的时间间隔,若相邻的两个谱峰对应的时间间隔小于预设阈值δ 3,则舍弃该相邻的两个谱峰中自相关值较小的峰值,只保留该相邻的两个谱峰中自相关值较大的峰值;
重复所述虚假谱峰抑制步骤,直到所有谱峰间隔均大于预设阈值δ 3
进一步地,将重复所述虚假谱峰抑制步骤后谱峰的总个数记为L,所述只保留该相邻的两个谱峰中自相关值较大的峰值后,且所述计算每对相邻的两个谱峰对应的自相关值的差值前还包括步骤:
只保留自相关值最大的谱峰、其左边相邻的[L/6]个谱峰及其右边相邻的[L/6]个谱峰,[]表示取整。
进一步地,所述计算每对相邻的两个谱峰对应的自相关值的差值时,每个谱峰只参与一次差值计算。
按照本发明的第二个方面,还提供了一种变换域通信系统信号的检测系统,包括:
自相关函数求解模块,用于按设定的采样率对接收的变换域通信系统信号进行采样,获取信号采样序列,求解所述信号采样序列的自相关函数,所述自相关函数的自变量为延迟时间,因变量为自相关值;
谱峰判断模块,用于将所述自相关函数的最大自相关值与比例因子相乘的结果作为谱峰判断阈值δ 1,将大于所述谱峰判断阈值δ 1的自相关值作为谱峰;
周期计算模块,用于计算每对相邻的两个谱峰对应的延迟时间的差值,将所述差值作为周期;
输出模块,用于统计每个周期出现的频次比,若频次比最大的周期占比超过预设阈值δ 2,则判断所述接收信号包括变换域通信系统信号,并将频次比最大的周期作为所述变换域通信系统信号的周期。
总体而言,通过本发明所构思的以上技术方案与现有技术相比,能够取得下列有益效果:
(1)可以检测TDCS信号并对其基础调制波形的周期给予估计,特别是对于低信噪比的TDCS信号具有较好的检测精度。并且该方法在工程上容易实现。
(2)算法效果好且复杂度低,通过去除虚假谱峰解决虚假谱峰对周期判断的影响,通过减小谱峰搜索范围既减少了数据量提高了算法速度,又抛除了噪声影响较大的点,提高了算法准确度。
附图说明
为了更清楚地说明本申请实施例中的技术方案,下面将对实施例中所需使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本申请实施例提供的变换域通信系统信号的检测方法的流程示意图;
图2为本申请实施例提供的TDCS信号基函数时域图;
图3为本申请实施例提供的TDCS信号基函数s 0(n)自相关函数示意图;
图4为本申请实施例提供的TDCS信号基函数s 0(n)、s 1(n)互相关函数示意图;
图5为本申请实施例提供的自相关函数示意图;
图6为本申请实施例提供的比例因子在不同信噪比下检测正确率示意图;
图7为本申请实施例提供的自相关函数的左右侧局部放大示意图;
图8为本申请另一实施例提供的变换域通信系统信号的检测方法的流程示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
本申请的说明书和权利要求书及上述附图中包含了一系列步骤或单元的过程、方法、系统、产品或设备没有限定于已列出的步骤或单元,而是可选地还包括没有列出的步骤或单元,或可选地还包括对于这些过程、方法、产品或设备固有的其他步骤或单元。
如图1所示,本发明实施例的一种变换域通信系统信号的检测方法,包括:
S101,按设定的采样率对接收信号进行采样,获取信号采样序列,求解信号采样序列的自相关函数,自相关函数的自变量为延迟时间,因变量为自相关值。
本申请要解决的问题是检测接收信号中是否包含变换域通信系统信号,即TDCS信号,并且若存在TDCS信号,判断TDCS信号的周期。检测并估计TDCS信号的周期,具有重要意义。对于干扰方,侦察TDCS信号并对其周期等特征参数做出一定估计,能够有效地对TDCS实施干扰。其次,对于接收方,虽然采用TDCS信号协作通信,但是由于收发双方所处电磁环境的不同,会导致基础调制波形的差异,降低TDCS信号通信的性能,因而,这就要求接收方实现对TDCS信号的检测,并对接收的TDCS信号基函数周期等特征参数给予估计,从而对本地产生的基础调制波形做一定修正,提高相关解调的准确度,降低系统的误码率。
本申请中检测TDCS信号是基于自相关函数进行的。下面具体说明其原理。
TDCS信号的基函数时域图如图2所示。TDCS信号的码元0和码元1是由基础调制波形s 0(n)、s 1(n)表示:
Figure 812948DEST_PATH_IMAGE003
(1)
Figure 676998DEST_PATH_IMAGE004
(2)
N为基础调制波形的长度。A k 频谱感知的结果,有干扰的频点A k =0,空闲的频点A k =1。e为自然常数,j为虚数单位,n表示第n个采样点,φ(k)是伪随机相位。C为功率归一化因子,N A A k 中1的总数,即可用频点的个数,ε s 为一个码元的能量。
Figure 225791DEST_PATH_IMAGE005
(5)
由循环码移键控(Cyclic Code Shift Keying,CCSK)调制方式可知,s 1(n)是由s 0(n)的循环移位得来,故两者最大自相关值相等。
为了便于公式推导和公式简洁,令
Figure 782675DEST_PATH_IMAGE006
n为第n个采样点,φ(p)为p点的伪随机相位值,φ(q)为q点的伪随机相位值。为了简化分析,假定A k 均为1,则当τ=0时:
Figure 518550DEST_PATH_IMAGE007
(6)
Figure 451870DEST_PATH_IMAGE008
表示s 0(n)在τ=0处的自相关值。p,q表示第pq个采样点。考察函数Gj,函数周期T=N/|p-q|,因为p, q均为整数值,所以当采样点数为N时,包含了|p-q|个整数周期,故在一个周期N/|p-q|内的采样累加和为0,即公式(6)后两项均为0,又因为s 1(n)由s 0(n)循环移位而来,故两者的最大自相关值相等,继而得到:
Figure 855170DEST_PATH_IMAGE009
(9)
Figure 839348DEST_PATH_IMAGE010
表示s 1(n)在τ=0处的自相关值。
当τ=N/2时,因为s 1(n)是由s 0(n)的循环移位N/2得来,很显然:
Figure 593677DEST_PATH_IMAGE011
(10)
当τ≠0且τ≠N/2时,如图3和图4所示,无论是自相关值还是互相关值,相较于τ=0时的自相关值和τ=N/2互相关值来说,可以忽略不计。
综上,在只分析τ≥0情况时,多码元TDCS信号的自相关函数R ss (τ)可表示为:
Figure 65110DEST_PATH_IMAGE012
(11)
式中,R ss (τ)表示TDCS信号的自相关函数,num表示TDCS信号共有num个码元,R p+k-1,p (N/2)表示第p+k-1个码元与第p个码元在τ=N/2处的互相关值,R p,p+k (N/2) 表示第p个码元与第p+k个码元在τ=N/2处的互相关值,R p,p+k (0) 表示第p个码元与第p+k个码元在τ=0处的互相关值。
由如上分析我们可以得到这样的结论:一、TDCS信号的自相关函数在N/2的整数倍处出现谱峰且谱峰的大小与信息码元的分布有关;二、随着τ的增大,加数项减少,因而峰值幅度总体趋势是减小的。
基于上述分析,因此可以基于自相关函数进行TDCS信号检测。
将信号采样序列记为x(n),将x(n)延时τ后的信号记为x(n-τ),则信号采样序列的自相关函数R(τ)公式为:
Figure 322916DEST_PATH_IMAGE013
(1)
如图5所示,变换域系统信号的自相关函数在基础调制波形周期的整数倍处出现谱峰。
S102,将自相关函数的最大自相关值与比例因子相乘的结果作为谱峰判断阈值δ 1,将大于谱峰判断阈值δ 1的自相关值作为谱峰。
因为接收信号的最大自相关值的大小与其信号本身有关且最大相关值出现的位置固定,所以门限值设计以最大相关值R(0)作为参考。以自相关函数的最大值R(0)作为标准,乘上比例因子μ后,作为初步判断谱峰的阈值。
进一步地,比例因子是与信噪比相关的函数,比例因子的函数采用仿真方法获得,具体包括:
设定信噪比和比例因子的变化范围,在信噪比和比例因子的变化范围内取M 1个信噪比和M 2个比例因子候选值,每个信噪比下仿真产生M 3个变换域通信系统信号;
对每个信噪比下的M 3个变换域通信系统信号,分别采用该M 2个比例因子候选值和的变换域通信系统信号的检测方法进行检测,判断检测的通信系统信号的周期是否正确,从该M 2个比例因子候选值中筛选周期检测正确的概率大于预设阈值δ 2的比例因子,获得M 1个信噪比下对应的M 1组比例因子;
根据M 1个信噪比下对应的M 1组比例因子拟合比例因子的函数。
在一个实施例中,设定信噪比范围为-20dB至10dB,比例因子的变化范围μ:0.003-0.090,步进0.001。每个信噪比下1000个TDCS信号,信号的不同源于每次加性高斯噪声的不同。对每组信噪比下的1000个TDCS信号分别做信号检测,将这1000次检测得到的正确的周期占比取平均值。
如图6所示,取δ 2为0.5,黑色线以内白色区域正确周期占比均大于0.5,符合后续处理中阈值的要求。故只需在白色区域内找一条曲线连接信噪比最小处至最大处。显然,这种曲线有无数条,为了便于曲线的数学表达,取一条直线作为最简路径,如图中的黑色直线,其数学表达即比例因子μ的经验式:
Figure 487181DEST_PATH_IMAGE014
(2)
SNR为信噪比,snr(fs, np)为与信号采样率fs和采样点数np相关的函数,表示求在采样率为fs和采样点数为np的条件下,使得上述M 1×M 2×M 3次信号周期检测准确的概率大于预设阈值δ 2的最小信噪比,μ snr 为最小信噪比下的比例因子。
fs=1000,np≥256000时,比例因子的函数为:
Figure 463227DEST_PATH_IMAGE015
S103,执行虚假谱峰筛选步骤,所述虚假谱峰筛选步骤包括:取所有谱峰对应的时间间隔,若相邻的两个谱峰对应的时间间隔小于预设阈值δ 3,则舍弃该相邻的两个谱峰中自相关值较小的峰值,只保留该相邻的两个谱峰中自相关值较大的峰值。重复所述虚假谱峰筛选步骤,直到所有谱峰间隔均大于预设阈值δ 3
依靠单一的谱峰判断阈值δ 1确定谱峰,无法抑制虚假谱峰的现象。即一定区域内可能出现连续的峰值,这是不符合TDCS信号特点的,因而需要设定阈值δ 3,若连续的间隔δ 3内同时出现多个谱峰,只保留该区域内最大峰值的坐标。
S104,将只保留该相邻的两个谱峰中自相关值较大的峰值后谱峰的总个数记为L,只保留自相关值最大的谱峰、其左边相邻的[L/6]个谱峰及其右边相邻的[L/6]个谱峰,[]表示取整。
如图7所示,划定坐标范围原因有两个:一是随着延迟点数的增大,时域相关值减小,受噪声影响大。二是这样可以减小数据量,提高算法的效率。
S105,计算每对相邻的两个谱峰对应的延迟时间的差值,将差值作为周期。
在划定的坐标范围内,每两个谱峰取一次差值,作为TDCS信号基函数周期的估计,每个谱峰不重复参与计算,即第1个谱峰与第2个谱峰计算差值,第3个谱峰与第4个谱峰计算差值,而不是第1个谱峰与第2个谱峰计算差值,第2个谱峰与第3个谱峰计算差值。从而避免一个峰值的偏移,影响两次周期估计结果。
S106,统计每个周期出现的频次比,若频次比最大的周期占比超过预设阈值δ 2,则判断接收信号包括变换域通信系统信号,并将频次比最大的周期作为变换域通信系统信号的周期。
对计算出来的周期进行统计,得到每个周期的占比,也即检测正确的概率。频次比最大的周期占比将最大频次比与阈值比较,阈值δ 2设定0.5。大于阈值δ 2时,认为存在变换域通信系统信号,最大频次比对应的周期即为变换域通信系统信号基础调制波形的周期,小于阈值δ 2时,不存在变换域通信系统信号。
上述步骤S103、S104并非是必须的步骤,但是包含步骤S103、S104会使得本发明的算法效果更好,通过去除虚假谱峰解决虚假谱峰对周期判断的影响,通过减小谱峰搜索范围既减少了数据量提高了算法速度,又抛除了噪声影响较大的点,提高了算法准确度。
图8是本发明另一实施例提供的变换域通信系统信号的检测方法的流程示意图,如上述实施例的不同在于,进行虚拟峰值去除后,还需要先计算峰值数量,若峰值数量不大于预设阈值,则判断不存在TDCS信号,只有峰值数量大于预设阈值,则进行后续周期检测等步骤。
本发明实施例的一种变换域通信系统信号的检测系统,包括:
自相关函数求解模块,用于按设定的采样率对接收的变换域通信系统信号进行采样,获取信号采样序列,求解信号采样序列的自相关函数,自相关函数的自变量为延迟时间,因变量为自相关值;
谱峰判断模块,用于将自相关函数的最大自相关值与比例因子相乘的结果作为谱峰判断阈值δ 1,将大于谱峰判断阈值的自相关值作为谱峰;
周期计算模块,用于计算每对相邻的两个谱峰对应的延迟时间的差值,将差值作为周期;
输出模块,用于统计每个周期出现的频次比,若频次比最大的周期占比超过预设阈值δ 2,则判断接收信号包括变换域通信系统信号,并将频次比最大的周期作为变换域通信系统信号的周期。
进一步地,比例因子是与信噪比相关的函数,比例因子的函数采用仿真方法获得,具体包括:
设定信噪比和比例因子的变化范围,在信噪比和比例因子的变化范围内取M 1个信噪比和M 2个比例因子候选值,每个信噪比下仿真产生M 3个变换域通信系统信号;
对每个信噪比下的M 3个变换域通信系统信号,分别采用该M 2个比例因子候选值和的变换域通信系统信号的检测方法进行检测,判断检测的通信系统信号的周期是否正确,从该M 2个比例因子候选值中筛选周期检测正确的概率大于预设阈值δ 2的比例因子,获得M 1个信噪比下对应的M 1组比例因子;
根据M 1个信噪比下对应的M 1组比例因子拟合比例因子的函数。
进一步地,变换域通信系统信号的检测系统还包括谱峰筛选模块,用于在计算每对相邻的两个谱峰对应的自相关值的差值前,执行步骤:
执行虚假谱峰抑制步骤,虚假谱峰抑制步骤包括:取所有谱峰对应的时间间隔,若相邻的两个谱峰对应的时间间隔小于预设阈值δ 3,则舍弃该相邻的两个谱峰中自相关值较小的峰值,只保留该相邻的两个谱峰中自相关值较大的峰值;
重复虚假谱峰抑制步骤,直到所有谱峰间隔均大于预设阈值δ 3
将重复虚假谱峰抑制步骤后谱峰的总个数记为L,只保留该相邻的两个谱峰中自相关值较大的峰值后,且计算每对相邻的两个谱峰对应的自相关值的差值前还包括步骤:
只保留自相关值最大的谱峰、其左边相邻的[L/6]个谱峰及其右边相邻的[L/6]个谱峰,[]表示取整。
上述系统的实现原理和上述方法相同,此处不再赘述。
以上所述者,仅为本公开的示例性实施例,不能以此限定本公开的范围。即但凡依本公开教导所作的等效变化与修饰,皆仍属本公开涵盖的范围内。本领域技术人员在考虑说明书及实践这里的公开后,将容易想到本公开的其实施方案。本申请旨在涵盖本公开的任何变型、用途或者适应性变化,这些变型、用途或者适应性变化遵循本公开的一般性原理并包括本公开未记载的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的,本公开的范围和精神由权利要求限定。
以上实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (10)

1.一种变换域通信系统信号的检测方法,其特征在于,包括:
按设定的采样率对接收信号进行采样,获取信号采样序列,求解所述信号采样序列的自相关函数,所述自相关函数的自变量为延迟时间,因变量为自相关值;
将所述自相关函数的最大自相关值与比例因子相乘的结果作为谱峰判断阈值δ 1,将大于所述谱峰判断阈值δ 1的自相关值作为谱峰;
计算每对相邻的两个谱峰对应的延迟时间的差值,将所述差值作为周期;
统计每个周期出现的频次比,若频次比最大的周期占比超过预设阈值δ 2,则判断所述接收信号包括变换域通信系统信号,并将频次比最大的周期作为所述变换域通信系统信号的周期。
2.如权利要求1所述的变换域通信系统信号的检测方法,其特征在于,所述比例因子是与信噪比相关的函数,所述比例因子的函数采用仿真方法获得,具体包括:
设定信噪比和比例因子的变化范围,在信噪比和比例因子的变化范围内取M 1个信噪比和M 2个比例因子候选值,每个信噪比下仿真产生M 3个变换域通信系统信号;
对每个信噪比下的M 3个变换域通信系统信号,分别采用该M 2个比例因子候选值和所述的变换域通信系统信号的检测方法进行检测,判断检测的通信系统信号的周期是否正确,从该M 2个比例因子候选值中筛选周期检测正确的概率大于预设阈值δ 2的比例因子,获得M 1个信噪比下对应的M 1组比例因子;
根据M 1个信噪比下对应的M 1组比例因子拟合比例因子的函数。
3.如权利要求2所述的变换域通信系统信号的检测方法,其特征在于,所述比例因子的函数为:
Figure 137722DEST_PATH_IMAGE001
其中,μ为比例因子,SNR为信噪比,snr(fs, np)为与信号采样率fs和采样点数np相关的函数,表示求在采样率为fs和采样点数为np的条件下,使得上述M 1×M 2×M 3次信号周期检测准确的概率大于预设阈值δ 2的最小信噪比,μ snr 为最小信噪比下的比例因子。
4.如权利要求3所述的变换域通信系统信号的检测方法,其特征在于,当fs=1000,np≥256000时,所述比例因子的函数为:
Figure 789283DEST_PATH_IMAGE002
5.如权利要求1所述的变换域通信系统信号的检测方法,其特征在于,所述计算每对相邻的两个谱峰对应的自相关值的差值前还包括:
执行虚假谱峰抑制步骤,所述虚假谱峰抑制步骤包括:取所有谱峰对应的时间间隔,若相邻的两个谱峰对应的时间间隔小于预设阈值δ 3,则舍弃该相邻的两个谱峰中自相关值较小的峰值,只保留该相邻的两个谱峰中自相关值较大的峰值;
重复所述虚假谱峰抑制步骤,直到所有谱峰间隔均大于预设阈值δ 3
6.如权利要求5所述的变换域通信系统信号的检测方法,其特征在于,将重复所述虚假谱峰抑制步骤后谱峰的总个数记为L,所述只保留该相邻的两个谱峰中自相关值较大的峰值后,且所述计算每对相邻的两个谱峰对应的自相关值的差值前还包括步骤:
只保留自相关值最大的谱峰、其左边相邻的[L/6]个谱峰及其右边相邻的[L/6]个谱峰,[]表示取整。
7.如权利要求1所述的变换域通信系统信号的检测方法,其特征在于,所述计算每对相邻的两个谱峰对应的自相关值的差值时,每个谱峰只参与一次差值计算。
8.一种变换域通信系统信号的检测系统,其特征在于,包括:
自相关函数求解模块,用于按设定的采样率对接收的变换域通信系统信号进行采样,获取信号采样序列,求解所述信号采样序列的自相关函数,所述自相关函数的自变量为延迟时间,因变量为自相关值;
谱峰判断模块,用于将所述自相关函数的最大自相关值与比例因子相乘的结果作为谱峰判断阈值δ 1,将大于所述谱峰判断阈值δ 1的自相关值作为谱峰;
周期计算模块,用于计算每对相邻的两个谱峰对应的延迟时间的差值,将所述差值作为周期;
输出模块,用于统计每个周期出现的频次比,若频次比最大的周期占比超过预设阈值δ 2,则判断所述接收信号包括变换域通信系统信号,并将频次比最大的周期作为所述变换域通信系统信号的周期。
9.如权利要求8所述的变换域通信系统信号的检测系统,其特征在于,所述比例因子是与信噪比相关的函数,所述比例因子的函数采用仿真方法获得,具体包括:
设定信噪比和比例因子的变化范围,在信噪比和比例因子的变化范围内取M 1个信噪比和M 2个比例因子候选值,每个信噪比下仿真产生M 3个变换域通信系统信号;
对每个信噪比下的M 3个变换域通信系统信号,分别采用该M 2个比例因子候选值和所述的变换域通信系统信号的检测方法进行检测,判断检测的通信系统信号的周期是否正确,从该M 2个比例因子候选值中筛选周期检测正确的概率大于预设阈值δ 2的比例因子,获得M 1个信噪比下对应的M 1组比例因子;
根据M 1个信噪比下对应的M 1组比例因子拟合比例因子的函数。
10.如权利要求8所述的变换域通信系统信号的检测系统,其特征在于,还包括谱峰筛选模块,用于在所述计算每对相邻的两个谱峰对应的自相关值的差值前,执行步骤:
执行虚假谱峰抑制步骤,所述虚假谱峰抑制步骤包括:取所有谱峰对应的时间间隔,若相邻的两个谱峰对应的时间间隔小于预设阈值δ 3,则舍弃该相邻的两个谱峰中自相关值较小的峰值,只保留该相邻的两个谱峰中自相关值较大的峰值;
重复所述虚假谱峰抑制步骤,直到所有谱峰间隔均大于预设阈值δ 3
将重复所述虚假谱峰抑制步骤后谱峰的总个数记为L,所述只保留该相邻的两个谱峰中自相关值较大的峰值后,且所述计算每对相邻的两个谱峰对应的自相关值的差值前还包括步骤:
只保留自相关值最大的谱峰、其左边相邻的[L/6]个谱峰及其右边相邻的[L/6]个谱峰,[]表示取整。
CN202210517701.7A 2022-05-13 2022-05-13 一种变换域通信系统信号的检测方法及系统 Active CN114650108B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210517701.7A CN114650108B (zh) 2022-05-13 2022-05-13 一种变换域通信系统信号的检测方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210517701.7A CN114650108B (zh) 2022-05-13 2022-05-13 一种变换域通信系统信号的检测方法及系统

Publications (2)

Publication Number Publication Date
CN114650108A CN114650108A (zh) 2022-06-21
CN114650108B true CN114650108B (zh) 2022-08-05

Family

ID=81997173

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210517701.7A Active CN114650108B (zh) 2022-05-13 2022-05-13 一种变换域通信系统信号的检测方法及系统

Country Status (1)

Country Link
CN (1) CN114650108B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116032703B (zh) * 2023-03-29 2023-06-27 中国人民解放军海军工程大学 一种变换域通信系统信号码元数量的估计方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002221546A (ja) * 2001-01-26 2002-08-09 Nippon Telegr & Teleph Corp <Ntt> ノイズ中の時間周期性微弱信号検出方法、その装置、そのプログラム及びその記録媒体
CN102841248A (zh) * 2012-09-04 2012-12-26 西安石油大学 一种任意频率任意幅度微弱信号检测方法
CN104301272A (zh) * 2013-07-17 2015-01-21 上海无线通信研究中心 基于循环自相关函数的统计谱域传输信号的检测方法
CN106100762A (zh) * 2016-08-23 2016-11-09 桂林电子科技大学 一种循环平稳谱分析的弱通信信号检测方法
CN114039679A (zh) * 2022-01-10 2022-02-11 中国人民解放军海军工程大学 一种低频正交天线信号检测方法及系统

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002221546A (ja) * 2001-01-26 2002-08-09 Nippon Telegr & Teleph Corp <Ntt> ノイズ中の時間周期性微弱信号検出方法、その装置、そのプログラム及びその記録媒体
CN102841248A (zh) * 2012-09-04 2012-12-26 西安石油大学 一种任意频率任意幅度微弱信号检测方法
CN104301272A (zh) * 2013-07-17 2015-01-21 上海无线通信研究中心 基于循环自相关函数的统计谱域传输信号的检测方法
CN106100762A (zh) * 2016-08-23 2016-11-09 桂林电子科技大学 一种循环平稳谱分析的弱通信信号检测方法
CN114039679A (zh) * 2022-01-10 2022-02-11 中国人民解放军海军工程大学 一种低频正交天线信号检测方法及系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Autocorrelation Based Detection of DSSS Signal for Cognitive Radio System;Zhipeng Deng,etc.;《2011 International Conference on Wireless Communications and Signal Processing (WCSP)》;20111208;全篇 *
基于循环谱的DSSS信号检测方法研究;周钱,等;《通信技术》;20171110;全篇 *

Also Published As

Publication number Publication date
CN114650108A (zh) 2022-06-21

Similar Documents

Publication Publication Date Title
Yang et al. Cyclostationary feature detection based spectrum sensing algorithm under complicated electromagnetic environment in cognitive radio networks
CN108512791B (zh) 基于定时频偏补偿的星载ais解调方法
CN101662305B (zh) 直接序列扩频系统伪随机码估计方法
CN109802912A (zh) 宽带无线通信系统的同步方法、装置、设备及存储介质
CN109412644B (zh) 一种直扩msk信号的多普勒频率估计方法
CN112187316A (zh) 信号处理方法、装置、接收机及存储介质
CN114650108B (zh) 一种变换域通信系统信号的检测方法及系统
CN101588191B (zh) 无线电信号认知方法及设备
CN110289926B (zh) 基于调制信号循环自相关函数对称峰值的频谱感知方法
CN104253659B (zh) 一种频谱检测方法及其装置
Baradkar et al. Implementation of energy detection method for spectrum sensing in cognitive radio based embedded wireless sensor network node
CN108900445B (zh) 一种信号符号率估计的方法及装置
CN105024966B (zh) Ofdm系统的符号定时同步方法
CN111539323B (zh) 一种循环前缀线性调频信号的频率估计方法与装置
CN109104215A (zh) 一种基于小波变换的跳频信号码元速率盲估计方法
CN105656511B (zh) 一种适应于有频偏和低信噪比环境下的差分相关捕获方法
CN115166650B (zh) 一种雷达信号识别与参数估计方法和系统
JP2022529799A (ja) 変調信号の周期的自己相関関数の対称なピーク値に基づくスペクトル検知方法
CN102111228A (zh) 一种基于循环对称性的认知无线电频谱感知方法
CN115378776A (zh) 一种基于循环谱参数的mfsk调制识别方法
CN109361632B (zh) 基于谱相关的多径boc信号的参数盲估计方法
CN102130862B (zh) 一种降低通信系统信道估计开销的方法
CN105959035B (zh) 一种直扩信号截获检测方法
Varma et al. Cognitive wideband sensing using correlation of inverted spectrum segments
CN109004998B (zh) 一种基于Gram-Charlier级数的OFDM信号频谱感知方法

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