CN102608594B - 一种调频步进信号合成高分辨像栅瓣抑制方法 - Google Patents

Abstract

本发明公开了一种调频步进信号合成高分辨像栅瓣抑制方法，首先对采样离散化后的调频步进信号的回波信号在快时间维进行加权失配滤波处理，加权函数满足

的条件，Δf为调频步进信号的步进间隔，B₁为调频步进信号的子脉冲带宽，K为加权函数相对于矩形窗的主瓣展宽系数；对加权失配滤波处理后结果在慢时间维进行合成子距离像处理，对子距离像进行合成全程一维距离像处理。使用本发明能够在不对***参数设计提出苛刻要求情况下，有效抑制调频步进信号合成高分辨像的时间栅瓣，提高频率步进雷达的成像及检测性能。

Description

一种调频步进信号合成高分辨像栅瓣抑制方法

技术领域

本发明涉及一种频率步进信号合成高分辨像栅瓣抑制方法，更具体地说，涉及一种调频步进信号合成高分辨像栅瓣抑制方法。

背景技术

调频步进信号是一种具有高分辨能力的雷达信号，它发射一组脉间载频线性跳变、脉内线性调频的雷达脉冲，通过对各个子脉冲回波进行脉内脉冲压缩和脉间逆傅立叶变换(IDFT)处理获得目标的距离高分辨像。调频步进信号用线性调频子脉冲代替了频率步进信号中的简单脉冲，因此克服了雷达作用距离和单脉冲距离分辨力之间的矛盾，又能够在保持距离分辨力不变的条件下减少脉冲数，提高数据率，降低对***处理带宽和采样率的要求。

为提高***分辨率需要增加频率步进信号的合成带宽即增大步进间隔Δf。但当子脉冲时宽T与步进间隔乘积大于1(TΔf＞1)时，等间隔跳频频率步进信号存在严重的时间栅瓣问题。频率步进信号的周期性时间栅瓣会造成测距模糊、目标遮挡并增加***的虚警概率对成像及检测带来不利影响，因此必须对频率步进信号的时间栅瓣进行有效抑制。

目前，频率步进信号的栅瓣抑制方法很多。在2003年IEEE Transactions onAerospace and Electronic Systems第39卷第二期第694-703页上的Nullifying ACFgrating lobes in stepped-frequency train of LFM pulses(消除调频步进脉冲串自相关栅瓣)一文中，介绍了通过调整子脉冲时宽、步进间隔、子脉冲带宽三者关系将栅瓣位置处置零，从而全部或部分消除调频步进信号栅瓣。但当子脉冲时宽带宽积较大时子脉冲之间的频带存在较大交叠，该方法并不实用。在2005年IEEE International Radar Conference第371-376页上的Suppression of gratinglobes in stepped-frequency train(频率步进脉冲串栅瓣抑制)一文中，对上述方法做出了改进，其思想为不使栅瓣位置为0，而是小于一个较小的数值。该方法在一定程度上放宽了对频谱交叠的要求。但上述两种方法均要求保证严格的位置关系，因此在实际应用到成像中时均需要较高的采样率(通常为子脉冲带宽的10倍以上)以达到期望的栅瓣抑制效果。在2008年Proceedings of SPIE，theInternational Society for Optical Engineering第6947卷第1-9页上的Low sidelobenonlinear stepped-freuency waveforms(低副瓣频率步进波形)一文中，提出根据某种非线性法则改变不同脉冲间的频率步进量以打破步进项时间周期性，从而实现频率步进信号的栅瓣抑制，该方法实现过于复杂难以实际应用。本发明将加权失配滤波技术用于调频步进信号合成高分辨像算法中，有效抑制了频率步进信号的时间栅瓣，提高了频率步进雷达体制的成像及检测性能，尤其适用于调频步进信号的低时间栅瓣一维成像。

发明内容

有鉴于此，本发明的目的在于提供一种调频步进信号合成高分辨像栅瓣抑制方法，在不对***参数设计提出苛刻要求情况下，有效抑制调频步进信号合成高分辨像的时间栅瓣，提高频率步进雷达的成像及检测性能。

为了达到上述目的，本发明提供一种基于加权失配滤波的调频步进信号合成高分辨像栅瓣抑制方法，所述方法包括：

