CN101295019A - 一种线性调频连续波合成孔径雷达视频信号的仿真方法 - Google Patents

Abstract

本发明提供了一种线性调频连续波合成孔径雷达视频信号的仿真方法，包括以下步骤：设置并读取仿真的合成孔径雷达系统参数信息；设置目标的信息参数；仿真时间的离散化处理；分别计算雷达平台与目标的瞬时位置，根据几何位置关系得到相位中心与目标的距离矢量；结合合成孔径雷达回波信号数学模型与雷达系统参数，得到一维的回波信号；一维信号按脉冲重复周期划分，格式成二维信号，存储完全部数据，得到了一个Na×Nr的二维信号，完成线性调频连续波合成孔径雷达视频信号的精细仿真。利用本发明可得到更准确的回波仿真信号，特别适用高速运动、远距离、雷达处于大斜视观测条件下的调频连续波信号的仿真方法，实现效率更优，具有更强的适应性和兼容性。

Description

一种线性调频连续波合成孔径雷达视频信号的仿真方法

技术领域

本发明属于雷达技术领域，特别涉及一种线性调频连续波(LFMCW)合成孔径雷达(SAR)视频回波信号的仿真方法。

背景技术

合成孔径雷达(SAR，Synthetic Aperture Radar)是一种高分辨率成像雷达，SAR主动发射电磁波(从微波波段到毫米波波段)，并接收目标反射的回波信号通过复杂的信号处理过程实现成像，具有全天时、全天候工作的能力。线性调频是合成孔径雷达常用的信号波形，传统的SAR系统采用脉冲多普勒体制的线性调频信号，其回波模型的仿真与成像算法相对比较成熟。

线性调频连续波SAR的概念最先由英国伦敦大学提出并首先应用与飞机高度计中。线性调频连续波SAR相对于传统脉冲多普勒体制SAR具有结构简单、体积小、重量轻、能耗小、成本低、截获率低、抗干扰能力强等特点。

LFMCW SAR具有上述的诸多优点，因此受到国内外越来越多的重视。对这一新体制雷达的研究，在研究初期，对成像方式的原理及相应处理算法的研究需要符合条件的SAR回波数据，应用合成孔径雷达的模拟技术可以降低雷达开发成本，缩短了开发的周期，使雷达各个单元之间的设计和测试可以同时进行。特别是对于高速运动平台，雷达在超大前斜视条件下成像的情况，由于客观条件的限制，国内外没有相关的真实数据，此时雷达信号的模拟可以在实际雷达系统前端不具备的条件下对雷达系统后期信号处理部分进行调试和测试，为雷达系统的设计、分析和性能测试提供了有效的工具。

雷达信号的仿真主要分为功能仿真和视频信号仿真，功能仿真通常是通过对目标和载体运动方程的计算，预测出目标会在某个时刻某个位置被检测到，得到该时刻目标的距离和方位数据；目标视频信号仿真包含了信号的相位信息，比功能仿真更复杂，逼真度更高。

传统体制的脉冲多普勒雷达采用线性调频脉冲信号作为发射信号，由于脉冲信号的持续时间非常短，占空比一般低于10％，可以忽略雷达从发射脉冲信号到接收到目标的回波信号这段时间内雷达平台的运动，即认为雷达在发射与接收信号时，平台是静止不动的，这就是常用的“STOP AND GO”模型的假设。该模型对于脉冲体制雷达具有一定的合理性，但线性调频连续波体制雷达有其特殊性。一方面，采用线性调频连续波技术的新体制雷达不同于脉冲体制雷达在工作时会反复的开关，它始终在不断发射信号，信号的持续时间占据了整个脉冲重复周期，而且线性调频连续波的重复周期比较长，必须考虑脉冲周期内雷达平台的运动；另一方面，为避免回波信号与发射信号的混叠，需要将发射天线与接收天线分开，使雷达系统在发射的同时接收，而信号的发射与接收又是充分隔离的。因此传统脉冲体制雷达的“STOPAND GO”模型不再适用，需要建立一种新的适用与线性调频连续波合成孔径雷达的回波仿真方法。

