发明名称 一种低采样率线扫频脉冲时延估计方法
摘要 本发明涉及一种低采样率线扫频脉冲时延估计方法,属于雷达信号处理领域。本发明对线扫频脉冲回波做简化分数阶傅里叶变换,通过判断峰值点周边相邻点的平均相位差所在的区间,实现脉冲时延解模糊,进而实现脉冲时延的估计。本发明提出的线扫频脉冲时延估计方法解决了低采样率条件下信号时延估计模糊问题,有效降低了接收信号的采样率和后续信号处理的运算量,并且可以通过快速傅里叶变换算法实现,计算复杂度低,为舰载雷达、合成孔径雷达等线扫频脉冲体制雷达信号处理提供了有效的工具。
申请公布号 CN102226839B 申请公布日期 2012.10.03
申请号 CN201110150808.4 申请日期 2011.06.07
申请人 北京理工大学 发明人 陶然;孟祥意;李元硕;郇浩
分类号 G01S7/285(2006.01)I 主分类号 G01S7/285(2006.01)I
代理机构 代理人
主权项 1.一种低采样率线扫频脉冲时延估计方法,其特征在于,包括如下步骤:步骤一、对接收到的线扫频雷达回波信号以△t为间隔进行时域采样,得到采样序列x(n),其中<img file="FDA00001783917200011.GIF" wi="259" he="68" />-a是线扫频脉冲雷达信号的调频率,a&gt;0;由脉冲重复时间T和脉冲宽度△T得到该采样序列x(n)的序列长度为N=T/Δt、脉宽长度为M=△T/Δt;步骤二、由步骤一确定的x(n)的序列长度N,对步骤一得到的采样序列x(n)进行p阶N点简化分数阶傅里叶变换,得到X<sub>p</sub>(m),即<maths num="0001"><![CDATA[<math><mrow><msub><mi>X</mi><mi>p</mi></msub><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><msup><mi>e</mi><mrow><mi>j</mi><mo>&CenterDot;</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>&CenterDot;</mo><mi>cot</mi><mrow><mo>(</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mi>p&pi;</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><msup><mi>n</mi><mn>2</mn></msup><mo>&CenterDot;</mo><mi>&Delta;</mi><msup><mi>t</mi><mn>2</mn></msup><mo>-</mo><mi>j</mi><mo>&CenterDot;</mo><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><mi>N</mi></mfrac><mo>&CenterDot;</mo><mi>mn</mi></mrow></msup><mi>x</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>其中,p=2arccot(a)/π;步骤三、搜索步骤二得到的X<sub>p</sub>(m)中幅值即|X<sub>p</sub>(m)|最大的点,并得到该点的坐标m<sub>0</sub>,由所采用的简化分数阶傅里叶变换阶次p和步骤一确定的序列长度N、时域采样间隔△t,按照下式获得线扫频脉冲的无模糊时延n<sub>0</sub>,即<maths num="0002"><![CDATA[<math><mrow><msub><mi>n</mi><mn>0</mn></msub><mo>=</mo><mfrac><mrow><mn>2</mn><mi>&pi;</mi><mo>&CenterDot;</mo><msub><mi>m</mi><mn>0</mn></msub></mrow><mrow><mi>N</mi><mo>&CenterDot;</mo><mi>cot</mi><mrow><mo>(</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mi>p&pi;</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>&Delta;</mi><msup><mi>t</mi><mn>2</mn></msup></mrow></mfrac><mo>;</mo></mrow></math>]]></maths>步骤四、在步骤三得到的X<sub>p</sub>(m<sub>0</sub>)左右两边各取最接近X<sub>p</sub>(m<sub>0</sub>)的b个点,然后求X<sub>p</sub>(m<sub>0</sub>-b)、…、X<sub>p</sub>(m<sub>0</sub>-1)、X<sub>p</sub>(m<sub>0</sub>)、X<sub>p</sub>(m<sub>0</sub>+1)、…、X<sub>p</sub>(m<sub>0</sub>+b)这2b+1个点中每两个相邻的点之间的相位差,然后对这2b个相位差求平均值<img file="FDA00001783917200014.GIF" wi="127" he="57" />步骤五、由步骤一得到的序列长度N、脉宽长度M、步骤四得到的相位差平均值<img file="FDA00001783917200015.GIF" wi="131" he="57" />步骤三得到的无模糊时延n<sub>0</sub>按照下式计算修正的平均相位差:<img file="FDA00001783917200016.GIF" wi="949" he="138" />其中,[·]<sub>mod2π</sub>表示关于2π的求余运算;步骤六、由所采用的简化分数阶傅里叶变换阶次p和步骤一确定的时域采样间隔△t确定在该时域采样间隔△t条件下的模糊区间长度△N,即<maths num="0003"><![CDATA[<math><mrow><mi>&Delta;N</mi><mo>=</mo><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><mrow><mi>cot</mi><mrow><mo>(</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mi>p&pi;</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>&Delta;</mi><msup><mi>t</mi><mn>2</mn></msup></mrow></mfrac></mrow></math>]]></maths>步骤七、由步骤一确定的序列长度N、步骤五得到的修正平均相位差<img file="FDA00001783917200018.GIF" wi="127" he="59" />步骤六确定的模糊区间长度△N,计算脉冲时延所在的模糊区间序号A,即<img file="FDA00001783917200021.GIF" wi="755" he="142" />其中,<img file="FDA00001783917200022.GIF" wi="56" he="67" />表示下取整运算;步骤八、由所采用的简化分数阶傅里叶变换阶次p、步骤一确定的时域采样间隔Δt、步骤三得到的无模糊时延n<sub>0</sub>、步骤七得到的模糊区间序号A得到脉冲的时延为τ=n<sub>0</sub>·Δt+A·Δτ其中Δτ=2π/[cot(pπ/2)△t]。
地址 100081 北京市海淀区中关村南大街5号