CN103957066B - 一种基于非参数核函数的短时突发通信信号检测方法 - Google Patents
一种基于非参数核函数的短时突发通信信号检测方法 Download PDFInfo
- Publication number
- CN103957066B CN103957066B CN201410215812.8A CN201410215812A CN103957066B CN 103957066 B CN103957066 B CN 103957066B CN 201410215812 A CN201410215812 A CN 201410215812A CN 103957066 B CN103957066 B CN 103957066B
- Authority
- CN
- China
- Prior art keywords
- function
- probability density
- prime
- matrix
- signal
- 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.)
- Expired - Fee Related
Links
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明属于通信信号盲检测技术领域,尤其涉及一种基于非参数核函数的短时突发通信信号检测。本发明通过高斯核函数来拟合概率密度函数比,采取最小化估计概率密度比与真实概率密度比的误差函数方法,将概率密度比的估计问题转化为一个凸优化问题,从而实现概率密度函数比的精确估计。利用得到的估计概率密度比函数构造Pearson距离得到信号的检测函数曲线,对检测函数曲线采用峰值搜索方法得到突发信号的起始点和终止点。本发明针对突发模式下的通信信号进行检测,能够精确地判断信号的起始点和终止点,为后续的信号处理提供了更多的先验信息,在非协作模式下对于信号的检测提供了更有效的处理手段,特别是突发模式通信通常处于非平稳信道。
Description
技术领域
本发明属于通信信号盲检测技术领域,尤其涉及一种基于非参数核函数的短时突发通信信号检测方法。
背景技术
传统的信号检测针对的主要是连续平稳信号,在协作模式下通常我们能够获取一些信号的先验信息,比如信号的调制方式、载频等,甚至信号可能带有前导头来作为信号检测的先验知识。但是随着社会的发展现有的通信方式已经不能满足时代的需求,我们对于非协作模式下的通信方式的需求越来越大,通常非协作模式能够提高信道的利用率。比如,认知无线电中非注册用户检测到注册用户已经下线,那么就可以使用信道来传输消息,这种智能化的传输方式由于能提高信道利用率得到越来越多的关注。同时,由于是非协作模式,信号的保密性得到保障,在军事领域中突发的通信信号因为其强抗干扰的能力在现代军事通信中也使用的越来越多。当然在民用领域中,突发模式通信也得到广泛的应用,比如卫星TDMA信号就可以看成是一种突发信号。综上所述,突发通信信号的应用领域越来越多,但是由于突发通信信号通常具有短时以及短间隔等特点,对于信号的起始点以及终止点的检测变得愈加困难,所以发展这种短时突发信号检测方法也显得理所当然。
现有的检测短时突发通信信号方法主要有基于时域和频域的以及时频结合的方法。基于时域的方法主要有能量检测方法,但是这种方法对噪声太敏感,当信噪比太低时的检测性能不尽令人满意。基于频域的方法相对于时域来看能够降低噪声的影响,但是对于起始点以及终止点的检测的准确性依然有待提高。另外,相对于协作模式下的信号传输,非协作模式下的突发信号传输信道通常是时变的,也就是说是非平稳的,这样对于传统的检测方法提出了更高的要求。这里我们提出的非参数核方法对于信道的平稳性没有要求,在非平稳信道中依然可以实现对信号起始点和终止点的有效检测。
发明内容
本发明的目的在于提供一种基于非参数核函数的短时突发通信信号盲检测方法,用以克服信道非平稳带来的性能缺陷,实现非协作模式下的突发通信信号起始点终止点的准确检测。
本发明的目的通过如下步骤实现:
S1、对接收机接收到的模拟信号进行采样,得到信号时间序列y[n],其中,n是不为零的自然数;
S2、用滑窗在S1所得y[n]上截取样本,构造样本向量Y(n),所述滑窗长度为k,其中,k≥1,k是自然数;
S3、根据样本向量Y(n)构造两个Hankel矩阵,具体如下:
S4、用高斯模型对S3所述两个Hankel矩阵进行估计,得到两个Hankel矩阵的样本的概率密度比其中,l在[1,2,3...,N]中遍历取值,为高斯核函数,σ为高斯宽度,p(Y)为样本矩阵Ψ(n)的概率密度函数,p'(Y)为样本矩阵Ψ(n+N)的概率密度函数,θ=(θ1,θ2,...,θN)T为模型参数,上标T表示矩阵转置;
S5、最小化真实密度比函数与估计的密度比函数的二阶误差 即得到凸优化问题的解同时得到密度比函数的估计其中,是为了规范化的目的,参数λ≥0,是一个N维向量,的第l个元素是是一个N×N矩阵,第(l,l')个元素是 Yi,Yj表示Hankel矩阵的不同列向量;
S6、根据交叉验证算法得到σ,λ,的最优值,以此求得最优概率密度比函数估计
S7、计算样本矩阵Ψ(n)和样本矩阵Ψ(n+N)的概率密度比函数的Pearson距离 其中, 是Ψ(n)作为概率密度比的分子,Ψ(n+N)作为概率密度比的分母时,是Ψ(n+N)作为概率密度比的分子,Ψ(n)作为概率密度比的分母时;
S8、通过对的峰值搜索得到突发信号的起始点和终止点,其中,是作为短时突发信号的检测函数,具体如下:
S81、对标准化后进行峰值搜索,滤除中的毛刺,得到L个最大峰值点,设置峰值滤波器门限为
S82、设置峰值间隔门限τ,当S81所述最大峰值点中相邻的点间隔小于τ,即|nm-nm-1|<τ,滤除第m个峰值点nm,滤除峰值处理后剩下的峰值点即为突发信号的起始点和终止点,其中,m∈[1,L],30≤τ≤50。
进一步地,S2所述滑窗长度满足:10≤k≤20。
进一步地,所述交叉验证算法得到σ,λ,的最优值,具体如下:
步骤1、设置候选参数高斯宽度σ=0.6dmed,0.8dmed,dmed,1.2dmed,1.4dmed,其中,dmed为Hankel矩阵的样本的中间距离,即矩阵[Ψ(n) Ψ(n+N)]各列向量之间的距离的中间值
步骤2、设置候选参数λ=10-3,10-2,10-1,100,10,设置交叉次数fold;
步骤3、将矩阵Ψ(n)的列平均分成不相交的fold段,每段n'个列向量,将Ψ(n+N)的列平均分成不相交的fold段,每段n'个列向量;
步骤4、将步骤3所述Ψ(n)的数据段中选取第r段数据,将步骤3所述Ψ(n+N)的数据段中选取第r段数据,记作r,剩余数据段记作R,其中,r∈{1,2,...,fold};
步骤5、根据步骤4所述R计算出并得到概率密度比函数估计
步骤6、根据步骤5所述和步骤4所述r计算概率密度函数估计误差函数 其中,表示矩阵Ψ(n)的第r段第i个列向量;
步骤7、将步骤6重复进行fold次,每次得到一个估计误差函数,求出总的平均误差函数
步骤8、选出最小做对应的σ,λ为最优参数。
进一步地,所述交叉次数fold=5。
本发明的有益效果是:
本发明通过高斯核函数来拟合概率密度函数比,采取最小化估计概率密度比与真实概率密度比的误差函数方法,将概率密度比的估计问题转化为一个凸优化问题,从而实现概率密度函数比的精确估计。利用得到的估计概率密度比函数构造Pearson距离得到信号的检测函数曲线,对检测函数曲线采用峰值搜索方法得到突发信号的起始点和终止点。本发明针对突发模式下的通信信号进行检测,能够精确地判断信号的起始点和终止点,为后续的信号处理提供了更多的先验信息,在非协作模式下对于信号的检测提供了更有效的处理手段,特别是突发模式通信通常处于非平稳信道。本发明在非平稳的条件下仍然能够提供很好的检测性能,因此尤其适合于非平稳信道下的突发信号检测。
附图说明
图1是本发明流程图。
图2是QPSK信号检测概率随信噪比变化曲线。
图3是QPSK信号虚警概率随信噪比变化曲线。
具体实施方式
下面结合附图对本发明的具体实施方式进行描述。
图1是本发明基于非参数核函数的短时突发通信信号检测方法的一种具体实施方式流程图。如图1所示,以QPSK信号为例,码速率2500波特率,每段突发信号段长度20ms,噪声段与信号段长度相同,采样率50kHz,载频10kHz,总共5段突发信号,突发信号与噪声信号相互间隔,总共进行200次蒙特卡洛实验。本发明基于非参数核函数的短时突发信号检测方法包括以下步骤:
S1、对接收机接收到的模拟信号进行采样,得到信号时间序列y[n],其中,n是不为零的自然数;
S2、用滑窗在S1所得y[n]上截取样本,构造样本向量Y(n),所述滑窗长度为k,其中,k=10;
S3、根据样本向量Y(n)构造两个Hankel矩阵,具体如下:
S4、用高斯模型对S3所述两个Hankel矩阵进行估计,得到两个Hankel矩阵的样本的概率密度比其中,l在[1,2,3...,N]中遍历取值,为高斯核函数,σ为高斯宽度,p(Y)为样本矩阵Ψ(n)的概率密度函数,p'(Y)为样本矩阵Ψ(n+N)的概率密度函数,θ=(θ1,θ2,...,θN)T为模型参数,上标T表示矩阵转置;
S5、最小化真实密度比函数与估计的密度比函数的二阶误差 即得到凸优化问题的解同时得到密度比函数的估计其中,是为了规范化的目的,参数λ≥0,是一个N维向量,的第l个元素是是一个N×N矩阵,第(l,l')个元素是 Yi,Yj表示Hankel矩阵的不同列向量;
S6、根据交叉验证算法得到σ,λ,的最优值,以此求得最优概率密度比函数估计
S7、计算样本矩阵Ψ(n)和样本矩阵Ψ(n+N)的概率密度比函数的Pearson距离 其中, 是Ψ(n)作为概率密度比的分子,Ψ(n+N)作为概率密度比的分母时,是Ψ(n+N)作为概率密度比的分子,Ψ(n)作为概率密度比的分母时;
S8、通过对的峰值搜索得到突发信号的起始点和终止点,其中,是作为短时突发信号的检测函数,具体如下:
S81、对标准化后进行峰值搜索,滤除中的毛刺,得到L个最大峰值点,设置峰值滤波器门限为
S82、设置峰值间隔门限τ,当S81所述最大峰值点中相邻的点间隔小于τ,即|nm-nm-1|<τ,滤除第m个峰值点nm,其中,m∈[1,L],τ=40;
S83、当|nm-n0|<20,则判定nm为检测到一个目标点,虚警概率定义为Pf=1-Pd,其中,n0为目标信号点,检测概率Pd为搜索到的信号点个数与L的比值,滤除峰值处理后剩下的峰值点即为突发信号的起始点和终止点。
其中,所述交叉验证算法得到σ,λ,的最优值,具体如下:
步骤1、设置候选参数高斯宽度σ=0.6dmed,0.8dmed,dmed,1.2dmed,1.4dmed,其中,dmed为Hankel矩阵的样本的中间距离,即矩阵[Ψ(n) Ψ(n+N)]各列向量之间的距离的中间值
步骤2、设置候选参数λ=10-3,10-2,10-1,100,10,设置交叉次数fold=5;
步骤3、将矩阵Ψ(n)的列平均分成不相交的fold段,每段n'个列向量,将Ψ(n+N)的列平均分成不相交的fold段,每段n'个列向量;
步骤4、将步骤3所述Ψ(n)的数据段中选取第r段数据,将步骤3所述Ψ(n+N)的数据段中选取第r段数据,记作r,剩余数据段记作R,其中,r∈{1,2,...,fold};
步骤5、根据步骤4所述R计算出并得到概率密度比函数估计
步骤6、根据步骤5所述和步骤4所述r计算概率密度函数估计误差函数 其中,表示矩阵Ψ(n)的第r段第i个列向量;
步骤7、将步骤6重复进行fold次,每次得到一个估计误差函数,求出总的平均误差函数
步骤8、选出最小做对应的σ,λ为最优参数。
图2为QPSK信号的检测概率随信噪比变化曲线,从图中可以看出,检测概率随信噪比增加而增加,到15dB检测概率几乎可以达到100%;图3为QPSK信号的虚警概率随信噪比变化曲线,从中可以看出虚检概率随信噪比增加而减小,大约在7dB左右虚警概率就可以达到0.1。
Claims (4)
1.一种基于非参数核函数的短时突发通信信号检测方法,其特征在于,包括如下步骤:
S1、对接收机接收到的模拟信号进行采样,得到信号时间序列y[n],其中,n是不为零的自然数;
S2、用滑窗在S1所得y[n]上截取样本,构造样本向量Y(n),所述滑窗长度为k,其中,k≥1,k是自然数;
S3、根据样本向量Y(n)构造两个Hankel矩阵,具体如下:
S4、用高斯模型对S3所述两个Hankel矩阵进行估计,得到两个Hankel矩阵的样本的概率密度比其中,l在[1,2,3...,N]中遍历取值,为高斯核函数,σ为高斯宽度,p(Y)为样本矩阵Ψ(n)的概率密度函数,p'(Y)为样本矩阵Ψ(n+N)的概率密度函数,θ=(θ1,θ2,...,θN)T为模型参数,上标T表示矩阵转置;
S5、最小化真实密度比函数与估计的密度比函数的二阶误差 即 得到凸优化问题的解同时得到密度比函数的估计其中,是为了规范化的目的,参数λ≥0,是一个N维向量,的第l个元素是 是一个N×N矩阵,第(l,l')个元素是 Yi,Yj表示Hankel矩阵的不同列向量;
S6、根据交叉验证算法得到的最优值,以此求得最优概率密度比函数估计
S7、计算样本矩阵Ψ(n)和样本矩阵Ψ(n+N)的概率密度比函数的Pearson距离 其中, 是Ψ(n)作为概率密度比的分子,Ψ(n+N)作为概率密度比的分母时,是Ψ(n+N)作为概率密度比的分子,Ψ(n)作为概率密度比的分母时;
S8、通过对的峰值搜索得到突发信号的起始点和终止点,其中,是作为短时突发信号的检测函数,具体如下:
S81、对标准化后进行峰值搜索,滤除中的毛刺,得到L个最大峰值点,设置峰值滤波器门限为
S82、设置峰值间隔门限τ,当S81所述最大峰值点中相邻的点间隔小于τ,即|nm-nm-1|<τ,滤除第m个峰值点nm,滤除峰值处理后剩下的峰值点即为突发信号的起始点和终止点,其中,m∈[1,L],30≤τ≤50。
2.根据权利要求1所述的一种基于非参数核函数的短时突发通信信号检测方法,其特征在于:S2所述滑窗长度满足:10≤k≤20。
3.根据权利要求1所述的一种基于非参数核函数的短时突发通信信号检测方法,其特征在于:所述交叉验证算法得到的最优值,具体如下:
步骤1、设置候选参数高斯宽度σ=0.6dmed,0.8dmed,dmed,1.2dmed,1.4dmed,其中,dmed为Hankel矩阵的样本的中间距离,即矩阵[Ψ(n)Ψ(n+N)]各列向量之间的距离的中间值
步骤2、设置候选参数λ=10-3,10-2,10-1,100,10,设置交叉次数fold;
步骤3、将矩阵Ψ(n)的列平均分成不相交的fold段,每段n'个列向量,将Ψ(n+N)的列平均分成不相交的fold段,每段n'个列向量;
步骤4、将步骤3所述Ψ(n)的数据段中选取第r段数据,将步骤3所述Ψ(n+N)的数据段中选取第r段数据,记作r,剩余数据段记作R,其中,r∈{1,2,...,fold};
步骤5、根据步骤4所述R计算出并得到概率密度比函数估计
步骤6、根据步骤5所述和步骤4所述r计算概率密度函数估计误差函数 其中,表示矩阵Ψ(n)的第r段第i个列向量;
步骤7、将步骤6重复进行fold次,每次得到一个估计误差函数,求出总的平均误差函数
步骤8、选出最小做对应的σ,λ为最优参数。
4.根据权利要求3所述的一种基于非参数核函数的短时突发通信信号检测方法,其特征在于:步骤2所述交叉次数fold=5。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410215812.8A CN103957066B (zh) | 2014-05-21 | 2014-05-21 | 一种基于非参数核函数的短时突发通信信号检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410215812.8A CN103957066B (zh) | 2014-05-21 | 2014-05-21 | 一种基于非参数核函数的短时突发通信信号检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103957066A CN103957066A (zh) | 2014-07-30 |
CN103957066B true CN103957066B (zh) | 2015-09-09 |
Family
ID=51334304
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410215812.8A Expired - Fee Related CN103957066B (zh) | 2014-05-21 | 2014-05-21 | 一种基于非参数核函数的短时突发通信信号检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103957066B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105553619B (zh) * | 2015-12-24 | 2018-10-12 | 成都国恒空间技术工程有限公司 | 基于调制识别辅助的突发信号检测方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1404309A (zh) * | 2002-10-24 | 2003-03-19 | 上海交通大学 | 正交频分复用系统的突发信号检测和自动增益控制方法 |
CN101848185A (zh) * | 2010-06-12 | 2010-09-29 | 电子科技大学 | 突发通信的帧同步信号设计技术 |
CN101931506A (zh) * | 2009-06-23 | 2010-12-29 | 大唐移动通信设备有限公司 | 一种基于概率数据辅助检测信号的方法及设备 |
CN102386984A (zh) * | 2011-10-19 | 2012-03-21 | 上海师范大学 | 一种认知无线电宽带频谱检测方法 |
CN103117964A (zh) * | 2013-01-09 | 2013-05-22 | 北京邮电大学 | 一种60GHz毫米波通信系统下的信号检测方法与装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE50111853D1 (de) * | 2001-07-17 | 2007-02-22 | Cit Alcatel | Überwachungseinheit für optische Burst-Signale |
US8705608B2 (en) * | 2012-04-23 | 2014-04-22 | Micrel, Inc. | Noise discriminator for enhanced noise detection in a passive optical network burst mode receiver |
-
2014
- 2014-05-21 CN CN201410215812.8A patent/CN103957066B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1404309A (zh) * | 2002-10-24 | 2003-03-19 | 上海交通大学 | 正交频分复用系统的突发信号检测和自动增益控制方法 |
CN101931506A (zh) * | 2009-06-23 | 2010-12-29 | 大唐移动通信设备有限公司 | 一种基于概率数据辅助检测信号的方法及设备 |
CN101848185A (zh) * | 2010-06-12 | 2010-09-29 | 电子科技大学 | 突发通信的帧同步信号设计技术 |
CN102386984A (zh) * | 2011-10-19 | 2012-03-21 | 上海师范大学 | 一种认知无线电宽带频谱检测方法 |
CN103117964A (zh) * | 2013-01-09 | 2013-05-22 | 北京邮电大学 | 一种60GHz毫米波通信系统下的信号检测方法与装置 |
Non-Patent Citations (4)
Title |
---|
Stochastic Behavior Analysis of the Gaussian Kernel Least-Mean-Square Algorithm;W.D.Parreira;《IEEE Transactions on Signal Processing》;20120531;第60卷(第5期);2208-2222 * |
基于Labview的短时突发信号的捕获与检测系统;钟骥;《通信技术》;20080331;第41卷(第3期);25-27 * |
复杂信号环境下非合作突发通信信号的实时检测;胡亚;《系统工程与电子技术》;20120731;第34卷(第7期);1463-1468 * |
突发信号宽带检测概率分析;罗胜恩;《计算机工程与应用》;20110930;第47卷(第25期);137-145 * |
Also Published As
Publication number | Publication date |
---|---|
CN103957066A (zh) | 2014-07-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103457890A (zh) | 一种有效识别非高斯噪声下数字调制信号的方法 | |
CN103051367A (zh) | 一种基于聚类的同步正交跳频信号盲源分离方法 | |
CN103346845B (zh) | 基于快速傅里叶变换的盲频谱感知方法及装置 | |
CN103607361B (zh) | Alpha稳定分布噪声下时频重叠信号的参数预估方法 | |
CN101577564A (zh) | 基于判决门限自适应的信号频谱感知与检测方法 | |
CN104502894A (zh) | 基于物理层信息的运动物体被动检测方法 | |
CN104618277B (zh) | 一种在窄带干扰下的短波ofdm信道估计方法 | |
CN101494508A (zh) | 基于特征循环频率的频谱检测方法 | |
CN102571230A (zh) | 基于高阶统计量与信噪比盲估计的分布式协同信号识别方法 | |
CN102798870B (zh) | 一种卫星导航系统兼容性分析方法 | |
CN106357574A (zh) | 基于顺序统计量的bpsk/qpsk信号调制盲识别方法 | |
CN101986583A (zh) | 基于协方差匹配的多天线频谱感知方法 | |
CN103616661A (zh) | 一种稳健的远场窄带信号源个数估计方法 | |
CN102404063A (zh) | 一种基于过采样的glrt检测方法 | |
CN103023831A (zh) | 一种适用于突发波形的载波频偏估计方法 | |
CN106899376A (zh) | 物理层网络编码连续相位调制信号的非相干检测方法 | |
CN103957066B (zh) | 一种基于非参数核函数的短时突发通信信号检测方法 | |
CN102025424B (zh) | 一种基于矢量传感器的ofdm水声通信方法 | |
CN103888201B (zh) | 一种利用空间分集的协作频谱感知方法 | |
CN105721083A (zh) | 一种基于自相关能量的频谱检测方法 | |
CN104931925A (zh) | 一种基于参考标签rfid定位的自适应方法 | |
CN109861793A (zh) | 一种基于二阶统计量的空频分组码信号盲识别方法 | |
WO2010075374A3 (en) | Extracting information from positioning pilot channel symbols in forward link only system | |
CN102594747B (zh) | 一种具有信噪比约束的无线传感器网络滚动时域信噪比估计方法 | |
CN102087313A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150909 Termination date: 20160521 |
|
CF01 | Termination of patent right due to non-payment of annual fee |