目前对于LFMCW SAR的回波模拟方法基本上可分为两类，一种是延用“STOPANDGO”模型，但之后会对雷达平台连续运动带来的影响进行补偿，另一种是考虑了雷达平台在脉冲重复周期内的运动，但在处理过程中，经过大量近似达到化简的目的。这两种方法对于慢速平台对近距离观测(无人机与小型机平台)的情况是可行的，但应用于高速运动平台，进行大前斜视远距观测时，会出现较大的误差；虽然在后续处理中可以进行补偿，增加了处理过程的复杂度，也限制了仿真方法的通用性。

发明内容

本发明提供了一种线性调频连续波合成孔径雷达视频信号的仿真方法，它具有精确、通用性强的特点，可以方便的得到各种参数条件下雷达回波信号的仿真数据。特别是在高速飞行、超大前斜视条件下，本发明仍能较快速得到准确的仿真数据。

本发明提供了一种高速运动、远距离、雷达处于大斜视观测条件下的调频连续波信号的仿真方法，它包括以下几个实现步骤：

(1)设置并读取仿真的SAR系统参数信息，包括仿真信号、雷达运动平台和雷达波束天线的参数信息等。信号参数选择为：信号所处的波段为Ka波段，脉冲带宽B_w为百兆量级，脉冲重复周期在毫秒的数量级。雷达平台做高速运动，即速度V＞1000m/s，雷达作用距离较远，即雷达与目标的距离为1-10km。雷达处于大斜视观测条件即是选取雷达天线波束中心与航向夹角为15°～30°。

(2)设置目标的信息参数，包括目标的个数、运动状态、位置等。多点目标可以看成是多个单点目标仿真信号的叠加。

(3)仿真时间的离散化处理。将模拟信号转换成离散的数字信号是通过采样来实现的，对模拟信号进行采样的频率即为信号的采样率。将整个仿真时间划分成若干等长的时间段，每段的时间长度即采样时间的确定决定着仿真的精度，分得越细越精确，但运算量也会随之增加。鉴于调频连续波的特殊性，对信号的采样可以分为两种：快采样和慢采样。以脉冲重复周期为周期对信号进行采样称为慢采样，即方位向上的采样率，数值上与脉冲重复频率f_p相等；对同一个周期内的信号进行采样称为快采样，即距离向上的采样率，用f_s表示。根据抽样定理，采样率要大于信号带宽，才能从抽样信号中恢复原连续信号。另一方面，选取的采样率越大，整个仿真过程中的数据率与计算量也会随之增加。SAR回波信号多为复数信号，因此采样率约为信号带宽的1.1倍就能满足信号不失真的条件。信号先经过慢采样，设得到Na个长度为T的时间段，再对每个时间段进行快采样，则每个时间段T又被分为Nr＝f_s·T个时间段，即整个仿真时间离散分为N＝Na·Nr个长度为的时间段，则瞬时时刻t_i可以表示为： $t_i=-\frac{N}{2f_s}+\frac{i}{f_s},$ i＝1、2、3…N。

(4)分别计算雷达平台与目标的瞬时位置，根据其几何位置关系得到相位中心与目标的距离矢量即雷达天线中心与目标连线的长度和视线夹角。设平台的初始位置为P_m0＝[X_m0 Y_m0 Z_m0]^T，目标的初始位置为P_T0＝[X_T0 Y_T0 Z_T0]^T，平台速度和它在各个方向的分量分别为V_T＝[V_mX V_mY V_mZ]^T，目标速度和它在各个方向的分量分别为V_T＝[V_TX V_TY V_TZ]^T，设t_i时刻为与快采样时间有关的瞬时时刻，即整个仿真时间被离散化为N个时长为Δt的时间段， $t_{i+1}-t_i=\mathrm{\Δt}=\frac{1}{f_s},$ 则t_i，i＝1、2……N时刻平台的位置为

$P_{\mathrm{mi}}=P_{m0}+V_m{t}_i=\left[\begin{array}{c}{X}_{m0}+V_{\mathrm{mX}}{t}_i\\ Y_{m0}+V_{\mathrm{mY}}{t}_i\\ Z_{m0}+V_{\mathrm{mZ}}{t}_i\end{array}\right]=\left[\begin{array}{c}{x}_{\mathrm{mi}}\\ y_{\mathrm{mi}}\\ z_{\mathrm{mi}}\end{array}\right]$ t_i，i＝1、2…N时刻目标的位置为

