CN113126087B - 一种星载干涉成像高度计天线 - Google Patents

一种星载干涉成像高度计天线 Download PDF

Info

Publication number
CN113126087B
CN113126087B CN202110260165.2A CN202110260165A CN113126087B CN 113126087 B CN113126087 B CN 113126087B CN 202110260165 A CN202110260165 A CN 202110260165A CN 113126087 B CN113126087 B CN 113126087B
Authority
CN
China
Prior art keywords
antenna
array
excitation
shaped
directional diagram
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
Application number
CN202110260165.2A
Other languages
English (en)
Other versions
CN113126087A (zh
Inventor
赵晓雯
张云华
董晓
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
National Space Science Center of CAS
Original Assignee
National Space Science Center of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by National Space Science Center of CAS filed Critical National Space Science Center of CAS
Priority to CN202110260165.2A priority Critical patent/CN113126087B/zh
Publication of CN113126087A publication Critical patent/CN113126087A/zh
Application granted granted Critical
Publication of CN113126087B publication Critical patent/CN113126087B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/882Radar or analogous systems specially adapted for specific applications for altimeters
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9021SAR image post-processing techniques
    • G01S13/9023SAR image post-processing techniques combined with interferometric techniques
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00

Abstract

本发明公开了一种星载干涉成像高度计天线,该天线为基于物理尺寸给定的赋形波束平面阵天线,包括若干个子阵,每个子阵包括若干个阵元,每个子阵中的阵元激励等幅同相,天线的距离向采用赋形波束设计,由入射角范围内的赋形曲线、副瓣电平、天线阵元以及幅相约束指标确定赋形波束及激励;天线的方位向设计采用子阵划分技术,通过子阵级幅度加权实现期望的泰勒加权方向图特性的笔形波束,根据可分离平面阵形成原理,由距离向和方位向的设计结果确定天线的激励参数。本发明的天线对天线距离向进行赋形波束综合,有利于提高干涉测高精度,对方位向进行子阵划分,不仅有利于提高天线带宽,而且简化整个天线系统结构,降低硬件实现成本和工程难度。

Description

一种星载干涉成像高度计天线
技术领域
本发明涉及无线通信领域,尤其涉及一种星载干涉成像高度计天线。
背景技术
星载干涉成像高度计是一种利用小角度入射、短基线干涉来实现对海陆表面高精度观测的微波遥感器,是海洋动力环境观测的主要遥感器之一(文献[1]:Y.Zhang,J.Jiang,X.Zhang,K.Xu,J.Yan,C.Jiang,et al,“Design and preliminary experimentof china imaging altimeter,”Microwave Remote Sensing of the Atmosphere andEnvironment III,pp.190-199,2003)。2016年,国际上第一个星载干涉成像高度计随天宫二号空间实验室发射升空,突破了传统高度计只能进行星下点沿飞行方向一维线观测、刈幅只有数公里的局限,除了可对海洋和陆地宽刈幅、高精度的观测外,还可对三维海面形态、海洋内波以及海面风速和风向进行观测和测量(参考文献[2]:X.Dong,Y.Zhang,andW.Zhai,“Design and algorithms of the Tiangong-2interferometric imaging radaaltimeter processor,”2017Progress In Electromagnetics Research Symposium-Spring(PIERS 2017),pp.3802-3803,2017;参考文献[3]:L.Ren,J.Yang,X.Dong,Y.Zhang,and Y.Jia,“Preliminary Evaluation and Correction of Sea Surface Height fromChinese Tiangong-2 Interferometric Imaging Radar Altimeter”,Remote Sensing,vol.12,no.15,2496,2020),对海洋科学、海洋环境观测和预报、全球气候变化研究具有重要的作用与意义。
天线作为星载干涉成像高度计重要组成部分,用来发射和接收电磁波,起着至关重要的作用。当前星载干涉成像高度计天线距离向采用笔形波束设计,然而,结合电磁场理论与星载数据分析可知,海面后向散射系数随入射角的增大而减小,可见回波能量在观测范围内分布不均匀,影响干涉测高精度,这不利于对中尺度海洋现象的高精度观测。
同时,为了获得更远的作用距离和更高的空间分辨率,星载干涉成像高度计天线必然也要朝着大型化的方向发展,但天线口径的增大,阵元数的增多必然导致馈电网络极其复杂,也影响着天线带宽,可见,当前采用的阵元级加权方法已不适用于大型天线,很难满足下一代星载干涉成像高度计宽幅宽、高分辨的系统需求。
发明内容
当前星载干涉成像高度计天线距离向观测不一致,馈电网络复杂,加工难度大,工作带宽窄,不利于对中尺度海洋现象高精度观测的问题。本发明的目的在于克服上述技术缺陷,提出了一种星载干涉成像高度计天线。
为了实现上述目的,本发明提出了一种星载干涉成像高度计天线,所述天线为基于物理尺寸给定的赋形波束平面阵天线,包括若干个子阵,每个子阵包括若干个阵元,每个子阵中的阵元激励等幅同相,所述天线的距离向采用赋形波束设计,由入射角范围内的赋形曲线、副瓣电平、天线阵元以及幅相约束指标确定赋形波束及激励;所述天线的方位向设计采用子阵划分技术,通过子阵级幅度加权实现期望的泰勒加权方向图特性的笔形波束,根据可分离平面阵形成原理,由距离向和方位向的设计结果确定天线的激励参数。
作为上述天线的一种改进,所述赋形波束设计包括以下步骤实现:
根据动态海面的观测数据,提取随入射角变化的后向散射系数的相对值;
由后向散射系数的相对值得到在入射角范围内的最大值和最小值,采用线性拟合得到入射角范围内的赋形方向图包络;
结合入射角范围内的赋形方向图包络、副瓣电平、天线阵元以及幅相约束指标,确定天线距离向的激励参数与赋形波束。
作为上述天线的一种改进,所述根据动态海面的观测数据,提取随入射角变化的后向散射系数的相对值;具体为:
对星载干涉成像高度计采集的动态海平面观测数据,进行合成孔径成像处理,得到图像功率值Pc
由天线距离向方向图G(θ)和观测斜距R,得到在入射角θ范围内,海面后向散射系数的相对值σ0为:
Figure GDA0003284940100000021
作为上述天线的一种改进,所述由后向散射系数的相对值得到在入射角范围内的最大值和最小值,采用线性拟合得到入射角范围内的赋形方向图包络;具体为:
由后向散射系数的相对值σ0得到入射角范围内的最大值G1和最小值G0
采用线性拟合得到入射角范围内的赋形方向图包络SB(θ)为:
Figure GDA0003284940100000022
其中,[θ01]为入射角的范围。
作为上述天线的一种改进,所述结合入射角范围内的赋形方向图包络、副瓣电平、天线阵元以及幅相约束指标,确定天线距离向的激励参数与赋形波束;具体为:
根据入射角范围内的赋形方向图包络SB(θ),期望的副瓣电平SLLd,阵元数M以及阵元均匀间隔d的指标要求,采用凸优化算法快速优化复激励,在无幅相约束下,得到期望赋形方向图Fd(θ),由此提取出主波束的两个零点θNULL1和θNULL2以及期望的赋形主波束方向图Fdn),θn为赋形主波束方向图第n个入射角,θNULL1≤θn≤θNULL2
将阵元激励和幅度作为优化变量,根据幅相约束指标,以对期望的主波束赋形方向图Fdn)和期望的副瓣电平SLLd的拟合为适应度函数,采用入侵杂草优化算法进行优化,计算出阵元激励,并由此得到天线距离向的激励参数与赋形波束。
作为上述天线的一种改进,所述适应度函数fit(θ)为:
Figure GDA0003284940100000031
其中,N为赋形主波束方向图的入射角个数,F(θn)为赋形方向图,SLL为副瓣电平:
SLL=max(|F(χt)|),χt∈[-90,θNULL1)∪(θNULL2,90]。
作为上述天线的一种改进,所述方位向采用子阵划分技术,通过子阵级幅度加权实现期望的泰勒加权方向图特性的笔形波束;具体为:
将天线方位向分布的P个均匀分布阵元划分成Q个子阵,并仅在Q个子阵端口进行激励加权,每个子阵中的阵元激励等幅同相,相应的阵因子F2(θ)表示为:
Figure GDA0003284940100000032
其中,dp表示第p个阵元的位置,wp表示第p个阵元的激励值,p=1,…,P,aq表示第q个子阵端口的激励权值,q=1,…,Q,j为虚部;
根据子阵加权矢量Wsub=[a1,a2,…,aQ]T,阵元激励矢量Wele=[w1,w2,…,wP]T,具体表示为:
Figure GDA0003284940100000033
其中,Nq是第q个子阵所包含的阵元数,用于表征相应子阵的尺寸,并且ΣNq=P;根据信号处理理论可知,阵元激励矢量Wele是一个可压缩矢量,可以用一个稀疏矢量来表示,由矢量Wele的元素特征可知,相应稀疏矢量的稀疏度为Q,即仅包含Q个非零元,稀疏矢量S具体表示为:
Figure GDA0003284940100000041
由此得出稀疏矢量S与阵元激励矢量Wele之间的数值关系为:
Wele=TS
其中,转换矩阵T为:
Figure GDA0003284940100000042
辐射特性Fy(θ)的矢量形式为:
Fy(θ)=ψ(θ)S
ψ(θ)=A(θ)T
其中,导向矢量A(θ)为:
Figure GDA0003284940100000043
根据方向图指标要求建立最小化子阵数的最优子阵划分模型为:
Figure GDA0003284940100000044
Figure GDA0003284940100000045
其中,FTaylor(θ)表示由方向图指标所得的泰勒加权方向图;
采用凸优化算法进行求解,并从稀疏解S中提取出子阵数、子阵尺寸和子阵激励,从而在获得最优子阵划分的同时,实现了与泰勒加权阵相当的辐射特性Fy(θ)。
作为上述天线的一种改进,所述根据可分离平面阵形成原理,由距离向和方位向的设计结果确定平面阵天线的激励参数;具体为:
根据可分离平面阵的形成原理,由距离向赋形方向图Fx(θ)和方位向辐射特性Fy(θ)得到平面阵方向图:
Figure GDA0003284940100000051
其中,
Figure GDA0003284940100000052
为俯仰角,在阵列结构确定的情况下,每个阵元的电流值等于方位向所得结果和距离向所得结果的乘积,由此即可确定平面阵天线的阵元激励。
与现有技术相比,本发明的优势在于:
本发明的天线对天线距离向进行赋形波束综合,实现入射角范围内回波能量一致性观测,有利于提高干涉测高精度,对方位向进行子阵划分,不仅有利于提高天线带宽,而且极大地降低数字处理通道数和馈电网络的复杂度,大大简化整个天线系统结构,降低硬件实现成本和工程难度,并且通道数远小于阵元数,实现了接收信号的降维处理。
附图说明
图1是本发明的星载干涉成像高度计天线的结构图;
图2是本发明的星载干涉成像高度计天线的设计方法流程图;
图3是仿真实例的天线距离向方向图特性及赋形包络;
图4(a)是仿真实例在天线距离向口径为10.5λ上设计的22个阵元位置与激励幅度分布图;
图4(b)是仿真实例在天线距离向口径为10.5λ上设计的22个阵元位置与激励相位分布图;
图5是仿真实例在天线方位向口径为109.5λ上220个阵元的子阵划分结果;
图6是仿真实例在天线方位向的方向图特性。
具体实施方式
本发明公开了一种星载干涉成像高度计天线,具体设计方法包括:由传统星载干涉成像高度计对动态海面的观测数据,提取后向散射系数随入射角的变化特征;计算出入射角范围内后向散射系数的变化量,并对变化曲线进行反加权,确定入射角范围内天线的赋形包络;结合副瓣电平、天线阵元、间隔以及幅相约束指标,确定天线距离向的激励参数与赋形波束;然后利用子阵划分技术完成天线方位向的设计,通过子阵级幅度加权实现与阵元级泰勒加权相当的辐射特性;最终根据可分离平面阵形成原理,由距离向和方位向的设计结果形成星载干涉成像高度计的平面阵天线。
下面结合附图对本发明的技术方案进行详细的说明。
现结合附图对本发明作进一步的描述。
如图1所示,为本发明的星载干涉成像高度计天线的一种示意,基于物理尺寸给定的赋形波束平面阵天线,包括若干个子阵,每个子阵包括若干个阵元,每个子阵中的阵元激励等幅同相,所述天线的距离向采用赋形波束设计,由入射角范围内的赋形曲线、副瓣电平、天线阵元以及幅相约束指标确定赋形波束及激励;所述天线的方位向设计采用子阵划分技术,通过子阵级幅度加权实现期望的泰勒加权方向图特性的笔形波束,根据可分离平面阵形成原理,由距离向和方位向的设计结果确定天线的激励参数。需要说明的是,图1仅仅是示意,以子阵个数8为例。
具体设计方法包括以下步骤:如图2所示。
S01:提取后向散射系数随入射角的变化特征;
根据动态海面的观测数据,提取后向散射系数随入射角的变化特征。
步骤S01具体包括如下过程:
对星载干涉成像高度计接收通道的动态海平面观测数据,进行合成孔径成像处理,得到图像功率值Pc,考虑到噪声功率比较小,可直接忽略,然后校正掉天线距离向方向图G(θ)和观测斜距R后,得到在入射角θ范围内,海面后向散射系数的相对值,即为
Figure GDA0003284940100000061
S02:确定入射角范围内的赋形方向图包络;
计算入射角范围内后向散射系数的变化量,对变化曲线进行反加权,确定入射角范围内天线方向图的赋形包络。
步骤S02具体包括如下步骤:
由后向散射系数σ0变化特征可得其在入射角范围内的最大值G1和最小值G0,那么为了实现距离向一致性观测,有效提高测高精度,赋形波束SB(θ)在入射角范围内应满足如下准则,即:
Figure GDA0003284940100000071
其中[θ01]为入射范围,据此,采用线性拟合可得入射角范围内的赋形方向图包络为
Figure GDA0003284940100000072
S03:天线距离向波束赋形;
结合入射角范围内的赋形包络、副瓣电平、天线阵元以及幅相约束指标,确定天线距离向的激励参数与赋形波束。
步骤S03具体包括如下步骤:
步骤S03-1)天线距离向阵元数为M,根据入射角范围内的赋形包络SB(θ),期望的副瓣电平SLLd,阵元数M以及阵元均匀间隔d的指标要求,可得距离向的赋形波束优化模型为:
Figure GDA0003284940100000073
其中ε为赋形误差,w为M个阵元的复激励,矩阵ψ(θ)可具体表示为:
Figure GDA0003284940100000074
为了对该问题进行高效求解,分别在副瓣区域SL和主瓣区域SB中对方向图进行采样,采样点位置分别为θk(k=1,2,…,K)和θn(n=1,2,…,N),综合问题可以进一步描述为
Figure GDA0003284940100000075
上式中的第一个约束函数等价于
Figure GDA0003284940100000076
分析可知,如果已知wr与wl中的任意一个,则优化问题相对于另一个而言则退化为凸优化问题,假如已知wl,相应的优化问题则转化为一个凸问题,具体表示为
Figure GDA0003284940100000081
易知,wr=wl是上述凸问题与原始综合问题相互等效的前提条件,这就需要在每次迭代过程中,根据计算的wr与当前迭代中的wl进行平滑操作来更新wl,并将其作为下次迭代的初值,直到max|wr-wl|≤10-3则迭代终止。这样就可以利用连续凸优化算法快速优化复激励,在无幅相约束下,获得期望赋形方向图Fd(θ),然后由此提取出主波束的两个零点θNULL1和θNULL2以及赋形主波束方向图Fdn)(θNULL1≤θn≤θNULL2);
步骤S03-2)以阵元激励幅度和相位作为优化变量,并根据实际需求约束二者的动态范围,结合期望的副瓣电平指标SLLd和主波束赋形方向图Fdn)(θNULL1≤θn≤θNULL2),建立优化过程中的代价函数为:
Figure GDA0003284940100000082
其中副瓣电平SLL可由下式确定,即
SLL=max(|F(χt)|),χt∈[-90,θNULL1)∪(θNULL2,90] (10)
利用入侵杂草优化算法对上述综合问题进行求解,计算出阵元激励,并由此得到天线距离向赋形方向图Fx(θ)。
S04:天线方位向子阵划分;
根据天线方位向方向图的指标要求,将P个均匀分布阵元划分成Q个子阵,并仅在Q个子阵端口进行激励加权,每个子阵中的阵元激励等幅同相,建立最小化子阵数的最优子阵划分模型,利于凸优化算法对其高效求解,实现与期望的泰勒加权方向图相当的辐射特性Fy(θ)。
步骤S04具体包括如下步骤:
将天线方位向分布的P个均匀分布阵元划分成Q个子阵,并仅在Q个子阵端口进行激励加权,每个子阵中的阵元激励等幅同相,那么相应的阵因子F2(θ)可表示为:
Figure GDA0003284940100000091
其中dp表示第p个阵元的位置,wp表示第p个阵元的激励值,其中p=1,…,P,aq表示第q个子阵端口的激励权值,其中q=1,…,Q。由此可知,子阵加权矢量Wsub=[a1,a2,…,aQ]T,阵元激励矢量Wele=[w1,w2,…,wP]T,由于阵列天线仅在子阵级加权,因此,阵元激励矢量Wele可以用子阵权值具体表示为
Figure GDA0003284940100000092
其中Nq是第q个子阵所包含的阵元数,用来表征相应子阵的尺寸,并且ΣNq=P。由信号处理理论可知,阵元激励矢量Wele是一个可压缩矢量,可以用一个稀疏矢量来表示,由矢量Wele的元素特征可知,相应稀疏矢量的稀疏度为Q,即仅包含Q个非零元,稀疏矢量S具体表示为:
Figure GDA0003284940100000093
由此可推导出稀疏矢量S与阵元激励矢量Wele之间的数值关系为:
Wele=TS (14)
其中转换矩阵T是一个下三角0-1矩阵,具体为
Figure GDA0003284940100000094
由此可知,稀疏矢量S是阵元激励矢量Wele的一种稀疏表示,结合稀疏矢量S和阵因子表达式,可以推导出阵因子的矢量形式:
Fy(θ)=ψ(θ)S (16)
其中
ψ(θ)=A(θ)T (17)
上式中导向矢量A(θ)为
Figure GDA0003284940100000101
根据方向图指标要求,子阵划分问题就可以等效为一个凸优化模型即:
Figure GDA0003284940100000102
其中,FTaylor(θ)表示由方向图指标所得的泰勒加权方向图,这就可以利用凸优化算法对上述子阵划分问题进行求解,并从稀疏解S中提取出子阵划分的所有参数,即子阵数、子阵尺寸和子阵激励,最终在获得最优子阵划分的同时,实现了与泰勒加权阵相当的辐射特性Fy(θ)。
S05:形成星载干涉成像高度计的平面阵天线;
根据可分离平面阵形成原理,由距离向和方位向的设计结果形成星载干涉成像高度计的平面天线。
步骤S05具体包括如下步骤:
根据可分离平面阵的形成原理,由天线距离向和方位向的综合结果可得平面阵方向图:
Figure GDA0003284940100000103
其中
Figure GDA0003284940100000104
为俯仰角,在阵列结构确定的情况下,每个阵元的电流值等于方位向所得结果和距离向所得结果的乘积,由此可以确定平面阵天线的激励参数。
本发明所提出的一种星载干涉成像高度计天线设计方法,可以进一步通过以下具体的仿真实例来验证并说明。
仿真实例:
本实例在给定的阵列口面Lx×Ly=10.5λ×109.5λ内半波长均匀排布天线单元,沿距离向和方位向分别实现一个期望的赋形波束和笔形波束,其中要求距离向的入射角区域为[0°,7°],由海面后向散射系数可推导出赋形区域的指标要求,即波束最大指向为7°,且随着入射角的增大,赋形波束功率方向图增大7dB,赋形区域内波纹抖动不超过0.5dB,并要求赋形波束左右两边的副瓣电平SLLd分别为-22dB和-18dB,各阵元相位差不超过300°,幅度动态范围不小于0.2;要求方位向实现-25dB的泰勒加权方向图特性。采用本发明提出的一种星载干涉成像高度计天线设计方法进行设计,距离向波束仿真结果如图3所示,由此可知,本发明方法综合的赋形波束极大地保证了距离向一致性观测,利于提高干涉测高精度,相应的阵元位置与激励分布如图4(a)、图4(b)所示,根据泰勒加权方向图特性,将天线方位向按左右对称进行子阵划分,结果如图5所示,8个子阵中的阵元数分别为32、16、21、25、25、21、16和32,子阵划分具有结构对称特性,仅在子阵级激励加权,而每个子阵中阵元的激励等幅同相,这样便于设计与加工,相应的方向图特性如图6所示,显然实现了期望的笔形波束,同时将通道数由阵元级的220个显著降低为8个,有效降低了馈电网络的复杂度,也利于展开天线带宽,这具有非常显著的工程意义和应用价值。由可分离平面阵的形成原理可知,由上述距离向和方位向的综合结果,即可得到平面阵天线辐射特性。
最后所应说明的是,以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,对本发明的技术方案进行修改或者等同替换,都不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。

