CN112684444B - 一种基于天线方向图合成抑制距离模糊的方法及装置 - Google Patents
一种基于天线方向图合成抑制距离模糊的方法及装置 Download PDFInfo
- Publication number
- CN112684444B CN112684444B CN202011378766.5A CN202011378766A CN112684444B CN 112684444 B CN112684444 B CN 112684444B CN 202011378766 A CN202011378766 A CN 202011378766A CN 112684444 B CN112684444 B CN 112684444B
- Authority
- CN
- China
- Prior art keywords
- determining
- target wave
- signal
- weight
- parameters
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 24
- 230000015572 biosynthetic process Effects 0.000 title claims abstract description 16
- 238000003786 synthesis reaction Methods 0.000 title claims abstract description 16
- 239000011159 matrix material Substances 0.000 claims abstract description 67
- 238000005457 optimization Methods 0.000 claims abstract description 38
- 238000005259 measurement Methods 0.000 claims abstract description 34
- 230000008054 signal transmission Effects 0.000 claims abstract description 13
- 239000011295 pitch Substances 0.000 claims description 12
- 239000013598 vector Substances 0.000 claims description 12
- 238000005070 sampling Methods 0.000 claims description 6
- 238000011022 operating instruction Methods 0.000 claims description 4
- 230000000007 visual effect Effects 0.000 claims 1
- 230000002194 synthesizing effect Effects 0.000 abstract description 3
- 230000002401 inhibitory effect Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 11
- 238000004364 calculation method Methods 0.000 description 3
- 230000001629 suppression Effects 0.000 description 3
- 230000033228 biological regulation Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000002592 echocardiography Methods 0.000 description 2
- 238000009434 installation Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000004580 weight loss Effects 0.000 description 1
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本申请实施例公开了一种基于天线方向图合成抑制距离模糊的方法,包括:获取目标波位工作指令;根据所述目标波位工作指令确定目标波位信号参数;初始化测量参数,其中,所述测量参数包括所述目标波位信号参数、SAR系统参数和自定义参数;根据所述测量参数确定与所述目标波位信号对应的信号发射权值;根据所述信号发射权值和所述测量参数确定所述目标波位信号的信号相关矩阵及距离模糊相关矩阵;根据所述信号相关矩阵和所述距离模糊相关矩阵确定初始阵列权值;根据所述目标波位工作指令确定凸优化问题;根据所述初始阵列权值和所述凸优化问题确定对应所述目标波位信号的接收权值。本申请实施例还公开了一种基于天线方向图合成抑制距离模糊的装置。
Description
技术领域
本申请涉及合成孔径雷达技术领域,尤其涉及一种基于天线方向图合成抑制距离模糊的方法及装置。
背景技术
合成孔径雷达(Synthetic Aperture Radar,SAR)是一种主动式微波成像设备,相比光学雷达具有更强的穿透性,可以实现全天时、全天候对地观测,在遥感领域有着广泛的应用。SAR通过处理宽带脉冲信号及方位向多普勒信号取得高分辨率图像,由于天线俯仰向方向图不可避免的存在旁瓣,在接收回波时同样会接收到处于测绘带外部的回波,影响最终的图像质量,这种干扰被称为距离模糊。为了提升SAR的性能,必须对距离模糊进行抑制。
基于天线方向图合成的距离模糊抑制方法根据距离模糊分布寻找在对应位置存在零陷的天线方向图,从而降低接收的距离模糊能量水平。这种方法不提升系统的复杂性,可以在保证天线增益、波束宽度、旁瓣水平满足要求的情况下有效抑制距离模糊。
判断这类距离模糊抑制方法性能的标准是计算资源占用量、算法速度、对方向图指标的控制能力,同时兼顾这些要求十分困难,当前已有的算法均不能兼顾这些要求,因此无法满足星载SAR的任务要求。
发明内容
本申请实施例提供一种基于天线方向图合成抑制距离模糊的方法,所述方法包括:
获取目标波位工作指令;
根据所述目标波位工作指令确定目标波位信号参数;
初始化测量参数,其中,所述测量参数包括所述目标波位信号参数、合成孔径雷达SAR系统参数和自定义参数;
根据所述测量参数确定与所述目标波位信号对应的信号发射权值;
根据所述信号发射权值和所述测量参数确定所述目标波位信号的信号相关矩阵及距离模糊相关矩阵;
根据所述信号相关矩阵和所述距离模糊相关矩阵确定初始阵列权值;
根据所述目标波位工作指令确定凸优化问题;
根据所述初始阵列权值和所述凸优化问题确定对应所述目标波位信号的接收权值。
上述实施例中,所述根据所述信号发射权值和所述测量参数确定所述目标波位信号的信号相关矩阵及距离模糊相关矩阵,包括:
在参考坐标系中,确定指定波位下视角的变化范围[αmin,αmax],在该范围内对下视角进行N点均匀采样获得序列α10,…,αN0;
根据公式获得采样的下视角序列对应的斜距R10,…,RN0,其中Re为地球半径;
寻找斜距对应的所有模糊斜距并编号,对于序列第i个元素Ri0,记其对应的模糊斜距为Rij(j=1,…,M),下视角为αij(j=1,…,M);
对远场球坐标系(θ,φ)进行如下变换,k0为波数
u=k0 sinθcosφ
v=k0 sinθsinφ
并计算αij(j=0,…,M)对应的(uij,vij);
引入向量
并记
其中Ek(u,v)是第k个阵元的远场方向图,(xk,yk)是第k个阵元在参考坐标系中的坐标;
分别计算信号矩阵与模糊矩阵
其中,
Cr(i,j)=V(i,j)V(i,j)H
Gt(i,j)为发射方向图在αij方向的增益,ηij为Rij对应的入射角,为对应ηij的反射率,p1,…,p6为经验参数。
上述实施例中,所述根据所述信号相关矩阵和所述距离模糊相关矩阵确定初始阵列权值,包括:
根据下式计算对角加载后的相关矩阵Dr
Dr=Ar+σ2I
其中σ0为非负的经验常数;
所述初始阵列权值w0为两矩阵(S,Dr)的最大广义特征值对应的特征向量。
上述实施例中,所述根据所述目标波位工作指令确定凸优化问题,包括:
根据实际需求确定可容忍的最大动态范围比ρ;
从α10,…,αN0中选出需要被特别抑制的角度,记为α1,…,αp;
At为对应α1,…,αp的相关矩阵,nt为α1,…,αp在原序列中对应的序号数;
所述凸优化问题可以写作
其中,wl为预先给定的权值,μt为用于加权的常数。
上述实施例中,所述根据所述初始阵列权值和所述凸优化问题确定对应所述目标波位信号的接收权值,包括:
S1,初始化wl=w0,设定最大容忍差值g,用于更新迭代的参数α、δ;
S2,根据所述凸优化问题获得w;
S3,更新权值:wl=(0.5-α)wl+(0.5+α)w;
S4,更新参数α=α-δ;
S5,如果|wl-w|≥g|wl|且α≥0,返回S2;
S6,得到w;
此时获得的w即为满足条件的所述接收权值,对应天线方向图可由公式F(u,v)=VHw获得。
本申请实施例提供一种基于天线方向图合成抑制距离模糊的装置,所述装置包括:
控制单元,用于获取目标波位工作指令;根据所述目标波位工作指令确定目标波位信号参数;初始化测量参数,其中,所述测量参数包括所述目标波位信号参数、合成孔径雷达SAR系统参数和自定义参数;根据所述测量参数确定与所述目标波位信号对应的信号发射权值;
运算单元,用于根据所述测量参数确定与所述目标波位信号对应的信号发射权值;根据所述信号发射权值和所述测量参数确定所述目标波位信号的信号相关矩阵及距离模糊相关矩阵;根据所述信号相关矩阵和所述距离模糊相关矩阵确定初始阵列权值;根据所述目标波位工作指令确定凸优化问题;根据所述初始阵列权值和所述凸优化问题确定对应所述目标波位信号的接收权值。
上述实施例中,所述运算单元,具体用于:
在参考坐标系中,确定指定波位下视角的变化范围[αmin,αmax],在该范围内对下视角进行N点均匀采样获得序列α10,…,αN0;
根据公式获得采样的下视角序列对应的斜距R10,…,RN0,其中Re为地球半径;
寻找斜距对应的所有模糊斜距并编号,对于序列第i个元素Ri0,记其对应的模糊斜距为Rij(j=1,…,M),下视角为αij(j=1,…,M);
对远场球坐标系(θ,φ)进行如下变换,k0为波数
u=k0sinθcosφ
v=k0sinθsinφ
并计算αij(j=0,…,M)对应的(uij,vij);
引入向量
并记
其中Ek(u,v)是第k个阵元的远场方向图,(xk,yk)是第k个阵元在参考坐标系中的坐标;
分别计算信号矩阵与模糊矩阵
其中,
Cr(i,j)=V(i,j)V(i,j)H
Gt(i,j)为发射方向图在αij方向的增益,ηij为Rij对应的入射角,为对应ηij的反射率,p1,…,p6为经验参数。
上述实施例中,所述运算单元,还具体用于:
根据下式计算对角加载后的相关矩阵Dr
Dr=Ar+σ2I
其中σ0为非负的经验常数;
所述初始阵列权值w0为两矩阵(S,Dr)的最大广义特征值对应的特征向量。
上述实施例中,所述运算单元,还具体用于:
根据实际需求确定可容忍的最大动态范围比ρ;
从α10,…,αN0中选出需要被特别抑制的角度,记为α1,…,αp;
At为对应α1,…,αp的相关矩阵,nt为α1,…,αp在原序列中对应的序号数;
所述凸优化问题可以写作
其中,wl为预先给定的权值,μt为用于加权的常数。
上述实施例中,所述运算单元,还具体用于:
S1,初始化wl=w0,设定最大容忍差值g,用于更新迭代的参数α、δ;
S2,根据所述凸优化问题获得w;
S3,更新权值:wl=(0.5-α)wl+(0.5+α)w;
S4,更新参数α=α-δ;
S5,如果|wl-w|≥g|wl|且α≥0,返回S2;
S6,得到w;
此时获得的w即为满足条件的所述接收权值,对应天线方向图可由公式F(u,v)=VHw获得。
本申请实施例提供的一种基于天线方向图合成抑制距离模糊的方法相对已有技术具备以下优点:其一,本发明在未使用全局优化算法的情况下将待解的非凸问题松弛为易解的凸优化问题,并通过更新迭代保证了结果的准确性,减少了占用计算资源并提升了运算速度;其二,本发明使用了一种有效快捷的迭代初值确定方法,保证了算法的收敛性并进一步的提升了算法的运行速度;其三,通过对约束条件的调整,本发明的优化结果的指标可以被灵活的控制。
附图说明
附图以示例而非限制的方式大体示出了本文中所讨论的各个实施例。
图1为本申请实施例一种基于天线方向图合成抑制距离模糊的方法流程示意图;
图2为本申请一实施例的发射天线方向示意图;
图3为本申请一实施例的接收方向图与发射方向图相同时的天线双程方向图以及距离模糊区域分布示意图;
图4为本申请一实施例的距离模糊水平示意图;
图5为本申请一实施例的初始权值w0与优化权值w的对比示意图;
图6为本申请一实施例的w0及w对应的方向图;
图7为本申请一实施例的优化后的天线双程方向图以及距离模糊区域分布示意图;
图8为本申请一实施例的w0及w对应的距离模糊水平示意图;
图9为本申请实例二与实例一优化模糊特性的对比示意图;
图10为本申请实施例一种基于天线方向图合成抑制距离模糊的装置结构示意图。
具体实施方式
为了能够更加详尽地了解本申请实施例的特点与技术内容,下面结合附图对本申请实施例的实现进行详细阐述,所附附图仅供参考说明之用,并非用来限定本申请实施例。
本申请的主要目的在于提出了一种基于天线方向图合成抑制距离模糊的方法,该方法占用计算资源及算法速度以及方向图指标的控制能力整体优于已被发表的算法,有效的抑制了距离模糊。
为达到上述目的,本申请的技术方案是这样实现的:
图1为本申请实施例一种基于天线方向图合成抑制距离模糊的方法流程示意图,如图1所示,本申请实施例的基于天线方向图合成抑制距离模糊的方法包括以下步骤:
步骤101,获取目标波位工作指令。
具体地,在一些实施例中,各个波位的发射和接收是一个整体的过程,即一次工作指令至少包括指定波位信号的发射和接收。
步骤102,根据所述目标波位工作指令确定目标波位信号参数。
在一些实施例中,目标波位工作指令可以通过在SAR系统中输入波位工作指令而生成。
步骤103,初始化测量参数,其中,所述测量参数包括所述目标波位信号参数、合成孔径雷达SAR系统参数和自定义参数。
在一些实施例中,SAR系统参数和自定义参数不随选取的目标波位变化;SAR系统参数和自定义参数可以是预先或根据用户要求设定。
步骤104,根据所述测量参数确定与所述目标波位信号对应的信号发射权值。
步骤105,根据所述信号发射权值和所述测量参数确定所述目标波位信号的信号相关矩阵及距离模糊相关矩阵。
在一些实施例中,所述根据所述信号发射权值和所述测量参数确定所述目标波位信号的信号相关矩阵及距离模糊相关矩阵,包括:
在参考坐标系中,确定指定波位下视角的变化范围[αmin,αmax],在该范围内对下视角进行N点均匀采样获得序列α10,…,αN0;
根据公式获得采样的下视角序列对应的斜距R10,…,RN0,其中Re为地球半径;
寻找斜距对应的所有模糊斜距并编号,对于序列第i个元素Ri0,记其对应的模糊斜距为Rij(j=1,…,M),下视角为αij(j=1,…,M);
对远场球坐标系(θ,φ)进行如下变换,k0为波数
u=k0sinθcosφ
v=k0sinθsinφ
并计算αij(j=0,…,M)对应的(uij,vij);
引入向量
并记
其中Ek(u,v)是第k个阵元的远场方向图,(xk,yk)是第k个阵元在参考坐标系中的坐标;
分别计算信号矩阵与模糊矩阵
其中,
Cr(i,j)=V(i,j)V(i,j)H
Gt(i,j)为发射方向图在αij方向的增益,ηij为Rij对应的入射角,为对应ηij的反射率,p1,…,p6为经验参数。
步骤106,根据所述信号相关矩阵和所述距离模糊相关矩阵确定初始阵列权值。
在一些实施例中,所述根据所述信号相关矩阵和所述距离模糊相关矩阵确定初始阵列权值,包括:
根据下式计算对角加载后的相关矩阵Dr
Dr=Ar+σ2I
其中σ0为非负的经验常数;
所述初始阵列权值w0为两矩阵(S,Dr)的最大广义特征值对应的特征向量。
步骤107,根据所述目标波位工作指令确定凸优化问题。
在一些实施例中,所述根据所述目标波位工作指令确定凸优化问题,包括:
根据实际需求确定可容忍的最大动态范围比ρ;
从α10,…,αN0中选出需要被特别抑制的角度,记为α1,…,αp;
At为对应α1,…,αp的相关矩阵,nt为α1,…,αp在原序列中对应的序号数;
所述凸优化问题可以写作
其中,wl为预先给定的权值,μt为用于加权的常数。
步骤108,根据所述初始阵列权值和所述凸优化问题确定对应所述目标波位信号的接收权值。
在一些实施例中,所述根据所述初始阵列权值和所述凸优化问题确定对应所述目标波位信号的接收权值,包括:
S1,初始化wl=w0,设定最大容忍差值g,用于更新迭代的参数α、δ;
S2,根据所述凸优化问题获得w;
S3,更新权值:wl=(0.5-α)wl+(0.5+α)w;
S4,更新参数α=α-δ;
S5,如果|wl-w|≥g|wl|且α≥0,返回S2;
S6,得到w;
此时获得的w即为满足条件的所述接收权值,对应天线方向图可由公式F(u,v)=VHw获得。
下面结合实施例对本申请再作进一步详细的描述。
实施例一
本实施例通过具体的抑制效果对本申请的技术方案进行进一步的描述与说明。
本实施例中,设定卫星天线高程向阵元数为22,工作于L波段。选定波位PRF为3366Hz,天线安装角30.5°,观测波束近端与远端视角分别为25.1°与27.3°,轨道高度607km。
图2展示了发射天线方向图,权值仅附加了使波束改变指向的相位。
图3展示了接收方向图与发射方向图相同时的天线双程方向图以及距离模糊区域分布。
图4展示了此时的距离模糊水平,可容忍的最大值一般为-20dB,当前距离模糊水平远远超过可容忍值。
根据上述实施例中步骤101至106:
将本实施例的波位参数及发射天线方向图代入获得信号相关矩阵及距离模糊相关矩阵;
设定σ2=0.3,计算Dr,根据(S,Dr)求得初始权值w0;
凸优化问题构建如下,设定p=1,μ1=0,可容忍的比值为
设定g=0.02,α=0.2,δ=0.01并进行迭代计算获得w。
图5展示了初始权值w0与优化权值w的对比。归一化后w最小元素值被控制在-7dB以内,调控目标实现。
图6展示了根据公式F(u,v)=VHw得出的w0及w对应的方向图。二者相似度很高,证明了本申请确定初始权值方法的有效性。
图7展示了优化后的天线双程方向图以及距离模糊区域分布,同图3对比可发现在距离模糊区域出现了明显的零陷。
图8展示了此时w0及w对应的距离模糊水平。首先,同图4相比,RASR被显著抑制,证明了本申请抑制距离模糊的有效性。其次,优化后权值的抑制效果相对初始权值损失极小,证明了本申请迭代策略的有效性。
实施例二
本实施例通过重点抑制特定视角的距离模糊进一步描述说明本申请的技术方案。同样设定卫星天线高程向阵元数为22,工作于L波段。选定波位PRF为3366Hz,天线安装角30.5°,观测波束近端与远端视角分别为25.1°与27.3°,轨道高度607km,并重点抑制远端视角处的距离模糊。
在步骤101至104与实施例一相同的情况下,根据步骤104,凸优化问题构建如下,设定p=1,μ1=0.1,A1为对应最远端视角的相关矩阵,
图9展示了实例二与实例一优化模糊特性的对比,远端视角处的模糊被有效抑制,本申请对结果的灵活调控能力被进一步证明。
从上面的描述中可以看出,采用本申请提供的方法,能够在兼顾计算资源占用量、算法速度、对关键指标的控制能力的情况下有效抑制距离模糊。
图10为本申请实施例一种基于天线方向图合成抑制距离模糊的装置结构示意图,如图10所示,为实现上述基于天线方向图合成抑制距离模糊的方法,本申请实施例提供一种基于天线方向图合成抑制距离模糊的装置,所述装置包括:控制单元21和运算单元22;其中,
控制单元21,用于获取目标波位工作指令;根据所述目标波位工作指令确定目标波位信号参数;初始化测量参数,其中,所述测量参数包括所述目标波位信号参数、合成孔径雷达SAR系统参数和自定义参数;根据所述测量参数确定与所述目标波位信号对应的信号发射权值。
运算单元22,用于根据所述测量参数确定与所述目标波位信号对应的信号发射权值;根据所述信号发射权值和所述测量参数确定所述目标波位信号的信号相关矩阵及距离模糊相关矩阵;根据所述信号相关矩阵和所述距离模糊相关矩阵确定初始阵列权值;根据所述目标波位工作指令确定凸优化问题;根据所述初始阵列权值和所述凸优化问题确定对应所述目标波位信号的接收权值。
上述实施例中,所述运算单元22,具体用于:
在参考坐标系中,确定指定波位下视角的变化范围[αmin,αmax],在该范围内对下视角进行N点均匀采样获得序列α10,…,αN0;
根据公式获得采样的下视角序列对应的斜距R10,…,RN0,其中Re为地球半径;
寻找斜距对应的所有模糊斜距并编号,对于序列第i个元素Ri0,记其对应的模糊斜距为Rij(j=1,…,M),下视角为αij(j=1,…,M);
对远场球坐标系(θ,φ)进行如下变换,k0为波数
u=k0 sinθcosφ
v=k0 sinθsinφ
并计算αij(j=0,…,M)对应的(uij,vij);
引入向量
并记
其中Ek(u,v)是第k个阵元的远场方向图,(xk,yk)是第k个阵元在参考坐标系中的坐标;
分别计算信号矩阵与模糊矩阵
其中,
Cr(i,j)=V(i,j)V(i,j)H
Gt(i,j)为发射方向图在αij方向的增益,ηij为Rij对应的入射角,为对应ηij的反射率,p1,…,p6为经验参数。
上述实施例中,所述运算单元22,还具体用于:
根据下式计算对角加载后的相关矩阵Dr
Dr=Ar+σ2I
其中σ0为非负的经验常数;
所述初始阵列权值w0为两矩阵(S,Dr)的最大广义特征值对应的特征向量。
上述实施例中,所述运算单元,还具体用于:
根据实际需求确定可容忍的最大动态范围比ρ;
从α10,…,αN0中选出需要被特别抑制的角度,记为α1,…,αp;
At为对应α1,…,αp的相关矩阵,nt为α1,…,αp在原序列中对应的序号数;
所述凸优化问题可以写作
其中,wl为预先给定的权值,μt为用于加权的常数。
上述实施例中,所述运算单元22,还具体用于:
S1,初始化wl=w0,设定最大容忍差值g,用于更新迭代的参数α、δ;
S2,根据所述凸优化问题获得w;
S3,更新权值:wl=(0.5-α)wl+(0.5+α)w;
S4,更新参数α=α-δ;
S5,如果|wl-w|≥g|wl|且α≥0,返回S2;
S6,得到w;
此时获得的w即为满足条件的所述接收权值,对应天线方向图可由公式F(u,v)=VHw获得。
本申请实施例所记载的技术方案之间,在不冲突的情况下,可以任意组合。
以上所述,仅为本申请的具体实施方式,但本申请的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本申请揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本申请的保护范围之内。因此,本申请的保护范围应以所述权利要求的保护范围为准。
Claims (10)
1.一种基于天线方向图合成抑制距离模糊的方法,其特征在于,所述方法包括:
获取目标波位工作指令;
根据所述目标波位工作指令确定目标波位信号参数;
初始化测量参数,其中,所述测量参数包括所述目标波位信号参数、合成孔径雷达SAR系统参数和自定义参数;
根据所述测量参数确定与所述目标波位信号对应的信号发射权值;
根据所述信号发射权值和所述测量参数确定所述目标波位信号的信号相关矩阵及距离模糊相关矩阵;
根据所述信号相关矩阵和所述距离模糊相关矩阵确定初始阵列权值;
根据所述目标波位工作指令确定凸优化问题;
根据所述初始阵列权值和所述凸优化问题确定对应所述目标波位信号的接收权值。
2.根据权利要求1所述的方法,其特征在于,所述根据所述信号发射权值和所述测量参数确定所述目标波位信号的信号相关矩阵及距离模糊相关矩阵,包括:
在参考坐标系中,确定指定波位下视角的变化范围[αmin,αmax],在该范围内对下视角进行N点均匀采样获得序列Δ10,…,ΔN0;
根据公式获得采样的下视角序列对应的斜距R10,…,RN0,其中Re为地球半径,H为轨道高度;
寻找斜距对应的所有模糊斜距并编号,对于序列第i个元素Ri0,记其对应的模糊斜距为Rij(j=1,…,M),下视角为Δij(j=1,…,M);
对远场球坐标系(θ,φ)进行如下变换,k0为波数
u=k0sinθcosφ
v=k0sinθsinφ
并计算Δij(j=0,…,M)对应的(uij,vij);
引入向量
并记
其中Ek(u,v)是第k个阵元的远场方向图,(xk,yk)是第k个阵元在参考坐标系中的坐标;
分别计算信号矩阵与模糊矩阵
其中,
Cr(i,j)=V(i,j)V(i,j)H
Gt(i,j)为发射方向图在Δij方向的增益,ηij为Rij对应的入射角,为对应ηij的反射率,p1,…,p6为经验参数。
3.根据权利要求2所述的方法,其特征在于,所述根据所述信号相关矩阵和所述距离模糊相关矩阵确定初始阵列权值,包括:
根据下式计算对角加载后的相关矩阵Dr
Dr=Ar+σ2I
其中σ0为非负的经验常数;
所述初始阵列权值w0为两矩阵(S,Dr)的最大广义特征值对应的特征向量。
4.根据权利要求3所述的方法,其特征在于,所述根据所述目标波位工作指令确定凸优化问题,包括:
根据实际需求确定可容忍的最大动态范围比ρ;
从α10,…,αN0中选出需要被特别抑制的角度,记为Δ1,…,Δp;
At为对应Δ1,…,αp的相关矩阵,nt为Δ1,…,αp在原序列中对应的序号数;
所述凸优化问题可以写作
minimizew g
wH(Dr+μtAt)w≤g,t=1…p
ξ≥0
其中,wl为预先给定的权值,μt为用于加权的常数。
5.根据权利要求4所述的方法,其特征在于,所述根据所述初始阵列权值和所述凸优化问题确定对应所述目标波位信号的接收权值,包括:
S1,初始化wl=w0,设定最大容忍差值g,用于更新迭代的参数α、δ;
S2,根据所述凸优化问题获得w;
S3,更新权值:wl=(0.5-α)wl+(0.5+α)w;
S4,更新参数α=α-δ;
S5,如果|wl-w|≥g|wl|且α≥0,返回S2;
S6,得到w;
此时获得的w即为满足条件的所述接收权值,对应天线方向图可由公式F(u,v)=VHw获得。
6.一种基于天线方向图合成抑制距离模糊的装置,其特征在于,所述装置包括:
控制单元,用于获取目标波位工作指令;根据所述目标波位工作指令确定目标波位信号参数;初始化测量参数,其中,所述测量参数包括所述目标波位信号参数、合成孔径雷达SAR系统参数和自定义参数;根据所述测量参数确定与所述目标波位信号对应的信号发射权值;
运算单元,用于根据所述测量参数确定与所述目标波位信号对应的信号发射权值;根据所述信号发射权值和所述测量参数确定所述目标波位信号的信号相关矩阵及距离模糊相关矩阵;根据所述信号相关矩阵和所述距离模糊相关矩阵确定初始阵列权值;根据所述目标波位工作指令确定凸优化问题;根据所述初始阵列权值和所述凸优化问题确定对应所述目标波位信号的接收权值。
7.根据权利要求6所述的装置,其特征在于,所述运算单元,具体用于:
在参考坐标系中,确定指定波位下视角的变化范围[αmin,αmax],在该范围内对下视角进行N点均匀采样获得序列α10,…,αN0;
根据公式获得采样的下视角序列对应的斜距R10,…,RN0,其中Re为地球半径,H为轨道高度;
寻找斜距对应的所有模糊斜距并编号,对于序列第i个元素Ri0,记其对应的模糊斜距为Rij(j=1,…,M),下视角为αij(j=1,…,M);
对远场球坐标系(θ,φ)进行如下变换,k0为波数
u=k0sinθcosφ
v=k0sinθsinφ
并计算αij(j=0,…,M)对应的(uij,vij);
引入向量
并记
其中Ek(u,v)是第k个阵元的远场方向图,(xk,yk)是第k个阵元在参考坐标系中的坐标;
分别计算信号矩阵与模糊矩阵
其中,
Cr(i,j)=V(i,j)V(i,j)H
Gt(i,j)为发射方向图在αij方向的增益,ηij为Rij对应的入射角,为对应ηij的反射率,p1,…,p6为经验参数。
8.根据权利要求7所述的装置,其特征在于,所述运算单元,还具体用于:
根据下式计算对角加载后的相关矩阵Dr
Dr=Ar+σ2I
其中σ0为非负的经验常数;
所述初始阵列权值w0为两矩阵(S,Dr)的最大广义特征值对应的特征向量。
9.根据权利要求8所述的装置,其特征在于,所述运算单元,还具体用于:
根据实际需求确定可容忍的最大动态范围比ρ;
从α10,…,αN0中选出需要被特别抑制的角度,记为α1,…,ap;
At为对应a1,…,ap的相关矩阵,nt为a1,…,ap在原序列中对应的序号数;
所述凸优化问题可以写作
minimizew g
wH(Dr+μtAt)w≤g,t=1…p
ξ≥0
其中,wl为预先给定的权值,μt为用于加权的常数。
10.根据权利要求9所述的装置,其特征在于,所述运算单元,还具体用于:
S1,初始化wl=w0,设定最大容忍差值g,用于更新迭代的参数a、δ;
S2,根据所述凸优化问题获得w;
S3,更新权值:wl=(0.5-a)wl+(0.5+α)w;
S4,更新参数α=α-δ;
S5,如果|wl-w|≥g|wl|且α≥0,返回S2;
S6,得到w;
此时获得的w即为满足条件的所述接收权值,对应天线方向图可由公式F(u,v)=VHw获得。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011378766.5A CN112684444B (zh) | 2020-11-30 | 2020-11-30 | 一种基于天线方向图合成抑制距离模糊的方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011378766.5A CN112684444B (zh) | 2020-11-30 | 2020-11-30 | 一种基于天线方向图合成抑制距离模糊的方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112684444A CN112684444A (zh) | 2021-04-20 |
CN112684444B true CN112684444B (zh) | 2023-09-22 |
Family
ID=75447072
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011378766.5A Active CN112684444B (zh) | 2020-11-30 | 2020-11-30 | 一种基于天线方向图合成抑制距离模糊的方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112684444B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113409223B (zh) * | 2021-06-30 | 2023-06-23 | 中国科学院空天信息创新研究院 | 基于粒子群算法与投影方法的距离模糊抑制方法和装置 |
CN113466802B (zh) * | 2021-06-30 | 2023-06-09 | 中国科学院空天信息创新研究院 | 基于粒子群算法与最优波束合成方法的距离模糊抑制方法 |
CN113376635B (zh) * | 2021-06-30 | 2023-06-09 | 中国科学院空天信息创新研究院 | 一种基于谱分解的距离模糊抑制方法、装置及存储介质 |
CN116500550B (zh) * | 2023-06-27 | 2023-08-25 | 中国科学院空天信息创新研究院 | 一种星载sar距离模糊抑制方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102393514A (zh) * | 2011-10-25 | 2012-03-28 | 北京航空航天大学 | 一种合成孔径雷达卫星的自适应波位设计方法和系统 |
CN102736073A (zh) * | 2012-06-19 | 2012-10-17 | 北京航空航天大学 | 一种通用模式下星载sar距离向模糊度的计算方法 |
CN106886656A (zh) * | 2017-03-15 | 2017-06-23 | 南京航空航天大学 | 一种基于改进mopso和凸优化算法的三维阵列天线方向图旁瓣抑制方法 |
CN107561535A (zh) * | 2017-09-07 | 2018-01-09 | 中国科学院电子学研究所 | 一种合成孔径雷达距离模糊抑制方法和装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012181052A (ja) * | 2011-02-28 | 2012-09-20 | Toshiba Corp | 相関抑圧フィルタ、ウェイト算出方法、ウェイト算出装置、アダプティブアレーアンテナ及びレーダ装置 |
-
2020
- 2020-11-30 CN CN202011378766.5A patent/CN112684444B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102393514A (zh) * | 2011-10-25 | 2012-03-28 | 北京航空航天大学 | 一种合成孔径雷达卫星的自适应波位设计方法和系统 |
CN102736073A (zh) * | 2012-06-19 | 2012-10-17 | 北京航空航天大学 | 一种通用模式下星载sar距离向模糊度的计算方法 |
CN106886656A (zh) * | 2017-03-15 | 2017-06-23 | 南京航空航天大学 | 一种基于改进mopso和凸优化算法的三维阵列天线方向图旁瓣抑制方法 |
CN107561535A (zh) * | 2017-09-07 | 2018-01-09 | 中国科学院电子学研究所 | 一种合成孔径雷达距离模糊抑制方法和装置 |
Non-Patent Citations (1)
Title |
---|
基于模糊粒子群算法的有向传感器网络路径覆盖策略;张聚伟 等;《模式识别与人工智能PR&AI》;第30卷(第2期);183-191 * |
Also Published As
Publication number | Publication date |
---|---|
CN112684444A (zh) | 2021-04-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112684444B (zh) | 一种基于天线方向图合成抑制距离模糊的方法及装置 | |
US6661366B2 (en) | Adaptive digital sub-array beamforming and deterministic sum and difference beamforming, with jamming cancellation and monopulse ratio preservation | |
US6697009B2 (en) | Adaptive digital beamforming architecture for target detection and angle estimation in multiple mainlobe and sidelobe jamming | |
US6531976B1 (en) | Adaptive digital beamforming radar technique for creating high resolution range profile for target in motion in the presence of jamming | |
US6720910B2 (en) | Pri-staggered post-doppler adaptive monopulse processing for detection and location of a moving target in ground clutter | |
Kang et al. | Efficient synthesis of antenna pattern using improved PSO for spaceborne SAR performance and imaging in presence of element failure | |
CN108020817A (zh) | 基于配准的机载前视阵雷达杂波抑制方法 | |
CN104076337B (zh) | 基于阵元幅相误差校正的机载雷达杂波抑制方法 | |
CN109375179A (zh) | 一种基于稀疏表示的机载雷达近程强杂波抑制方法 | |
JP2005512435A (ja) | 自動較正される縮小ランク適応プロセッサのためのシステムおよび方法 | |
CN111007503A (zh) | 基于频率谱精确定位的运动目标聚焦和定位方法及系统 | |
CN111352077B (zh) | 基于多比例分式规划的低截获频控阵mimo雷达系统的设计方法 | |
US20030184473A1 (en) | Adaptive digital sub-array beamforming and deterministic sum and difference beamforming, with jamming cancellation and monopulse ratio preservation | |
CN106291547B (zh) | 基于天线方向图辅助的多普勒模糊分量自适应抑制方法 | |
CN113466797A (zh) | 基于杂波脊匹配稀疏恢复的双基sar空-时杂波抑制方法 | |
CN105353386A (zh) | 一种利用惯性导航设备的导航接收机抗干扰的方法和装置 | |
CN116679306A (zh) | 一种基于系数修正的扩展孔径声纳成像方法 | |
CN113376635B (zh) | 一种基于谱分解的距离模糊抑制方法、装置及存储介质 | |
CN115436940A (zh) | 一种稀疏滑动聚束sar成像模式实现方法及装置 | |
CN115343683B (zh) | 基于等多普勒样本和等锥角样本结合的杂波抑制方法 | |
CN112904297B (zh) | 一种分维自适应单脉冲波束形成与角度估计方法 | |
CN118068272B (zh) | 最大化主瓣增益的稳健低副瓣发射波束形成方法 | |
Recchia et al. | Impact of the antenna stability on the Doppler Centroid frequency | |
Scheiber et al. | Comparison of digital beamforming techniques for enhanced ice sounding radar data processing | |
CN114779198B (zh) | 共形阵机载雷达空时杂波谱自适应补偿与杂波抑制方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |