CN113126087A - 一种星载干涉成像高度计天线 - Google Patents
一种星载干涉成像高度计天线 Download PDFInfo
- Publication number
- CN113126087A CN113126087A CN202110260165.2A CN202110260165A CN113126087A CN 113126087 A CN113126087 A CN 113126087A CN 202110260165 A CN202110260165 A CN 202110260165A CN 113126087 A CN113126087 A CN 113126087A
- Authority
- CN
- China
- Prior art keywords
- antenna
- array
- excitation
- shaped
- sub
- 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.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Systems 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/88—Radar or analogous systems specially adapted for specific applications
- G01S13/882—Radar or analogous systems specially adapted for specific applications for altimeters
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Systems 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/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9021—SAR image post-processing techniques
- G01S13/9023—SAR image post-processing techniques combined with interferometric techniques
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
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,”2017 Progress 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 SurfaceHeight from Chinese Tiangong-2 Interferometric Imaging Radar Altimeter”,Remote Sensing,vol.12,no.15,2496,2020),对海洋科学、海洋环境观测和预报、全球气候变化研究具有重要的作用与意义。
天线作为星载干涉成像高度计重要组成部分,用来发射和接收电磁波,起着至关重要的作用。当前星载干涉成像高度计天线距离向采用笔形波束设计,然而,结合电磁场理论与星载数据分析可知,海面后向散射系数随入射角的增大而减小,可见回波能量在观测范围内分布不均匀,影响干涉测高精度,这不利于对中尺度海洋现象的高精度观测。
同时,为了获得更远的作用距离和更高的空间分辨率,星载干涉成像高度计天线必然也要朝着大型化的方向发展,但天线口径的增大,阵元数的增多必然导致馈电网络极其复杂,也影响着天线带宽,可见,当前采用的阵元级加权方法已不适用于大型天线,很难满足下一代星载干涉成像高度计宽幅宽、高分辨的系统需求。
发明内容
当前星载干涉成像高度计天线距离向观测不一致,馈电网络复杂,加工难度大,工作带宽窄,不利于对中尺度海洋现象高精度观测的问题。本发明的目的在于克服上述技术缺陷,提出了一种星载干涉成像高度计天线。
为了实现上述目的,本发明提出了一种星载干涉成像高度计天线,所述天线为基于物理尺寸给定的赋形波束平面阵天线,包括若干个子阵,每个子阵包括若干个阵元,每个子阵中的阵元激励等幅同相,所述天线的距离向采用赋形波束设计,由入射角范围内的赋形曲线、副瓣电平、天线阵元以及幅相约束指标确定赋形波束及激励;所述天线的方位向设计采用子阵划分技术,通过子阵级幅度加权实现期望的泰勒加权方向图特性的笔形波束,根据可分离平面阵形成原理,由距离向和方位向的设计结果确定天线的激励参数。
作为上述天线的一种改进,所述赋形波束设计包括以下步骤实现:
根据动态海面的观测数据,提取随入射角变化的后向散射系数的相对值;
由后向散射系数的相对值得到在入射角范围内的最大值和最小值,采用线性拟合得到入射角范围内的赋形方向图包略;
结合入射角范围内的赋形方向图包略、副瓣电平、天线阵元以及幅相约束指标,确定天线距离向的激励参数与赋形波束。
作为上述天线的一种改进,所述根据动态海面的观测数据,提取随入射角变化的后向散射系数的相对值;具体为:
对星载干涉成像高度计采集的动态海平面观测数据,进行合成孔径成像处理,得到图像功率值Pc;
由天线距离向方向图G(θ)和观测斜距R,得到在入射角θ范围内,海面后向散射系数的相对值σ0为:
作为上述天线的一种改进,所述由后向散射系数的相对值得到在入射角范围内的最大值和最小值,采用线性拟合得到入射角范围内的赋形方向图包略;具体为:
由后向散射系数的相对值σ0得到入射角范围内的最大值G1和最小值G2,
采用线性拟合得到入射角范围内的赋形方向图包略SB(θ)为:
其中,[θ0,θ1]为入射角的范围。
作为上述天线的一种改进,所述结合入射角范围内的赋形方向图包略、副瓣电平、天线阵元以及幅相约束指标,确定天线距离向的激励参数与赋形波束;具体为:
根据入射角范围内的赋形方向图包略SB(θ),期望的副瓣电平SLLd,阵元数M以及阵元均匀间隔d的指标要求,采用凸优化算法快速优化复激励,在无幅相约束下,得到期望赋形方向图Fd(θ),由此提取出主波束的两个零点θNULL1和θNULL2以及期望的赋形主波束方向图Fd(θn),θn为赋形主波束方向图第n个入射角,θNULL1≤θn≤θNULL2;
将阵元激励和幅度作为优化变量,根据幅相约束指标,以对期望的主波束赋形方向图Fd(θn)和期望的副瓣电平SLLd的拟合为适应度函数,采用入侵杂草优化算法进行优化,计算出阵元激励,并由此得到天线距离向的激励参数与赋形波束。
作为上述天线的一种改进,所述适应度函数fit(θ)为:
其中,N为赋形主波束方向图的入射角个数,F(θn)为赋形方向图,SLL为副瓣电平:
SLL=max(|F(θ)|),θ∈[-90,θNULL1)∪(θNULL2,90]。
作为上述天线的一种改进,所述方位向采用子阵划分技术,通过子阵级幅度加权实现期望的泰勒加权方向图特性的笔形波束;具体为:
将天线方位向分布的P个均匀分布阵元划分成Q个子阵,并仅在Q个子阵端口进行激励加权,每个子阵中的阵元激励等幅同相,相应的阵因子F2(θ)表示为:
其中,dp表示第p个阵元的位置,wp表示第p个阵元的激励值,p=1,...,P,aq表示第q个子阵端口的激励权值,q=1,...,Q,j为虚部;
根据子阵加权矢量Wsub=[a1,a2,...,aQ]T,阵元激励矢量Wele=[w1,w2,...,wP]T,具体表示为:
其中,Nq是第q个子阵所包含的阵元数,用于表征相应子阵的尺寸,并且∑Nq=P;根据信号处理理论可知,阵元激励矢量Wele是一个可压缩矢量,可以用一个稀疏矢量来表示,由矢量Wele的元素特征可知,相应稀疏矢量的稀疏度为Q,即仅包含Q个非零元,稀疏矢量S具体表示为:
由此得出稀疏矢量S与阵元激励矢量Wele之间的数值关系为:
Wele=TS
其中,转换矩阵T为:
辐射特性Fy(θ)的矢量形式为:
Fy(θ)=ψ(θ)S
ψ(θ)=A(θ)T
其中,导向矢量A(θ)为:
根据方向图指标要求建立最小化子阵数的最优子阵划分模型为:
其中,FTaylor(θ)表示由方向图指标所得的泰勒加权方向图;
采用凸优化算法进行求解,并从稀疏解S中提取出子阵数、子阵尺寸和子阵激励,从而在获得最优子阵划分的同时,实现了与泰勒加权阵相当的辐射特性Fy(θ)。
作为上述天线的一种改进,所述根据可分离平面阵形成原理,由距离向和方位向的设计结果确定平面阵天线的激励参数;具体为:
根据可分离平面阵的形成原理,由距离向赋形方向图Fx(θ)和方位向辐射特性Fy(θ)得到平面阵方向图:
与现有技术相比,本发明的优势在于:
本发明的天线对天线距离向进行赋形波束综合,实现入射角范围内回波能量一致性观测,有利于提高干涉测高精度,对方位向进行子阵划分,不仅有利于提高天线带宽,而且极大地降低数字处理通道数和馈电网络的复杂度,大大简化整个天线系统结构,降低硬件实现成本和工程难度,并且通道数远小于阵元数,实现了接收信号的降维处理。
附图说明
图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后,得到在入射角θ范围内,海面后向散射系数的相对值,即为
S02:确定入射角范围内的赋形方向图包略;
计算入射角范围内后向散射系数的变化量,对变化曲线进行反加权,确定入射角范围内天线方向图的赋形包略。
步骤S02具体包括如下步骤:
由后向散射系数σ0变化特征可得其在入射角范围内的最大值G1和最小值G2,那么为了实现距离向一致性观测,有效提高测高精度,赋形波束SB(θ)在入射角范围内应满足如下准则,即:
其中[θ0,θ1]为入射范围,据此,采用线性拟合可得入射角范围内的赋形方向图包略为
S03:天线距离向波束赋形;
结合入射角范围内的赋形包略、副瓣电平、天线阵元以及幅相约束指标,确定天线距离向的激励参数与赋形波束。
步骤S03具体包括如下步骤:
步骤S03-1)天线距离向阵元数为M,根据入射角范围内的赋形包略SB(θ),期望的副瓣电平SLLd,阵元数M以及阵元均匀间隔d的指标要求,可得距离向的赋形波束优化模型为:
其中ε为赋形误差,w为M个阵元的复激励,矩阵ψ(θ)可具体表示为:
为了对该问题进行高效求解,分别在副瓣区域SL和主瓣区域SB中对方向图进行采样,采样点位置分别为θk(k=l,2,…,K)和θn(n=l,2,…,N),综合问题可以进一步描述为
上式中的第一个约束函数等价于
分析可知,如果已知wr与wl中的任意一个,则优化问题相对于另一个而言则退化为凸优化问题,假如已知wl,相应的优化问题则转化为一个凸问题,具体表示为
易知,wr=wl是上述凸问题与原始综合问题相互等效的前提条件,这就需要在每次迭代过程中,根据计算的wr与当前迭代中的wl进行平滑操作来更新wl,并将其作为下次迭代的初值,直到max|wr-wl|≤10-3则迭代终止。这样就可以利用连续凸优化算法快速优化复激励,在无幅相约束下,获得期望赋形方向图Fd(θ),然后由此提取出主波束的两个零点θNLL1和θNULL2以及赋形主波束方向图Fd(θn)(θNULL1≤θn≤θNULL2);
步骤S03-2)以阵元激励幅度和相位作为优化变量,并根据实际需求约束二者的动态范围,结合期望的副瓣电平指标SLLd和主波束赋形方向图Fd(θn)(θNULL1≤θn≤θNULL2),建立优化过程中的代价函数为:
其中副瓣电平SLL可由下式确定,即
SLL=max(|F(θ)|),θ∈[-90,θNULL1)∪(θNULL2,90] (10)
利用入侵杂草优化算法对上述综合问题进行求解,计算出阵元激励,并由此得到天线距离向赋形方向图Fx(θ)。
S04:天线方位向子阵划分;
根据天线方位向方向图的指标要求,将P个均匀分布阵元划分成Q个子阵,并仅在Q个子阵端口进行激励加权,每个子阵中的阵元激励等幅同相,建立最小化子阵数的最优子阵划分模型,利于凸优化算法对其高效求解,实现与期望的泰勒加权方向图相当的辐射特性Fy(θ)。
步骤S04具体包括如下步骤:
将天线方位向分布的P个均匀分布阵元划分成Q个子阵,并仅在Q个子阵端口进行激励加权,每个子阵中的阵元激励等幅同相,那么相应的阵因子F2(θ)可表示为:
其中dp表示第p个阵元的位置,wp表示第p个阵元的激励值,其中p=1,...,P,aq表示第q个子阵端口的激励权值,其中q=1,...,Q。由此可知,子阵加权矢量Wsub=[a1,a2,...,aQ]T,阵元激励矢量Wele=[w1,w2,...,wP]T,由于阵列天线仅在子阵级加权,因此,阵元激励矢量Wele可以用子阵权值具体表示为
其中Nq是第q个子阵所包含的阵元数,用来表征相应子阵的尺寸,并且∑Nq=P。由信号处理理论可知,阵元激励矢量Wele是一个可压缩矢量,可以用一个稀疏矢量来表示,由矢量Wele的元素特征可知,相应稀疏矢量的稀疏度为Q,即仅包含Q个非零元,稀疏矢量S具体表示为:
由此可推导出稀疏矢量S与阵元激励矢量Wele之间的数值关系为:
Wele=TS (14)
其中转换矩阵T是一个下三角0-1矩阵,具体为
由此可知,稀疏矢量S是阵元激励矢量Wele的一种稀疏表示,结合稀疏矢量S和阵因子表达式,可以推导出阵因子的矢量形式:
Fy(θ)=ψ(θ)S (16)
其中
ψ(θ)=A(θ)T (17)
上式中导向矢量A(θ)为
根据方向图指标要求,子阵划分问题就可以等效为一个凸优化模型即:
其中,FTaylor(θ)表示由方向图指标所得的泰勒加权方向图,这就可以利用凸优化算法对上述子阵划分问题进行求解,并从稀疏解S中提取出子阵划分的所有参数,即子阵数、子阵尺寸和子阵激励,最终在获得最优子阵划分的同时,实现了与泰勒加权阵相当的辐射特性Fy(θ)。
S05:形成星载干涉成像高度计的平面阵天线;
根据可分离平面阵形成原理,由距离向和方位向的设计结果形成星载干涉成像高度计的平面天线。
步骤S05具体包括如下步骤:
根据可分离平面阵的形成原理,由天线距离向和方位向的综合结果可得平面阵方向图:
本发明所提出的一种星载干涉成像高度计天线设计方法,可以进一步通过以下具体的仿真实例来验证并说明。
仿真实例:
本实例在给定的阵列口面Lx×Ly=10.5λ×109.5λ内半波长均匀排布天线单元,沿距离向和方位向分别实现一个期望的赋形波束和笔形波束,其中要求距离向的入射角区域为[0°,7°],由海面后向散射系数可推导出赋形区域的指标要求,即波束最大指向为7°,且随着入射角的增大,赋形波束功率方向图增大7dB,赋形区域内波纹抖动不超过0.5dB,并要求赋形波束左右两边的副瓣电平SLLd分别为-22dB和-18dB,各阵元相位差不超过300°,幅度动态范围不小于0.2;要求方位向实现-25dB的泰勒加权方向图特性。采用本发明提出的一种星载干涉成像高度计天线设计方法进行设计,距离向波束仿真结果如图3所示,由此可知,本发明方法综合的赋形波束极大地保证了距离向一致性观测,利于提高干涉测高精度,相应的阵元位置与激励分布如图4(a)、(b)所示,根据泰勒加权方向图特性,将天线方位向按左右对称进行子阵划分,结果如图5所示,8个子阵中的阵元数分别为32、16、21、25、25、21、16和32,子阵划分具有结构对称特性,仅在子阵级激励加权,而每个子阵中阵元的激励等幅同相,这样便于设计与加工,相应的方向图特性如图6所示,显然实现了期望的笔形波束,同时将通道数由阵元级的220个显著降低为8个,有效降低了馈电网络的复杂度,也利于展开天线带宽,这具有非常显著的工程意义和应用价值。由可分离平面阵的形成原理可知,由上述距离向和方位向的综合结果,即可得到平面阵天线辐射特性。
最后所应说明的是,以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,对本发明的技术方案进行修改或者等同替换,都不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (8)
1.一种星载干涉成像高度计天线,其特征在于,所述天线为基于物理尺寸给定的赋形波束平面阵天线,包括若干个子阵,每个子阵包括若干个阵元,每个子阵中的阵元激励等幅同相,所述天线的距离向采用赋形波束设计,由入射角范围内的赋形曲线、副瓣电平、天线阵元以及幅相约束指标确定赋形波束及激励;所述天线的方位向设计采用子阵划分技术,通过子阵级幅度加权实现期望的泰勒加权方向图特性的笔形波束,根据可分离平面阵形成原理,由距离向和方位向的设计结果确定天线的激励参数。
2.根据权利要求1所述的星载干涉成像高度计天线,其特征在于,所述赋形波束设计包括以下步骤实现:
根据动态海面的观测数据,提取随入射角变化的后向散射系数的相对值;
由后向散射系数的相对值得到在入射角范围内的最大值和最小值,采用线性拟合得到入射角范围内的赋形方向图包略;
结合入射角范围内的赋形方向图包略、副瓣电平、天线阵元以及幅相约束指标,确定天线距离向的激励参数与赋形波束。
5.根据权利要求4所述的星载干涉成像高度计天线,其特征在于,所述结合入射角范围内的赋形方向图包略、副瓣电平、天线阵元以及幅相约束指标,确定天线距离向的激励参数与赋形波束;具体为:
根据入射角范围内的赋形方向图包略SB(θ),期望的副瓣电平SLLd,阵元数M以及阵元均匀间隔d的指标要求,采用凸优化算法快速优化复激励,在无幅相约束下,得到期望赋形方向图Fd(θ),由此提取出主波束的两个零点θNULL1和θNULL2以及期望的赋形主波束方向图Fd(θn),θn为赋形主波束方向图第n个入射角,θNULL1≤θn≤θNULL2;
将阵元激励和幅度作为优化变量,根据幅相约束指标,以对期望的主波束赋形方向图Fd(θn)和期望的副瓣电平SLLd的拟合为适应度函数,采用入侵杂草优化算法进行优化,计算出阵元激励,并由此得到天线距离向的激励参数与赋形波束。
7.根据权利要求6所述的星载干涉成像高度计天线,其特征在于,所述方位向采用子阵划分技术,通过子阵级幅度加权实现期望的泰勒加权方向图特性的笔形波束;具体为:
将天线方位向分布的P个均匀分布阵元划分成Q个子阵,并仅在Q个子阵端口进行激励加权,每个子阵中的阵元激励等幅同相,相应的阵因子F2(θ)表示为:
其中,dp表示第p个阵元的位置,wp表示第p个阵元的激励值,p=1,...,P,aq表示第q个子阵端口的激励权值,q=1,...,Q,j为虚部;
根据子阵加权矢量Wsub=[a1,a2,...,aQ]T,阵元激励矢量Wele=[w1,w2,...,wP]T,具体表示为:
其中,Nq是第q个子阵所包含的阵元数,用于表征相应子阵的尺寸,并且∑Nq=P;根据信号处理理论可知,阵元激励矢量Wele是一个可压缩矢量,可以用一个稀疏矢量来表示,由矢量Wele的元素特征可知,相应稀疏矢量的稀疏度为Q,即仅包含Q个非零元,稀疏矢量S具体表示为:
由此得出稀疏矢量S与阵元激励矢量Wele之间的数值关系为:
Wele=TS
其中,转换矩阵T为:
辐射特性Fy(θ)的矢量形式为:
Fy(θ)=ψ(θ)S
ψ(θ)=A(θ)T
其中,导向矢量A(θ)为:
根据方向图指标要求建立最小化子阵数的最优子阵划分模型为:
其中,FTaylor(θ)表示由方向图指标所得的泰勒加权方向图;
采用凸优化算法进行求解,并从稀疏解S中提取出子阵数、子阵尺寸和子阵激励,从而在获得最优子阵划分的同时,实现了与泰勒加权阵相当的辐射特性Fy(θ)。
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 true CN113126087A (zh) | 2021-07-16 |
CN113126087B 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) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113794062A (zh) * | 2021-11-17 | 2021-12-14 | 南京天朗防务科技有限公司 | 和差波束天线的子阵划分和电压分布方法、装置 |
CN114201891A (zh) * | 2022-02-17 | 2022-03-18 | 中国科学院空天信息创新研究院 | 一种天线方位多通道间基线和相位中心的配置方法 |
CN115755051A (zh) * | 2022-11-18 | 2023-03-07 | 北京卫星信息工程研究所 | 双星高分宽幅sar信号分布式在轨处理方法及系统 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102110883A (zh) * | 2010-12-01 | 2011-06-29 | 西安空间无线电技术研究所 | 一种赋形可变波束阵列天线的波束赋形方法 |
CN102299421A (zh) * | 2011-05-31 | 2011-12-28 | 西安空间无线电技术研究所 | 一种幅相加权的窄边波导缝隙阵列天线 |
WO2013011023A1 (de) * | 2011-07-20 | 2013-01-24 | Deutsches Zentrum für Luft- und Raumfahrt e.V. | Reflektorantenne für ein radar mit synthetischer apertur |
CN104035094A (zh) * | 2014-06-25 | 2014-09-10 | 内蒙古工业大学 | 接收距离向多孔径宽幅星载sar回波的方法和装置 |
CN104993251A (zh) * | 2015-06-26 | 2015-10-21 | 中国船舶重工集团公司第七二四研究所 | 一种大型平面阵列天线方向图综合级联优化方法 |
CN106981728A (zh) * | 2017-03-14 | 2017-07-25 | 中国科学院国家空间科学中心 | 一种方向图可重构的稀布直线阵综合方法 |
CN108987941A (zh) * | 2018-05-22 | 2018-12-11 | 中国科学院国家空间科学中心 | 一种基于压缩感知的一维子阵划分方法 |
DE102018208366B3 (de) * | 2018-05-28 | 2019-11-28 | Deutsches Zentrum für Luft- und Raumfahrt e.V. | Synthetik-Apertur-Radarverfahren und Synthetik-Apertur-Radarvorrichtung |
CN111262612A (zh) * | 2020-01-16 | 2020-06-09 | 电子科技大学 | 一种基于笔形波束的阵列赋形波束解析合成方法 |
-
2021
- 2021-03-10 CN CN202110260165.2A patent/CN113126087B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102110883A (zh) * | 2010-12-01 | 2011-06-29 | 西安空间无线电技术研究所 | 一种赋形可变波束阵列天线的波束赋形方法 |
CN102299421A (zh) * | 2011-05-31 | 2011-12-28 | 西安空间无线电技术研究所 | 一种幅相加权的窄边波导缝隙阵列天线 |
WO2013011023A1 (de) * | 2011-07-20 | 2013-01-24 | Deutsches Zentrum für Luft- und Raumfahrt e.V. | Reflektorantenne für ein radar mit synthetischer apertur |
CN104035094A (zh) * | 2014-06-25 | 2014-09-10 | 内蒙古工业大学 | 接收距离向多孔径宽幅星载sar回波的方法和装置 |
CN104993251A (zh) * | 2015-06-26 | 2015-10-21 | 中国船舶重工集团公司第七二四研究所 | 一种大型平面阵列天线方向图综合级联优化方法 |
CN106981728A (zh) * | 2017-03-14 | 2017-07-25 | 中国科学院国家空间科学中心 | 一种方向图可重构的稀布直线阵综合方法 |
CN108987941A (zh) * | 2018-05-22 | 2018-12-11 | 中国科学院国家空间科学中心 | 一种基于压缩感知的一维子阵划分方法 |
DE102018208366B3 (de) * | 2018-05-28 | 2019-11-28 | Deutsches Zentrum für Luft- und Raumfahrt e.V. | Synthetik-Apertur-Radarverfahren und Synthetik-Apertur-Radarvorrichtung |
CN111262612A (zh) * | 2020-01-16 | 2020-06-09 | 电子科技大学 | 一种基于笔形波束的阵列赋形波束解析合成方法 |
Non-Patent Citations (6)
Title |
---|
BRIAN D. POLLARD等: "The Wide Swath Ocean Altimeter: Radar Interferometry for Global Ocean Mapping with Centimetric Accuracy", 《PROCEEDINGS, IEEE AEROSPACE CONFERENCE》 * |
刘燕 等: "混合入侵杂草算法用于阵列天线波束赋形", 《西安电子科技大学学报(自然科学版)》 * |
梁广 等: "基于子阵列的低轨星载多波束相控阵天线的设计与实现", 《电子与信息学报》 * |
肖疆: "合成孔径雷达天线技术的若干关键问题研究", 《万方数据知识服务平台 学位论文》 * |
袁诚 等: "一种新型亚毫米波反射面天线设计", 《计算机仿真》 * |
赵晓雯: "稀布阵列天线的压缩感知和入侵杂草优化算法研究", 《中国优秀博硕士学位论文全文数据库(博士)信息科技辑》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113794062A (zh) * | 2021-11-17 | 2021-12-14 | 南京天朗防务科技有限公司 | 和差波束天线的子阵划分和电压分布方法、装置 |
CN114201891A (zh) * | 2022-02-17 | 2022-03-18 | 中国科学院空天信息创新研究院 | 一种天线方位多通道间基线和相位中心的配置方法 |
CN115755051A (zh) * | 2022-11-18 | 2023-03-07 | 北京卫星信息工程研究所 | 双星高分宽幅sar信号分布式在轨处理方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN113126087B (zh) | 2021-11-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113126087B (zh) | 一种星载干涉成像高度计天线 | |
CN109212526B (zh) | 用于高频地波雷达的分布式阵列目标角度测量方法 | |
CN109870687B (zh) | 一种用于互质采样星载sar的成像处理方法 | |
CN108459307B (zh) | 基于杂波的mimo雷达收发阵列幅相误差校正方法 | |
CN109581352B (zh) | 一种基于毫米波雷达的超分辨测角系统 | |
CN103399315B (zh) | 实孔径相控阵雷达高分辨探测成像方法 | |
CN113314832B (zh) | 一种毫米波车载mimo雷达天线阵列装置及设计方法 | |
CN107942295B (zh) | 一种前视阵列sar系统的稀疏天线 | |
CN109765521B (zh) | 一种基于子阵划分的波束域成像方法 | |
CN111257878B (zh) | 基于俯仰维频率脉内扫描高分宽幅sar的波形设计方法 | |
CN109031226B (zh) | 基于机电耦合的变形有源相控阵雷达探测性能快速评估方法 | |
CN111537997B (zh) | 一种基于mimo和压缩感知技术的三维雷达成像方法 | |
CN111190151A (zh) | 扫描模式下多模态小卫星sar的系统参数设计及发射功率优化方法 | |
CN113671485B (zh) | 基于admm的米波面阵雷达二维doa估计方法 | |
CN106707284A (zh) | 一种用于双基地逆合成孔径雷达的成像方法 | |
CN111090094B (zh) | 脉冲多普勒雷达的双波束角度测量方法、系统及存储介质 | |
CN111239742A (zh) | 海杂波背景下基于mbmc的低空风切变风速估计方法 | |
CN113359196B (zh) | 基于子空间法和dbf的多目标生命体征探测方法 | |
RU2752878C2 (ru) | Способ пеленгации широкополосных сигналов с повышенной разрешающей способностью | |
CN116184343A (zh) | 基于相控阵雷达的三维空间蜂群目标检测及信息估计方法 | |
CN108594235B (zh) | 一种提高方位向雷达反射率系数分辨率的方法和系统 | |
CN110618403B (zh) | 一种基于双波束雷达的着陆飞行器参数测量方法 | |
CN112180368A (zh) | 一种数据处理方法、设备、系统及存储介质 | |
CN113671477B (zh) | 一种基于图信号处理的雷达目标距离估计方法 | |
CN113671494B (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 |