自适应环境变化的恒虚警检测方法
【技术领域】
[0001] 本发明设及一种恒虚警检测方法，尤其设及一种自适应环境变化的恒虚警检测方 法。
【背景技术】
[0002] 现有的恒虚警检测算法很多，如单元平均、平均选大、平均选小、有序统计、杂波图 法等，它们都有各自的适用范围。单元平均对数据特性要求较高，主要基于W下两个假设： a.目标是独立的；b.参考窗内的所有干扰数据样本是独立同分布的，且和包含目标的单元 内的干扰同分布。而有序统计主要目的在于抑制遮蔽效应所引起的性能恶化，但有序统计 虚警率较高。杂波图法需要对目标所处环境的杂波特性具有先验知识。综上，在实际应用 中应根据实际情况合理选择恒虚警检测算法。
[0003] 单元平均在非均匀背景中性能严重下降，平均选大在多目标环境中检测性能下降 非常严重，平均选小在均匀杂波背景中检测性能和杂波边缘中虚警性能都很差。
【发明内容】
[0004] 为解决现有技术存在的上述问题，本发明提供了一种自适应环境变化的恒虚警检 巧巧法，利用了空间和时间上的相关信息，既能在均匀背景中利用单元平均恒虚警使得性 能最优，又能在非均匀杂波背景中利用有序统计方法，使得检测性能和虚警性能显著提高。
[0005] 为实现上述目的，本发明所采用的技术方案是：自适应环境变化的恒虚警检测方 法，具体步骤如下：
[0006] 对干扰存储器进行初始化；
[0007] S1 :将第一个周期FFT数据模值共Q点复制M次，分别存储到M个长度为Q的干扰 存储器中；
[000引 S2:利用当前周期数据的频谱信息和前面M个周期的干扰频谱信息来估计当前周 期数据的干扰频谱；
[0009] S3:将当前周期数据的频谱求平方得到当前周期数据的功率谱，对计算出的干扰 频谱求平方得到干扰功率谱；
[0010] S4 :当前周期数据的预判口限的乘积因子为a，则预判口限值为：
[0011]
[0012] 该数值P将作为干扰参考窗内是否被目标信号污染的预判依据，D,b为计算出的干 扰频谱，k为干扰存储器里数据的序号；
[0013] S5 :检测参考窗长L当前数据点所在位置为i，数据长度为Q:
[0014] A、如ia，即当前需要判决的数据点前窗数据长度小于以则选择其对应干扰功率 谱数组中第i个点位置的后窗的L个数据计算判决口限值，如果后窗的数据均值大于P，则 采用有序统计CFAR方法，反之采用单元平均CFAR方法；
[0015]B、如L< =K=Q-L即当前需要判决的数据点前后窗的数据长度均大于等于L则选择其对应干扰功率谱数组中第i个点位置的前后窗共化个数据用于计算其判决口限 值；如前后窗的数据均值均大于P，则采用有序统计CFAR;如只有前（后）窗的数据均值大 于P，则采用后（前）窗数据进行单元平均CFAR;如前后窗的数据均值均小于等于P，则采 用前后窗数据进行单元平均CFAR;
[0016] C、如i〉Q-L即当前需要判决的数据点后窗数据长度小于L则选择其对应干扰功 率谱数组中第i个点位置的前窗的L个数据计算判决口限值，如果前窗的数据均值大于P， 则采用有序统计CFAR方法，反之采用单元平均CFAR方法；
[0017] S6:当前周期数据功率谱的数据点与对应的口限值做比较，选取超过口限的峰值 点作为目标：如果仅有一个数据超过口限，则该数据直接判为目标，如果有多个连续数据 超过口限，则选取其峰值点作为目标，记录并输出检测出的目标数目、目标在功率谱中的位 置、目标对应功率谱的数据大小；
[0018] S7 :判决出目标后对当前数据进行更准确的干扰估计，作为下一个周期数据口限 检测的依据，做法如下：如果当前周期的功率谱数据Pi超过判决口限值K1倍，则将该点的 数据W及得出的口限值进行加权。
[0019] 进一步的，所述步骤S1中的Q点是FFT的变换点数，为2的指数次幕。
[0020] 进一步的，所述步骤S1中的干扰存储器存放的干扰信号包括杂波和噪声。
[0021] 进一步的，所述步骤S2中估计当前数据的干扰频谱的公式为：
[0022]
[002引其中数组Di内存储当前周期的前Q点FFT数据模值，CM为干扰存储器，Dzb为计算 出的干扰频谱、k为干扰存储器里数据的序号；
[0024] 进一步的，在步骤S5中的单元平均CFAR算法过程为：将所选窗中的数据排序，分 别去除掉10% -20%的最大值和10% -20%的最小值后，做单元平均检测，输出口限值，进 行单元平均检测的公式如下： 阳0巧]
[0026] 其中，Ta,i为口限值，L'为做单元平均的数据长度，P为口限因子，Pm为虚警概 率.
[0027] 进一步的，在步骤S5中的有序统计CFAR算法过程为：分别将前后窗中数据进行升 序排列，然后分别选取第3/4位置的数据作为前窗口限值和后窗口限值，在前后窗口限值 中选取较小的口限值作为当前数据点的判决口限值；
[0028] 作为更进一步的，在步骤S7中的加权公式为：
[0029]
[0030] 将得到的干扰估计值Di'替换掉原数据中的Di，完成当前周期的干扰信号频谱估 计，并将当前周期的干扰信号频谱作为一组新的干扰参考单元Cm,替换掉干扰参考单元中 第一个周期的干扰序列Cl，而原来的干扰参考单元C冲的数据放入干扰参考单元Ck冲，所 述k= 2~M，完成干扰参考单元M个周期序列的更新。
[0031] 本发明专利由于采用W上技术方案，能够取得如下的技术效果：由于利用了空间 和时间上的相关信息，既能在均匀背景中利用单元平均恒虚警使得性能最优，又能在非均 匀杂波背景中利用有序统计方法，使得检测性能和虚警性能显著提高。通过对所处环境的 干扰进行估计并不断根据当前数据对干扰信息更新，能够自适应所处干扰环境的变化。
【附图说明】
[0032] 本发明共有附图2幅：
[0033] 图1是本发明的工作流程图；
[0034] 图2是步骤S5的工作流程图。
【具体实施方式】
[0035] 下面通过具体实施例，并结合附图，对本发明的技术方案作进一步的解释说明。
[0036] 如图1所示，本发明提供的是：自适应环境变化的恒虚警检测方法，具体步骤如 下：
[0037] 对干扰存储器进行初始化；
[0038] S1 :将第一个周期FFT数据模值的共Q点复制M次，分别存储到M个长度为Q的干 扰存储器中；Q点是FFT的变换点数，为2的指数次幕；本设计中杂波RAM数量默认为5,即 M= 5,Q= 256 ;干扰存储器存放的干扰信号包括杂波和噪声；
[0039] S2:利用当前周期数据的频谱信息和前面M个周期的干扰频谱信息来估计当前周 期数据的干扰频谱；估计过程如式：
[0040]
[0041] 其中数组Di内存储当前周期的前256点FFT数据模值，Cm为干扰存储器，D,b为计 算出的干扰频谱、k为干扰存储器里数据的序号；
[0042] 本设计中M= 5,则干扰的频谱估计公式为：
[0043]
[0044] S3:将当前周期数据的频谱求平方得到当前周期数据的功率谱，对计算出的干扰 频谱求平方得到干扰功率谱；
[0045] S4 :当前周期数据的预判口限的乘积因子为a= 1，则预判口限
[0046] 值为
[0047] 该数值将作为干扰参考窗内是否被目标信号污染的预判依据；
[0048] S5 :检测参考窗长L=50,当前数据点所在位置为i，数据长度Q:
[0049] A、如ia，即当前需要判决的数据点前窗数据长度小于以则选择其对应干扰功率 谱数组中第i个点位置的后窗的L个数据计算判决口限值，如果后窗的数据均值大于P，则 采用有