1)对采样离散化后的调频步进信号的回波信号在快时间维(采样维)进行加权失配滤波处理；其中，加权失配滤波处理所采用的加权函数需要满足如下条件：

$\frac{1}{\mathrm{\Δf}}\≥\frac{1}{B_1}(K+1)$

其中，Δf为调频步进信号的步进间隔，B₁为调频步进信号的子脉冲带宽，K为加权失配滤波时选取的加权函数相对于矩形窗的主瓣展宽系数；

2)对步骤1)处理后的信号在慢时间维(PRT维)进行合成子距离像处理；

3)对子距离像进行合成全程一维距离像处理。

较佳地，所述步骤2)具体包括：

步骤S301：对加权失配滤波后结果在慢时间维进行零插值处理；在参数设计时保证插值后点数N_zp与调频步进信号的步进间隔Δf二者乘积为***采样率的整数倍；

步骤S302，对零插值后的信号在慢时间维上做IFFT处理，得到不同采样点上的子距离像。优选地，本步骤S302进行IFFT处理后得到M个距离长度为δR_c＝c/2Δf的子距离像，子距离像中心位于R_k；参数设计时满足δR_c＞δR_g的条件；

其中，δR_g＝cT_g/2为进行快时间维处理时的粗分辨单元，T_g为调频步进信号的采样间隔，c为光速；M为调频步进信号的采样点数；R_k＝R_m+k·δR_g为成像距离窗口内第k个采样点对应的距离，R_m为成像距离窗口起始位置。

所述步骤3)具体包括：

步骤S401：舍弃子距离像中距离段[-δR_g/2+R_k，δR_g/2+R_k]以外部分，仅抽取子距离像高分辨谱线中以R_k为中心位置的

个点，其中

为单采样点所代表距离对应的高分辨谱线个数；

步骤S402：将不同子距离像的抽取结果顺次拼接并进行取模操作，得到全程一维距离像。

本发明与现有的合成高分辨像栅瓣抑制方法相比，具有以下有益效果：

(1)本发明将加权失配滤波技术用于调频步进信号合成高分辨像算法中，并设计加权函数的主瓣展宽系数K满足条件

使得步进项栅瓣不落入调频项主瓣中，有效抑制了频率步进信号的时间栅瓣，提高了频率步进雷达体制的成像及检测性能，尤其适用于调频步进信号的低时间栅瓣一维成像。本发明不需要精密设计子脉冲时宽、步进间隔、子脉冲带宽之间的关系，当调频步进信号的子脉冲时宽带宽积较大时，仍可以有效降低子脉冲间的频谱重叠和***采样率从而放宽对参数设计的限制；

(2)本发明在合成子距离像的插值过程中，设计插值后点数N_zp与调频步进信号的步进间隔Δf二者乘积为***采样率的整数倍，从而保证最终抽取结果的正确性

(3)本发明在合成子距离像时，通过参数的设计(δR_c＞δR_g)保证相邻子距离像间不出现距离空隙；在抽取时以目标实际位置作为中心进行设定点数的抽取，在保证不同子距离像的抽取结果不会重叠的基础上，使得合成高分辨像时不破坏调频步进信号模糊函数所限定的信号内部结构关系，从而保证最终合成结果的正确性，可以有效降低调频步进信号一维成像结果的时间栅瓣；

(4)实现简单、具有较强的实用性。

附图说明

图1调频步进信号模糊函数主峰零多普勒切面图。

图2通过参数设计抑制调频步进信号模糊栅瓣示例图。

图3加权失配滤波抑制调频步进信号模糊栅瓣示例图。

图4基于加权失配滤波的合成高分辨处理原理框图。

图5快时间维加权失配滤波处理框图。

图6慢时间维合成子距离像处理框图。

图7合成全程一维距离像处理框图。

图8实施例中***参数下实测数据合成高分辨像栅瓣抑制结果(左侧为常规合成高分辨像算法结果，右侧为本发明栅瓣抑制后合成高分辨像结果)

具体实施方式

首先，结合附图1至附图3阐明本发明的基本思路。

附图1、附图2、附图3均是调频步进信号模糊函数主峰零多普勒切面图，附图1是一般参数设计下的调频步进信号模糊函数主峰零多普勒切面图，附图2是通过参数设计方法抑制调频步进信号模糊栅瓣示例图，附图3是加权失配滤波抑制调频步进信号模糊栅瓣示例图。下面结合附图1、附图2、附图3对本发明的基本思路进行描述。

