[发明专利]一种基于抛物线拟合的单脉冲鉴相方法

权利要求书

1.一种基于抛物线拟合的单脉冲鉴相方法，其特征在于，包括以下步骤：

步骤1，利用数字信道化测量主天线中频信号的脉冲标记和中频在信道化中的通道位置，利用数字信道化输出的通道号与FFT输出的序号的对应关系，计算出通道位置对应在512点FFT时的序号，然后将序号左右各扩展3-5个点后作为FFT最大值搜索范围；

步骤2，对FPGA输入主天线中频信号做匹配延时，与数字信道化输出脉冲标记对齐，并将该中频输入信号转换为四相，当该中频有效信号不足512点时，将不足部分置零后作为测频的中频输入；

步骤3，利用两组128点复数FFT运算，对输入的四相中频信号做FFT，将两路复数FFT运算结果合成512点FFT，并分两项将前256点有效结果输出；

步骤4，对并行两路的FFT结果串形求取其模，然后利用串形比较求取最大值和左右次最大值以及对应的序号，最后在两路中并行求取最大值和左右次最大值以及对应的序号；

步骤5，利用最大值和左右次最大值以及对应的序号对512点FFT结果进行第一次拟合，得出结果大于或等于-0.5且小于或等于0.5的FFT序号修正值Delta，然后对Delta附近的两点做单点DFT运算；

步骤6，利用DFT运算结果对512点FFT结果进行第二次拟合，得到第二次修正值Deltb，利用Deltb引导对干涉仪四个天线对应的中频进行512点的单点DFT运算；

步骤7，对上述四路DFT运算结果进行反正切运算，分别求出四路中频的相位，根据干涉仪基线关系求出对应天线间的相位差，最后，利用最大值序号、Deltb、采样率和FFT点数求出频率值。

2.根据权利要求1所述的方法，其特征在于，步骤1数字信道化通道号与 FFT序号的对应关系计算方法为：

步骤11，定义数字信道化实通道数为FFT点数为中频在数字信道化中的通道号为i，M₁和N₁为自然数；

步骤12，数字信道化求出的通道号i是第二奈奎斯特区，FFT序号是第一奈奎斯特区，中频在数字信道化中的通道号i对应于FFT中的序号的

3.根据权利要求1所述的方法，其特征在于，步骤3两路复数128点FFT实现512点实数FFT具体包括以下步骤：

步骤31，中频信号x(n)，四相中频信号的1、2(x(4n)、x(4n+1))两相和3、4两相(x(4n+2)、x(4n+3))分别接入两路128点复数FFT，计算出两组复数128点FFT结果，分别为Y₁(k)和Y₂(k)，j为复数量，公式如下所示：

n为自然数，N₂为FFT点数且取值为512，k为自然数且取值范围为0到

步骤32，将复数128点FFT输出与其共轭反序结果相加得到第一相中频的128点FFT结果，复数128点FFT输出与其共轭反序结果相减得到第二相中频的128点FFT结果，进而求出的四相128点四路FFT结果，第一相为X₁₁(k)， 第二相为X₁₂(k)， 第三相为X₂₁(k)， 第四相为X₂₂(k)，

步骤33，四相中频分别求取128点的FFT结果分别乘上四组不同的旋转因子，便得到512点的FFT，512点FFT只需要输出前256点的FFT结果表示为X(k)和结果如下所示：

为旋转因子。