$P_{\mathrm{Ti}}=P_{T0}+V_T{t}_i=\left[\begin{array}{c}{X}_{T0}+V_{\mathrm{TX}}{t}_i\\ Y_{T0}+V_{\mathrm{TY}}{t}_i\\ Z_{T0}+V_{\mathrm{TZ}}{t}_i\end{array}\right]=\left[\begin{array}{c}{X}_{\mathrm{Ti}}\\ Y_{\mathrm{Ti}}\\ Z_{\mathrm{Ti}}\end{array}\right]$ P_mi、P_Ti为雷达平台与目标的瞬时位置，t_i

时刻等效斜距为 $r_i=\sqrt{{(x_{\mathrm{mi}}-X_{\mathrm{Ti}})}^2+{(y_{\mathrm{mi}}-Y_{\mathrm{Ti}})}^2+{(z_{\mathrm{mi}}-Z_{\mathrm{Ti}})}^2},$

$=\sqrt{r_0^2-{2r}_0{\mathrm{Vt}}_i\cos \mathrm{\φ}+{\left({\mathrm{Vt}}_i\right)}^2}$

其中r₀为t_i＝0时刻雷达平台与目标之间的距离，V为雷达平台相对于目标的相对运动速率。φ为斜视角，即r₀与V之间的夹角。斜距的计算不再是传统的雷达仿真方法中仅随慢采样时间即方位向时间变化，而是在采用了“瞬时停走”的假设以后，在每个快采样时间都会有相应变化。因此，本发明所提供的雷达信号仿真方法比现有的雷达信号仿真方法更加精确。

(5)结合合成孔径雷达回波信号数学模型与步骤(1)中的雷达系统参数，得到一维的回波信号。步骤(3)时间的离散化中对采样率的选择可以看出，信号的带宽越大，采样率也会随之增大。雷达信号的带宽多在百兆的量级，步骤(4)中对雷达平台与目标位置及等效斜距的计算都直接与采样率相关，若直接采用传统的雷达信号模型，会带来很大的运算量，也影响了整个仿真过程的速度，因此，本发明采用经过去调频(dechirping)处理的回波信号模型，下面详细说明去调频处理的过程。去调频是用一个时间固定，而频率、调频率相同的线性调频信号作为参考信号，用它与回波做差频处理，设接收信号为 $s_r=\mathrm{\α}(t_q-\frac{{2r}_i}{c})\exp \left\{j2\mathrm{\π}(f_c(t-\frac{{2r}_i}{c})+\frac{1}{2}b{(t_q-\frac{{2r}_i}{c})}^2)\right\},$ 其中b为发射线性调频信号的调频率；α(t)为发射脉冲包络，r_i为点目标到雷达的斜距瞬时变化，c为光速，f_c为信号载频，t_q为距离向上的时间， $t_q=\mathrm{mod}(t,T)-\frac{T}{2},$ mod(t，T)为t除以T的余数，，接收信号的延时时间 $\mathrm{\τ}=\frac{{2r}_i}{c};$ 设参考距离为r_ref，则参考信号为： $s_{\mathrm{ref}}=\mathrm{\α}(t_q-\frac{{2r}_{\mathrm{ref}}}{c})\exp \left\{j2\mathrm{\π}(f_c(t-\frac{{2r}_{\mathrm{ref}}}{c})+\frac{1}{2}b{(t_q-\frac{{2r}_{\mathrm{ref}}}{c})}^2)\right\},$ 混频信号 $s_{\mathrm{if}}=s_r\·s_{\mathrm{ref}}^{*},$ 上述处理为去调频过程，s_if为去调频处理后的信号。经过去调频处理的雷达信号数学模型表达式为：

$s_{\mathrm{if}}=(t,t_q,r)=\mathrm{\σ}{W}_a\left(t\right)\mathrm{\α}(t_q-\frac{{2r}_i}{c})\mathrm{\α}(t_q-\frac{{2r}_{\mathrm{ref}}}{c})\exp \{-\mathrm{jb}\frac{4\mathrm{\π}}{c}{t}_q(r_i-R_{\mathrm{ref}})\}$

                                                     ，σ为点目标后向散

