[發明專利]分段式徒動補償方法

在審

• 申請號: 201410539267.8
• 申請日: 2014-09-27
• 公開（公告）號: CN104281781A
• 公開（公告）日: 2015-01-14
• 發明（設計）人: 鄭敏
• 申請（專利權）人: 鄭敏
• 主分類號: G06F19/00
• 分類號: G06F19/00;G01S7/40
• 代理公司: 暫無信息
• 代理人: 暫無信息

1.一種分段式徒動補償方法，其特征在于，所述方法包括：

步驟1、分別將直達波信號幀s(n)和目標回波信號幀y(n)劃分為N₁個子幀，每個子幀的積累點數為N₂，且N₂滿足2v_maxN₂/f_s＜c/B及其中，v_max為目標最大可能速度，β_0，max為最大的時延變化率速度，B為發射信號帶寬，f_s為采樣率，并將每個子幀的序號記為n₁＝0，...，N₁-1，每個子幀積累點的序號記為n₂＝0，...，N₂-1，每幀信號采樣點的序號為n＝N₂n₁+n₂；

步驟2、令N＝N₁×N₂，則n＝0，...，N-1，設一幀的積累時間為T，則N＝T×f_s，則將直達波信號幀和目標回波信號幀進行p階分數階傅里葉域變換后輸出的頻譜序列的序列號為k＝N₁k₂+k₁，則計算直達波信號幀s(n)和目標回波信號幀y(n)的p階分數階互模糊函數為：

χp(m,k2,k1)=Ap(N)ejπ((N1k2+k1)Δu)2cotαΣn1=0N1-1Σn2=0N2-1ΩWN2n2k2WN1n1k1,]]>其中，

Ap(N)=(sinα+jcosα)/N,α=Δpπ/2,]]>

Δt和Δu分別是時域和p階分數階傅里葉域的采樣間隔，有Δu×Δt＝sinα/N；

從而得到Ω=s(n1,n2)y(n1,n2)ejπ((N2n1+n2)Δt)2cotαWNn2k1,]]>

步驟3、由于滿足則忽略子幀時間內的距離徒動，則：y(n1,n2)≈A0x(n1,n2-(m0-β0N2n1))ej2πkdc(N2n1+n2)ejπfdr(N2n1+n2)2Δt2+ω2(n1,n2)]]>其中，A₀為回波信號幅值，m₀為離散時延，k_dc＝f_dcΔt，其中，f_dc為多普勒頻率，f_dr為多普勒調頻率；

計算y(n₁，n₂)的階分數階互模擬函數得到

χp^(m,k2,k1)=A0Ap(N)ej2πN1k2k1Δu2cotα^·]]>

Σn1=0N1-1{Σn2=0N2-1[x*(n1,n2-m)x(n1,n2-m0+β0N2n1)ej2πkdc(N2n1+n2)WN2n2k1]WN2n2k2}WN1n1k1]]>

令

z(n1,m,k2,k1)=Σn2=0N2-1[x*(n1,n2-m)x(n1,n2-m0+β0N2n1)ej2πkdc(N2n1+n2)WN2n2k1]WN2n2k2]]>

步驟4、對z(n₁，m，k₂，k₁)中的m進行離散時間傅里葉變換，得到

Z(n1,l,k2,k1)=X*(n1,l)X(n1,lM-ξ)ej2πn1N2β0(kc+lM)ej2πm0N2(ξ-N2lM)ej2πn1ξβ0]]>其中，M為時延總點數；β₀為時延變化率；k_c＝f_cΔt；f_c為雷達信號載波頻率；

步驟5、對Z(n₁，l，k₂，k₁)的n₁進行keystone變換，并令有Z(n^1,l,k2,k1)=X·(n^1,l)X(n^1,N2lM-ξ)ej2πkdn^1N2ej2πm0N2(ξ-N2lM)ej2πβ0kcn^1ξl/M+kc]]>

步驟6、再將進行第二次FFT變換回時域

χp^(m,k2,k1)≈A0Ap^(N)ej2πN1k2k1Δu2cotα^·]]>

Σn^1=0N1-1{Σn2=0N2-1[x*(n^1,n2-m)x(n^1,n2-m0)ej2πkdc(N2n^1+n2)ej2πn^1ξβ0WN2n2k1]WN2n2k2}WN1n^1k1]]>從而使得距離徒動得到補償。