主权项 |
1.一种低采样率线扫频脉冲时延估计方法,其特征在于,包括如下步骤:步骤一、对接收到的线扫频雷达回波信号以△t为间隔进行时域采样,得到采样序列x(n),其中<img file="FDA00001783917200011.GIF" wi="259" he="68" />-a是线扫频脉冲雷达信号的调频率,a>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>Σ</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>·</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>·</mo><mi>cot</mi><mrow><mo>(</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mi>pπ</mi><mo>)</mo></mrow><mo>·</mo><msup><mi>n</mi><mn>2</mn></msup><mo>·</mo><mi>Δ</mi><msup><mi>t</mi><mn>2</mn></msup><mo>-</mo><mi>j</mi><mo>·</mo><mfrac><mrow><mn>2</mn><mi>π</mi></mrow><mi>N</mi></mfrac><mo>·</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>π</mi><mo>·</mo><msub><mi>m</mi><mn>0</mn></msub></mrow><mrow><mi>N</mi><mo>·</mo><mi>cot</mi><mrow><mo>(</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mi>pπ</mi><mo>)</mo></mrow><mo>·</mo><mi>Δ</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>ΔN</mi><mo>=</mo><mfrac><mrow><mn>2</mn><mi>π</mi></mrow><mrow><mi>cot</mi><mrow><mo>(</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mi>pπ</mi><mo>)</mo></mrow><mo>·</mo><mi>Δ</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]。 |