$\exp \left\{\mathrm{jb}\frac{4\mathrm{\π}}{c^2}(r_i^2-R_{\mathrm{ref}}^2)\right\}\exp \{-j\frac{4\mathrm{\π}}{\mathrm{\λ}}(r_i-R_{\mathrm{ref}})\}$

射特性，W_a(t)为方位向天线方向性函数，α(t)为发射脉冲包络，t为瞬时时刻，上式为本发明雷达回波信号数学模型，将步骤(4)中求出的瞬时斜距与步骤(1)中的雷达系统参数代入，即可得到随瞬时时间变化的长度为N的一维回波信号，将信号依次存入储存器中。本发明还可以加入误差模型，注入各种误差，即s＝s_if+ψ，ψ为误差表达式，可以得到存在误差的雷达回波信号的仿真数据。

(6)一维信号按脉冲重复周期划分，格式成二维信号。雷达图像处理需要二维的回波信号，将步骤(5)中得到的长度为N的一维信号按脉冲重复周期进行划分，设每个周期内的离散信号为Nr个，有Na个仿真周期，即N＝Nr·Na，将每一个重复周期时间内的信号叠放于上一个重复周期的数据后，得到了一个Na×Nr的二维信号，完成线性调频连续波合成孔径雷达视频信号的仿真。

本发明提供的线性调频连续波合成孔径雷达视频信号的精细仿真方法的优点：

(1)本发明提供的方法解决了雷达仿真处理在高速平台、大前斜视条件下使用“STOPAND GO”模型所引入的误差问题，使得雷达每个快采样时间内(瞬时时刻)的斜距信息更精确，进而得到更准确的回波仿真信号。

(2)本发明提供的方法吸收了对线性调频信号进行去调频的方式处理的方法，避免了雷达发射信号与接收信号的混叠，同时可以减小信号的带宽，降低了采样率，提高了处理效率。

(3)本发明提供的方法中的斜距计算，直接应用了点对点的空间几何关系来确定斜距，是一种比较准确的雷达距离模型，适用于正侧视、大斜视等情况下的斜距计算。

(4)具有较好的灵活性与可控性。可以根据应用场合的不同，将仿真时间按需进行划分，从而在满足应用需求的条件下，实现效率最优。

(5)强大的适应性，该平台可通过设置不同的参数完成机载、弹载等雷达平台的回波仿真运算。

(6)强大的兼容性。通过改变回波信号的数学模型，可以完成LFMCW信号在不同调制形式下的仿真如相位编码等技术。

(7)强大完善性。留有误差接口，可根据应用场合的不同注入各种误差。

附图说明

图1是本发明LFMCW雷达视频信号的仿真方法流程图；

图2是本发明中关于斜距计算时雷达平台与目标空间位置的示意图；

图3是本发明中将一维信号格式化为二维信号的示意图。

具体实施方式

下面结合附图对本发明进行详细说明。

本发明提供的一种线性调频连续波合成孔径雷达视频信号的仿真方法是通过以下步骤来实现的，如图1所示：

步骤一、读取仿真的SAR系统参数信息，包括仿真信号、雷达运动平台和雷达波束天线的参数信息等。

(a)仿真信号的参数：信号的波段(载频和波长λ)、脉冲带宽B_w、脉冲重复频率f_p(脉冲扫频周期T＝1/f_p)等，这类参数由仿真的具体需求决定，本发明可以很灵活的进行参数设置。根据调频连续波信号(FMCW)的特殊性，本发明信号参数选择为：信号所处的波段为Ka波段，脉冲带宽B_w为百兆量级，脉冲重复周期在毫秒的数量级。

(b)雷达运动平台参数：平台的初始位置、平台的运动状态等，这类参数决定了雷达运动的航迹。本发明可针对高速运动的雷达平台，即速度V＞1000m/s，雷达作用距离较远，即雷达与目标的距离为1-10km。