Claims (6)

1.一种星载干涉成像高度计天线,其特征在于,所述天线为基于物理尺寸给定的赋形波束平面阵天线,包括若干个子阵,每个子阵包括若干个阵元,每个子阵中的阵元激励等幅同相,所述天线的距离向采用赋形波束设计,由入射角范围内的赋形曲线、副瓣电平、天线阵元以及幅相约束指标确定赋形波束及激励;所述天线的方位向设计采用子阵划分技术,通过子阵级幅度加权实现期望的泰勒加权方向图特性的笔形波束,根据可分离平面阵形成原理,由距离向和方位向的设计结果确定天线的激励参数;所述赋形波束设计包括以下步骤实现:
根据动态海面的观测数据,提取随入射角变化的后向散射系数的相对值;
由后向散射系数的相对值得到天线方向图在入射角范围内的最大值和最小值,采用线性拟合得到入射角范围内的赋形方向图包络;
结合入射角范围内的赋形方向图包络、副瓣电平、天线阵元以及幅相约束指标,确定天线距离向的激励参数与赋形波束;
所述结合入射角范围内的赋形方向图包络、副瓣电平、天线阵元以及幅相约束指标,确定天线距离向的激励参数与赋形波束;具体为:
根据入射角范围内的赋形方向图包络SB(θ)、期望的副瓣电平SLLd、天线距离向阵元数M以及阵元均匀间隔d的指标要求,采用凸优化算法快速优化复激励,在无幅相约束下,得到期望赋形方向图Fd(θ),由此提取出主波束的两个零点θNULL1和θNULL2以及期望的赋形主波束方向图Fdn),θn为赋形主波束方向图第n个入射角,θNULL1≤θn≤θNULL2;将阵元激励和幅度作为优化变量,根据幅相约束指标,以对期望的赋形主波束方向图Fdn)和期望的副瓣电平SLLd的拟合为适应度函数,采用入侵杂草优化算法进行优化,计算出阵元激励,并由此得到天线距离向的激励参数与赋形波束。
2.根据权利要求1所述的星载干涉成像高度计天线,其特征在于,所述根据动态海面的观测数据,提取随入射角变化的后向散射系数的相对值;具体为:
对星载干涉成像高度计采集的动态海平面观测数据,进行合成孔径成像处理,得到图像功率值Pc
由天线距离向方向图G(θ)和观测斜距R,得到在入射角θ范围内,海面后向散射系数的相对值σ0为:
Figure FDA0003284940090000021
3.根据权利要求2所述的星载干涉成像高度计天线,其特征在于,所述由后向散射系数的相对值得到天线方向图在入射角范围内的最大值和最小值,采用线性拟合得到入射角范围内的赋形方向图包络;具体为:
由后向散射系数的相对值σ0得到天线方向图在入射角范围内的最大值G1和最小值G0
采用线性拟合得到入射角范围内的赋形方向图包络SB(θ)为:
Figure FDA0003284940090000022
其中,[θ01]为入射角的范围。
4.根据权利要求3所述的星载干涉成像高度计天线,其特征在于,所述适应度函数fit(θ)为:
Figure FDA0003284940090000023
其中,N为赋形主波束方向图的入射角个数,F(θn)为赋形方向图,SLL为副瓣电平:
SLL=max(|F(χt)|),χt∈[-90,θNULL1)∪(θNULL2,90]。
5.根据权利要求4所述的星载干涉成像高度计天线,其特征在于,所述方位向设计采用子阵划分技术,通过子阵级幅度加权实现期望的泰勒加权方向图特性的笔形波束;具体包括:
将天线方位向分布的P个均匀分布阵元划分成Q个子阵,并仅在Q个子阵端口进行激励加权,每个子阵中的阵元激励等幅同相,相应的阵因子F2(θ)表示为:
Figure FDA0003284940090000024
其中,dp表示第p个阵元的位置,wp表示第p个阵元的激励值,p=1,…,P,aq表示第q个子阵端口的激励权值,q=1,…,Q,j为虚部;
根据子阵加权矢量Wsub=[a1,a2,…,aQ]T,将阵元激励矢量Wele=[w1,w2,…,wP]T,具体表示为:
Figure FDA0003284940090000031
其中,Nq是第q个子阵所包含的阵元数,用于表征相应子阵的尺寸,并且ΣNq=P;根据信号处理理论可知,阵元激励矢量Wele是一个可压缩矢量,用一个稀疏矢量来表示,由矢量Wele的元素特征可知,相应稀疏矢量的稀疏度为Q,即仅包含Q个非零元,稀疏矢量S具体表示为:
Figure FDA0003284940090000032
由此得出稀疏矢量S与阵元激励矢量Wele之间的数值关系为:
Wele=TS
其中,转换矩阵T为:
Figure FDA0003284940090000033
辐射特性Fy(θ)的矢量形式为:
Fy(θ)=ψ(θ)S
ψ(θ)=A(θ)T
其中,导向矢量A(θ)为:
Figure FDA0003284940090000034
根据方向图指标要求建立最小化子阵数的最优子阵划分模型为:
Figure FDA0003284940090000035
Figure FDA0003284940090000036
其中,ε为赋形误差,FTaylor(θ)表示由方向图指标所得的泰勒加权方向图;
采用凸优化算法进行求解,并从稀疏矢量S中提取出子阵数、子阵尺寸和子阵激励,从而在获得最优子阵划分的同时,实现了与泰勒加权阵相当的辐射特性Fy(θ)。
6.根据权利要求5所述的星载干涉成像高度计天线,其特征在于,所述根据可分离平面阵形成原理,由距离向和方位向的设计结果确定天线的激励参数;具体为:
根据可分离平面阵的形成原理,由距离向赋形方向图Fx(θ)和方位向辐射特性Fy(θ)得到平面阵方向图:
Figure FDA0003284940090000041
其中,
Figure FDA0003284940090000042
Figure FDA0003284940090000043
为俯仰角,在阵列结构确定的情况下,每个阵元的电流值等于方位向所得结果和距离向所得结果的乘积,由此即可确定平面阵天线的阵元激励。
CN202110260165.2A 2021-03-10 2021-03-10 一种星载干涉成像高度计天线 Active CN113126087B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110260165.2A CN113126087B (zh) 2021-03-10 2021-03-10 一种星载干涉成像高度计天线

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110260165.2A CN113126087B (zh) 2021-03-10 2021-03-10 一种星载干涉成像高度计天线

