[发明专利]一种Loran-C定时和定位终端中连续波干扰检测方法

权利要求书

1.一种Loran-C定时和定位终端中连续波干扰检测方法，其特征在于，包括以下步骤：

步骤1，对Loran-C台站的脉冲组信号进行轮询捕获，确定出能够捕获到的Loran-C脉冲组位置并在时域进行标记，确定脉冲组的起始和结束概位；没有捕获到脉冲组信号则整个时域为静默区；

步骤2，在静默区内提取M组长度为N的连续数据序列x_i(n)，i表示第i组数据，i＝1…M；n表示每组数据序列中的第n个数据，n＝1…N；

步骤3，对M组数据中每组长度为N的数据序列x(n)加窗后进行快速傅里叶变换并取绝对值，得到数据序列x_i(n)的幅度谱X_i(n)；

步骤4，对M组幅度谱数据序列X_i(n)对应点数据进行累加取平均，得到非相干累积结果

步骤5，对数据X′(n)滑动中值滤波，输出数据序列为X″(n)；

步骤6，采用数据拟合方法对数据序列X″(n)进行拟合，得到一条光滑的趋势项拟合曲线

步骤7，利用非相干累积结果X′(n)减去趋势项拟合结果获得的噪低平稳的谱线

步骤8，对谱线数据Y(n)进行自适应门限计算；

步骤9，利用门限T_h与谱线数据Y(n)进行比较对干扰谱线脉冲检测，大于门限T_h则认为存在干扰谱，取一段大于门限数据的中间值作为连续波干扰频点，进行连续波干扰抑制。

2.根据权利要求1所述的Loran-C定时和定位终端中连续波干扰检测方法，其特征在于，所述的步骤2中长度N根据实际分辨率要求选取，分辨率计算公式为f_s/N，f_s为AD采样率。

3.根据权利要求1所述的Loran-C定时和定位终端中连续波干扰检测方法，其特征在于，所述的步骤5对数据X′(n)滑动取2m+1数据，m为阶数；滑动中值滤波按照正序/逆序方式排列数据，取中点数据作为滤波输出，输出数据序列为X″(n)，

X″(n)＝median(|X′(n-m)|,|X′(n-m+1)|,…|X′(n)|,…,|X′(n+m-1)|,|X′(n+m)|)

其中，median(*)表示对2m+1各数据滑动中值滤波。

4.根据权利要求1所述的Loran-C定时和定位终端中连续波干扰检测方法，其特征在于，所述的步骤6采用的拟合方法为最小二乘法拟合或多项式拟合。

5.根据权利要求1所述的Loran-C定时和定位终端中连续波干扰检测方法，其特征在于，所述的步骤8对谱线数据Y(n)进行升序排列，取前N/2点的谱线数据通过计算方差设置门限，门限式中为N/2点谱线的方差；c为门限系数。

6.根据权利要求5所述的Loran-C定时和定位终端中连续波干扰检测方法，其特征在于，所述的门限系数c取2～10。

7.根据权利要求1所述的Loran-C定时和定位终端中连续波干扰检测方法，其特征在于，所述的步骤9选取级联数字IIR或自适应LMS陷波器进行连续波干扰抑制。

8.根据权利要求1所述的Loran-C定时和定位终端中连续波干扰检测方法，其特征在于，在Loran定时或定位终端中连续重复进行步骤1～9，将连续波干扰抑制后的数字信号继续进行捕获，找出存在的脉冲组信号后重新确定静默区，并对静默区的AD数据进行连续波干扰检测，检测后进行抑制，并将抑制后的信号送入下一级Loran-C信号处理模块使用。