(c)雷达波束天线的参数：天线位置、天线尺寸、波束的中心视角等，天线尺寸决定了天线的观测范围，本发明可以对雷达在大斜视条件下进行观测的情况进行信号仿真，大斜视条件即天线中心指向与雷达平台航线夹角很小的情况，如夹角仅为15°～30°。

步骤二、设置目标的信息参数，包括目标的个数、运动状态、位置等。多点目标可以看成是多个单点目标仿真信号的叠加。

步骤三、仿真时间的离散化处理，即将整个仿真时间划分成若干等长的时间段，对整个仿真时间进行快采样。将模拟信号转换成离散的数字信号通过采样来实现，对模拟信号进行采样的频率即为信号的采样率。鉴于调频连续波的特殊性，对信号的采样可以分为两种：快采样和慢采样。以脉冲重复周期为周期进行采样称为慢采样，即方位向上的采样率，数值上与f_p相等；对同一个周期内的信号进行采样称为快采样，即距离向上的采样率，用f_s表示，根据抽样定理，采样率要大于信号带宽，才能从抽样信号中恢复原连续信号。另一方面，选取的采样率越大，整个仿真过程中的数据率与计算量也会随之增加，因此需要权衡准确性与计算量来选择采样率。SAR回波信号多为复数信号，因此采样率约为信号带宽的1.1倍就能满足信号不失真的条件。由于本发明是仿真的经过去调频后的回波信号，信号的带宽比实际信号的带宽窄得多，一般不会超过MHZ的量级。

设脉冲重复周期T，仿真时间为Na个脉冲重复周期，即仿真时间为T·Na，慢采样 $f_p=\frac{1}{T},$ 可看做将整个仿真时间按脉冲重复周期进行划分，快采样速率为f_s，即前面文字中的采样率，则整个仿真时间分为N＝Na·T·f_s个长度为

的时间段。则瞬时时刻t_i可以表示为： $t_i=-\frac{N}{2f_s}+\frac{i}{f_s},$ i＝1、2、3…N。

步骤四、分别计算雷达平台与目标的瞬时位置，根据几何位置关系得到相位中心与目标的距离矢量即雷达天线中心与目标连线的长度和视线夹角。步骤一和步骤二已经提供了雷达平台与目标的初始位置和运动参数，只需要将步骤三中离散化后的时间t_i代入P_mi＝P_m0+V_mt_i，P_Ti＝P_T0+V_Tt_i，式中P_m0、P_T0分别为雷达平台与目标的初始位置坐标，V_m、V_T分别为雷达平台与目标的运动速度，t_i为瞬时时刻，就可以得到相应时间雷达平台与目标的位置信息P_mi、P_Ti，由空间几何关系，可以得到雷达的斜距信息，并将斜距的值r_i顺序存入计算机存储器中。

下面结合附图2进行详细说明。设P_Ti＝1、2…N，P_mi，i＝1，2…N分别为目标和雷达平台在时刻t_i时刻的位置，平台飞行速度为V，则有 $P_{\mathrm{mi}}=P_{m0}+V_m{t}_i=\left[\begin{array}{c}{X}_{m0}+V_{\mathrm{mX}}{t}_i\\ Y_{m0}+V_{\mathrm{mY}}{t}_i\\ Z_{m0}+V_{\mathrm{mZ}}{t}_i\end{array}\right]=\left[\begin{array}{c}{x}_{\mathrm{mi}}\\ y_{\mathrm{mi}}\\ z_{\mathrm{mi}}\end{array}\right],$ $P_{\mathrm{Ti}}=P_{T0}+V_T{t}_i=\left[\begin{array}{c}{X}_{T0}+V_{\mathrm{TX}}{t}_i\\ Y_{T0}+V_{\mathrm{TY}}{t}_i\\ Z_{T0}+V_{\mathrm{TZ}}{t}_i\end{array}\right]=\left[\begin{array}{c}{X}_{\mathrm{Ti}}\\ Y_{\mathrm{Ti}}\\ Z_{\mathrm{Ti}}\end{array}\right]$

斜距

$r_i=\sqrt{{(x_{\mathrm{mi}}-X_{\mathrm{Ti}})}^2+{(y_{\mathrm{mi}}-Y_{\mathrm{Ti}})}^2+{(z_{\mathrm{mi}}-Z_{\mathrm{Ti}})}^2}.$