Publications (2)

Publication Number Publication Date
CN113126087A CN113126087A (zh) 2021-07-16
CN113126087B true CN113126087B (zh) 2021-11-09

Family

ID=76772899

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110260165.2A Active CN113126087B (zh) 2021-03-10 2021-03-10 一种星载干涉成像高度计天线

Country Status (1)

Country Link
CN (1) CN113126087B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113794062B (zh) * 2021-11-17 2022-03-29 南京天朗防务科技有限公司 和差波束天线的子阵划分和电压分布方法、装置
CN114201891B (zh) * 2022-02-17 2022-05-06 中国科学院空天信息创新研究院 一种天线方位多通道间基线和相位中心的配置方法
CN115755051B (zh) * 2022-11-18 2023-07-25 北京卫星信息工程研究所 双星高分宽幅sar信号分布式在轨处理方法及系统

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102110883B (zh) * 2010-12-01 2013-06-19 西安空间无线电技术研究所 一种赋形可变波束阵列天线的波束赋形方法
CN102299421B (zh) * 2011-05-31 2014-11-19 西安空间无线电技术研究所 一种幅相加权的窄边波导缝隙阵列天线
US9531081B2 (en) * 2011-07-20 2016-12-27 Deutsches Zentrum für Luft- und Raumfahrt e.V. Reflector antenna for a synthetic aperture radar
CN104035094B (zh) * 2014-06-25 2016-10-26 内蒙古工业大学 接收距离向多孔径宽幅星载sar回波的方法和装置
CN104993251B (zh) * 2015-06-26 2017-07-14 中国船舶重工集团公司第七二四研究所 一种大型平面阵列天线方向图综合级联优化方法
CN106981728B (zh) * 2017-03-14 2019-12-24 中国科学院国家空间科学中心 一种方向图可重构的稀布直线阵综合方法
CN108987941B (zh) * 2018-05-22 2020-07-03 中国科学院国家空间科学中心 一种基于压缩感知的一维子阵划分方法
DE102018208366B3 (de) * 2018-05-28 2019-11-28 Deutsches Zentrum für Luft- und Raumfahrt e.V. Synthetik-Apertur-Radarverfahren und Synthetik-Apertur-Radarvorrichtung
CN111262612B (zh) * 2020-01-16 2022-11-22 电子科技大学 一种基于笔形波束的阵列赋形波束解析合成方法

