一种基于高频振荡共振原理的信号检测方法与流程

本发明涉及通讯信号检测
技术领域：
，具体涉及其中的高频振荡共振调谐方法。
背景技术：
：在通讯信号传输、雷达信号检测等信息处理过程中，传输信号不仅受到电路和设备的固有噪声的污染，而且还受到不同调制方式、非线性效应、码间串扰等引起的周期性高频振荡干扰。为提高在噪声和高频振荡组合形成的混合干扰中弱信号检测概率，对于高频振荡、噪声以及抖动的建模和分析方面取得了一些理论进展，也出现了用于减缓高频振荡、噪声以及抖动的链路结构和电路。在这些传统的检测器设计方法中，检测器结构都是依据输入噪声和干扰的概率分布进行设计的。但是，这种设计方法在噪声和干扰的统计分布特性未知时将无法进行，同时，对于某些复杂的非高斯噪声概率分布，检测器结构设计异常复杂，难以实现。为了应对上述检测器设计缺陷，次优性能的广义相关检测器和非参数检测器由于其易于实现和鲁棒性强的特点被广泛采用。为提高广义相关检测器和非参数检测器的检测性能，多通道分集技术和自适应参数调节技术都得到广泛关注和发展，但是这些方法面临着无法逼近检测性能极限的局限性，特别是在非高斯噪声背景下的弱信号检测性能还不明确，没有一种能够简易实现的检测方法。同时，传统的检测方法都是将高频振荡、噪声以及抖动视为有害因素，认为其是信号检测的主要障碍，通过系统结构、算法和电路设计等对其进行消除。然而，Landa与McClintock[P.S.Landa,P.V.E.McClintock,Vibrationalresonance,JournalofPhysicsA：MathematicalandGeneral，33(45),L433–L438,2000]证实了高频振荡对系统输出信噪比、神经元的信息传输等具有有益的作用，即振荡共振理论。实际操作中高频振荡非常易于实现，利用高频振荡提高检测器的性能是本发明的主要技术思路。弱信号检测的基础模型是二元假设检验观测模型，基于奈曼-皮尔逊原理的恒虚警检测基础技术背景如下。二元假设检验观测模型为H1:xn＝θsn+wn,H0:xn＝wn这里H1表示信号存在的假设，弱信号强度θ＞0且θ→0，离散信号向量S＝[s1,s2,…,sN]T，观测向量为X＝[x1,x2,…,xN]T，数据长度为N，H0表示只有噪声的假设，白噪声向量W＝[w1,w2,…,wN]T，这里噪声概率密度模型不局限为高斯分布模型。利用广义相关检测器T(X)=Σn=1Ng(xn)sn>γ]]>判断为假设H1，否则判断为H0。我们假设转换函数g具有零均值E[g(x)]＝0和有限二阶矩E[g2(x)]＜∞。由于统计量T(X)是由大量同分布随机变量的和，因此T(X)的分布可以看成高斯分布。那么在假设H0是真时，统计量T(X)的均值为E[T|H0]＝0，方差为var(T|H0)≈PsE[g2(x)]，这里信号能量而在假设H1是真时，统计量T(X)的均值为E[T|H0]≈θPsE[g'(x)]，方差为var(T|H0)＝var(T|H1)。那么相关检测器的检测门限为γ=var[T|H0]Q-1(PF)]]>检测概率可以计算为PD=Q{Q-1(PF)-θPsE2[g′(x)]/E[g2(x)]}]]>这里右尾概率函数恒虚警概率Pf下相关检测器的检测性能与检测效率ξ＝E2[g'(x)]/E[g2(x)]成正比，或者与检测器输出信噪比SNR＝θ2Psξ＝Esξ成正比。高频振荡共振是类似于随机共振的一种非线性信息处理技术，通过调谐高频振荡的幅值和频率达到优化检测器性能的目的。高频振荡干扰不仅易于实现，而且还对于外部噪声的强度变化和类型改变具有较好的鲁棒性，不仅适用于确定性系统，也适用于随机系统，基于高频振荡共振原理的信号检测方法是一种非常适合实际应用的新检测技术。技术实现要素：本发明要解决的技术问题是提供一种基于高频振荡共振原理的信号检测方法。为了解决上述技术问题，本发明采用的技术方案是，一种基于高频振荡共振原理的信号检测方法，包括以下步骤：步骤1：初始化参数，包括接收数据X并记录其长度N，设定标准信号波形S＝{sn}，选择转移函数g，高频振荡频率范围[fL,fU]且fL＞＞fs，fs是信号频率上限，给定虚警概率Pf，假设H1表示信号存在，而假设H0表示仅有噪声，Es为实际信号能量；步骤2：将接收数据X复制M份，每份加入不同频率的高频振荡之后馈入转移函数g，全部函数的输出组合之后计算统计量均值以及两两不同函数输出组合之后求M(M-1)组点积互相关数据gmngpn，m≠p，n＝1,2,…,N，从而建立检测统计量步骤3：给定条件下，检测性能与通道M成正比，对于M＝∞的极限情况，可以由任意一个通道的输出均值求出统计量T(X)在假设H1下的均值E[T|H1]，依据均方极限定理，在假设H0下方差var[T|H0]可以转化为任意两个不同通道的点积期望值E[gmngpn]；步骤4：利用步骤3的结果可以计算结果计算检测效率ξ=E2[T|H1]var[T|H0]]]>则通过∂ξ∂Aη=0]]>求出最优振荡幅值调谐值从而完成自适应高频振荡共振调谐规律；步骤5：根据此最优值产生的高频振荡干扰加入接收数据X，更新数据为计算假设H1下的均值和方差计算H0假设下的均值(对称性)和方差得到最大检测效率依据给定的虚警概率Pf，可以计算出门限值γopt=var[T(X~)|H0]Q-1(Pf)]]>步骤6：将统计量和门限值进行比较，若T(X)=Σn=1Ng‾nsn>γopt]]>判断为假设H1，否则判断为H0；最后按上式获得检测概率。本发明的有益效果是：利用了高频振荡干扰简单易于实现的特点，将其馈入不同分集通道，利用不同通道的输出进行互相关计算后再分集的技术，形成自适应高频振荡共振调谐方法，依据能量累积和中心极限定理，该方法能够达到了检测器的Chernoff上界，通过高频振荡共振原理，本发明极大地提高了弱信号情况下检测器的检测概率。附图说明下面结合附图和具体实施方式对本发明作进一步详细的说明。图1是本发明基于高频振荡共振原理的信号检测方法实施例的高频振荡共振检测器示意图。图2是本发明基于高频振荡共振原理的信号检测方法实施例的检测概率的振荡共振图。图3是本发明基于高频振荡共振原理的信号检测方法实施例的不同输入信噪比下振荡共振曲线。图4是本发明基于高频振荡共振原理的信号检测方法实施例的不同输入信噪比下检测概率。图5是本发明基于高频振荡共振原理的信号检测方法实施例的不同噪声类型下最大检测概率maxPd(A)和最优振荡幅值(B)。具体实施方式图1是一种高频振荡共振检测器，以下是基于该高频振荡共振检测器的信号检测方法的具体步骤：(1)初始化参数，给定虚警概率转移函数g，信号强度θ，参考信号向量S，设信号最大频率fs，白噪声向量W以及其概率密度函数fw，观测数据向量X,数据长度为N，高频振荡频率范围[fL,fU]且fL＞＞fs，信号能量(2)取M组高频振荡ηmn＝Aηsin(2πfmn/N)加入到观测数据向量X中，m＝1,2,…,M，n＝1,2,…,N，将混合数据通过转移函数g得到各通道的输出ymn＝g(xn+ηmn)(1)将两两通道输出相互点积计算得到相关输出zn＝g(xn+ηin)g(xn+ηjn)，i,j＝1,3,5…M，i≠j(2)(3)各通道的分集和为g‾n=1MΣm=1Mymn---(3)]]>建立统计量T(xn)=Σn=1Ng‾nsn---(4)]]>(4)计算H0假设下统计量的输出均值E[T|H0]=Σn=1N1MΣm=1M∫g(x+ηmn)fw(x)dxsn=0---(5)]]>统计量的方差var[T|H0]=EsM2Σn=1N{Σm=1M∫g2(x+ηmn)fw(x)dx+Σi=1MΣj=1M∫g(x+ηin)g(x+ηjn)fw(x)dx},i≠j---(6)]]>(5)计算H1假设下统计量的输出均值E[T|H1]=Σn=1N1MΣm=1M∫g(x+sn+ηmn)fw(x)dxsn---(7)]]>统计量的方差var[T|H0]≈var[T|H1](8)(6)计算检测器检测功效ξ=E2[T|H1]Esvar[T|H0]---(9)]]>(7)给定系统、噪声、振荡幅值，检测功效是检测器数目M的单调递增函数，为得到M＝∞的检测极限值，我们将检测统计量的特征进行近似E[T|H1]=limM→∞Σn=1N1MΣm=1M∫g(x+sn+ηmn)fw(x)dxsn=Σn=1N∫g(x+sn+ηmn)fw(x)dxsn,∀m∈M---(10)]]>var[T|H0]=limM→∞EsM2Σn=1N{Σm=1M∫g2(x+ηmn)fw(x)dx+Σi=1MΣj=1M∫g(x+ηin)g(x+ηjn)fw(x)dx}=EsΣn=1N∫g(x+ηin)g(x+ηjn)fw(x)dx,∀i,j∈M,i≠j---(11)]]>(8)那么可以利用输出信噪比的表达式，得到最优振荡幅值的解满足的超越方程Σn=1N∫∂g(x+sn+ηmn)∂Aηfw(x)dxΣn=1N∫g(x+ηin)g(x+ηjn)fw(x)dx=Σn=1N∫g(x+sn+ηmn)fw(x)dxΣn=1N∫∂g(x+ηin)∂Aηg(x+ηjn)fw(x)dx---(12)]]>获得最优振荡幅值调谐值(9)将具有最优振荡幅值调谐值的高频振荡重新馈入传递函数，得到更新的均值和方差E[T|H1]=Σn=1N∫g(x+sn+ηmnopt)fw(x)dxsn,∀m∈M---(13)]]>var[T|H0]=EsΣn=1N∫g(x+ηinopt)g(x+ηjnopt)fw(x)dx,∀i,j∈M,i≠j---(14)]]>此时的检测器功效计算为ξopt=(Σn=1N∫g(x+sn+ηmnopt)fw(x)dx)2Σn=1N∫g(x+ηinopt)g(x+ηjnopt)fw(x)dx---(15)]]>(10)最终检测参数为：检测门限γopt=EsΣn=1N∫g(x+ηinopt)g(x+ηjnopt)fw(x)dxQ-1(PF)---(16)]]>检测器判定方法T(X)=Σn=1Ng‾nsn>γopt---(17)]]>检测概率可以计算为PD=Q[Q-1(PF)-Esξ]---(18)]]>由于式(13)与(14)是M＝∞情况下检测统计量T(X)，其检测概率式(18)就是给定检测器的Chernoff界。实验结果(1)在给定的噪声强度离散输入信号sn＝0.1sin(2πfsn/N)，fs为任意的，输入信噪比为-22dB，虚警概率Pf＝0.01，外部噪声的概率密度函数为fw(x)=α2βΓ(α-1)exp(-|xβ|α)]]>这里α是形状参数，β是尺度参数，此例中α＝8。检测器的转换函数为硬限幅器g(x)=1,x>θ0,-θ≤x≤θ-1,x<-θ]]>这里θ为阈值参数。高频振荡ηmn＝Aηsin(2πfmn/N)，fm＝(20+m)fs，我们模拟无穷多并联检测器，取fi＝21fs和fj＝22fs(m可以是任意的整数，i≠j)。图2为检测概率式(25)中PD与阈值θ和高频振荡幅值Aη三维图。我们可以清晰地看出，存在一个最优高频振荡幅值使得PD能够取得最大值。(2)给定阈值θ＝2，图3给出了三种输入信噪比下振荡共振曲线，其他参数同图3。我们取其最大检测概率，对比不同输入信噪比的性能曲线在图4中给出，可以看出即使输入信噪比的恶化为-26dB，检测概率依然能够高于80％。(3)给定输入信噪比-22dB，对于不同类型的噪声，由α参数控制，我们在图5中给出最大检测概率maxPd(图5A)和相应的最优振荡幅值(图5B)。由图5(A)可以看出，同样输入信噪比和检测器情况下，高斯噪声是最差的噪声监测环境，非高斯噪声环境下都能够得到更高的检测概率。以上所述的本发明实施方式，并不构成对本发明保护范围的限定。任何在本发明的精神和原则之内所作的修改、等同替换和改进等，均应包含在本发明的权利要求保护范围之内。当前第1页1 2 3