[发明专利]一种窄带干扰频域抑制方法

权利要求书

1.一种窄带干扰频域抑制方法，其特征在于，包括如下步骤：

(1)获取准时支路滤波信号和延时支路滤波信号：

(1a)对接收端接收的时域模拟信号进行离散采样，总采样点数为N，得到时域原始采样信号，其中N＝2ⁿ，n为大于1的正整数；

(1b)将两个长度为N的Blackman-Harris窗函数序列中的一个与时域原始采样信号相乘，得到准时支路滤波信号，将另一个窗函数序列的前N/2个点并接到后N/2个点末尾，得到延时窗函数序列，再将延时窗函数序列与时域原始采样信号相乘，得到延时支路滤波信号；

(2)获取准时支路滤波信号的频谱及幅度谱和延时支路滤波信号的频谱及幅度谱：

对准时支路滤波信号和延时支路滤波信号分别进行N点快速傅里叶变换，得到准时支路频谱R₀(f)和延时支路频谱R₁(f)，并根据R₀(f)计算准时支路幅度谱φ₀(f)，根据R₁(f)计算延时支路幅度谱φ₁(f)，其中f表示R₀(f)和R₁(f)的信号频率；

(3)估计时域原始采样信号的载波频率f_c：

根据R₀(f)或R₁(f)计算边带相关值序列F(f)，并选取F(f)中最大值所对应的频率作为时域原始采样信号的载波频率f_c；

(4)获取准时支路幅度谱阈值和延时支路幅度谱阈值

计算φ₀(f)对应的准时支路幅度谱均值和方差以及φ₁(f)对应的延时支路幅度谱均值和方差并根据和计算准时支路幅度谱阈值根据和计算延时支路幅度谱阈值

(5)对准时支路频谱R₀(f)和延时支路频谱R₁(f)进行干扰抑制：

(5a)将φ₀(f)中关于f_c对称的幅度谱φ₀(f_c-σ)和φ₀(f_c+σ)逐对与准时支路幅度谱阈值进行比较，并根据比较结果对准时支路频谱R₀(f)中的对应位置R₀(f_c-σ)和R₀(f_c+σ)进行处理：

当且时，采用限幅值rR₀(f_c-σ)替换R₀(f_c-σ)，用限幅值rR₀(f_c+σ)替换R₀(f_c+σ)；

当且时，保留R₀(f_c+σ)的原值，并采用限幅值对R₀(f_c-σ)进行替换；

当且时，保留R₀(f_c-σ)的原值，并采用限幅值对R₀(f_c+σ)进行替换；

当且时，保留R₀(f_c-σ)和R₀(f_c+σ)的原值；

将经过处理的准时支路频谱R₀(f)记为R′₀(f)，并将R′₀(f)作为干扰抑制后的准时支路频谱，其中σ表示频率间隔，0≤σ≤Rs，Rs为原始采样信号码元速率，r表示限幅因子，0≤r＜1；

(5b)将φ₁(f)中关于f_c对称的幅度谱φ₁(f_c-σ)和φ₁(f_c+σ)逐对与延时支路幅度谱阈值进行比较，并根据比较结果对延时支路频谱R₁(f)中的对应位置R₁(f_c-σ)和R₁(f_c+σ)进行处理：

当且时，采用限幅值rR₁(f_c-σ)替换R₁(f_c-σ)，用限幅值rR₁(f_c+σ)替换R₁(f_c+σ)；

当且时，保留R₁(f_c+σ)的原值，并采用限幅值对R₁(f_c-σ)进行替换；

当且时，保留R₁(f_c-σ)的原值，并采用限幅值对R₁(f_c+σ)进行替换；

当且时，保留R₁(f_c-σ)和R₁(f_c+σ)的原值；

将经过处理的延时支路频谱R₁(f)记为R′₁(f)，并将R′₁(f)作为干扰抑制后的延时支路频谱，其中σ表示频率间隔，0≤σ≤Rs，Rs为原始采样信号码元速率，r表示限幅因子，0≤r＜1；

(6)计算干扰抑制后的准时支路幅度谱φ′₀(f)和延时支路幅度谱φ′₁(f)：

根据干扰抑制后的准时支路频谱R′₀(f)计算干扰抑制后的准时支路幅度谱φ′₀(f)，同时根据干扰抑制后的延时支路频谱R′₁(f)计算干扰抑制后的延时支路幅度谱φ′₁(f)；

(7)判断干扰抑制后的各支路幅度谱与对应阈值的大小：

分别比较干扰抑制后的准时支路幅度谱φ′₀(f)与准时支路幅度谱阈值以及延时支路幅度谱φ′₁(f)与延时支路幅度谱阈值的大小，当φ′₀(f)中每一点均小于等于且φ′₁(f)中每一点均小于等于执行步骤(8)，否则令φ₀(f)＝φ′₀(f)，同时令φ₁(f)＝φ′₁(f)，并执行步骤(4)；

(8)获取干扰抑制后的时域中频信号：

对R′₀(f)与R′₁(f)的和进行傅里叶反变换，得到干扰抑制后的时域中频信号。