Also Published As

Publication number Publication date
CN113126087A (zh) 2021-07-16

Similar Documents

Publication Publication Date Title
CN113126087B (zh) 一种星载干涉成像高度计天线
CN109212526B (zh) 用于高频地波雷达的分布式阵列目标角度测量方法
CN109870687B (zh) 一种用于互质采样星载sar的成像处理方法
CN108459307B (zh) 基于杂波的mimo雷达收发阵列幅相误差校正方法
CN109061642B (zh) 一种贝叶斯迭代重加权稀疏自聚焦阵列sar成像方法
CN108549059B (zh) 一种复杂地形条件下的低空目标仰角估计方法
CN103399315B (zh) 实孔径相控阵雷达高分辨探测成像方法
CN109581352B (zh) 一种基于毫米波雷达的超分辨测角系统
CN108693511B (zh) 时分复用mimo雷达的运动目标角度计算方法
CN109765521B (zh) 一种基于子阵划分的波束域成像方法
CN110346794B (zh) 一种资源优化配置的分布式雷达成像方法
CN113314832B (zh) 一种毫米波车载mimo雷达天线阵列装置及设计方法
CN107942295B (zh) 一种前视阵列sar系统的稀疏天线
CN111537997B (zh) 一种基于mimo和压缩感知技术的三维雷达成像方法
CN108828504B (zh) 基于部分相关波形的mimo雷达目标方向快速估计方法
CN110646765B (zh) 一种基于黎曼距离的广义旁瓣对消算法
CN106707284A (zh) 一种用于双基地逆合成孔径雷达的成像方法
CN108828586B (zh) 一种基于波束域的双基地mimo雷达测角优化方法
Hélal et al. Radar imaging and high-resolution array processing applied to a classical VHF-ST profiler
CN111090094B (zh) 脉冲多普勒雷达的双波束角度测量方法、系统及存储介质
CN113671485B (zh) 基于admm的米波面阵雷达二维doa估计方法
CN111812612A (zh) 基于子阵正交lfm信号的mimo雷达部分相关波形设计方法
CN113359196B (zh) 基于子空间法和dbf的多目标生命体征探测方法
CN116184343A (zh) 基于相控阵雷达的三维空间蜂群目标检测及信息估计方法
RU2752878C2 (ru) Способ пеленгации широкополосных сигналов с повышенной разрешающей способностью

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