假定频率步进信号起始载频为f₀，步进阶梯为Δf，步进点数为N，脉冲重复周期为T_R，脉冲重复频率为PRF，子脉冲脉宽为T，子脉冲带宽为B₁，子脉冲调频斜率为k＝B₁/T，第n个子脉冲的载频为f_n＝f₀+nΔf，n＝0，1，...，N-1；目标初始距离为R₀，目标径向速度V_R；采样间隔T_g，采样点数为M，成像窗口起点τ_m，插值后的步进点数为N_zp，每组IFFT结果的抽取点数为

调频步进发射信号模型：

$S_T\left(t\right)=\underset{n=0}{\overset{N-1}{\mathrm{\Σ}}}{x}_{\mathrm{LFM}}(t-{\mathrm{nT}}_R)\exp \left[j2\mathrm{\π}(f_0+\mathrm{n\Δf})(t-{\mathrm{nT}}_R)\right]---\left(1\right)$

其中，

x_LFM(t)＝exp(jπkt²)rect(t/T)    (2)

$\mathrm{rect}\left(\frac{t}{T}\right)=\left\{\begin{array}{c}1,\left|t\right|\≤T/2\\ 0,\mathrm{others}\end{array}\right.----\left(3\right)$

调频步进信号模糊函数主峰零多普勒切面的表达式如下：

$\left|R\left(\mathrm{\τ}\right)\right|=|1-\frac{\left|\mathrm{\τ}\right|}{T}|\left|\sin c\right[B_1\mathrm{\τ}(1-\frac{\left|\mathrm{\τ}\right|}{T})\left]\right|\left|\frac{\sin \left(\mathrm{N\π\τ\Δf}\right)}{N\sin \left(\mathrm{\π\τ\Δf}\right)}\right|---\left(4\right)$

其中，τ为两目标间的相对时延。式(4)中第一项由矩形包络产生，第二项由线性调频项产生，而第三项由频率步进项产生。前两项对最终模糊函数结果产生包络调制，而第三项引起栅瓣，栅瓣位置出现在：

${\mathrm{\τ}}_n=\frac{n}{\mathrm{\Δf}},n=\±1,\±2,K,\mathrm{int}\left[\mathrm{T\Δf}\right]---\left(5\right)$

一般参数设计下的调频步进信号模糊函数主峰零多普勒切面图如附图1所示，其左图为公式(4)中三项的分解图(即分别将三项画于图中)，右图实线给出了最终模糊图结果。

已有方法抑制调频步进信号模糊栅瓣的一种基本思路为：通过调整子脉冲时宽T、步进间隔Δf、子脉冲带宽B₁三者关系将栅瓣位置处置零或使其小于一个较小的数值，如附图2所示。但这要求式(4)中的三项保证严格的位置关系，因此在实际应用到成像中时需要较高的采样率(通常为子脉冲带宽的10倍以上)以达到期望的栅瓣抑制效果，这使得工程实现代价太大，并不实用。另一种已有方法提出根据某种非线性法则改变不同脉冲间的频率步进量以打破式(4)中步进项时间周期性，从而实现频率步进信号的栅瓣抑制，但该方法实现过于复杂难以实际应用。

由附图1可以很直观的看出，在保证式(4)中的步进项栅瓣不落入调频项主瓣中(即B₁＞Δf)的前提下，只要能够压低调频项的旁瓣便可以有效的抑制一维成像结果的栅瓣，并且不需要保证严格的位置关系(不需要高采样率)。基于上述思想，本发明抑制调频步进信号模糊栅瓣的基本思路为：通过加权失配滤波压低式(4)中调频子脉冲脉压结果旁瓣，进而最终抑制频率步进信号合成高分辨成像结果栅瓣。为更直观理解上述思想，附图3给出一组参数下的仿真结果。

其次，结合附图4至附图8阐明本发明的具体实施步骤。

附图4是本发明在工程上实际实施的信号处理总框图，附图5、附图6、附图7为附图4的具体实现细化框图，附图8为本实施例中***参数下的实测数据高分辨像栅瓣抑制结果图。参照附图，对本发明进行详细的描述。

如图4所示，本发明算法以调频步进信号视频回波为输入，先后经过：快时间维加权失配滤波处理、慢时间维合成子距离像处理、合成全程一维距离像处理3个大步骤，最终获得低栅瓣全程一维距离像。下面分步骤将本发明的原理和实现方法描述如下：

步骤S101，慢时间维对视频回波进行加权失配滤波处理。

雷达发射机发射一组调频步进脉冲串信号，该信号经目标反射后，得到相对发射信号延迟为τ(t)的回波信号，将该回波信号与相参本振信号进行下变频处理，获得视频回波。

假定成像距离窗口为[R_m，R_m+ΔR]，与其对应的成像时间窗口为[τ_m，τ_m+Δτ]，其中R_m为起始距离，ΔR为成像窗口大小，τ_m＝2R_m/c，Δτ＝2ΔR/c；粗分辨单元为δR_g＝cT_g/2。

每个PRT在成像窗口内以采样间隔T_g对回波信号进行采样，距离窗口内的不同采样点对应的距离为R_k＝R_m+kδR_g，其中，k＝0，1，...，M-1并且M＝int(ΔR/δR_g)，int(·)代表取整操作。采样后每帧频率步进信号形成N×M维回波矩阵。分别对N个不同PRT的回波信号在快时间维上(采样维上)做加权失配滤波处理。

下面结合图5对该步骤的具体实现流程进行描述。

步骤S201，由式(1)经下变频处理后可得调频步进信号的视频回波表达式如下：

$S_R(t,{\mathrm{nT}}_R)=\underset{n=0}{\overset{N-1}{\mathrm{\Σ}}}{x}_{\mathrm{LFM}}(t-\mathrm{\τ}\left(t\right)-{\mathrm{nT}}_R)\exp \left[j2\mathrm{\πn\Δf}(t-\mathrm{\τ}\left(t\right)-{\mathrm{nT}}_R)\right]\exp [-j2\mathrm{\π}{f}_0\mathrm{\τ}\left(t\right)]\exp \left[j2\mathrm{\π}{f}_{\mathrm{dn}}t\right]---\left(6\right)$

其中，τ(t)和对应f_n的多普勒频率f_dn分别为：

$\mathrm{\τ}\left(t\right)=\frac{2R\left(t\right)}{c}\≈\frac{2R_0}{c}+\frac{2V_{R}t}{c}={\mathrm{\τ}}_0+\frac{2V_{R}t}{c}---\left(7\right)$

$f_{\mathrm{dn}}=-\frac{2V_R}{c}(f_0+\mathrm{n\Δf})---\left(8\right)$

当V_R＜＜c时，式(8)可化简为：

$f_{\mathrm{dn}}\≈f_d=-\frac{2V_R}{c}{f}_0---\left(9\right)$

f_d为对应f₀的多普勒频率。

步骤S202，分别对N个不同PRT的回波信号在脉内做快速傅里叶变换(FFT)。

步骤S203，第n个PRT加权失配滤波处理的参考信号表达式如下：

$S_{\mathrm{Ref}}^0(t,n{T}_R)=x_{\mathrm{LFM}}\left(t\right)\exp \left[j2\mathrm{\πn\Δft}\right]---\left(10\right)$

步骤S204，对加权失配滤波处理的参考信号做加权处理：

S_Ref(t，nT_R)＝x_WLFM(t)exp[j2πnΔft]     (11)

其中，

x_WLFM(t)＝x_LFM(t)×W(t)，W(t)为加权函数。(12)

加权失配滤波处理在压低脉压结果旁瓣的同时会造其主瓣展宽及一定的峰值幅度损失，但不会降低最终自相关结果的分辨能力。此处重点考虑失配滤波后子脉冲脉压结果主瓣展宽对频率步进信号栅瓣的影响。当展宽后的主瓣超出1/Δf时，步进项的第一栅瓣落入调频项的主瓣内，此时不但不能有效压低栅瓣，有时反而会造成栅瓣抬高。为得到较好的栅瓣抑制效果，参数设计时应保证步进项的栅瓣不落入加权失配滤波后调频项的第一旁瓣以内，这里还留出了一个栅瓣位置的裕量，因此，在加权函数W(t)选取和参数Δf、B₁设计应折中考虑，使其遵守以下原则：

$\frac{1}{\mathrm{\Δf}}\≥\frac{1}{B_1}(K+1)---\left(13\right)$

其中，K为加权失配滤波时选取的加权函数相对于矩形窗的主瓣展宽系数，即K＝τ₂/τ₁(τ₁为矩形窗主瓣宽度，τ₂为所选加权函数的主瓣宽度)。加权函数通常可使用窗函数，加权函数确定后K值确定。在实际中，式(13)所限定的加权函数选取及参数设计原则是保证本发明最终合成高分辨像的栅瓣抑制效果的一个非常重要的准则。

步骤S205，对加权处理后的参考信号在脉内做FFT。

步骤S206，对步骤S205结果取共轭。

步骤S207，将步骤S206结果分别与步骤S202各PRT回波信号的脉内FFT结果相乘。

步骤S208，将步骤S207各PRT结果分别在脉内做快速傅里叶反变换(IFFT)，IFFT后表达式如下：

(14)

其中，

$x_{\mathrm{UMF}}(t-{\mathrm{\τ}}_0-{\mathrm{nT}}_R,f_{\mathrm{dn}})=\left[x_{\mathrm{LFM}}(t-{\mathrm{\τ}}_0-n{T}_R)\exp \left(j2\mathrm{\π}{f}_{\mathrm{dn}}(t-{\mathrm{\τ}}_0-{\mathrm{nT}}_R)\right)\right]\⊗{x^{*}}_{\mathrm{WLFM}}(-t)---\left(15\right)$

x_UMF(t-τ₀-nT_R，f_dn)即为调频步进信号的线性调频子脉冲加权失配滤波结果。

现将式(14)给出的公式变为离散化形式。记t＝t_k＝kT_g+nT_R+τ_m代入式(14)并整理得：

(16)

k＝0，1，2，...，M-1；n＝0，1，2，...，N-1；

其中，A＝exp[-j2πf₀τ₀]exp[j2πf_dτ₀]。当V_RN/(PRF·δ_R)＜2时，式(14)简化为：

经过以上八个步骤，便完成了调频步进信号视频回波的加权失配滤波处理。

步骤S102，对加权失配滤波结果进行慢时间维合成子距离像处理。

下面结合图6对该步骤的具体实现流程进行描述。

步骤S301，对加权失配滤波后形成的M个长度为N的向量分别零插值到N_zp点。在参数设计时必须保证N_zpΔf为f_s＝1/T_g的整数倍，否则会造成最终抽取结果不正确，从而保证合成高分辨像时不能破坏调频步进信号模糊函数所限定的信号内部结构关系。

步骤S302，对零插值后的M个长度为N_zp的向量分别在慢时间维上做IFFT处理，得到M个距离长度为δR_c＝c/2Δf的子距离像，子距离像中心位于R_k且其最小分辨单元为δr＝c/(2N_zpΔf)。为了保证相邻子距离像间不出现距离空隙，参数设计时必须保证δR_c＞δR_g，即T_g＜1/Δf。IFFT处理后的表达式为：

(18)

经过以上两个步骤，便完成了调频步进信号慢时间维合成子距离像处理。

步骤S103，对子距离像进行合成全程一维距离像处理。

下面结合图7对该步骤的具体实现流程进行描述。

步骤S401，舍弃子距离像中距离段[-δR_g/2+R_k，δR_g/2+R_k]以外部分，即仅子距离像高分辨谱线中以R_k为中心位置的

个点，从而实现以目标实际位置作为中心进行抽取，其中

为单采样点所代表距离对应的高分辨谱线个数。最终，每个子距离像中所抽取的高分辨谱线位于R_k+pδr处，

使得不同子距离像的抽取结果不会重叠，从而保证本发明合成高分辨像时不能破坏调频步进信号模糊函数所限定的信号内部结构关系。

抽取后表达式为

${N^{\left(c\right)}}_{\mathrm{zp}}=T_g{N}_{\mathrm{zp}}\mathrm{\Δf}$

步骤S402，将不同子距离像的抽取结果顺次拼接，并进行取模操作便得到了最终的全程一维距离像，其表达式如下：

$p_r\left(\frac{l}{\mathrm{\Δf}{N}_{\mathrm{zp}}}\right)=\underset{k=0}{\overset{M-1}{\mathrm{\Σ}}}\left|x_{\mathrm{UMF}}({\mathrm{kT}}_g-{\mathrm{\τ}}_0+{\mathrm{\τ}}_m,f_d)\right|$

$\×\left|\frac{\sin \left(\mathrm{\πN\Δf}(\frac{l-{\mathrm{kN}}_{\mathrm{zp}}^{\left(c\right)}}{\mathrm{\Δ}{\mathrm{fN}}_{\mathrm{zp}}}+({\mathrm{kT}}_g+{\mathrm{\τ}}_m-{\mathrm{\τ}}_0)+\frac{f_d{T}_R}{\mathrm{\Δf}})\right)}{\sin \left(\mathrm{\π\Δf}(\frac{l-N_{\mathrm{zp}}^{\left(c\right)}}{\mathrm{\Δf}{N}_{\mathrm{zp}}}+({\mathrm{kT}}_g+{\mathrm{\τ}}_m-{\mathrm{\τ}}_0)+\frac{f_d{T}_R}{\mathrm{\Δf}})\right)}\right|\×{W^{\left(c\right)}}_{N_{\mathrm{zp}}}(l-{\mathrm{kN}}_{\mathrm{zp}}^{\left(c\right)})---\left(20\right)$

$\frac{N_{\mathrm{zp}}^{\left(c\right)}}{2}\≤l\≤\frac{N_{\mathrm{zp}}^{\left(c\right)}}{2}-1+(M-1)N_{\mathrm{zp}}^{\left(c\right)}$ $N_{\mathrm{zp}}^{\left(c\right)}=T_g{N}_{\mathrm{zp}}\mathrm{\Δf}$

经过以上两个步骤，便完成了调频步进信号合成全程一维距离像处理，从而本发明的整个处理流程结束。

通过以上三个大步骤对线性调频步进信号视频回波的合成高分辨处理就可以得到超低时间栅瓣的合成高分辨像，此时频率步进信号的时间栅瓣对后续信号处理的影响基本可以忽略。这便可以有效消除频率步进信号时间栅瓣带来的测距模糊、目标遮挡，大大提升频率步进信号的成像及检测性能，进而提升后续目标跟踪性能及识别效果。

附图8为上述实施例的一个静止点目标实测数据高分辨像栅瓣抑制结果图，左图为常规合成高分辨算法成像结果，右图为本发明合成高分辨算法成像结果，其***参数如下：20120320

表1***仿真参数

以上所述的具体描述，对发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (1)

1.一种调频步进信号合成高分辨像栅瓣抑制方法，其特征在于，其具体包括如下步骤： 

1）对采样离散化后的调频步进信号的回波信号在快时间维进行加权失配滤波处理；所述加权失配滤波处理所采用的加权函数满足如下条件： 

其中，Δf为调频步进信号的步进间隔，B₁为调频步进信号的子脉冲带宽，K为所述加权函数相对于矩形窗的主瓣展宽系数； 

2）对加权失配滤波处理后结果在慢时间维进行合成子距离像处理，具体包括： 

步骤S301：对加权失配滤波后结果在慢时间维进行零插值处理；在参数设计时保证插值后点数N_zp与调频步进信号的步进间隔Δf二者乘积为***采样率的整数倍； 

步骤S302，对零插值后的信号在慢时间维上做IFFT处理，得到M个距离长度为δR_c＝c/2Δf的子距离像，子距离像中心位于R_k；参数设计时满足δR_c＞δR_g的条件； 

其中，δR_g＝cT_g/2为进行快时间维处理时的粗分辨单元，T_g为调频步进信号的采样间隔，c为光速；M为调频步进信号的采样点数；R_k＝R_m+k·δR_g为成像距离窗口内第k个采样点对应的距离，R_m为成像距离窗口起始位置； 

3）对子距离像进行合成全程一维距离像处理，具体包括： 

步骤S401：舍弃子距离像中距离段[-δR_g/2+R_k,δR_g/2+R_k]以外部分，仅抽取子距离像高分辨谱线中以R_k为中心位置的

个点，其中为单采样点所代表距离对应的高分辨谱线个数； 

步骤S402：将不同子距离像的抽取结果顺次拼接并进行取模操作，得到全程一维距离像。 

Images