$=\sqrt{r_0^2-{2r}_0{\mathrm{Vt}}_i\cos \mathrm{\φ}+{\left({\mathrm{Vt}}_i\right)}^2}$

步骤五、结合合成孔径雷达回波信号数学模型与步骤一中的雷达系统参数，得到一维的回波信号。雷达平台沿其航线飞行，平台与目标间的距离随着时间推进也在不断变化。这种距离上的变化，在雷达信号中则反映的是信号的延时的不同，反映了雷达距离向的变化。

前面的步骤可以看到，每个快采样时间都记录下斜距值会带来很大的运算量，步骤三中已经说明了采样率应该如何选择。本发明中采用了去调频后的雷达视频信号，直接模拟经过去调频处理后的信号。它包括了信号的所有信息，同时也大大降低了信号快采样的采样率，减少了计算机的运算量。

(a)雷达延时的计算雷达发射的信号是电磁波，速度为光速c，发射的信号到达目标后再次被接收，这段时间可以表示为τ＝(r_t+r_r)/c，r_t为发射信号时刻的斜距，r_r为接收信号时刻的斜距。本发明采用“瞬时停走”的假设，认为信号在发射与接收时刻的斜距是相等的，即每个快采样的时间间隔内方位向位置的变化引起的斜距变化可以忽略，此时信号延时表达式简化为τ＝2r_i/c，r_i为快采样时刻的斜距。

(b)经过去调频处理的雷达信号数学模型表达式为：

$s_{\mathrm{if}}(t,t_q,r)=\mathrm{\σ}{W}_a\left(t\right)\mathrm{\α}(t_q-\frac{{2r}_i}{c})\mathrm{\α}(t_q-\frac{{2r}_{\mathrm{ref}}}{c})\exp \{-\mathrm{jb}\frac{4\mathrm{\π}}{c}{t}_q(r_i-R_{\mathrm{ref}})\}$

                                 ，σ为点目标后向散射特

$\exp \left\{\mathrm{jb}\frac{4\mathrm{\π}}{c^2}(r_i^2-R_{\mathrm{ref}}^2)\right\}\exp \{-j\frac{4\mathrm{\π}}{\mathrm{\λ}}(r_i-R_{\mathrm{ref}})\}$

性，W_a(t)为方位向天线方向性函数，α(t)为发射脉冲包络，b为发射线性调频信号的调频率，r_i为点目标到雷达的斜距瞬时变化，c为光速，λ为信号波长，t_q为距离向上的时间， $t_q=\mathrm{mod}(t,T)-\frac{T}{2},$ mod(t，T)为t除以T的余数，R_ref为去调频处理时选取的参考斜距，t为瞬时时刻。此处斜距r_i不再是传统的雷达仿真方法中仅随慢采样时间即方位向时间变化，而是在采用了“瞬时停走”的假设以后，在每个快采样时间都会有相应变化。因此，本发明所提供的雷达信号仿真方法比现有的雷达信号仿真方法更加精确。

(c)将步骤四中得到的斜距信息r_i代入上面的经过去调频处理的雷达信号数学模型表达式，结合与步骤一中雷达参数信息，得到一维的雷达信号，将其存入寄存器。

(d)初始化存储器，将一维寄存器中的回波信号与存储器中数据进行叠加。

(e)判断是否完成全部点目标的仿真，如果判断为“否”，则重复步骤二至步骤五的处理过程，得到一维的回波信号；如果判断为“是”，则转入步骤六。

步骤六、一维信号按脉冲重复周期划分，格式成二维信号。

下面结合附图3进行详细说明。将步骤五中得到的随快采样变化的信号按时间的先后顺次存入储存器，得到长度为N的一维信号，Na为仿真的脉冲周期个数，即为方位向的采样数；Nr为对一个周期内信号进行快采样后的信号数，即为距离向的采样点数，Nr＝T·f_s。一维信号按照脉冲周期(慢采样时间)进行划分，每个脉冲周期内有Nr段数据，从起始时刻起，前Nr个数据存为一行，第Nr+1个数据存到第二行第一个位置，第Nr+2个数据存到第二行第二个位置，第2Nr个数据存到第二行第Nr个位置，第2Nr+1个数据存到第三行第一个位置，以此类推……得到一个Na×Nr的二维信号完成了线性调频连续波合成孔径雷达视频信号的仿真。

