本发明属于雷达技术领域，它特别涉及了合成孔径雷达(sar)成像技术领域。
背景技术：
作为一种工作在微波波段的有源雷达，合成孔径雷达(syntheticapertureradar,sar)具有全天时、全天候的成像能力，即无论是白天或黑夜、晴天还是雷雨风雪天气，都可以随时随地成像，克服了光学和红外系统不能在晚上和复杂天气条件进行成像的缺点。传统的sar成像一般只具有二维成像分辨率，在一些起伏比较大的地方比如陡峭的山峰、峡谷以及城市中矗立挺拔的高楼时，传统sar成像存在的失真(阴影遮挡效应、空间模糊、顶底倒置等)导致空间的一些重要信息(比如高度)丢失，所以能对目标进行三维成像是非常有必要的，为了适应这种需求，目前常见的三维成像技术有圆周sar(circularsar)三维成像、层析sar(tomographysar)三维成像、线阵sar(arraysar，asar)三维成像。
线阵sar(lasar)三维成像的基本原理是在切行迹向添加阵列天线，通过沿航迹向平台的飞行形成虚拟的面阵进而获得二维分辨率，距离向再通过脉冲压缩技术获得第三维的分辨率。相比于圆周sar三维成像，线阵sar三维成像不需要圆周运动的轨迹；相比于层析sar三维成像需要航过多次，线阵sar三维成像只需一次航过，所以线阵sar三维成像相对于层析sar和圆周sar三维成像有更强的灵活性。目前线阵sar三维成像技术在地形测绘、城市测绘、灾难救援、军事探测等领域发挥着重要的作用。
使用匹配滤波器算法进行lasar三维成像时，lasar的阵列成像分辨率受线性阵列的长度限制。为了实现高质量的成像，lasar必须满足以下两个要求:
(1)线性阵列中相邻元素之间的间距必须满足nyquist采样定理，以避免光栅波瓣。
(2)必须采用整个线性阵列的回波信号，尽量避免旁瓣干扰。
以上两个要求使得线性阵列中的元素数量非常庞大，在真实的硬件系统中实现lasar是非常复杂和昂贵的。
压缩感知(cs)算法可以利用随机采样信号恢复原始稀疏信号，根据lasar成像场景的稀疏性，已成功地将其引入lasar成像。采用不同cs算法进行lasar稀疏成像时，提高了lasar成像分辨率。为了减少用于稀疏成像的阵元，通常用随机采样阵代替线性阵元。
在进行稀疏成像之前，一旦预设参数设置不当，cs算法可能会受到严重的旁瓣干扰，甚至不能正确估计散射系数。它们的计算复杂度可能显著增加。在不同的lasar数据下，这个预设参数不能总是满足高质量稀疏成像的要求，需要反复调试。由于预设参数通常不需要在不同的lasar数据下重新调试，sbrim算法在cs算法中获得了较好的成像结果。然而，由于高维矩阵的运算(如矩阵的求逆)，其计算复杂度非常巨大，硬件系统必须有很强的计算能力才能完成高维矩阵的运算。这一缺点严重限制了sbrim算法在三维稀疏成像中的应用。因此，我们需要研究如何在保证成像质量的前提下提高sbrim算法的计算效率。
由于三维成像场景的稀疏性强，目标区域通常只占整个成像场景的一小部分，稀疏成像的测量矩阵通常包含大量与目标无关的元素。因此，通过对测量矩阵进行简化，可以有效地提高稀疏成像的计算效率。因此，我们需要研究一种新的快速cs算法用于lasar三维成像。为了提高基于cs算法的lasar成像的运算效率，本文发明提出了基于相关向量机的快速贝叶斯压缩感知(fbcs-rvm)算法。
技术实现要素：
为了提高基于cs算法的lasar成像的运算效率，本文发明提出了基于相关向量机的快速贝叶斯压缩感知(fbcs-rvm)方法，该方法在相关向量机(rvm)框架下计算最大似然边际函数，获得最优超参数，提取非零最优超参数对应的散射单元作为成像场景中的目标区域。然后，根据目标区域，简化测量矩阵并进行稀疏成像。此外，在低信噪比(snr)或高稀疏度的情况下，不能总是准确地提取目标区域，其中可能包含多个元素，与其他元素相比，这些元素的散射系数太小或更接近于0。使对角矩阵奇异不可逆，散射系数无法正确估计。为了解决这一问题，将奇异矩阵的逆矩阵替换为截断的奇异值分解(tsvd)算法得到的广义逆矩阵，以正确估计散射系数。根据奇异矩阵的秩，提取并消除目标区域中散射系数太小的元素，得到更精确的目标区域。
为了方便描述本发明的内容，首先作以下术语定义：
定义1、合成孔径雷达(sar)
合成孔径雷达是将雷达固定于载荷运动平台上，结合平台的运动以合成等效阵列以实现阵列向的分辨率，再利用雷达波束向回波延时实现距离一维成像，从而实现对观测目标二维成像的一种合成孔径雷达技术。
定义2、标准合成孔径雷达回波数据距离向脉冲压缩
标准合成孔径雷达回波数据距离向脉冲压缩是指利用合成孔径雷达发射信号参数，采用匹配滤波技术对合成孔径雷达的距离向信号进行信号聚焦成像的过程。详见文献“雷达成像技术”，保铮，邢孟道，王彤，电子工业出版社，2005。
定义3、范数
设x是复数域上线性空间，其中表示复数域，若它满足如下性质：||x||≥0，且当||x||＝0仅有x＝0；||ax||＝|a|||x||，其中a为任意常数；||x1 x2||≤||x1|| ||x2||，则称||x||为x空间上的范数(norm)，其中x1和x2为x空间上的任意两个值。对于定义1中的n×1维离散信号向量x＝[x1,x2,…,xn]t，向量x的lp范数表达式为其中xi为向量x的第i个元素，∑|·|表示绝对值求和运算符号，向量x的l1范数表达式为向量x的l2范数表达式为向量x的l0范数表达式为且xi≠0。详见文献“矩阵理论”，黄廷祝等编著，高等教育出版社出版。
定义4、方位向、距离向
将雷达平台运动的方向叫做方位向，将垂直于方位向的方向叫做距离向。
定义5、压缩感知稀疏重构理论
如果一个信号是稀疏的或可压缩的，那么该信号就可以用远低于奈奎斯特采样定理所要求的采样率来无失真的重构出该信号。如果信号稀疏，并且测量矩阵满足不相干和rip属性，使用压缩感知恢复的信号稀疏重建可以通过解决以下最优化问题来实现：
其中，α是估计信号，y是测量信号，θ是测量矩阵，ε是噪声门限。详见文献“阵列三维合成孔径雷达稀疏成像技术研究”韦顺军，2013。
定义6、迭代最小化稀疏贝叶斯重构算法，简称sbrim方法
迭代最小化稀疏贝叶斯重构算法(sparseautofocusbayesianrecoveryviaiterativeminimum)由电子科技大学的韦顺军副教授于2013年提出,详见文献“阵列三维合成孔径雷达稀疏成像技术研究”韦顺军，2013。
定义7、合成孔径雷达原始回波仿真方法
合成孔径雷达原始回波仿真方法是指基于合成孔径雷达成像原理仿真出一定系统参数条件下具有合成孔径雷达回波信号特性的原始信号的方法，详见文献“张朋，合成孔径雷达回波信号仿真研究，西北工业大学博士论文，2004”。
定义8、线阵sar的快时刻和慢时刻
线阵sar运动平台飞过一个方位向合成孔径长度所需要的时间称为慢时间，雷达系统以一定时间长度的重复周期发射接收脉冲，因此慢时间可以表示为一个以脉冲重复周期为步长的离散化时间变量，其中每一个脉冲重复周期离散时间变量值为一个慢时刻。快时刻是指在一个脉冲重复周期内，距离向采样回波信号的时间间隔变量。详见文献“合成孔径雷达成像原理”，皮一鸣等编著，电子科技大学出版社出版。
定义9、信号线性测量模型
对于一个数字信号测量系统，假设n×1维离散信号向量x＝[x1,x2,…,xn]t为该测量系统需要测量的信号，向量y＝[y1,y2,…,ym]t为该测量系统输出的m维离散信号向量，其中t为转置运算符号，y1为向量y中的第一个元素，y2表示向量y中的第二个元素，ym表示向量y中的第m个元素，信号的线性测量模型是指测量信号y和被测量信号x的关系可以表示为y＝ax，其中a为m×n矩阵，矩阵a为线性测量模型中信号x的测量矩阵。详见文献“阵列三维合成孔径雷达稀疏成像技术研究，韦顺军，2013”。
定义10、传统理论成像分辨率
线阵sar成像传统理论分辨率是指利用经典匹配滤波理论成像算法得到线阵sar系统在距离向、方位向和切航迹向的成像分辨率。对于收发共用天线，线阵sar距离向的分辨率记为ρr，近似表达式为其中c为电磁波在空气中传播的速度，br为线阵sar发射信号带宽；方位向分辨率记为ρa，近似表达为其中da为天线在方位向的真实孔径；切航迹向的分辨率记为其中λ为线阵sar雷达载波波长，r0为线阵sar平台到成像场景中心的参考斜距，l为阵列天线长度。
定义11、标准的降序排列方法
在正常的数值型数据排序中，升序排列是按照数据从高到低排列，即在要排序的一组数中，选出最大的一个数与第一个位置的数交换，然后在剩下的数当中再找最大的数与第二个位置的数交换，如此循环到倒数第二个数和最后一个数比较为止。
定义12、奇异值分解
设aha的特征值为λ1≥λ2≥…≥λr＞λr 1＝…＝λn＝0，则称为矩阵a的正奇异值；设σ1,σ2,…,σr，是a的r个正奇异值，则存在m阶酉矩阵u和存在n阶酉矩阵v，使得为矩阵a的奇异值分解，其中d＝diag(σ1,σ2,…,σr)。
本发明提出的一种基于相关向量机的线性阵列sar三维成像快速贝叶斯压缩感知方法，它包括如下步骤：
步骤1、初始化sar系统参数：
初始化sar系统参数包括：平台速度矢量记为阵列天线各阵元初始位置矢量，记做其中l为天线各阵元序号，nl为阵列天线的阵元总数；阵列天线长度，记做l；雷达发射信号载频为fc；雷达发射信号的调频斜率为fdr；脉冲重复时间记为pri；雷达系统的脉冲重复频率为prf；雷达发射信号带宽记做br；电磁波在空气中的传播速度记做c；距离向快时间记做n，n＝1,…,nr，nr为距离向快时刻总数，方位向慢时刻记做k，k＝1,…,k，k为方位向慢时刻总数；上述参数均为sar系统标准参数，其中雷达信号载频fc，雷达发射信号的调频斜率fdr，脉冲重复时间pri，雷达系统的脉冲重复频率prf，雷达发射信号带宽br，阵列天线的阵元总数nl，阵列天线长度l在线阵sar系统设计过程中已经确定；平台速度矢量阵列天线各阵元初始位置矢量在sar观测方案设计中已经确定；根据sar成像系统方案和观测方案，sar成像方法需要的初始化成像系统参数均为已知；
步骤2、划分sar的成像场景空间：
以雷达波束照射场区域地平面和垂直于该地平面向上的单位向量所构成的空间直角坐标系作为线阵sar的成像场景目标空间ω，其中水平横向和水平纵向构成阵列维成像空间；初始化水平横向成像场空间长度为lx，水平纵向成像场空间长度为ly；将成像场景目标空间ω均匀划分为大小相等的立体单元格，成像场景空间在水平横向、水平纵向单元格数分别为mx，my；
采用公式计算得到水平横向、水平纵向的单元格大小，分别记为dx和dy；成像场景空间高度向单元格总数为nr，高度向单元格的大小为线阵sar成像系统距离向分辨率，记为dz；
采用公式计算得到划分后的成像场景空间ω中第n个等距离单元格的阵列维成像空间中第mx个水平横向单元格第my个水平纵向单元格所对应的元素的位置，记为其中mx＝1，…,mx，my＝1，…,my；
采用公式得到成像场景空间散射系数矩阵，记做δ，其中m＝(my-1)mx mx＝1，…,m，n＝1,…,nr，nr为距离向快时刻总数，为划分后的成像场景空间ω中第n个等距离单元格的阵列维成像空间中第m个元素的散射系数，m＝mx·my为第n个等距离单元格的阵列维成像空间中的等效单元格总数；
步骤3、建立线阵sar(lineararraysar,lasar)的测量矩阵：
步骤3.1、在实际线阵sar成像中，原始回波数据由数据接收机提供；线阵sar中二维等效阵列中第i个天线阵元在第n个距离向快时刻的原始回波数据记做s(n,l,k)；采用标准合成孔径雷达距离向脉冲压缩方法对s(n,l,k)进行距离向脉冲压缩，得到距离向压缩后的线阵sar数据，记做sac(n,l,k)；
采用公式sn＝{sn(i)＝sac(n,l,k)}，i＝k×nl l，k＝1,…,k，l＝1,…,nl，计算得到第n个等距离单元格回波信号，记为sn，n＝1,…,nr，nl为步骤1中初始化得到的阵列天线阵元总数，k为步骤1中初始化得到的方位向慢时刻总数，其中nr为步骤1中初始化得到的距离向快时刻总数，na＝k×nl为回波信号的行数，即线阵sar中二维等效阵列中的元素的总数；
步骤3.2、采用公式计算得到第l个阵列天线在第k个方位向慢时刻的位置矢量，记为其中nl为步骤1中初始化得到的阵列天线阵元总数，k为步骤1中初始化得到的方位向慢时刻总数，为步骤1中初始化得到的阵列天线各阵元初始位置，为步骤1中初始化得到的平台速度，prf为步骤1中初始化得到的雷达系统的脉冲重复频率；
采用公式计算得到在第k个方位向慢时刻线阵sar成像场景目标空间ω中第n个等距离单元格中阵列维成像空间中第m个元素到第l个天线阵元的时间延时，记为
采用公式计算得到回波信号sn与散射系数矩阵δ之间的测量矩阵，记做θn；其中m为步骤2中得到的第n个等距离单元格的阵列维成像空间中的等效单元格总数，||·||2表示定义3中定义的向量l2范数，为步骤2中得到的划分后的成像场景目标空间ω中第n个等距离单元格的阵列维成像空间中第m个元素的位置，c为步骤1中初始化得到的电磁波在空气中的传播速度，sn为步骤3.1中得到的脉冲压缩后第n个等距离单元格回波信号，n＝1,…,nr，nr为距离向快时刻总数；
步骤4、对测量矩阵、回波信号、散射系数进行实数域表示：
采用公式和计算得到实数域回波信号和散射系数其中na为步骤3.1中得到的回波信号sn的行数，m为步骤2中得到的第n个等距离单元格的阵列维成像空间中的等效单元格总数，re(·)和im(·)分别表示复数域数据的实部和虚部，为划分后的成像场景空间ω中第n个等距离单元格的阵列维成像空间中第m个元素的散射系数；采用公式计算得到实数域测量矩阵，记为sn为步骤3.1中得到的第n个等距离单元格回波信号，θn为步骤3.2中得到的第n个等距离单元格的复数域测量矩阵；
步骤5、初始化目标区域：
步骤5.1、初始化噪声方差β0；
步骤5.2、采用公式计算得到的第m列向量，记为其中为步骤4中的实数域测量矩阵；采用公式计算得到第n个等距离单元格中第m个散射系数，记为其中为步骤4中第n个等距离单元格的实数域回波信号，||·||为定义3中的范数，∑(·)为求和运算符号，(·)t为矩阵转置的运算符号；
采用公式初始化第n个等距离单元格回波信号的散射系数，记为其中m为步骤2中得到的第n个等距离单元格的阵列维成像空间中的等效单元格总数；采用公式初始化超参数，记为其中为第m个散射系数对应的超参数；
步骤5.3、遍历步骤4中的散射系数得到最大元素对应的散射单元，记为id0；并根据公式得到初始化目标区域，记为
采用公式计算得到散射系数中的最大值对应的超参数，记为其中为的第id0列元素，为步骤4中实数域测量矩阵，为步骤4中第n个等距离单元格的实数域回波信号，为步骤5.2中第n个等距离单元格回波信号的散射系数的初始值，||·||为定义3中定义的范数，β0为步骤5.1初始化得到的噪声方差；
步骤6、计算最优超参数来提取目标区域：
步骤6.1、初始化迭代终止门限εe，初始化第i个超参数的迭代误差其中i＝1,…,2m，m为步骤2中得到的第n个等距离单元格的阵列维成像空间中的等效单元格总数；采用公式初始化超参数的迭代误差，记为e0；
在第t次迭代中，若t＝1，则为初始目标区域；若t＞1，为经过t-1迭代得到的目标区域；其中为步骤5.3中得到的初始目标区域；
步骤6.2、采用公式计算得到第t-1次迭代的超参数的对角矩阵，记为其中为第t-1次迭代中第i个散射系数对应的超参数，m为步骤2中得到的第n个等距离单元格的阵列维成像空间中的等效单元格总数；
采用公式计算得到协方差矩阵，记为ct-1；其中ι为单位矩阵，β为信号噪声的方差，为步骤4中得到的实数域测量矩阵，na为步骤3.1中得到的回波信号sn的行数，(·)t为步骤5.2中矩阵转置的运算符号；
采用公式更新第i个超参数对应的稀疏因子、质量因子，记为fit和其中为的第i列元素，为步骤4中第n个等距离单元格的实数域回波信号；
步骤6.3、利用步骤6.2中的fit和采用公式
采用公式更新第t次迭代后的第i个超参数的估计误差，记为并采用公式计算得到所有超参数的迭代误差，记为et；其中m为步骤2中得到的第n个等距离单元格的阵列维成像空间中的等效单元格总数；
步骤6.4、如果则迭代终止，采用公式得到最终目标区域，记为其中e0为步骤6.1中所有超参数的迭代误差的初始值，||·||为定义3中定义的范数，εe为步骤6.1中的迭代终止门限；
否则，选取第t次迭代的eit中最大误差项对应的散射单元，记为idt；
采用公式更新步骤6.1中的目标区域记为其中为步骤6.2中第t次迭代中第idt个散射单元的为步骤6.2中第t次迭代中第idt个散射单元的fit，为第t次迭代中第idt个散射单元的超参数，∪为添加元素运算符号，/为删除元素运算符号；
步骤6.5、采用公式计算得到第t次迭代后的目标区域对应的测量矩阵，记为φt；其中为步骤4中实数域测量矩阵的第idi列元素，i＝1,…,t；
步骤7、目标区域转换，得到复数域目标区域；
采用公式计算得到复数域目标区域，记为gn；其中nn为目标区域中的元素个数，ωr为目标区域gn中第r个元素在成像场景中的序号，为实值目标区域中第r个元素，round(·)为四舍五入取整运算符号；
步骤8、利用目标区域简化测量矩阵和散射系数：
采用公式计算得到θn的第ωr列向量，记为其中θn为步骤4中的复数域测量矩阵，na是步骤3.1中得到的回波信号sn的行数，ωr为步骤7中的目标区域gn中第r个元素在成像场景中的序号；
采用公式计算得到目标区域的简化散射系数向量，记为αω；其中为序号ωr对应的散射系数；采用公式计算得到目标区域的简化测量矩阵，记为θω；其中nn为步骤7中目标区域gn中的元素个数；采用公式初始化目标区域的散射系数估计值记为其中(·)h为hermite矩阵运算符号，其中sn为步骤3.1中得到的第n个等距离单元格回波信号；
步骤9、估计目标区域的散射系数：
步骤9.1、初始化正则化参数λ、平滑因子η、范数项p、迭代次数t、最大迭代次数nmax；误差终止迭代门限ε0；
步骤9.2、采用公式计算得到第t次迭代的对角矩阵，记为ωr为步骤7中目标区域gn中第r个元素在成像场景中的序号，为步骤8中测量矩阵θω的第r列向量，为第t-1次迭代后信号噪声方差的估计值，为第t-1次迭代后估计值的第r个元素，na是步骤3.1中得到的回波信号sn的行数，nn为步骤7中目标区域gn中的元素个数，diag{·}为对角矩阵符号，(·)h为步骤8中hermite矩阵运算符号，λ为步骤9.1中初始化得到的正则化参数，p为步骤9.1中初始化得到范数项；η为步骤9.1中初始化得到的平滑因子；
步骤9.3、采用公式计算得到的秩，记为k；rank(·)为矩阵的秩运算符号，为步骤9.2中第t次迭代的对角矩阵。
如果满足k＝nn；采用公式计算得到第t次迭代的散射系数和信号噪声方差的估计值，分别记为和θω为步骤8中的简化测量矩阵，na是步骤3.1中得到的回波信号sn的行数，sn为步骤3.1中得到的第n个等距离单元格回波信号，(·)h为步骤8中hermite矩阵运算符号，||·||为定义3中定义的范数，其中k为步骤9.2中得到对角矩阵的秩；；
否则，采用公式对对角矩阵进行定义12中的奇异值分解，其中是的第i个奇异值，并且满足和为对应于的奇异向量，nn为步骤7中目标区域gn中的元素个数，∑(·)为步骤5.2中求和运算符号；采用公式计算得到广义逆矩阵，记为其中k为步骤9.2中的秩；采用公式计算得到第t次迭代的散射系数的估计值，记为其中θω为步骤8中的简化测量矩阵，(·)h为步骤8中hermite矩阵运算符号；采用定义11中的降序排列对按照幅值以降序排列，并得到其中第k个元素，记为采用公式重新提取目标区域，记为gn；其中第t次迭代中序号为ωr的散射系数的估计值，ωr为步骤7中目标区域gn中第r个元素在成像场景中的序号；
步骤10、判断是否终止迭代：
如果且t＜nmax，则继续执行步骤9，且t＝t 1，其中为步骤9.3中得到的t次迭代后第n个等距离单元格的阵列维成像空间的成像结果，nmax为步骤9.1中初始化得到的高分辨率成像的最大迭代次数，ε0为步骤9.1中初始化得到的误差迭代终止门限；
若不满足和n≥nmax任一条件，算法迭代终止，则输出得到的高分辨率成像中第n个等距离单元格的散射系数向量估计值
步骤11、全场景三维成像：
采用公式将各等距离单元格散射系数向量排列为三维矩阵形式，得到线阵sar成像场景空间的三维高分辨率的成像结果，记为其中nr为步骤1中初始化得到的距离向快时刻总数，其中为步骤10中得到的高分辨率成像中第n个等距离单元格散射系数估计结果；
至此，我们全场景线阵sar的三维成像结果，整个重构方法结束。
本发明的创新点是：针对线阵sar三维成像的运算量巨大的问题，本发明在定义6中的sbrim算法的基础上，结合相关向量机，首先利用相关向量机提取成像空间中的目标区域，然后根据目标区域简化测量矩阵，并利用简化的测量矩阵估计目标区域的散射系数，并利用截断奇异值分解算法去除目标区域中与其他元素相比太小的元素，以保证散射系数正确估计。
本发明优点在于针对线阵sar三维成像的运算量巨大的问题，在定义6中的sbrim算法的基础上，结合了相关向量机，通过利用相关向量机实现目标与背景的分类进而提取成像场景空间中目标可能存在的区域，并根据目标可能存在的区域代替全成像场景空间简化测量矩阵，并估计目标区域的散射系数，该算法成功的避免了线阵sar三维成像中的高维度矩阵运算，极大地提高算法的运算效率；同时测量矩阵更好的表征了成像场景空间中的目标特性，更好的抑制了虚假目标、旁瓣干扰对于高质量成像的影响，成功的提高了线阵sar三维成像的成像质量。本算法具有重构精度高、运算效率较高的优势，本发明可适用于线阵合成孔径雷达三维成像等领域。
附图说明
图1为本发明流程图；
图2为系统参数表；
具体实施方式
本发明主要采用计算机仿真的方法进行验证，所有步骤、结论都在matlab-2017b上验证正确。具体实施步骤如下：
步骤1、初始化sar系统参数：
初始化sar系统参数包括：平台速度矢量记为阵列天线各阵元初始位置矢量，记做其中l为天线各阵元序号，nl＝64为阵列天线的阵元总数；阵列天线长度，记做l＝3m；雷达发射信号载频为fc＝37.5ghz；雷达发射信号的调频斜率为fdr＝4×1014hz/s；脉冲重复时间记为pri＝2μs；雷达系统的脉冲重复频率为prf＝0.5mhz；雷达发射信号带宽记做br＝0.8ghz；电磁波在空气中的传播速度记做c＝3×108m/s；距离向快时间记做n，n＝1,…,nr，nr＝512为距离向快时刻总数，方位向慢时刻记做k，k＝1,…,k，k＝64为方位向慢时刻总数；上述参数均为sar系统标准参数，其中雷达信号载频fc，雷达发射信号的调频斜率fdr，脉冲重复时间pri，雷达系统的脉冲重复频率prf，雷达发射信号带宽br，阵列天线的阵元总数nl，阵列天线长度l在线阵sar系统设计过程中已经确定；平台速度矢量阵列天线各阵元初始位置矢量在sar观测方案设计中已经确定；根据sar成像系统方案和观测方案，sar成像方法需要的初始化成像系统参数均为已知；
步骤2、划分sar的成像场景空间：
以雷达波束照射场区域地平面和垂直于该地平面向上的单位向量所构成的空间直角坐标系作为线阵sar的成像场景目标空间ω，其中水平横向和水平纵向构成阵列维成像空间；初始化水平横向成像场空间长度为lx＝50m，水平纵向成像场空间长度为ly＝70m；将成像场景目标空间ω均匀划分为大小相等的立体单元格，成像场景空间在水平横向、水平纵向单元格数分别为mx＝101，my＝101；根据公式计算得到水平横向、水平纵向的单元格大小，分别记为dx和dy；成像场景空间高度向单元格总数为nr＝512，高度向单元格的大小为线阵sar成像系统距离向分辨率，记为dz＝0.12m；根据公式得到划分后的成像场景空间ω中第n个等距离单元格的阵列维成像空间中第mx个水平横向单元格第my个水平纵向单元格所对应的元素的位置，记为其中mx＝1，…,101，my＝1，…,101；
根据公式初始化得到成像场景空间散射系数矩阵，记做δ，其中m＝(my-1)mx mx＝1，…,m，n＝1,…,nr，nr＝512为距离向快时刻总数，为划分后的成像场景空间ω中第n个等距离单元格的阵列维成像空间中第m个元素的散射系数，m＝mx·my＝10201为第n个等距离单元格的阵列维成像空间中的等效单元格总数；
步骤3、建立线阵sar(lineararraysar,lasar)的测量矩阵：
步骤3.1、在实际线阵sar成像中，原始回波数据由数据接收机提供；线阵sar中二维等效阵列中第i个天线阵元在第n个距离向快时刻的原始回波数据记做s(n,l,k)；采用标准合成孔径雷达距离向脉冲压缩方法对s(n,l,k)进行距离向脉冲压缩，得到距离向压缩后的线阵sar数据，记做sac(n,l,k)；
采用公式sn＝{sn(i)＝sac(n,l,k)}，i＝k×nl l，k＝1,…,k，l＝1,…,nl，计算得到第n个等距离单元格回波信号，记为sn，n＝1,…,nr，nl＝64为步骤1中初始化得到的阵列天线阵元总数，k＝64为步骤1中初始化得到的方位向慢时刻总数，其中nr＝512为步骤1中初始化得到的距离向快时刻总数，na＝k×nl＝4096为回波信号的行数，即线阵sar中二维等效阵列中的元素的总数；
步骤3.2、采用公式计算得到第l个阵列天线在第k个方位向慢时刻的位置矢量，记为其中nl＝64为步骤1中初始化得到的阵列天线阵元总数，k＝64为步骤1中初始化得到的方位向慢时刻总数，为步骤1中初始化得到的阵列天线各阵元初始位置，为步骤1中初始化得到的平台速度，prf＝0.5mhz为步骤1中初始化得到的雷达系统的脉冲重复频率；
采用公式计算得到在第k个方位向慢时刻线阵sar成像场景目标空间ω中第n个等距离单元格中阵列维成像空间中第m个元素到第l个天线阵元的时间延时，记为
采用公式计算得到回波信号sn与散射系数矩阵δ之间的测量矩阵，记做θn；其中m＝10201为步骤2中得到的第n个等距离单元格的阵列维成像空间中的等效单元格总数，||·||2表示定义3中定义的向量l2范数，为步骤2中得到的划分后的成像场景目标空间ω中第n个等距离单元格的阵列维成像空间中第m个元素的位置，c为步骤1中初始化得到的电磁波在空气中的传播速度；sn为步骤3.1中得到的脉冲压缩后第n个等距离单元格回波信号，n＝1,…,nr，nr＝512为距离向快时刻总数；
步骤4、对测量矩阵、回波信号、散射系数进行实数域表示：
采用公式和计算得到实数域回波信号和散射系数其中na＝4096为步骤3.1中得到的回波信号sn的行数，m＝10201为步骤2中第n个等距离单元格的阵列维成像空间中的等效单元格总数，re(·)和im(·)分别表示复数域数据的实部和虚部；为步骤2中得到的第n个等距离单元格散射单元中第m各散射单元的散射系数；采用公式计算得到实数域测量矩阵，记为sn为步骤3.1中得到的第n个等距离单元格回波信号，θn为步骤3.2中得到的第n个等距离单元格的复数域测量矩阵；
步骤5、初始化目标区域：
步骤5.1、初始化噪声方差β0＝37326；
步骤5.2、采用公式计算得到的第m列向量，记为其中为步骤4中的实数域测量矩阵；采用公式计算得到第n个等距离单元格中第m个散射系数，记为其中为步骤4中第n个等距离单元格的实数域回波信号，||·||为定义3中的范数，∑(·)为求和运算符号，(·)t为矩阵转置的运算符号；
采用公式初始化第n个等距离单元格回波信号的散射系数，记为其中m＝10201为步骤2中第n个等距离单元格的阵列维成像空间中的等效单元格总数；采用公式初始化超参数，记为其中为第m个散射系数对应的超参数；
步骤5.3、遍历步骤4中的散射系数得到最大元素对应的散射单元，记为id0；并根据公式得到初始化目标区域，记为
采用公式计算得到散射系数中的最大值对应的超参数，记为其中为的第id0列元素，为步骤4中实数域测量矩阵，为步骤4中第n个等距离单元格的实数域回波信号，为步骤5.2中第n个等距离单元格回波信号的散射系数的初始值，||·||为定义3中定义的范数，β0＝37326为步骤5.1初始化得到的噪声方差；
步骤6、计算最优超参数来提取目标区域：
步骤6.1、初始化迭代终止门限εe＝10-10，初始化第i个超参数的迭代误差其中i＝1,…,2m，m＝10201为步骤2中第n个等距离单元格的阵列维成像空间中的等效单元格总数；采用公式初始化超参数的迭代误差，记为e0；
在第t次迭代中，若t＝1，则为初始目标区域；若t＞1，为经过t-1迭代得到的目标区域；其中为步骤5.2中得到的初始目标区域；
步骤6.2、采用公式计算得到第t-1次迭代的超参数的对角矩阵，记为其中为第t-1次迭代中第i个散射系数对应的超参数，m＝10201为步骤2中第n个等距离单元格的阵列维成像空间中的等效单元格总数；
采用公式计算得到协方差矩阵，记为ct-1；其中ι为单位矩阵，β为信号噪声的方差，为步骤4中得到的实数域测量矩阵，(·)t为步骤5.2中矩阵转置的运算符号；
采用公式更新第i个超参数对应的稀疏因子、质量因子，记为fit和其中为的第i列元素，为步骤4中第n个等距离单元格的实数域回波信号；
步骤6.3、利用步骤6.2中的fit和采用公式
采用公式更新第t次迭代后的第i个超参数的估计误差，记为并采用公式计算得到所有超参数的迭代误差，记为et，其中m＝10201为步骤2中第n个等距离单元格的阵列维成像空间中的等效单元格总数；
步骤6.4、如果则迭代终止，采用公式得到最终目标区域，记为其中e0为步骤6.1中所有超参数的迭代误差的初始值，||·||为定义3中定义的范数，εe＝10-10为步骤6.1中的迭代终止门限；
否则，选取第t次迭代的中最大误差项对应的散射单元，记为idt；
采用公式更新步骤6.1中的目标区域记为其中为步骤6.2中第t次迭代中第idt个散射单元的为步骤6.2中第t次迭代中第idt个散射单元的fit，为第t次迭代中第idt个散射单元的超参数，∪为添加元素运算符号，/为删除元素运算符号；
步骤6.5、采用公式计算得到第t次迭代后的目标区域对应的测量矩阵，记为φt；其中为步骤4中实数域测量矩阵的第idi列元素，i＝1,…,t；
步骤7、目标区域转换，得到复数域目标区域；
采用公式计算得到复数域目标区域，记为gn；其中nn为目标区域中的元素个数，ωr为目标区域gn中第r个元素在成像场景中的序号，为实值目标区域中第r个元素，round(·)为四舍五入取整运算符号；
步骤8、利用目标区域简化测量矩阵和散射系数：
采用公式计算得到θn的第ωr列向量，记为其中θn为步骤4中的复数域测量矩阵，na＝4096为步骤3.1中得到的回波信号sn的行数，ωr为步骤7中的目标区域gn中第r个元素在成像场景中的序号；
采用公式计算得到目标区域的简化散射系数向量，记为αω；其中为序号ωr对应的散射系数；采用公式计算得到目标区域的简化测量矩阵，记为θω；其中nn为步骤7中目标区域gn中的元素个数；采用公式初始化目标区域的散射系数估计值记为其中(·)h为hermite矩阵运算符号，其中sn为步骤3.1中得到的第n个等距离单元格回波信号；
步骤9、估计目标区域的散射系数：
步骤9.1、初始化正则化参数λ＝1、平滑因子η＝10-5、范数项p＝0.8、迭代次数t＝1、最大迭代次数nmax＝20；误差终止迭代门限ε0＝10-10；
步骤9.2、采用公式计算得到第t次迭代的对角矩阵，记为ωr为步骤7中目标区域gn中第r个元素在成像场景中的序号，为步骤8中测量矩阵θω的第r列向量，为第t-1次迭代后信号噪声方差的估计值，为第t-1次迭代后估计值的第r个元素，na＝4096为步骤3.1中得到的回波信号sn的行数，nn为步骤7中目标区域gn中的元素个数，diag{·}为对角矩阵符号，(·)h为步骤8中hermite矩阵运算符号，λ＝1为步骤9.1中初始化得到的正则化参数，p＝0.8为步骤9.1中初始化得到范数项；η＝10-5为步骤9.1中初始化得到的平滑因子；
步骤9.3、采用公式计算得到的秩，记为k，rank(·)为矩阵的秩运算符号，为步骤9.2中第t次迭代的对角矩阵。
如果满足k＝nn；采用公式计算得到第t次迭代的散射系数和信号噪声方差的估计值，分别记为和θω为步骤8中的简化测量矩阵，na＝4096为步骤3.1中得到的回波信号sn的行数，sn为步骤3.1中得到的第n个等距离单元格回波信号，(·)h为步骤8中hermite矩阵运算符号，||·||为定义3中定义的范数，其中k为步骤9.2中得到对角矩阵的秩；
否则，采用公式对对角矩阵进行定义12中的奇异值分解，其中是的第i个奇异值，并且满足和为对应于的奇异向量，nn为步骤7中目标区域gn中的元素个数，∑(·)为步骤5.2中求和运算符号；采用公式计算得到广义逆矩阵，记为其中k为步骤9.2中的秩；采用公式计算得到第t次迭代的散射系数的估计值，记为其中θω为步骤8中的简化测量矩阵，(·)h为步骤8中hermite矩阵运算符号；采用定义11中的降序排列对按照幅值以降序排列，并得到其中第k个元素，记为采用公式重新提取目标区域，记为gn；其中第t次迭代中序号为ωr的散射系数的估计值，ωr为步骤7中目标区域gn中第r个元素在成像场景中的序号；
步骤10、判断是否终止迭代：
如果且t＜nmax，则继续执行步骤9，且t＝t 1，其中为步骤9.3中得到的t次迭代后第n个等距离单元格的阵列维成像空间的成像结果，nmax＝20为步骤9.1中初始化得到的高分辨率成像的最大迭代次数，ε0＝10-10为步骤9.1中初始化得到的误差迭代终止门限；
若不满足和n≥nmax任一条件，算法迭代终止，则输出得到的高分辨率成像中第n个等距离单元格的散射系数向量估计值
步骤11、全场景三维成像：
采用公式将各等距离单元格散射系数向量排列为三维矩阵形式，得到线阵sar成像场景空间的三维高分辨率的成像结果，记为其中nr＝512为步骤1中初始化得到的距离向快时刻总数，其中为步骤10中得到的高分辨率成像中第n个等距离单元格散射系数估计结果；
至此，我们全场景线阵sar的三维成像结果，整个重构方法结束。
经过计算机仿真及实测数据结果证明，本发明首先获得全成像空间中的目标区域，并利用目标区域简化测量矩阵并进行稀疏成像。本发明通过利用目标区域代替全成像场景空间进行高分辨率成像进而成功的避免了高维度的矩阵运算，同时极大地提升了线阵sar三维成像的运算效率。
技术特征：1.一种基于相关向量机的线性阵列sar三维成像快速贝叶斯压缩感知方法，其特征是它包括如下步骤：
步骤1、初始化sar系统参数：
初始化sar系统参数包括：平台速度矢量记为阵列天线各阵元初始位置矢量，记做其中l为天线各阵元序号，nl为阵列天线的阵元总数；阵列天线长度，记做l；雷达发射信号载频为fc；雷达发射信号的调频斜率为fdr；脉冲重复时间记为pri；雷达系统的脉冲重复频率为prf；雷达发射信号带宽记做br；电磁波在空气中的传播速度记做c；距离向快时间记做n，n＝1,…,nr，nr为距离向快时刻总数，方位向慢时刻记做k，k＝1,…,k，k为方位向慢时刻总数；上述参数均为sar系统标准参数，其中雷达信号载频fc，雷达发射信号的调频斜率fdr，脉冲重复时间pri，雷达系统的脉冲重复频率prf，雷达发射信号带宽br，阵列天线的阵元总数nl，阵列天线长度l在线阵sar系统设计过程中已经确定；平台速度矢量阵列天线各阵元初始位置矢量在sar观测方案设计中已经确定；根据sar成像系统方案和观测方案，sar成像方法需要的初始化成像系统参数均为已知；
步骤2、划分sar的成像场景空间：
以雷达波束照射场区域地平面和垂直于该地平面向上的单位向量所构成的空间直角坐标系作为线阵sar的成像场景目标空间ω，其中水平横向和水平纵向构成阵列维成像空间；初始化水平横向成像场空间长度为lx，水平纵向成像场空间长度为ly；将成像场景目标空间ω均匀划分为大小相等的立体单元格，成像场景空间在水平横向、水平纵向单元格数分别为mx，my；
采用公式计算得到水平横向、水平纵向的单元格大小，分别记为dx和dy；成像场景空间高度向单元格总数为nr，高度向单元格的大小为线阵sar成像系统距离向分辨率，记为dz；
采用公式计算得到划分后的成像场景空间ω中第n个等距离单元格的阵列维成像空间中第mx个水平横向单元格第my个水平纵向单元格所对应的元素的位置，记为其中mx＝1，…,mx，my＝1，…,my；
采用公式得到成像场景空间散射系数矩阵，记做δ，其中m＝(my-1)mx mx＝1，…,m，n＝1,…,nr，nr为距离向快时刻总数，为划分后的成像场景空间ω中第n个等距离单元格的阵列维成像空间中第m个元素的散射系数，m＝mx·my为第n个等距离单元格的阵列维成像空间中的等效单元格总数；
步骤3、建立线阵sar(lineararraysar,lasar)的测量矩阵：
步骤3.1、在实际线阵sar成像中，原始回波数据由数据接收机提供；线阵sar中二维等效阵列中第i个天线阵元在第n个距离向快时刻的原始回波数据记做s(n,l,k)；采用标准合成孔径雷达距离向脉冲压缩方法对s(n,l,k)进行距离向脉冲压缩，得到距离向压缩后的线阵sar数据，记做sac(n,l,k)；
采用公式sn＝{sn(i)＝sac(n,l,k)}，i＝k×nl l，k＝1,…,k，l＝1,…,nl，计算得到第n个等距离单元格回波信号，记为sn，n＝1,…,nr，nl为步骤1中初始化得到的阵列天线阵元总数，k为步骤1中初始化得到的方位向慢时刻总数，其中nr为步骤1中初始化得到的距离向快时刻总数，na＝k×nl为回波信号的行数，即线阵sar中二维等效阵列中的元素的总数；
步骤3.2、采用公式计算得到第l个阵列天线在第k个方位向慢时刻的位置矢量，记为其中nl为步骤1中初始化得到的阵列天线阵元总数，k为步骤1中初始化得到的方位向慢时刻总数，为步骤1中初始化得到的阵列天线各阵元初始位置，为步骤1中初始化得到的平台速度，prf为步骤1中初始化得到的雷达系统的脉冲重复频率；
采用公式计算得到在第k个方位向慢时刻线阵sar成像场景目标空间ω中第n个等距离单元格中阵列维成像空间中第m个元素到第l个天线阵元的时间延时，记为
采用公式i＝k×nl i，k＝1,…,k，l＝1,…,nl，m＝1,…,m，计算得到回波信号sn与散射系数矩阵δ之间的测量矩阵，记做θn；其中m为步骤2中得到的第n个等距离单元格的阵列维成像空间中的等效单元格总数，||·||2表示向量l2范数，为步骤2中得到的划分后的成像场景目标空间ω中第n个等距离单元格的阵列维成像空间中第m个元素的位置，c为步骤1中初始化得到的电磁波在空气中的传播速度，sn为步骤3.1中得到的脉冲压缩后第n个等距离单元格回波信号，n＝1,…,nr，nr为距离向快时刻总数；
步骤4、对测量矩阵、回波信号、散射系数进行实数域表示：
采用公式和i＝1,…,na，m＝1,…,m，计算得到实数域回波信号和散射系数其中na为步骤3.1中得到的回波信号sn的行数，m为步骤2中得到的第n个等距离单元格的阵列维成像空间中的等效单元格总数，re(·)和im(·)分别表示复数域数据的实部和虚部，为划分后的成像场景空间ω中第n个等距离单元格的阵列维成像空间中第m个元素的散射系数；采用公式计算得到实数域测量矩阵，记为sn为步骤3.1中得到的第n个等距离单元格回波信号，θn为步骤3.2中得到的第n个等距离单元格的复数域测量矩阵；
步骤5、初始化目标区域：
步骤5.1、初始化噪声方差β0；
步骤5.2、采用公式计算得到的第m列向量，记为其中为步骤4中的实数域测量矩阵；采用公式计算得到第n个等距离单元格中第m个散射系数，记为其中为步骤4中第n个等距离单元格的实数域回波信号，||·||为范数，∑(·)为求和运算符号，(·)t为矩阵转置的运算符号；
采用公式初始化第n个等距离单元格回波信号的散射系数，记为其中m为步骤2中得到的第n个等距离单元格的阵列维成像空间中的等效单元格总数；采用公式初始化超参数，记为其中为第m个散射系数对应的超参数；
步骤5.3、遍历步骤4中的散射系数得到最大元素对应的散射单元，记为id0；并根据公式得到初始化目标区域，记为
采用公式计算得到散射系数中的最大值对应的超参数，记为其中为的第id0列元素，为步骤4中实数域测量矩阵，为步骤4中第n个等距离单元格的实数域回波信号，为步骤5.2中第n个等距离单元格回波信号的散射系数的初始值，||·||为范数，β0为步骤5.1初始化得到的噪声方差；
步骤6、计算最优超参数来提取目标区域：
步骤6.1、初始化迭代终止门限εe，初始化第i个超参数的迭代误差其中i＝1,…,2m，m为步骤2中得到的第n个等距离单元格的阵列维成像空间中的等效单元格总数；采用公式初始化超参数的迭代误差，记为e0；
在第t次迭代中，若t＝1，则为初始目标区域；若t＞1，为经过t-1迭代得到的目标区域；其中为步骤5.3中得到的初始目标区域；
步骤6.2、采用公式计算得到第t-1次迭代的超参数的对角矩阵，记为其中为第t-1次迭代中第i个散射系数对应的超参数，m为步骤2中得到的第n个等距离单元格的阵列维成像空间中的等效单元格总数；
采用公式计算得到协方差矩阵，记为ct-1；其中ι为单位矩阵，β为信号噪声的方差，为步骤4中得到的实数域测量矩阵，na为步骤3.1中得到的回波信号sn的行数，(·)t为步骤5.2中矩阵转置的运算符号；
采用公式更新第i个超参数对应的稀疏因子、质量因子，记为fit和其中为的第i列元素，为步骤4中第n个等距离单元格的实数域回波信号；
步骤6.3、利用步骤6.2中的fit和采用公式更新第t次迭代后的第i个超参数，记为
采用公式更新第t次迭代后的第i个超参数的估计误差，记为并采用公式计算得到所有超参数的迭代误差，记为et；其中m为步骤2中得到的第n个等距离单元格的阵列维成像空间中的等效单元格总数；
步骤6.4、如果则迭代终止，采用公式得到最终目标区域，记为其中e0为步骤6.1中所有超参数的迭代误差的初始值，||·||为范数，εe为步骤6.1中的迭代终止门限；
否则，选取第t次迭代的中最大误差项对应的散射单元，记为idt；
采用公式更新步骤6.1中的目标区域记为其中为步骤6.2中第t次迭代中第idt个散射单元的为步骤6.2中第t次迭代中第idt个散射单元的fit，为第t次迭代中第idt个散射单元的超参数，∪为添加元素运算符号，/为删除元素运算符号；
步骤6.5、采用公式计算得到第t次迭代后的目标区域对应的测量矩阵，记为φt；其中为步骤4中实数域测量矩阵的第idi列元素，i＝1,…,t；
步骤7、目标区域转换，得到复数域目标区域；
采用公式计算得到复数域目标区域，记为gn；其中nn为目标区域中的元素个数，ωr为目标区域gn中第r个元素在成像场景中的序号，为实值目标区域中第r个元素，round(·)为四舍五入取整运算符号；
步骤8、利用目标区域简化测量矩阵和散射系数：
采用公式计算得到θn的第ωr列向量，记为其中θn为步骤4中的复数域测量矩阵，na是步骤3.1中得到的回波信号sn的行数，ωr为步骤7中的目标区域gn中第r个元素在成像场景中的序号；
采用公式计算得到目标区域的简化散射系数向量，记为αω；其中为序号ωr对应的散射系数；采用公式计算得到目标区域的简化测量矩阵，记为θω；其中nn为步骤7中目标区域gn中的元素个数；采用公式初始化目标区域的散射系数估计值记为其中(·)h为hermite矩阵运算符号，其中sn为步骤3.1中得到的第n个等距离单元格回波信号；
步骤9、估计目标区域的散射系数：
步骤9.1、初始化正则化参数λ、平滑因子η、范数项p、迭代次数t、最大迭代次数nmax；误差终止迭代门限ε0；
步骤9.2、采用公式计算得到第t次迭代的对角矩阵，记为ωr为步骤7中目标区域gn中第r个元素在成像场景中的序号，为步骤8中测量矩阵θω的第r列向量，为第t-1次迭代后信号噪声方差的估计值，为第t-1次迭代后估计值的第r个元素，na是步骤3.1中得到的回波信号sn的行数，nn为步骤7中目标区域gn中的元素个数，diag{·}为对角矩阵符号，(·)h为步骤8中hermite矩阵运算符号，λ为步骤9.1中初始化得到的正则化参数，p为步骤9.1中初始化得到范数项；η为步骤9.1中初始化得到的平滑因子；
步骤9.3、采用公式计算得到的秩，记为k；rank(·)为矩阵的秩运算符号，为步骤9.2中第t次迭代的对角矩阵；
如果满足k＝nn；采用公式计算得到第t次迭代的散射系数和信号噪声方差的估计值，分别记为和θω为步骤8中的简化测量矩阵，na是步骤3.1中得到的回波信号sn的行数，sn为步骤3.1中得到的第n个等距离单元格回波信号，(·)h为步骤8中hermite矩阵运算符号，||·||为范数，其中k为步骤9.2中得到对角矩阵的秩；；
否则，采用公式对对角矩阵进行奇异值分解，其中是的第i个奇异值，并且满足和为对应于的奇异向量，nn为步骤7中目标区域gn中的元素个数，∑(·)为步骤5.2中求和运算符号；采用公式计算得到广义逆矩阵，记为其中k为步骤9.2中的秩；采用公式计算得到第t次迭代的散射系数的估计值，记为其中θω为步骤8中的简化测量矩阵，(·)h为步骤8中hermite矩阵运算符号；采用标准的降序排列方法对按照幅值以降序排列，并得到其中第k个元素，记为采用公式重新提取目标区域，记为gn；其中第t次迭代中序号为ωr的散射系数的估计值，ωr为步骤7中目标区域gn中第r个元素在成像场景中的序号；
步骤10、判断是否终止迭代：
如果且t＜nmax，则继续执行步骤9，且t＝t 1，其中为步骤9.3中得到的t次迭代后第n个等距离单元格的阵列维成像空间的成像结果，nmax为步骤9.1中初始化得到的高分辨率成像的最大迭代次数，ε0为步骤9.1中初始化得到的误差迭代终止门限；
若不满足和n≥nmax任一条件，算法迭代终止，则输出得到的高分辨率成像中第n个等距离单元格的散射系数向量估计值步骤11、全场景三维成像：
采用公式将各等距离单元格散射系数向量排列为三维矩阵形式，得到线阵sar成像场景空间的三维高分辨率的成像结果，记为其中nr为步骤1中初始化得到的距离向快时刻总数，其中为步骤10中得到的高分辨率成像中第n个等距离单元格散射系数估计结果；
至此，我们全场景线阵sar的三维成像结果，整个重构方法结束。
技术总结本发明公开了一种基于相关向量机的快速贝叶斯压缩感知(FBCS‑RVM)方法，它是在SBRIM算法的基础上，结合了相关向量机，通过利用相关向量机实现目标与背景的分类进而提取成像场景空间中目标可能存在的区域，并根据目标可能存在的区域代替全成像场景空间简化测量矩阵，并估计目标区域的散射系数，本发明成功的避免了线阵SAR三维成像中的高维度矩阵运算，极大地提高算法的运算效率；同时测量矩阵更好的表征了成像场景空间中的目标特性，更好的抑制了虚假目标、旁瓣干扰对于高质量成像的影响，成功的提高了线阵SAR三维成像的成像质量。本发明具有重构精度高、运算效率较高的优势，本发明适用于线阵合成孔径雷达三维成像等领域。
技术研发人员：张晓玲;王楠;胥小我;田博坤;韦顺军;师君
受保护的技术使用者：电子科技大学
技术研发日：2021.04.30
技术公布日：2021.08.03
