[發明專利]一種線性調頻干擾對消方法

有效

• 申請號: 202010801764.6
• 申請日: 2020-08-11
• 公開（公告）號: CN111917676B
• 公開（公告）日: 2021-07-27
• 發明（設計）人: 申冰冰;程郁凡;賈榮波;吳垚忠;柏果;何金埔
• 申請（專利權）人: 電子科技大學
• 主分類號: H04L27/10
• 分類號: H04L27/10;H04L27/16
• 代理公司: 成都點睛專利代理事務所(普通合伙) 51232
• 代理人: 孫一峰

1.一種線性調頻干擾對消方法，假設發射信號s(t)經過高斯白噪聲信道和線性調頻干擾信號J(t)之后得到接收信號r(t)，r(t)＝J(t)+s(t)+w(t)，w(t)是零均值加性復高斯白噪聲，實部和虛部方差均為σ²；其特征在于，所述對消方法包括以下步驟：

S1、以采樣頻率f_s對r(t)進行均勻采樣，得到具有N個采樣點的離散信號r[n]，其中A為干擾幅度、f₀為干擾初始頻率、為初相、μ為調頻率，n＝0,1,2…N-1，w[n]～CN(0,2σ²)；

S2、根據離散信號r[n]得到調頻率μ的估計值具體包括：

S21、初始化變量迭代次數m＝0，中間變量r₁[n]＝r[n]，其中為向下取整；

S22、判斷m是否小于M，若是，轉S23；否則，轉S211；

S23、令v＝2^m，計算r₁[n]的自相關序列y_m[n]：

y_m[n]＝r₁[n+v]×r₁[n]^*

其中[·]^*為取共軛，n為正整數，0≤n≤N-2^m-1；

S24、判斷m是否等于0，若是，轉S25；否則轉S26；

S25、利用如下公式，

得到L_m點頻譜Y_m＝[Y_m[0],Y_m[1],…,Y_m[L_m-1]]，其中L_m＝N-2^m，m＝0,1,2,…M-1；

S26、利用如下公式：

得到整數k_m的估計值0≤k_mL_m；

S27、利用如下公式：

得到小數偏差δ_m的估計值其中為取實部，|δ_m|≤0.5；

S28、利用如下公式

得到μ_m的估計值

S29、利用如下公式

對r₁[n]更新；

S210、m＝m+1，回到步驟S22；

S211、利用如下公式：

得到μ_m的最終估計值其中v＝2^m；

S3、根據得到的對離散信號r[n]進行調頻率消除得到調頻率消除后的信號r₁′[n]：

S4、忽略信號、噪聲和殘余調頻率的影響，r₁′[n]為一個復正弦信號，初始頻率為f₀，初始相位為幅度為A，表示為：

S5、對r₁′[n]進行N點FFT得到R′[k]，找到最大譜線的位置利用基于傅里葉系數的迭代頻域插值算法計算最大頻譜的左右兩個與其相距半頻譜處譜線的幅值得到頻率偏差得到初始頻率估計值

S6、計算初相估計值和幅度估計值

S7、重構線性調頻信號，得到

S8、對采樣得到的離散信號r[n]進行干擾對消：

2.根據權利要求1所述的一種線性調頻干擾對消方法，其特征在于，步驟S5所述初始頻率估計值的具體計算步驟如下：

S51、對r₁′[n]進行N點FFT得到R′[k]；

S52、利用如下公式得到R′[k]幅值最大值位置

S53、然后利用如下公式得到小數倍頻偏位置

S54、利用如下公式得到初始頻率估計值