本发明主要针对线性调频连续波信号的特点，在“瞬时停走”模型的假设的基础上，提出了一种适用于高速运动平台，在大斜视条件下的线性调频连续波合成孔径雷达视频信号的精细仿真方法。仿真过程中逐个脉冲记下瞬时时刻雷达斜距变化，能够较准确地实现信号的仿真，同时该方法中选用的是与参考信号进行去调频处理的回波信号，能在很大程度上解决由于记录瞬时斜距带来的计算量大的问题。

Claims (3)

1、一种线性调频连续波合成孔径雷达视频信号的仿真方法，其特征在于：它包括以下几个实现步骤，

步骤一：设置并读取仿真的合成孔径雷达系统参数信息，包括仿真信号、雷达运动平台和雷达波束天线的参数信息；其中信号所处的波段为Ka波段，脉冲带宽B_w为百兆量级，脉冲重复周期在毫秒的数量级，雷达平台做高速运动，即速度V＞1000m/s，雷达作用距离即雷达与目标的距离为1-10km，雷达处于大斜视观测条件即是选取雷达天线波束中心与航向夹角为15°～30°；

步骤二：设置目标的信息参数，包括目标的个数、运动状态、位置，多点目标看作是多个单点目标仿真信号的叠加；

步骤三：仿真时间的离散化处理，对整个仿真时间进行快采样，将整个仿真时间划分成若干等长的时间段；

本步骤是将模拟信号转换成离散的数字信号，对模拟信号进行采样的频率即为信号的采样率，设脉冲重复周期T，仿真时间为Na个脉冲重复周期，即仿真时间为T·Na，慢采样 $f_p=\frac{1}{T},$ 将整个仿真时间按脉冲重复周期进行划分；快采样速率为f_s，即前面文字中的采样率，则整个仿真时间分为N＝Na·T·f_s个长度为

的时间段，则瞬时时刻t_i表示为： $t_i=-\frac{N}{2f_s}+\frac{i}{f_s},i=\mathrm{1,2}\·\·\·N;$

步骤四：分别计算雷达平台与目标的瞬时位置，根据几何位置关系得到相位中心与目标的距离矢量即雷达天线中心与目标连线的长度和视线夹角；步骤一和步骤二已经提供了雷达平台与目标的初始位置和运动参数，只需要将步骤三中离散化后的时间t_i代入P_mi＝P_m0+V_mt_i，P_Ti＝P_T0+V_Tt_i，式中P_m0、P_T0分别为雷达平台与目标的初始位置坐标，V_m、V_T分别为雷达平台与目标的运动速度，t_i为瞬时时刻，就得到相应时间雷达平台与目标的位置信息P_mi、P_Ti，由空间几何关系，得到雷达的斜距信息，并将斜距的值r_i顺序存入计算机存储器中；设P_Ti，i＝1、2…N P_mi，i＝1，2…N分别为目标和雷达平台在时刻t_i时刻的位置，平台飞行速度为V，则有

$P_{\mathrm{mi}}=P_{m0}+V_m{t}_i=\left[\begin{array}{c}{X}_{m0}+V_{\mathrm{mX}}{t}_i\\ Y_{m0}+V_{\mathrm{mY}}{t}_i\\ Z_{m0}+V_{\mathrm{mZ}}{t}_i\end{array}\right]=\left[\begin{array}{c}{x}_{\mathrm{mi}}\\ y_{\mathrm{mi}}\\ z_{\mathrm{mi}}\end{array}\right]$

$P_{\mathrm{Ti}}=P_{T0}+V_T{t}_i=\left[\begin{array}{c}{X}_{T0}+V_{\mathrm{TX}}{t}_i\\ Y_{T0}+V_{\mathrm{TY}}{t}_i\\ Z_{T0}+V_{\mathrm{TZ}}{t}_i\end{array}\right]=\left[\begin{array}{c}{x}_{\mathrm{Ti}}\\ Y_{\mathrm{Ti}}\\ Z_{\mathrm{Ti}}\end{array}\right]$

