提高大场景目标成像质量的MIMO雷达成像方法与流程

本发明涉及雷达
技术领域：
，尤其涉及一种提高大场景目标成像质量的mimo雷达成像方法。
背景技术：
：雷达成像(radarimaging)技术是雷达发展史上的一个飞跃，其拓展了雷达最初的检测(detection)和测距(ranging)的功能，使得雷达能利用获取的电磁散射信息，获得对场景的全景雷达图像。合成孔径雷达(sar)虽然具有高的方位向分辨率，但其重访周期长，在需要长时间连续凝视观测和成像的场景下，存在不可避免的先天劣势。传统的实孔径雷达凝视成像，其角度分辨率受实际天线阵列孔径限制，限制了其在实际中的应用。近年来，雷达领域中引入了一种新体制雷达——多输入多输出(mimo,multiple-inputmultiple-output)雷达，可以实现远距离、高分辨的凝视成像。mimo雷达在发射端和接收端均采用多天线结构，各个发射天线同时辐射周期性、相互正交的信号波形，接收端的每个天线接收所有发射信号并在后端进行信号分选，从而得到了远多于实际收、发阵元数目的观测通道和自由度。空间并存的多观测通道使得mimo雷达能够实时采集携带有目标不同幅度、时延或相位的回波信息，这种并行、多通道获取信息的能力正是mimo雷达的根本优势所在。源于这种体制优势，mimo雷达与传统的单/多基地雷达或是相控阵雷达相比极大地提高了雷达的总体性能。将mimo雷达成像应用于大场景成像时，由于mimo雷达成像能力与观测通道数有关，一方面，对于某一观测通道数有限的mimo雷达系统，其成像能力也是有限的，而大场景具有大量的目标信息，目标信息量会超出雷达系统成像能力极限，会导致成像质量下降；另一方面，并且由于mimo雷达成像机理，在推导空间谱分布的时候，应用到目标的远场近似处理，导致目标距离成像中心越远时，反演结果与真实值偏差越大；因此，当mimo雷达系统对于大场景目标成像时，采用传统方法成像效果较差。技术实现要素：本发明的目的是提供一种提高大场景目标成像质量的mimo雷达成像方法，可以极大的提高大场景目标成像质量。本发明的目的是通过以下技术方案实现的：一种提高大场景目标成像质量的mimo雷达成像方法，包括：mimo成像雷达阵列的各个辐射天线单元多通道且同步发射周期性且相互正交的窄脉冲信号；利用窄脉冲信号的雷达距离门特征，将斜视方向共同辐射的二维成像区域在距离向上划分为一定宽度的多个独立成像条带；对于每一成像条带，mimo成像雷达阵列均形成以相应成像条带中心为成像中心的新mimo雷达构型，并确定新mimo雷达构型下的空间谱支撑域；对于各接收机所接收的散射回波信号进行回波距离门划分和回波子带划分处理，得到对应于每一成像条带的有效散射回波信号；对每一成像条带的有效散射回波信号分别进行处理，得到每一成像条带每个观测通道下的空间谱值；利用每一成像条带每个观测通道下的空间谱值，以及相应的空间谱支撑域，进行反演成像，获得每一成像条带的成像结果；将所有成像条带的成像结果进行拼接，从而得到整个场景的成像结果。由上述本发明提供的技术方案可以看出，采用mimo条带划分方法和独立条带成像方式，有效解决了大场景高分辨mimo雷达成像条件下需要较多的收发阵列单元和巨大的多通道数据处理压力，有利于提高大场景目标成像质量。附图说明为了更清楚地说明本发明实施例的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域的普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他附图。图1为本发明实施例提供的提高大场景目标成像质量的mimo雷达成像方法的流程图；图2为本发明实施例提供的提高大场景目标成像质量的mimo雷达成像方法的成像场景示意图；图3为本发明实施例提供的散射回波信号的距离门与回波子带划分示意图；图4为本发明实施例提供的仿真目标示意图；图5为本发明实施例提供的利用传统mimo雷达成像方法在未进行成像条带划分条件下的仿真成像结果；图6为本发明实施例提供的采用本发明提供的mimo雷达成像方法在划分4个成像条带条件下的仿真成像结果。具体实施方式下面结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明的保护范围。本发明实施例提供一种提高大场景目标成像质量的mimo雷达成像方法，该方法采用了mimo条带划分和独立条带成像方式，该方法的流程如图1所示，相关成像场景如图2所示。如图1所示，相关成像方法主要包括如下步骤：步骤11、利用窄脉冲信号的雷达距离门特征，将斜视方向共同辐射的二维成像区域在距离向上划分为一定宽度的多个独立成像条带。本发明实施例中，mimo成像雷达阵列的各个辐射天线单元多通道且同步发射周期性且相互正交的窄脉冲信号，所述mimo成像雷达阵列包括一定分布构型下的p个辐射天线单元与q个接收机。mimo成像雷达阵列中各辐射天线单元相对二维成像区域s的斜视角为θ，二维成像区域s为尺寸为a×a的方形区域，mimo成像雷达阵列与二维成像区域s的垂直高度为h；以整个场景的二维成像区域s为xoy平面，高度向为z轴建立合适的空间极坐标系o，得到第p个辐射天线单元的位置矢量(p＝1,2,...,p)，得到第q个接收机的位置矢量(q＝1,2,...,q)；其中的r,θ,分别对应模值、俯仰角和方位角。各个辐射天线单元多通道且同步发射周期性且相互正交的窄脉冲信号，可以为但不限于随机跳频窄脉冲信号，则第p个辐射天线单元发射的第l个脉冲信号可以表示为：其中，j表示虚数；t为时刻；窄脉冲信号为矩形包络，脉冲宽度为τ，脉冲周期为t，fp,l为第p个辐射天线单元发射的第l个脉冲信号的跳频频率，在辐射天线单元带宽b内随机选择，每次脉冲的所有fp,l组成正交跳频图案。区别于传统的mimo雷达成像方法，本发明实施例中首先利用窄脉冲雷达距离门特征，将斜视方向共同辐射的二维成像区域s划分多个成像条带，具体过程如下：窄脉冲信号在二维成像区域s距离向上产生距离门，将其在距离向切分为n条成像条带，则每个成像条带宽度为：其中，c为光速，θ为辐射天线单元相对于二维成像区域s的斜视角，成像条带总数目n表示为n＝a/δl，a为二维成像区域s距离向宽度；第n个成像条带到所述二维成像区域s近端距离范围为[(n-1)δl,nδl],n＝(1,2,…,n)。步骤12、对于每一成像条带，mimo成像雷达阵列均形成以相应成像条带中心为成像中心的新mimo雷达构型，并确定新mimo雷达构型下的空间谱支撑域。区别于传统方法，本发明实施例中在分割成像条带后，以每个成像条带中心作为成像中心分别计算空间谱支撑域，则可相对于传统方法得到更多的场景构型以及空间谱填充，具体过程如下：确定第n个成像条带sn的条带中心on的位置矢量为：其中，为二维成像区域s中心的位置矢量，假定成像条带总数目n为偶数，所有成像条带对称地分布于场景中心的两侧，为相邻条带中心间距离向外指向的位置矢量，矢量长度等于条带宽度δl；本发明实施例中，确定成像条带中心后，需计算收、发信机对应每个条带的相对位置坐标。以成像条带sn的条带中心on为新的成像中心和坐标原点，建立新的空间极坐标系on'；所述mimo成像雷达阵列相对于成像条带sn形成新mimo雷达构型，第p个辐射天线单元与第q个接收机相对于新的空间极坐标系on'的位置矢量分别为：其中，分别对应第p个辐射天线单元、第q个接收机的位置矢量；中的三个元素分别为位置矢量在所述新的空间极坐标系on'中的模值、俯仰角和方位角。对应于成像条带sn的新mimo雷达构型，第p个辐射天线单元与第q个接收机观测通道的二维空间谱支撑域表示为：步骤13、对于各接收机所接收的散射回波信号进行回波距离门划分和回波子带划分处理，得到对应于每一成像条带的有效散射回波信号。本步骤所涉及的回波距离门划分、回波子带划分及选择可参见图3，该步骤的具体过程如下：1)对每一个成像条带进行二维空间离散化，划分为m个网格，m＝i×v，i为方位向分辨单元数，v为距离向分辨单元数；ρa、ρr分别为方位向和距离向的网格间距；第n个成像条带sn第m个网格中心点的位置记为m＝(1,2,..,m)，其目标后向散射系数为2)将第l个窄脉冲信号的散射回波波至点(l-1)t+(τp,q,n,m)min在时间轴左移(l-1)t，实现所有l个脉冲回波的对齐，与每个划分的成像条带相对应，将对齐后的每个窄脉冲回波按照脉冲宽度τ划分为多个回波距离门；其中，为第p个辐射天线单元发射的信号经位于的目标网格单元散射，由第q接收机所接收信号的双程传播时延，(τp,q,n,m)min为其中最小时延；式中的分别对应第p个辐射天线单元、第q个接收机的位置矢量，c为光速；3)由于每个成像条带在距离向是可分辨的，含有v个宽度为ρr的成像网格，则在每个回波距离门内，按照宽度δ划分为多个回波子带，其中δ＝2ρrsinθ/c，θ为辐射天线单元相对于二维成像区域s的斜视角；4)选取第l个脉冲散射回波的第n个距离门内的最后一个回波子带，作为对应于第n个成像条带sn的有效散射回波信号，其仅包含成像条带sn内所有目标分辨单元散射信息，而不包含其他成像条带的目标分辨单元散射信息。步骤14、对每一成像条带的有效散射回波信号分别进行处理，得到每一成像条带每个观测通道下的空间谱值。本步骤所述的处理主要包括匹配滤波、分选、归一化、去载频等相关处理，具体过程如下：1)由第q个接收机所接收到的第n个成像条带sn的散射回波信号表示为：其中，rect为矩形函数，t为时刻，ap,q为考虑发、收双程空间自由传播的衰减因子，fp,l为第p个辐射天线单元发射的第l个脉冲信号的跳频频率，辐射天线单元总数为p；2)对散射回波信号sq,n(t)经过回波距离门与回波子带划分，再进行匹配滤波分选与归一化处理，得到第p个辐射天线单元发射、第q个接收机所接收的，第l个窄脉冲信号经过第n个成像条带sn目标区域散射的有效散射回波信号，即观测通道(p,q,l,n)的回波信号为：3)对回波信号sp,q,l,n(t)进行去载波处理，解调出其中的不变量，得到信号sp,q,n,l：其中，τp,n，τq,n分别为第p个辐射天线单元、第q个接收机到第n个成像条带sn的条带中心on的电磁波自由传播固定时延，分别有：其中，为第n个成像条带sn的条带中心on的位置矢量；4)经过对信号sp,q,n,l的整理与变换处理后，获得所述观测通道(p,q,l,n)的空间谱值，表示为：其中，(kx,p,q,l,n，ky,p,q,l,n)为对应于第n个成像条带sn的新mimo雷达构型下的空间谱支撑域。步骤15、利用每一成像条带每个观测通道下的空间谱值，以及相应的空间谱支撑域，进行反演成像，获得每一成像条带的成像结果。本发明实施例中，基于mimo雷达成像的空间谱理论，利用第n个成像条带sn的空间谱值gp,q,n,l(kx,p,q,l,n,ky,p,q,l,n)及相应的空间谱支撑域(kx,p,q,l,n，ky,p,q,l,n)进行第n个成像条带sn的反演成像，反演成像过程表示为：其中，ψ为反演算法，常见方法有ifft算法、sf-rm算法、ola-rm算法、nufft算法等，对于不同场景可以选取适合的算法，本领域技术人员能够理解基于空间谱理论的mimo雷达目标散射系数反演求算过程，在此不做赘述。步骤16、将所有成像条带的成像结果进行拼接，从而得到整个场景的成像结果。所有成像条带成像结果的拼接过程表示如下：其中，为第n个成像条带sn的成像结果，n＝(1,2,…,n)。本发明实施例上述方案，采用mimo条带划分方法和独立条带成像方式，有效解决了大场景高分辨mimo雷达成像条件下需要较多的收发阵列单元和巨大的多通道数据处理压力，有利于提高大场景目标成像质量。为了便于说明本发明上述方案的有效性，下面结合具体成像场景，通过附图和具体实施方式对本发明的原理和效果作详细描述。下面以一个较为典型的具体实例为例，给出本发明的具体实施方式，但是，不应以此限定本发明的实际应用以及保护范围。仿真条件：仿真参数参数设置辐射源阵列口径大小20m×20m辐射天线单元个数41接收机个数41二维成像区域s与辐射雷达阵列的竖直距离h＝75m雷达阵列对二维成像区域s的斜下视角θ＝45°二维成像区域s范围a×a16m×16m成像区域网格划分个数101×101＝10201发射信号形式正交跳频脉冲信号脉冲周期t1us系统带宽b＝1000mhz成像场景可以参见图2，仿真目标为位于地面停泊的飞机模型，对飞机模型进行了相应处理，只有限的离散强散射点，如图4所示。为了说明本发明所述方法的有效性，以下仿真利用雷达图像的积分旁瓣比(integralsideloberatio，islr)来评价成像质量。通常islr的值越小，图像的对比度将越低，图像中将充满大量的阴影区域，同时，过小的islr将导致图像中含有虚假目标。当发射信号的脉冲宽度取τ＝100ns的条件下，在成像区域不进行成像条带的划分，利用传统mimo雷达成像方法对整个场景进行反演成像，成像结果如图5所示；当发射信号的脉冲宽度取τ＝18.856ns的条件下，在利用窄脉冲雷达的距离门特征，把整个成像区域划分为n＝4个成像条带，利用本发明提出的mimo雷达成像方法对整个场景进行反演成像，成像结果如图6所示；同时，在两种条件下，分别计算反演成像的恢复图像积分旁瓣比如下表1所示。显然，与不划分条带的反演成像结果相比，划分成像条带后，mimo雷达的反演成像结果看起来更为清晰，其积分旁瓣比islr也提高了约1.2db，成像时间更短，该仿真结果说明本发明提出的mimo雷达成像方法的有效性。通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到上述实施例可以通过软件实现，也可以借助软件加必要的通用硬件平台的方式来实现。基于这样的理解，上述实施例的技术方案可以以软件产品的形式体现出来，该软件产品可以存储在一个非易失性存储介质(可以是cd-rom，u盘，移动硬盘等)中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述的方法。以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本
技术领域：
的技术人员在本发明披露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求书的保护范围为准。当前第1页12