[发明专利]对密集频谱进行快速高精度细化校正的方法

权利要求书

1.一种对密集频谱进行快速高精度细化校正的方法，包括以下步骤：

(1)确定细化频带和输出点数；

(2)将细化频带转换为单位圆上的一段圆弧，确定线性调频Z变换的路径，确定起点、终点和间隔点的位置；

(3)在给定的路径上，对信号进行线性调频Z变换；

(4)根据细化频段内频率点位置和线性调频Z变换(CZT)的结果，得到相应的细化谱；

(5)构成一个N点的汉宁窗；

(6)汉宁窗对自己求卷积，得到2N-1点的卷积窗；

(7)全相位法预处理；

(8)对细化后的频谱采用全相位相位差法进行频谱校正。

2.如权利要求1所述的对密集频谱进行快速高精度细化校正的方法，其特征在于所述的步骤(2)包括如下步骤：

(21).已知信号x(n)，其Z变换的定义为：

X(z)=Σn=0∞x(n)z-n---(1)]]>

式中：x(n)为采集的信号序列；s为拉普拉斯变量，为实数，圆频率ω＝ΩT_s为一角度：

(22).令z_r＝AW^-r，式中则式中：A₀、W₀为任意的正实数。这样，给定A₀、W₀、θ₀、当r＝0，1，…，∞时，可得到在z平面上的点z₀，z₁，...，z_∞，取这些点上的Z变换，有：

X(zr)=Σn=0∞x(n)zr-n=Σn=0∞x(n)A-nWnr---(2)]]>

(23).由于nr＝[r²+n²-(r-n)²]/2，式(2)又可写成：

X(zr)=Σn=0∞x(n)A-nWr2/2Wn2/2W-(r-n)2/2=Wr2/2Σn=0∞g(n)h(r-n)=Wr2/2[g(r)*h(r)]=Wr2/2y(r)---(3)]]>

式中：g(n)=x(n)A-nWn2/2---(4)]]>

h(n)=W-n2/2---(5)]]>

y(r)＝g(r)*h(r)                           (6)

(24).当圆弧半径r＝0时，该点在z平面上的幅度为A₀，幅角为θ₀。当r＝1时，幅度变为幅角变为即在θ₀的基础上有增量不难想象，当输入x(n)为N点序列，输出X(z_r)为M点序列时，随着r的变化，点z₀,z₁,...,z_M-1就构成了CZT变换的路径，为一条螺旋线；

(25).频率轴上0～f_s/2对应单位圆上幅角0～πrad，因此假定欲细化的频带为0≤f₁＜f₂≤f_s/2，且有M条独立谱线，则对应单位圆上幅角范围为2πf₁/f_s～2πf₂/f_s，圆弧上也有M点取值；为此CZT的路径为单位圆上的起点为终点为zM-1=ej2πf2/fs,]]>间隔为Δz=ej2π(f2-f1)/[(M-1)fs]]]>的一段圆弧。