基于l_q正则化的偏置相位中心天线成像方法

基于l_q 正则化的偏置相位中心天线成像方法
【专利摘要】本发明提供了一种基于lq正则化的偏置相位中心天线成像方法，其构建了偏置相位中心天线成像雷达的特性矩阵和重建矩阵，依据成像要求，选择lq正则化模型及其对应的阈值函数；基于成像雷达的回波数据及其特性矩阵和重建矩阵，得到DPCA数据处理算子；将该算子引入到迭代阈值运算过程中，计算出观测场景的后向散射系数矩阵。本发明能够利用由非均匀采样获取的偏置相位中心回波数据，实现对观测场景的无模糊雷达成像；在完成对大空间尺度观测场景的高分辨率成像工作时，相比由矩阵?向量乘法运算实现的稀疏微波成像算法，花费的时间更短、占用的资源更少；能够有效地抑制因方位向非均匀采样产生的杂波干扰，同时，对回波数据中的加性噪声具有较强的健壮性。
【专利说明】
基于Iq正则化的偏置相位中心天线成像方法
技术领域
[0001] 本发明涉及稀疏微波成像技术领域，尤其涉及一种基于l<〇H则化的偏置相位中心 天线成像方法。
【背景技术】
[0002] 随着合成孔径雷达技术的发展，微波遥感在工程中的实际应用对现代成像雷达系 统的性能要求越来越高。作为一种能够实现高分辨率、宽测绘带微波成像的有效方法，偏置 相位中心天线(Displaced Phase CenterAntenna，DPCA)技术在高性能合成孔径雷达成像 领域获得了广泛使用。为了获得最优的成像性能，偏置相位中心天线成像雷达系统通常需 要满足在方位向上符合奈奎斯特采样定理要求的均匀采样条件，即运载平台在雷达系统发 射相邻信号脉冲的时间间隔内，它移动的距离等于天线沿顺轨方向总长度的一半。可是，对 于某些特定入射角度下的观测约束，根据上述条件求得的脉冲重复频率（Pulse Repetition Frequency，PRF)很可能不符合斑马图的要求，由此便会导致在合成孔径时间 内对方位向接收信号进行非均匀采样，致使一些经典的合成孔径雷达成像算法(例如:距离 多普勒成像算法)无法直接应用于上述回波数据的成像处理。
[0003] 目前，解决该问题的方法主要分为两种，即基于重建滤波器组的多普勒频谱重建 方法，以及基于压缩感知的偏置相位中心天线成像方法。基于重建滤波器组的多普勒频谱 重建方法将偏置相位中心天线成像雷达系统的数据采集过程近似地看成由单站合成孔径 雷达系统接收的回波信号经过一组线性滤波器处理后获取的数据。当回波数据在方位向上 满足广义的香农-奈奎斯特采样定理时，完整的无混叠多普勒频谱便可以利用由上述线性 滤波器组所决定的权值函数，处理每个接收通道的回波数据来重建获得。但是，当回波数据 中存在加性噪声时，方位向的非均匀采样会导致重建的雷达图像中出现大量的杂波干扰， 给图像的解译工作带来极大的困难。而基于压缩感知的偏置相位中心天线成像方法则是根 据成像雷达系统接收的回波数据与观测场景的后向散射系数之间的时域关系，来构建精确 的雷达观测模型。然后，利用稀疏重建算法对该模型进行求解，便可获得观测场景的无模糊 雷达图像。但是，当所要重建的观测场景空间维度较大时，因包含观测矩阵的矩阵-向量乘 法运算所产生的庞大计算量与内存损耗会使现有的数据处理器难以承受，从而限制了基于 压缩感知的偏置相位中心天线成像方法在高性能合成孔径雷达成像领域中的实际应用。

【发明内容】

