[发明专利]变周期三角波、恒频体制恒虚警检测及数据处理方法

权利要求书

1.变周期三角波、恒频体制恒虚警检测及数据处理方法，其特征在于，包括：

S1：对不同天线时域数据去直流；

S2：不同天线的时域数据作FFT变换，并进行频谱移位；

S3：频谱移位后对其中一个天线的数据进行门限检测，并记录目标数目及目标在频谱中的位置信息；

S4：设置系统参数，系统带宽B、工作频率f0、测距范围Rmin～Rmax、测速范围Vmin～Vmax，三角波周期、三角波的采样率fs及恒频波的采样率fs1、FFT变换点数N_FFT；

S5.通过三角波上、下扫频部分进行目标配对得到目标的距离及速度信息，并与恒频波部分检测得到的目标速度进行匹配，得到所需要的目标距离及速度，并记录该目标在频谱中的位置；

S6.对于静止目标，根据步骤S5得到的目标在三角波上扫频或者下扫频部分频谱中的对应位置，和根据比相测角，取不同天线的对应点进行测角；

对于运动目标，根据步骤S5得到的目标在恒频波部分频谱中的对应位置，和根据比相测角，取不同天线的对应点进行测角；

三角波为三个，三个三角波的周期，分别为2*T1、2*T2、2*T3；

步骤S3中门限检测，采用的方法是单元平均选小门限检测，具体步骤如下：

第一步.对进行FFT变换后的频谱取平方，得到功率谱信息；

第二步.设置系统参数：虚警概率Pfa，前后参考窗的长度N_qc、N_hc，保护窗的长度N_pro，根据测速范围计算恒频波部分检测点对应的最小位置kmin_hp和最大位置kmax_hp，根据测距及测速范围计算三角波部分检测点对应的最小位置kmin_up、kmin_down和最大位置kmax_up、kmax_down，设光速为c；

第三步.对于恒频波频谱，判断当前检测点是否大于kmin且小于kmax；

如果当前检测点不在kmin-kmax之间，则不进行门限检测；

如果当前检测点在kmin-kmax之间，则对前后窗数据求平均之后进行比较，选取平均值小的参考窗数据进行门限检测；

第四步.对第三步中选取的平均值乘以门限因子得到门限值；

其中，P_fa＝10^-10

其中，Bk为选取的参考窗内数据，L为参考窗长度，β为门限因子，Pfa为虚警概率，T_mx即为求得的门限值；

第五步.将该门限值与功率谱对应点进行比较，记录超过门限的检测点位置及数目；

第二步中恒频波的最小检测点位置为：

kmin_hp＝N_FFT/2+Vmin*N_FFT*f0/(fs1*c)；

最大检测点位置为：

kmax_hp＝N_FFT/2+Vmax*N_FFT*f0/(fs1*c)；

分别对上式中的kmin和kmax取整数部分，即得所需要的位置编号。

2.根据权利要求1所述的变周期三角波、恒频体制恒虚警检测及数据处理方法，其特征在于，第二步中第一个三角波信号周期为2*T1，

对应第一个三角波上扫频的最小检测点位置为:

kmin_up＝N_FFT/2-2*Rmax*B*N_FFT/(c*T1*fs)+2*f0*N_FFT*Vmin/(fs*c)

第一个三角波上扫频的最大检测点位置为：

kmax_up＝N_FFT/2-2*Rmin*B*N_FFT/(c*T1*fs)+2*f0*N_FFT*Vmax/(fs*c)

第一个三角波下扫频的最小检测点位置为:

kmin_down＝N_FFT/2+2*Rmin*B*N_FFT/(c*T1*fs)+2*f0*N_FFT*Vmin/(fs*c)

第一个三角波下扫频的最大检测点位置为：

kmax_down＝N_FFT/2+2*Rmax*B*N_FFT/(c*T1*fs)+2*f0*N_FFT*Vmax/(fs*c)

分别对上面所得到的位置进行取整处理，对于第二个三角波周期2*T2、第三个三角波周期2*T3，只需要将上述公式中第一个三角波中周期T1换成T2或者T3，即可得到第二个、第三个三角波上下扫频对应的最大、最小检测点位置。