斜距 $r_i=\sqrt{{(x_{\mathrm{mi}}-X_{\mathrm{Ti}})}^2+{(y_{\mathrm{mi}}-Y_{\mathrm{Ti}})}^2+{(z_{\mathrm{mi}}-Z_{\mathrm{Ti}})}^2}$

$=\sqrt{r_0^2-{2r}_0{\mathrm{Vt}}_i\cos \mathrm{\φ}+{\left({\mathrm{Vt}}_i\right)}^2};$

步骤五：结合合成孔径雷达回波信号数学模型与步骤一中的雷达系统参数，得到一维的回波信号；具体步骤如下：

(a)雷达延时的计算雷达发射的信号是电磁波，速度为光速c，假设信号在发射与接收时刻的斜距是相等的，即每个快采样的时间间隔内方位向位置的变化引起的斜距变化忽略不计，此时信号延时表达式简化为τ＝2r_i/c，r_i为快采样时刻的斜距；

(b)经过去调频处理的雷达信号数学模型表达式为：

$s_{\mathrm{if}}=(t,t_q,r)=\mathrm{\σ}{W}_a\left(t\right)\mathrm{\α}(t_q-\frac{{2r}_i}{c})\mathrm{\α}(t_q-\frac{{2r}_{\mathrm{ref}}}{c})\exp \{-\mathrm{jb}\frac{4\mathrm{\π}}{c}{t}_q(r_i-R_{\mathrm{ref}})\}$

$\exp \left\{\mathrm{jb}\frac{4\mathrm{\π}}{c^2}(r_i^2-R_{\mathrm{ref}}^2)\right\}\exp \{-j\frac{4\mathrm{\π}}{\mathrm{\λ}}(r_i-R_{\mathrm{ref}})\}.$

σ为点目标后向散射特性，W_a(t)为方位向天线方向性函数，α(t)为发射脉冲包络，b为发射线性调频信号的调频率，r_i为点目标到雷达的斜距瞬时变化，c为光速，λ为信号波长，t_q为距离向上的时间， $t_q=\mathrm{mod}(t,T)-\frac{T}{2},$ mod(t，T)为t除以T的余数，R_ref为去调频处理时选取的参考斜距，t为瞬时时刻；

(c)将步骤四中得到的斜距信息r_i代入上面的经过去调频处理的雷达信号数学模型表达式，结合与步骤一中雷达参数信息，得到一维的雷达信号，将其存入寄存器；

(d)初始化存储器，将一维寄存器中的回波信号与存储器中数据进行叠加；

(e)判断是否完成全部点目标的仿真，如果判断为“否”，则重复步骤二至步骤五的处理过程，得到一维的回波信号；如果判断为“是”，则转入步骤六；

步骤六：一维信号按脉冲重复周期划分，格式成二维信号；

将步骤五中得到的随快采样变化的信号按时间的先后顺次存入储存器，得到长度为N的一维信号，Na为仿真的脉冲周期个数，即为方位向的采样数；Nr为对一个周期内信号进行快采样后的信号数，即为距离向的采样点数，Nr＝T·f_s；一维信号按照脉冲周期即慢采样时间进行划分，每个脉冲周期内有Nr段数据，从起始时刻起，前Nr个数据存为一行，第Nr+1个数据存到第二行第一个位置，第Nr+2个数据存到第二行第二个位置，第2Nr个数据存到第二行第Nr个位置，第2Nr+1个数据存到第三行第一个位置，以此类推，得到一个Na×Nr的二维信号，完成线性调频连续波合成孔径雷达视频信号的仿真。

2、根据权利要求1所述的一种线性调频连续波合成孔径雷达视频信号的仿真方法，其特征在于，步骤三所述的慢采样是方位向上的采样率，以脉冲重复周期为周期进行采样，数值上与f_p相等。

3、根据权利要求1所述的一种线性调频连续波合成孔径雷达视频信号的仿真方法，其特征在于，步骤三所述的快采样是对同一个周期内的信号进行采样，即距离向上的采样率，用f_s表示。

Images