[0004] (一)要解决的技术问题
[0005] 为了解决现有技术问题，本发明提供了一种基于l<〇H则化的偏置相位中心天线成 像方法。
[0006] (二)技术方案
[0007] 本发明提供了一种基于l<〇H则化的偏置相位中心天线成像方法，包括:步骤A:根 据偏置相位中心天线成像雷达观测模型的二次相位近似，基于成像雷达系统参数，得到成 像雷达的特性矩阵和重建矩阵;步骤B:选择lq正则化模型及其对应的阈值函数;步骤C:基 于成像雷达的回波数据以及特性矩阵和重建矩阵，得到DPCA数据处理算子；以及步骤D:基 于成像雷达的回波数据和DPCA数据处理算子，利用l<〇H则化模型的阈值函数进行迭代阈值 运算，得到观测场景的后向散射系数矩阵。
[0008] (三)有益效果
[0009] 从上述技术方案可以看出，本发明的基于l<〇H则化的偏置相位中心天线成像方法 具有以下有益效果：
[0010] (1)当偏置相位中心天线成像雷达的脉冲重复频率不满足在方位向上的均匀采样 条件时，本发明能够利用由非均匀采样获取的偏置相位中心回波数据，重建无混叠的多普 勒频谱，从而实现对观测场景的无模糊雷达成像；
[0011] (2)由于利用DPCA数据处理算子代替了迭代阈值算法中观测矩阵的功能，因此，在 完成针对大空间尺度观测场景的高分辨率成像工作时，相比由矩阵_向量乘法运算实现的 稀疏微波成像算法，本发明所要花费的时间更短、占用的资源更少；
[0012] (3)当利用由非均匀采样获取的偏置相位中心回波数据重建观测场景的雷达图像 时，本发明能够有效地抑制因方位向非均匀采样产生的杂波干扰，同时，对回波数据中的加 性噪声具有较强的健壮性。
【附图说明】
[0013] 图1为根据本发明实施例的基于lq正则化的偏置相位中心天线成像方法流程图；
[0014] 图2(a)、图2(b)和图2(c)为当回波数据中无加性噪声且方位向均匀采样符合度等 于33%时，利用不同的成像方法对雷达图像中方位模糊的去除结果；
[0015]图3(a)、图3(b)、图3(c)、图3(d)、图3(e)和图3(f)为当信噪比等于20dB且方位向 均匀采样符合度分别取100%、33%和5%时，利用不同的成像方法对雷达图像中杂波干扰 的抑制结果；
[0016] 图4为不同的成像方法抗加性噪声干扰能力的测试结果。
【具体实施方式】
[0017] 稀疏信号处理作为近些年来信号与信息处理研究领域的前沿性课题，已引起数学 界与工程学界相关领域学者与专家的广泛关注。在稀疏信号处理理论中，lq正则化(〇<q< 1)是用来解决在稀疏约束下信号重建问题的有效方法。由lq正则化模型推导获得的迭代阈 值算法通常具有较快的收敛速度、较强的抗加性噪声干扰的能力和简明的实现方式，所以 这类算法在与稀疏信号处理有关的工程领域获得了广泛使用。而在合成孔径雷达成像应用 中，已有的研究成果表明，基于lq正则化的迭代阈值算法能够有效地抑制雷达图像中的旁 瓣与杂波干扰，并可以去除因频谱混叠造成的距离/方位模糊。
[0018] 为此，我们根据lq正则化理论，依照迭代阈值算法的算法框架，推导出一种新的偏 置相位中心天线成像方法。通过将基于广义的香农-奈奎斯特采样定理获得的DPCA数据处 理算子引入到该方法的实现过程中，以代替根据偏置相位中心天线成像雷达观测模型所构 建的观测矩阵的功能，新的成像方法可以在花费较短时间和较少资源的情况下，高效地完 成针对大空间尺度观测场景的高分辨率成像工作。仿真实验的实验结果表明，该种方法可 以利用由非均匀采样获取的偏置相位中心回波数据无模糊地重建出方位向图像，并能有效 地抑制因方位向非均匀采样产生的杂波干扰，同时，也对回波数据中的加性噪声具有较强 的健壮性。
[0019] 为使本发明的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照 附图，对本发明做进一步地详细说明。
[0020] 本发明的实施例提供了一种基于lq正则化的偏置相位中心天线成像方法，其具体 包括：
[0021] 步骤A:根据偏置相位中心天线成像雷达观测模型的二次相位近似，基于成像雷达 系统参数，得到成像雷达的特性矩阵和重建矩阵。
[0022]步骤A具体包括：
[0023]子步骤A1:根据偏置相位中心天线成像雷达观测模型的二次相位近似，得到成像 雷达接收通道的特性滤波器在多普勒频域的脉冲响应函数。
[0024]子步骤A1具体包括：
(1)
[0026] 其中，Hdf)表示成像雷达第i个接收通道的特性滤波器在多普勒频域的脉冲响应 函数，△11表示成像雷达第i个接收天线与成像雷达发射天线之间的距离；iG[l，2，...， N]，N为成像雷达的接收通道数;A表示成像雷达信号的载波波长;RQ表示成像雷达与观测目 标间的最短斜距;v表示成像雷达运载平台的运动速度，f为雷达信号的频率。
[0027] 子步骤A2:由成像雷达接收通道的特性滤波器在多普勒频域的脉冲响应函数，得 到成像雷达的特性矩阵。
[0028] 当成像雷达接收通道对方位向回波信号的采样频率等于脉冲重复频率PRF时，可 以利用由特性滤波器Hdf)构建的特性矩阵H(f)，描述具有N个接收通道的偏置相位中心天 线成像雷达的系统特性，从而，
[0029] 子步骤A2具体包括:特性矩阵可以表示为： 一 哺） 哺) … 哺、 ' , "(: PK.F) (f 4 FRF) … HJ f- FKF) _ ,
[0030] H(/> , 1 , ， (2J //;(/-;. (Ar. /人(Ar-F;FRF)
[0031] 其中，iMf)表示成像雷达第i个接收通道的特性滤波器在多普勒频域的脉冲响应 函数，i = l，2,. . .，N，N为成像雷达的接收通道数;PRF为脉冲重复频率。
[0032] 子步骤A3:由成像雷达的特性矩阵得到成像雷达的重建矩阵。
[0033] 由广义的香农-奈奎斯特采样定理可知，对于一个频带受限的信号，当使用1/N的 奈奎斯特采样频率对其N个独立表征进行采样时，原信号的非混叠频谱可以由上述采样数 据得到的N个混叠频谱经过相应的加权处理得到，因此，
[0034]子步骤A3具体包括： (3)
[0036] 其中，P(f)为成像雷达的重建矩阵，Pdf)为成像雷达第i个接收通道的重建滤波 器，该重建滤波器由N个带通滤波器P^(f)(j = l，…，N)组成，每个带通滤波器P^(f)的通带 范围是,由此可见，带通滤波器Pdf)构成该重建矩阵。
[0037] 步骤B:选择lq正则化模型及其对应的阈值函数。
[0038] 在得到iMfWPPdf)之后，需要根据成像雷达的具体要求，选定合适的lq正则化(0 <q<l)模型，由该模型推导出适用于迭代阈值运算的阈值函数。
[0039] 优选地，可以选择h正则化模型或11/2正则化模型及其对应的阈值函数。
[0040] 对于h正则化模型，其阈值函数表示为：
(4)
[0042]对于11/2正则化模型，其阈值函数表示为：
(5)
[0045] 其中，z表示复数域的输入量;k>〇,表示正则化参数;]i>〇表示迭代阈值运算收敛 性的控制参数。
[0046] 步骤C:基于成像雷达的回波数据以及重建矩阵和特性矩阵，得到DPCA数据处理算 子。
[0047] 为了提高迭代阈值算法的成像效率，根据雷达成像原理与逆成像回波仿真思想， 构建了能够用于处理由非均匀采样获取的偏置相位中心回波数据的算子，即DPCA数据处理 算子。该算子包括雷达图像重建项与回波数据生成项，分别对应雷达图像重建过程与回波 数据生成过程，其中，雷达图像重建过程可以被划分为四个主要的数据处理步骤，即距离向 脉冲压缩、无混叠多普勒频谱重建、距离徙动校正和方位向脉冲压缩。
[0048]步骤C具体包括：
[0049] 子步骤C1:基于成像雷达的回波数据以及重建矩阵，得到DPCA数据处理算子的雷 达图像重建项，该雷达图像重建项表达式为：
[0050] ：T(Y) = Ff 既 C)叫泣 ?[^[[[單和％]#^^ (6)
[0051] 其中，Yi表示由成像雷达第i个接收通道采样获取的回波数据组成的回波数据矩 阵;Y表示由成像雷达采样获取的回波数据组成的回波数据矩阵，且Y= U A;FjPFa分别表 示在距离向与方位向上用来处理数据的离散傅里叶变换矩阵，而则是它们的共辄 转置;MjPMa分别表示在距离向与方位向上用来处理数据的频域匹配滤波器矩阵;Pi表示由 重建矩阵中的重建滤波器Pdf)的解析式计算获得的重建滤波器矩阵表示用来完成 距离徙动校正的sine插值算子;?表示哈达马乘积运算。
[0052]由sine插值的基本原理可知，该图像重建过程是可逆的，因此，
[0053]子步骤C2:对雷达图像重建项取逆，获得DPCA数据处理算子的回波数据生成项，回 波数据生成项的表达式为：
[0054] 朵(X)= [[If[明If (7)
[0055] 其中，X表示观测场景的后向散射系数矩阵表示由特性矩阵中的特性滤波器出 (f)的解析式计算获得的特性滤波器矩阵;上角标*表示取矩阵元素的共辄，上角标-1表示 相应算子的逆过程。
[0056]步骤D:基于成像雷达的回波数据和DPCA数据处理算子，利用lq正则化模型的阈值 函数进行迭代阈值运算，得到观测场景的后向散射系数矩阵。
[0057]步骤D具体包括:将成像雷达的回波数据组成的回波数据矩阵、DPCA数据处理算子 引入到由lq正则化模型推得的阈值函数中，进行如下迭代阈值运算：
[0058] X"'山=' + 成(yx，j 08)
[0059] 其中，^为由成像雷达第i个接收通道采样获取的回波数据组成的回波数据矩阵， Y为由成像雷达采样获取的回波数据组成的回波数据矩阵，且Y=UA; k为正则化参数;y为 收敛性控制参数，并令XW = 0、设最大迭代次数为K;
[0060] 当迭代次数达到K或满足收敛精度，得到观测场景的后向散射系数矩阵的估计X。 [0061 ]其中，该收敛精度具体包括：
[0062] p(k+1) = | | x(k+1)-X(k) | 12/ | | X(k) | 12 (9)
[0063] P(k+1)小于1(T6,在(9)式中，|卜| |2表示取矩阵的二范数运算；即当P(k+1)小于1(T6 时结束迭代阈值运算，得到观测场景的后向散射系数矩阵的估计i。
[0064] 根据上述实施例，本发明利用仿真实验，对本实施例所述方法的有益效果进行验 证。在仿真实验中，偏置相位中心天线成像雷达的系统参数如表1所示。
[0065] 表 1
[0067]在仿真实验一中，观测场景内放置了一个面散射体和两个点散射体，重建雷达图 像所用回波数据无加性噪声干扰，且其方位向均匀采样符合度等于33%。图2(a)~图2(c) 为利用不同的成像方法对雷达图像中方位模糊的去除结果。其中，图2(a)是利用距离多普 勒成像算法处理单通道回波数据获得的对雷达图像中方位模糊的去除结果；图2(b)是利用 基于重建滤波器组的多普勒频谱重建方法对雷达图像中方位模糊的去除结果；图2(c)是利 用基于l q正则化的偏置相位中心天线成像方法对雷达图像中方位模糊的去除结果。比较图 2所示的仿真实验结果可知，相比距离多普勒成像算法，基于重建滤波器组的多普勒频谱重 建方法和基于l q正则化的偏置相位中心天线成像方法均能够有效地去除雷达图像中的方 位模糊。因为由成像雷达系统两个接收通道采样获取的回波数据并不是严格独立的，所以 由基于重建滤波器组的多普勒频谱重建方法重建的雷达图像中仍存在少量的方位模糊。而 由基于l q正则化的偏置相位中心天线成像方法重建的雷达图像中则完全没有方位模糊存 在。
[0068]在仿真实验二中，观测场景内放置了三个幅值相差10dB的面散射体，重建雷达图 像所用回波数据的信噪比为20dB。图3(a)~图3(f)为利用不同的成像方法对雷达图像中杂 波干扰的抑制结果。其中，图3(a)、图3(c)和图3(e)分别是在方位向均匀采样符合度取 100%、33%和5%的条件下，利用基于重建滤波器组的多普勒频谱重建方法对雷达图像中 杂波干扰的抑制结果；图3(b)、图3(d)和图3(f)分别是在方位向均匀采样符合度取100%、 33%和5%的条件下，利用基于l q正则化的偏置相位中心天线成像方法对雷达图像中杂波 干扰的抑制结果。由图3(a)、图3(c)和图3(e)可知，随着方位向均匀采样符合度的减小，SP 方位向采样间隔非均匀程度的增大，在由基于重建滤波器组的多普勒频谱重建方法获得的 雷达图像中，杂波干扰对图像质量的影响越来越大，甚至造成较弱的观测目标被杂波电平 完全淹没。而基于l q正则化的偏置相位中心天线成像方法则能够有效地抑制杂波干扰对重 建结果图像质量的影响（如图3(b)、图3(d)和图3(f)所示）。
[0069]在仿真实验三中，为了研究不同的成像方法抗加性噪声干扰的能力，我们在信噪 比分别取5dB、10dB和15dB的条件下，绘制出基于重建滤波器组的多普勒频谱重建方法和基 于lq正则化的偏置相位中心天线成像方法随方位向均匀采样符合度变化的雷达图像信杂 比曲线。值得注意的是，在图4的标注中，"DSR"表示基于重建滤波器组的多普勒频谱重建方 法，而"Lq-DIA"表示基于l q正则化的偏置相位中心天线快速成像方法。根据图4所示的仿真 结果可知，在同一横坐标处，对应信噪比为5dB的Lq-DIA信杂比曲线的取值几乎都要高于对 应信噪比为15dB的DSR信杂比曲线的取值。因此，相比基于重建滤波器组的多普勒频谱重建 方法，基于l q正则化的偏置相位中心天线成像方法对回波数据中的加性噪声具有更强的健 壮性。
[0070]至此，已经结合附图对本实施例进行了详细描述。依据以上描述，本领域技术人员 应当对本发明的基于lq正则化的偏置相位中心天线成像方法有了清楚的认识。
[0071]本发明的基于l<〇H则化的偏置相位中心天线成像方法，不但能够利用由非均匀采 样获取的偏置相位中心回波数据，高效地完成针对大空间尺度观测场景的高分辨率成像工 作，还可以有效地去除雷达图像中的方位模糊，并抑制因方位向非均匀采样产生的杂波干 扰，同时，也对回波数据中的加性噪声具有较强的健壮性。
[0072]需要说明的是，在附图或说明书正文中，未绘示或描述的实现方式，均为所属技术 领域中普通技术人员所知的形式，并未进行详细说明。此外，上述对各元件的定义并不仅限 于实施例中提到的各种方式，本领域普通技术人员可对其进行简单地更改或替换，例如：
[0073] (1)实施例中提到的方向用语，例如"上"、"下"、"前"、"后"、"左"、"右"等，仅是参 考附图的方向，并非用来限制本发明的保护范围；
[0074] (2)上述实施例可基于设计及可靠度的考虑，彼此混合搭配使用或与其他实施例 混合搭配使用，即不同实施例中的技术特征可以自由组合形成更多的实施例。
[0075] 以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详 细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限制本发明，凡 在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保 护范围之内。
【主权项】
1. 一种基于lq正则化的偏置相位中心天线成像方法，其特征在于，包括： 步骤A:根据偏置相位中心天线成像雷达观测模型的二次相位近似，基于成像雷达系统 参数，得到成像雷达的特性矩阵和重建矩阵； 步骤B:选择lq正则化模型及其对应的阈值函数； 步骤C:基于成像雷达的回波数据以及特性矩阵和重建矩阵，得到DPCA数据处理算子； 以及 步骤D:基于成像雷达的回波数据和DPCA数据处理算子，利用1<〇Η则化模型的阈值函数 进行迭代阈值运算，得到观测场景的后向散射系数矩阵。2. 如权利要求1所述的基于lq正则化的偏置相位中心天线成像方法，其特征在于，所述 步骤A具体包括： 子步骤A1:根据偏置相位中心天线成像雷达观测模型的二次相位近似，得到成像雷达 接收通道的特性滤波器在多普勒频域的脉冲响应函数； 子步骤A2:由成像雷达接收通道的特性滤波器在多普勒频域的脉冲响应函数，得到成 像雷达的特性矩阵；以及 子步骤A3:由成像雷达的特性矩阵得到成像雷达的重建矩阵。3. 如权利要求2所述的基于1<〇Η则化的偏置相位中心天线成像方法，其特征在于，所述 步骤C具体包括： 子步骤C1:基于成像雷达的回波数据以及重建矩阵，得到DPCA数据处理算子的雷达图 像重建项，所述雷达图像重建项表达式为：其中，h表示由成像雷达第i个接收通道采样获取的回波数据组成的回波数据矩阵;Y表 示由成像雷达采样获取的回波数据组成的回波数据矩阵，且分别表示在距 离向与方位向上用来处理数据的离散傅里叶变换矩阵，和if则是它们的共辄转置;M r和 Ma分别表示在距离向与方位向上用来处理数据的频域匹配滤波器矩阵;？1表示由重建矩阵 中的重建滤波器?:^)的解析式计算获得的重建滤波器矩阵；表示用来完成距离徙动 校正的sine插值算子；0表示哈达马乘积运算；以及 子步骤C2:对雷达图像重建项取逆，获得DPCA数据处理算子的回波数据生成项，回波数 据生成项的表达式为：其中，X表示观测场景的后向散射系数矩阵表示由特性矩阵中的特性滤波器Hdf)的 解析式计算获得的特性滤波器矩阵;上角标*表示取矩阵元素的共辄，上角标-1表示相应算 子的逆过程。4. 如权利要求3所述的基于1<〇Η则化的偏置相位中心天线成像方法，其特征在于，所述 步骤D具体包括： 将成像雷达的回波数据组成的回波数据矩阵、DPCA数据处理算子引入到由1<〇Η则化模 型推得的阈值函数中，进行如下迭代阈值运算：其中，Tku,q为iq正则化模型的阈值函数;YiS由成像雷达第i个接收通道采样获取的回 波数据组成的回波数据矩阵，γ为由成像雷达采样获取的回波数据组成的回波数据矩阵，且 Y= U iYi; κ为正则化参数;μ为收敛性控制参数;令χ(()) = 〇;设最大迭代次数为κ; 当迭代次数达到Κ或满足收敛精度，便可得到观测场景的后向散射系数矩阵的估计X。5. 如权利要求4所述的基于1<〇Η则化的偏置相位中心天线成像方法，其特征在于，所述 收敛精度具体包括： P(k+1)= I |x(k+i)_x(k)| |2，| |x(k)| |2 (9) p(k+1)小于10-6;其中，Μ · ||2表示取矩阵的二范数运算；即当P(k+1)小于1〇_6时结束迭代 阈值运算，得到观测场景的后向散射系数矩阵的估计X。6. 如权利要求2所述的基于1<〇Η则化的偏置相位中心天线成像方法，其特征在于，所述 子步骤A1具体包括：其中，HKf)表示成像雷达第i个接收通道的特性滤波器在多普勒频域的脉冲响应函数， Δ ^表示成像雷达第i个接收天线与成像雷达发射天线之间的距离；i e [ 1，2，. . .，N]，N为 成像雷达的接收通道数;λ表示成像雷达信号的载波波长;R〇表示成像雷达与观测目标间的 最短斜距;v表示成像雷达运载平台的运动速度，f为雷达信号的频率。7. 如权利要求2所述的基于1<〇Η则化的偏置相位中心天线成像方法，其特征在于，所述 子步骤A2具体包括： 所述特性矩阵为：其中，HKf)表示成像雷达第i个接收通道的特性滤波器在多普勒频域的脉冲响应函数， i = l，2,. . .，N，N为成像雷达的接收通道数;PRF为脉冲重复频率。8. 如权利要求2所述的基于1<〇Η则化的偏置相位中心天线成像方法，其特征在于，所述 子步骤A3具体包括： 所述重建矩阵为：其中，P(f)为成像雷达的重建矩阵，Pdf)为成像雷达第i个接收通道的重建滤波器，所 述重建滤波器由N个带通滤波器?^(〇(」=1，-_，《组成，每个带通滤波器?^(〇的通带范9. 如权利要求1所述的基于lq正则化的偏置相位中心天线成像方法，其特征在于： 所述步骤B具体包括:选择h正则化模型以及对应的阈值函数； 所述步骤D中，利用1:正则化模型的阈值函数进行迭代阈值运算。10. 如权利要求1所述的基于lq正则化的偏置相位中心天线成像方法，其特征在于： 所述步骤B具体包括:选择1 1/2正则化模型以及对应的阈值函数； 所述步骤D中，利用11/2正则化模型的阈值函数进行迭代阈值运算。
【文档编号】G01S13/90GK105929397SQ201610202747
【公开日】2016年9月7日
【申请日】2016年4月1日
【发明人】吴戎, 吴一戎, 全相印, 张冰尘, 张柘
【申请人】中国科学院电子学研究所