技术领域本发明涉及一种弱信号目标检测的优化方法，属于计算机算法技术领域。
背景技术：
相关检测技术是根据噪声与噪声、噪声与信号均不相关，而信号与信号则完全相关的特性，通过相关运算达到去除噪声的一种技术。现已普遍证明，它是从噪声中提取有用信号，提高输出信噪比的有效方法。它在自动控制、通信、雷达等领域都获得了广泛的应用。四重相关和三谱(四重相关函数的频谱函数)分析就是高阶矩范围内的信号处理方法，由于零均值的高斯平稳随机过程的四重相关等于零，以及四重相关具有位移和旋转不变性等特点，并且三谱富有丰富的冗余信息。被动式监控探测系统，例如红外监控探测系统较之于主动式探测系统具有难以比拟的优势。然而，被动式目标探测系统所要探测的信号通常迭加有强噪声。因此，从强噪声中将有用的信号提取出来，就成为目标识别的前提。四重相关检测技术可实现信号与噪声的分离，去除噪声。采用匹配滤波器，给出最大的信噪比。通常，对监控探测系统所探测到的目标信号的处理主要是利用二阶统计矩，即功率谱和二重相关函数，这就是二重相关匹配滤波技术。从频域角度来分析，二重相关匹配滤波器可以最大限度地吸收有用信号的能量，而最大限度地抑制信号频带以外的噪声。根据信号与噪声的不同频谱特点，采用线性滤波器消除噪声频谱。当线性滤波器传输函数为输入信号频谱函数的复共轭时，该滤波器能够给出最大的信噪比。然而，许多实验研究表明，无论二重相关滤波器设计得多么精细，信号频带之内的噪声仍然是难以抑制的。为了进一步提高信噪比，应该设法抑制信号频带之内的噪声。一个有效的解决办法就是采用高阶矩范围内的信号处理方法。四重相关和三谱(四重相关函数的频谱函数)分析就是高阶矩范围内的信号处理方法，由于零均值的高斯平稳随机过程的四重相关等于零，以及四重相关具有位移和旋转不变性等特点，并且三谱富有丰富的冗余信息，因此，四重相关匹配滤波技术较之二重相关匹配滤波技术可以进一步抑制噪声从而提高信噪比。随着这一技术研究的深入，基于四重相关的分析方法发展了许多新的信号处理方法，使得这一技术更加完善.并且在红外成像、信号的高阶谱分析等诸多领域获得了广泛的应用。在滤波方面，线性滤波器可以地滤除信号中的噪声，但滤波以后的信噪比不够高，当线性滤波器传输函数为输入信号频谱函数的复共轭时，该滤波器能够给出最大的信噪比。多正弦窗谱估计具有较小的偏差。同时，由于进行了多个特征谱的加权平均，对谱图会有一定的平滑效果，因此谱估计的方差性能较传统的周期图法会有显著的改善。同时能够保证一定的频率分辨能力。
技术实现要素：
本发明的目的在于提供一种弱信号目标检测的优化方法，以便更好地针对弱信号目标检测进行优化，采用更好的算法予以优化。为了实现上述目的，本发明的技术方案如下。一种弱信号目标检测的优化方法，其具体步骤如下：(1)读取离散化后的输入信号：x(n)n＝0，1，2……N-1，其中N为信号x(t)的采样点数：(2)对输入信号x(n)进行四重自相关，具体方法为，设信号x(n)包括有用信号s(n)和随机加性噪声u(n)，信号表达式为：x(n)＝s(n)+u(n)，首先对输入信号x(n)进行二重自相关：rxx(m)=limN→∞1NΣn=0N-1x(n)x(n+m)=limN→∞1NΣn=0N-1[s(n)+u(n)][s(n+m)+u(n+m)],m=0,1,2......N-1]]>根据信号与信号相关，信号与噪声的不相关性，有：rxx(m)=limN→∞1NΣn=0N-1s(n)s(n+m)]]>即：rxx(m)＝rss(m)，表明含噪信号x(n)的二重相关近似等于有用信号s(n)的二重自相关。信号x(n)的三重自相关可表示为：rxxx(m)=limN→∞1NΣn=0N-1rxx(n)x(n+m)=limN→∞1NΣn=0N-1rxx(n)[s(n+m)+u(n+m)]]]>根据rxx(n)与s(n)相关，rxx(n)与噪声的不相关性，有：rxxx(m)＝rsss(m)表明含噪信号x(n)的三重相关近似等于有用信号s(n)的三重自相关。同理，可得，信号x(n)的四重相关可表示为：rxxxx(m)＝rssss(m)即含噪信号x(n)的四重相关近似等于有用信号s(n)的四重自相关。(3)计算相关函数rxxxx(m)和输入信号x(n)的频谱：R(ωk)=Σm=0N-1rxxxx(m)exp(-j2πNkm),k=0,1,2......N-1]]>X(ωk)=Σn=0N-1x(n)exp(-j2πNkn),k=0,1,2......N-1]]>(4)由相关信号的频谱，根据信号的频谱与自相关函数频谱的关系，利用公式由四重相关信号频谱计算有用信号s(n)的频谱；(5)计算匹配滤波器的传输函数：当线性滤波器传输函数为输入信号频谱函数的复共轭时，该滤波器能够给出最大的信噪比，这种滤波器称为匹配滤波器；匹配滤波器的传输函数为：H(ωk)＝X*(ωk)；(6)求信号经过匹配滤波器的频谱：经过滤波后的信号x1(n)的频谱为：X1(ωk)＝S(ωk)H(ωk)(7)采用多正弦窗进行谱估计：窗函数取为：ωl(n)=2N+1sin[π(l+1)nN+1],n=0,1,2......N-1]]>设K为正弦窗数量，μk为第k个正弦窗的加权系数，则经过滤波后的信号x1(n)的多正弦窗谱估计为：Xmt(ωk)=Σl=1KμlXlmt(ωk)]]>其中为第k个特征谱。而Xlmt(ωk)=|Σn=1Nx1(n)ωl(n)exp(-j2πNkn)|2,k=0,]]>1，2……N-1所以多正弦窗谱估计可表示为：Xmt(ωk)=Σl=1KπμlN+1|X1(ωk+π(l+1)N+1)-X1(ωk-π(l+1)N+1)|2]]>(8)将信号频谱还原为时域信号：x2(n)=1NΣk=0N-1Xmt(ωk)exp(j2πNkn),n=0,1,2......N-1]]>该发明的有益效果在于：本发明方法中，利用四重相关和三谱(四重相关函数的频谱函数)分析，实现高阶矩范围内信号处理，由于零均值的高斯平稳随机过程的四重相关等于零，以及四重相关具有位移和旋转不变性等特点，并且三谱含有丰富的冗余信息，采用匹配滤波器进行滤波，该滤波器能够给出最大的信噪比。因此，四重相关匹配滤波技术较之二重相关匹配滤波技术，可以进一步抑制噪声，从而提高信噪比。多正弦窗谱估计在有效改进谱估计的方差性能的同时，能够保证一定的频率分辨能力。附图说明图1是本发明实施例中所使用算法流程图。具体实施方式下面结合附图对本发明的具体实施方式进行描述，以便更好的理解本发明。实施例本实施例中的弱信号目标检测的优化方法，其流程如图1所示，具体步骤如下：(1)读取离散化后的输入信号：x(n)n＝0，1，2……N-1，其中N为信号x(t)的采样点数；(2)对输入信号x(n)进行四重自相关，信号x(n)的四重相关可表示为：rxxxx(m)＝rssss(m)即含噪信号x(n)的四重相关近似等于有用信号s(n)的四重自相关。用Matlab语言实现输入信号x(n)二重相关：R2m＝autocorr(xn)；四重相关：R4m＝autocorr(R2m)；(3)计算相关函数rxxxx(m)和输入信号x(n)的频谱：R(ωk)=Σm=0N-1rxxxx(m)exp(-j2πNkm),k=0,1,2......N-1]]>X(ωk)=Σn=0N-1x(n)exp(-j2πNkn),k=0,1,2......N-1]]>用Matlab语言计算相关函数rxxxx(m)的频谱：n＝0：N-1；k＝0：N-1；％序列范围WN＝exp(-j＊2＊pi/N)；％设定wn因子nk＝n′＊k；％赋值nk转置为n行1列矩阵与k相乘WNnk＝WN.^nk；％矩阵点幂Rk＝r4m＊WNnk；％相关函数rxxxx(m)的频谱Xk＝xn＊WNnk；％输人信号x(n)的频谱(4)由相关信号的频谱，根据信号的频谱与自相关函数频谱的关系，利用公式由四重相关信号频谱计算有用信号s(n)的频谱；用Matlab语言计算有用信号s(n)的频谱：Rk_abs＝abs(Rk)；％R(ωk)的幅频Sk_abs＝Rk_abs.^(1/4)；％S(ωk)的幅频Sk_angle＝angle(Rk)/4；％S(ωk)的相频(5)计算匹配滤波器的传输函数：当线性滤波器传输函数为输入信号频谱函数的复共轭时，该滤波器能够给出最大的信噪比，这种滤波器称为匹配滤波器；匹配滤波器的传输函数为：H(ωk)＝X*(ωk)用Matlab语言计算匹配滤波器的传输函数：Hk＝conj(Xk)；(6)求信号经过匹配滤波器的频谱：经过滤波后的信号x1(n)的频谱为：X1(ωk)＝S(ωk)H(ωk)用Matlab语言计算经过滤波后的信号x1(n)的频谱为：X1k＝Sk.＊Hk；(7)采用多正弦窗进行谱估计：窗函数取为：ωl(n)=2N+1sin[π(l+1)nN+1],n=0,1,2......N-1]]>设K为正弦窗数量，μk为第k个正弦窗的加权系数，则经过滤波后的信号x1(n)的多正弦窗谱估计为：Xmt(ωk)=Σl=1KμlXlmt(ωk)]]>其中为第k个特征谱。而Xlmt(ωk)=|Σn=1Nx1(n)ωl(n)exp(-j2πNkn)|2,k=0,]]>1,2......N-1]]>所以多正弦窗谱估计可表示为：Xmt(ωk)=Σl=1KπμlN+1{X1(ωk+π(l+1)N+1)-X1(ωk-π(l+1)N+1)|2]]>用Matlab语言计算多正弦窗谱估计：K1＝0：N-1；％序列范围dw＝2＊pi/N；％频率分辨率dwk＝K1＊dw；％生成频率序列dwi＝ones(1，N)；％产生全为1的序列Fork＝1：NXmtk(k)＝0；Forl＝1：Kwi＝dwk(k)+pi＊(I+1)/(N+1)；％多弦窗估计的偏移后的频率Xmtk1_abs＝interp1(dwk，abs(X1k)，wi)；％插值计算的幅频Xmtk1_angle＝interp1(dwk，angle(X1k)，wi)；％的相频Xmtk2_abs＝interp1(dwk，abs(X1k)，wi)；％的幅频Xmtk2_angle＝interp1(dwk，angle(X1k)，wi)；％的相频Xmtk1＝Xmtk1_abs＊exp(j＊Xmtk1_angle)；％的复数表示Xmtk2＝Xmtk2_abs＊exp(j＊Xmtk2_angle)；％的复数表示Window1＝Xmtk1+Xmtk2；％X1(ωk+π(l+l)N+1)+X1(ωk-π(l+1)N+1)]]>Window1_abs1＝abs(Window1)^2；％|X1(ωk+π(l+l)N+1)+X1(ωk-π(l+1)N+1)|2]]>Window1_abs2＝window1_abs＊pi＊u(I)/(N+1)；％πμlN+1|X1(ωk+π(l+1)N+1)-X1(ωk-π(l+1)N+1)|]]>Xmtk(k)＝Xmtk(k)+Window1_abs2；％多正弦窗谱估计Xmt(ωk)=Σl=1KπμlN+1|X1(ωk+π(l+1)N+1)-X1(ωk-π(l+1)N+1)|2]]>EndEnd(8)将信号频谱还原为时域信号：x2(n)=1NΣk=0N-1Xmt(ωk)exp(j2πNkn),n=0,1,2......N-1]]>用Matlab语言将信号频谱还原为时域信号：n＝0：N-1；k＝0：N-1；％序列范围WN1＝exp(j＊2＊pi/N)；％设定wn因子kn＝k′＊n；％赋值nk转置为n行1列矩阵与k相乘WNnk1＝WN1.^kn；％矩阵点幂x2n＝Xmtk＊WNnk1；％求出时域信号x2(n)以上所述是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也视为本发明的保护范围。