CN114609581B - 一种比较扫描方向矢量投影极大值的阵列超分辨测向方法 - Google Patents
一种比较扫描方向矢量投影极大值的阵列超分辨测向方法 Download PDFInfo
- Publication number
- CN114609581B CN114609581B CN202210226222.XA CN202210226222A CN114609581B CN 114609581 B CN114609581 B CN 114609581B CN 202210226222 A CN202210226222 A CN 202210226222A CN 114609581 B CN114609581 B CN 114609581B
- Authority
- CN
- China
- Prior art keywords
- vector
- signal
- subspace
- projection
- angle
- 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
- 239000013598 vector Substances 0.000 title claims abstract description 269
- 238000000034 method Methods 0.000 title claims abstract description 24
- 239000011159 matrix material Substances 0.000 claims description 19
- 238000007781 pre-processing Methods 0.000 claims description 15
- 230000003595 spectral effect Effects 0.000 claims description 7
- 238000000354 decomposition reaction Methods 0.000 description 6
- 238000001228 spectrum Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
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
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/02—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
- G01S3/14—Systems for determining direction or deviation from predetermined direction
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明涉及一种比较扫描方向矢量投影极大值的阵列超分辨测向方法,在划分信号子空间和噪声子空间后,首先根据阵列结构在测向角度范围内生成每个角度的扫描方向矢量;第二,计算扫描方向矢量向信号子空间的投影矢量;第三,计算扫描方向矢量向噪声子空间的投影矢量;第四,选取扫描方向矢量向信号子空间投影矢量的极大值;第五,选取扫描方向矢量向噪声子空间投影矢量的极大值;第六,构建测向目标函数,形成空间谱图,估计信号入射方向角度。本发明综合利用扫描方向矢量向信号子空间和噪声子空间投影矢量的极大值关系,构建了新的测向目标函数,具有超高的信号方向测向分辨率。
Description
技术领域
本发明涉及阵列信号处理中的测向领域,尤其是通过综合利用扫描方向矢量向信号子空间和噪声子空间投影矢量的极大值关系,构建新的测向目标函数,能够获得超高的信号测向分辨率。
背景技术
测向是阵列信号处理的主要研究方向之一。早期雷达通过机械转动天线进行空中目标测向,接着相控阵雷达通过改变移相器相位进行空中目标测向,测向分辨率低;1959年Capon通过在主瓣方向约束增益不变的条件下使阵列输出平均功率极小化,提出了自适应波束形成器,显著提高了阵列测向分辨率,进入了高分辨测向时代。1986年Schmidt通过利用信号子空间和噪声子空间正交关系,提出了多重信号分类方法(MUSIC),进入了超分辨测向时代。此后,技术进展较小。
鉴于以上分析,有必要研究新的具有超高测向分辨率的新方法,在超分辨测向的基础上,进一步提高测向分辨率。
发明内容
本发明的技术解决问题:克服现有技术的不足,提供一种比较扫描方向矢量投影极大值的阵列超分辨测向方法,综合利用信号方向矢量与信号子空间等价关系、信号子空间和噪声子空间正交关系,分别计算扫描方向矢量向信号子空间的投影矢量、扫描方向矢量向噪声子空间的投影矢量,通过选取扫描方向矢量向信号子空间投影矢量的极大值、扫描方向矢量向噪声子空间投影矢量的极大值,构建新的测向目标函数,从形成的空间谱图中估计出信号入射方向角度,具有超分辨测向能力。
本发明的目的是通过如下技术方案实现的:
本发明一种比较扫描方向矢量投影极大值的阵列超分辨测向方法,对阵列接收数据进行预处理,获得信号子空间和噪声子空间,根据阵列结构在测向角度范围内生成每个角度的扫描方向矢量,并利用信号方向矢量与信号子空间等价关系、信号方向矢量与噪声子空间正交关系,分别计算扫描方向矢量向信号子空间的投影矢量和扫描方向矢量向噪声子空间的投影矢量,进一步选取扫描方向矢量向信号子空间投影矢量的极大值和扫描方向矢量向噪声子空间投影矢量的极大值,并构建新的测向目标函数,最后从形成的空间谱图中估计出信号入射方向角度,具有超分辨测向能力;
具体包括如下步骤:
步骤1、根据阵列结构在测向角度范围内生成每个角度的扫描方向矢量,当该角度与信号入射方向角度相同时,称扫描方向矢量为信号方向矢量,当该角度与信号入射方向角度不同时,称扫描方向矢量为其它扫描方向矢量;扫描方向矢量包括信号方向矢量和其它扫描方向矢量;
步骤2、对阵列接收数据预处理获得的信号子空间和噪声子空间,利用信号方向矢量属于信号子空间、其它扫描方向矢量不完全属于或不属于信号子空间的性质,计算扫描方向矢量向信号子空间的投影矢量;
步骤3、对阵列接收数据预处理获得的信号子空间和噪声子空间,利用信号方向矢量不属于噪声子空间、其它扫描方向矢量属于或不完全属于噪声子空间的性质,计算扫描方向矢量向噪声子空间的投影矢量;
步骤4、选取扫描方向矢量向信号子空间投影矢量的极大值;
步骤5、选取扫描方向矢量向噪声子空间投影矢量的极大值;
步骤6、根据步骤4的扫描方向矢量向信号子空间投影矢量的极大值和步骤5的扫描方向矢量向噪声子空间投影矢量的极大值,构建新的测向目标函数,在入射方向角度范围内,计算超分辨测向目标函数的值,形成空间谱图,空间谱图的谱峰所对应的角度即为信号入射方向角度的估计值。
所述步骤1包括:
步骤1、根据阵列结构生成对应入射方向角度的扫描方向矢量,以线阵为例,M元阵列接收空间不同入射方向角度θ1、θ2、...、θL的L个远场窄带信号,以下简称信号,L小于M,对于测向角度范围[-90°,90°]内的每个角度θ,生成扫描方向矢量a(θ),||||2是矢量l2范数算子;当该角度θ与信号入射方向角度θ1、θ2、...、θL相同时,扫描方向矢量分别为a(θ1)、a(θ2)、…、a(θL),扫描方向矢量就是信号方向矢量,当该角度θ与信号入射方向角度θ1、θ2、...、θL不同时,扫描方向矢量不是信号方向矢量,称为其它扫描方向矢量;
所述步骤2包括:
步骤2、对阵列接收数据预处理获得的信号子空间Es和噪声子空间En,计算扫描方向矢量a(θ)向信号子空间Es的投影矢量bs(θ):
bs(θ)=EsEs Ha(θ),-90°≤θ≤90°;
上标H是矩阵和矢量的共轭转置算子;对于信号入射方向角度θ1、θ2、...、θL,信号方向矢量a(θ1)、a(θ2)、…、a(θL)在信号子空间Es中,并与噪声子空间En正交,信号方向矢量a(θ1)、a(θ2)、…、a(θL)向信号子空间Es投影,投影矢量bs(θ1)、bs(θ2)、…、bs(θL)的元素绝对值较大;对于其它入射方向角度θ,扫描方向矢量a(θ)不完全在信号子空间Es中,甚至在噪声子空间En中,扫描方向矢量a(θ)向信号子空间Es的投影矢量bs(θ)的元素绝对值变小甚至为0。
所述步骤3包括:
步骤3、对阵列接收数据预处理获得的信号子空间Es和噪声子空间En,计算扫描方向矢量a(θ)向噪声子空间En的投影矢量bn(θ):
bn(θ)=EnEn Ha(θ),-90°≤θ≤90°;
对于信号入射方向角度θ1、θ2、...、θL,信号方向矢量a(θ1)、a(θ2)、…、a(θL)向噪声子空间En投影,投影矢量bn(θ1)、bn(θ2)、…、bn(θL)的元素绝对值为0;对于其它入射方向角度θ,扫描方向矢量a(θ)向噪声子空间En的投影矢量bn(θ)的元素绝对值变大。
所述步骤4包括:
步骤4、选取扫描方向矢量a(θ)对信号子空间Es的投影矢量bs(θ)的极大值:
vs(θ)=max|bs(θ)|=max|EsEs Ha(θ)|,-90°≤θ≤90°,
其中,||为将矢量元素取绝对值算子,max为取矢量元素极大值算子。
所述步骤5包括:
步骤5、选取扫描方向矢量a(θ)对噪声子空间En的投影矢量bn(θ)的极大值:
vn(θ)=max|bn(θ)|=max|EnEn Ha(θ)|,-90°≤θ≤90°。
所述步骤6包括如下步骤:
步骤61、根据步骤4选取的扫描方向矢量a(θ)对信号子空间Es的投影矢量bs(θ)的极大值vs(θ)和步骤5选取的扫描方向矢量a(θ)对噪声子空间En的投影矢量bn(θ)的极大值vn(θ),构建如下的新的测向目标函数:
或
步骤62、当对θ在[-90°,90°]范围内进行扫描,对于θ的每一个值,计算F(θ),画出空间谱图F(θ),-90°≤θ≤90°;空间谱图的L个峰值对应的角度是信号入射方向角度θ1、θ2、...、θL的估计值,记为是通过谱峰搜索获得的。
本发明与现有技术相比的优点在于:由上述本发明提供的技术方案可以看出,对阵列接收数据协方差矩阵进行特征分解,获得信号子空间和噪声子空间;利用信号方向矢量与信号子空间等价关系,计算扫描方向矢量向信号子空间的投影矢量;利用信号方向矢量和噪声子空间正交关系,计算扫描方向矢量向噪声子空间的投影矢量;选取扫描方向矢量向信号子空间投影矢量的极大值,选取扫描方向矢量向噪声子空间投影矢量的极大值;基于以上选取的两个极大值,构建新的测向目标函数;从形成的空间谱图中进行谱峰搜索来估计出信号入射方向角度,具有超分辨测向能力,降低了谱峰搜索的复杂性。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域的普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他附图。
图1为本发明实施例提供的比较扫描方向矢量投影极大值的阵列超分辨测向方法流程图;
图2为本发明实施例提供的阵列信号接收模型的一种示意图。为了简洁,图2只给出了M元线性阵列接收空间中1个窄带远场信号的示意图,信号的入射方向与阵列法线的角度为θ,并认为信号是以平面波的形式入射到各个阵元,右侧第1个阵元设为参考阵元,d1,d2,…,dM-1为其它阵元与参考阵元之间的间距。
具体实施方式
本发明实施例对阵列接收数据协方差矩阵进行特征分解,获得信号子空间和噪声子空间;计算扫描方向矢量向信号子空间的投影矢量;计算扫描方向矢量向噪声子空间的投影矢量;选取扫描方向矢量向信号子空间投影矢量的极大值;选取扫描方向矢量向噪声子空间投影矢量的极大值;构建新的测向目标函数,形成空间谱图,由谱峰估计出信号入射方向角度。本发明综合利用扫描方向矢量向信号子空间和噪声子空间投影矢量的极大值关系,构建新的测向目标函数,具有超高的信号方向测向分辨率,降低了谱峰搜索的复杂性。
如图1所示,本发明先对阵列接收数据进行预处理,包括:估计出阵列协方差矩阵,对该矩阵进行特征分解,获得信号子空间和噪声子空间。
主要包括如下步骤:
步骤1、根据阵列结构在测向角度范围内生成每个角度的扫描方向矢量,当该角度与信号入射方向角度相同时,称扫描方向矢量为信号方向矢量,当该角度与信号入射方向角度不同时,称扫描方向矢量为其它扫描方向矢量;
步骤2、对阵列接收数据预处理获得的信号子空间和噪声子空间,利用信号方向矢量属于信号子空间、其它扫描方向矢量不完全属于或不属于信号子空间的性质,计算扫描方向矢量向信号子空间的投影矢量;
步骤3、对阵列接收数据预处理获得的信号子空间和噪声子空间,利用信号方向矢量不属于噪声子空间、其它扫描方向矢量属于或不完全属于噪声子空间的性质,计算扫描方向矢量向噪声子空间的投影矢量;
步骤4、选取扫描方向矢量向信号子空间投影矢量的极大值;
步骤5、选取扫描方向矢量向噪声子空间投影矢量的极大值;
步骤6、根据步骤4的扫描方向矢量向信号子空间投影矢量的极大值和步骤5的扫描方向矢量向噪声子空间投影矢量的极大值,构建新的测向目标函数,在入射方向角度范围内,计算超分辨测向目标函数的值,形成空间谱图,空间谱图的谱峰所对应的角度即为信号入射方向角度的估计值。
本发明上述方案,相比较于已有的超分辨测向方法,对阵列接收数据协方差矩阵进行特征分解,获得信号子空间和噪声子空间;基于信号方向矢量与信号子空间等价关系,计算扫描方向矢量向信号子空间的投影矢量;基于信号方向矢量和噪声子空间正交关系,计算扫描方向矢量向噪声子空间的投影矢量;通过选取扫描方向矢量向信号子空间投影矢量的极大值和扫描方向矢量向噪声子空间投影矢量的极大值,构建新的测向目标函数,并从形成的空间谱图中估计出信号入射方向角度,具有超分辨测向能力。
为了便于理解,先介绍多重信号分类方法(MUSIC),然后介绍预处理,接着针对上述七个步骤做详细的说明。
本发明实例适用于任意类型的阵列形式,包括线阵、圆阵、共形阵等等,适用的波达方向包括一维方位角、一维俯仰角、二维方位角和俯仰角。为了计算方便起见,这里只针对图2给出的线阵进行讨论,具体的阵列信号模型如下:
考虑M元线阵接收空间不同入射方向角度θ1、θ2、...、θL的L(L<M)个远场窄带信号,右侧第1个阵元设为参考阵元,从右至左的其它阵元相对于参考阵元的间距分别为d1、d2、...、dM-1。由于各个信号入射方向角度不同,各个信号的平面波前就不同,到达各个阵元相对于到达参考阵元的时延也就不同,形成的信号方向矢量也就不同。则阵列在观测时间k的接收数据(称为阵列接收的第k个快拍数据)表示为:
x(k)=xs(k)+xn(k);
其中,xs(k)和xn(k)分别表示信号和噪声,sl(k)是第l个信号的波形,各信号sl(k)均为零均值且互不相关,al是第l个信号的方向矢量,/>l=1,2,…,L,xn(k)是加性独立同分布零均值白噪声,各信号sl(k)与各阵元噪声互不相关。
阵列接收信号矢量的M×M维协方差矩阵R为:
基于假设条件,Rs和A的秩均为L,因此ARsAH是秩为L的厄米特半正定矩阵,其L个非零正特征值按大小排列为μ1≥μ2≥…≥μL>0。R为厄米特正定矩阵,其M个非零正特征值按大小排列满足:
对应的特征矢量分别为e1,e2,…,eL,eL+1,…,eM,则
对于所有l'>L,由特征分解性质可得:
因此,
ARsAHel'=0,l'>L
这意味着:
aH(θl)el'=0,l=1,2,…,L,l'=L+1,L+2,…,M
上式说明最小特征值所对应的特征矢量与信号方向矢量正交。L个大特征值所对应的特征矢量e1,e2,…,eL张成一个子空间,由剩余的M-L相等的小特征值所对应的特征矢量eL+1,eL+2,…,eM张成另一个子空间。由于这两个子空间正交,故知前者是和信号有关的,称为信号子空间,记为Es;后者是信号子空间的补空间,称为噪声子空间,记为En。
建立如下函数:
当对θ进行扫描时,其L个峰值与信号入射方向角度对应。
在实际情况下,理想的阵列协方差矩阵R难以获取,只能用阵列样本协方差矩阵来代替R,阵列样本的M×M维协方差矩阵/>为
其中,K为阵列接收数据快拍数。
本发明的目的是:综合利用扫描方向矢量向信号子空间和噪声子空间投影矢量的极大值关系,构建新的测向目标函数,从形成的空间谱图中估计出信号入射方向角度,能够取得超高的测向分辨率。在预处理之后分如下七个步骤实施。
预处理:
利用M元阵列接收数据估计阵列的M×M维协方差矩阵对/>进行特征分解,其M个特征值从大到小排列为λ1,λ2,…,λL,λL+1,λL+2,…,λM,λ1,λ2,…,λL,λL+1,λL+2,…,λM对应的特征矢量分别为e1,e2,…,eL,eL+1,eL+2,…,eM,获得M×L维信号子空间Es=[e1,e2,…,eL]及其对应的L×L维信号特征值矩阵Λs=diag{λ1,λ2,…,λL},获得M×(M-L)维噪声子空间En=[eL+1,eL+2,…,eM]及其对应的(M-L)×(M-L)维噪声特征值矩阵Λn=diag{λL+1,λL+2,…,λM},Es和En是正交的,/>其中,下标s和n分别是区分信号和噪声的符号,diag{}为对角矩阵算子,上标H是矩阵的共轭转置算子。
步骤1:根据阵列结构生成对应入射方向角度的扫描方向矢量,M元阵列接收空间不同入射方向角度θ1、θ2、...、θL的L个远场窄带信号,以下简称信号,L小于M,对于测向角度范围[-90°,90°]内的每个角度θ,生成扫描方向矢量a(θ),|| ||2是矢量l2范数算子;当该角度θ与信号入射方向角度θ1、θ2、...、θL相同时,扫描方向矢量分别为a(θ1)、a(θ2)、…、a(θL),扫描方向矢量就是信号方向矢量,当该角度θ与信号入射方向角度θ1、θ2、...、θL不同时,扫描方向矢量不是信号方向矢量,称为其它扫描方向矢量;扫描方向矢量包括信号方向矢量和其它扫描方向矢量;
步骤2:对阵列接收数据预处理获得的信号子空间Es和噪声子空间En,计算扫描方向矢量a(θ)向信号子空间Es的投影矢量bs(θ):
bs(θ)=EsEs Ha(θ),-90°≤θ≤90°;
上标H是矩阵和矢量的共轭转置算子;对于信号入射方向角度θ1、θ2、...、θL,信号方向矢量a(θ1)、a(θ2)、…、a(θL)在信号子空间Es中,并与噪声子空间En正交,信号方向矢量a(θ1)、a(θ2)、…、a(θL)向信号子空间Es投影,投影矢量bs(θ1)、bs(θ2)、…、bs(θL)的元素绝对值较大;对于其它入射方向角度θ,扫描方向矢量a(θ)不完全在信号子空间Es中,甚至在噪声子空间En中,扫描方向矢量a(θ)向信号子空间Es的投影矢量bs(θ)的元素绝对值变小甚至为0。
步骤3、对阵列接收数据预处理获得的信号子空间Es和噪声子空间En,计算扫描方向矢量a(θ)向噪声子空间En的投影矢量bn(θ):
bn(θ)=EnEn Ha(θ),-90°≤θ≤90°;
对于信号入射方向角度θ1、θ2、...、θL,信号方向矢量a(θ1)、a(θ2)、…、a(θL)向噪声子空间En投影,投影矢量bn(θ1)、bn(θ2)、…、bn(θL)的元素绝对值为0;对于其它入射方向角度θ,扫描方向矢量a(θ)向噪声子空间En的投影矢量bn(θ)的元素绝对值变大。
步骤4:
选取扫描方向矢量a(θ)对信号子空间Es的投影矢量bs(θ)的极大值:
vs(θ)=max|bs(θ)|=max|EsEs Ha(θ)|,-90°≤θ≤90°,
其中,| |为将矢量元素取绝对值算子,max为取矢量元素极大值算子。
步骤5:
选取扫描方向矢量a(θ)对噪声子空间En的投影矢量bn(θ)的极大值:
vn(θ)=max|bn(θ)|=max|EnEn Ha(θ)|,-90°≤θ≤90°。
步骤6:
步骤61、根据步骤4选取的扫描方向矢量a(θ)对信号子空间Es的投影矢量bs(θ)的极大值vs(θ)和步骤5选取的扫描方向矢量a(θ)对噪声子空间En的投影矢量bn(θ)的极大值vn(θ),构建如下的新的测向目标函数:
或
步骤62、当对θ在[-90°,90°]范围内进行扫描,对于θ的每一个值,计算F(θ),画出空间谱图F(θ),-90°≤θ≤90°;空间谱图的L个峰值对应的角度是信号入射方向角度θ1、θ2、...、θL的估计值,记为是通过谱峰搜索获得的。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到上述实施例可以通过软件实现,也可以借助软件加必要的通用硬件平台的方式来实现。基于这样的理解,上述实施例的技术方案可以以软件产品的形式体现出来,该软件产品可以存储在一个非易失性存储介质(可以是CD-ROM,U盘,移动硬盘等)中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述的方法。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明披露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求书的保护范围为准。
Claims (7)
1.一种比较扫描方向矢量投影极大值的阵列超分辨测向方法,其特征在于:对阵列接收数据进行预处理,获得信号子空间和噪声子空间,根据阵列结构在测向角度范围内生成每个角度的扫描方向矢量,并利用信号方向矢量与信号子空间等价关系、信号方向矢量与噪声子空间正交关系,分别计算扫描方向矢量向信号子空间的投影矢量和扫描方向矢量向噪声子空间的投影矢量,进一步选取扫描方向矢量向信号子空间投影矢量的极大值和扫描方向矢量向噪声子空间投影矢量的极大值,并构建新的测向目标函数,最后从形成的空间谱图中估计出信号入射方向角度,具有超分辨测向能力;
具体包括如下步骤:
步骤1、根据阵列结构在测向角度范围内生成每个角度的扫描方向矢量,当该角度与信号入射方向角度相同时,称扫描方向矢量为信号方向矢量,当该角度与信号入射方向角度不同时,称扫描方向矢量为其它扫描方向矢量;
步骤2、对阵列接收数据预处理获得的信号子空间和噪声子空间,利用信号方向矢量属于信号子空间、其它扫描方向矢量不完全属于或不属于信号子空间的性质,计算扫描方向矢量向信号子空间的投影矢量;
步骤3、对阵列接收数据预处理获得的信号子空间和噪声子空间,利用信号方向矢量不属于噪声子空间、其它扫描方向矢量属于或不完全属于噪声子空间的性质,计算扫描方向矢量向噪声子空间的投影矢量;
步骤4、选取扫描方向矢量向信号子空间投影矢量的极大值;
步骤5、选取扫描方向矢量向噪声子空间投影矢量的极大值;
步骤6、根据步骤4的扫描方向矢量向信号子空间投影矢量的极大值和步骤5的扫描方向矢量向噪声子空间投影矢量的极大值,构建新的测向目标函数,在入射方向角度范围内,计算超分辨测向目标函数的值,形成空间谱图,空间谱图的谱峰所对应的角度即为信号入射方向角度的估计值。
2.根据权利要求1所述的比较扫描方向矢量投影极大值的阵列超分辨测向方法,其特征在于:所述步骤1包括:
步骤1、根据阵列结构生成对应入射方向角度的扫描方向矢量,M元阵列接收空间不同入射方向角度θ1、θ2、...、θL的L个远场窄带信号,以下简称信号,L小于M,对于测向角度范围[-90°,90°]内的每个角度θ,生成扫描方向矢量a(θ),|| ||2是矢量l2范数算子;当该角度θ与信号入射方向角度θ1、θ2、...、θL相同时,扫描方向矢量分别为a(θ1)、a(θ2)、…、a(θL),扫描方向矢量就是信号方向矢量,当该角度θ与信号入射方向角度θ1、θ2、...、θL不同时,扫描方向矢量不是信号方向矢量,称为其它扫描方向矢量;扫描方向矢量包括信号方向矢量和其它扫描方向矢量。
3.根据权利要求1所述的比较扫描方向矢量投影极大值的阵列超分辨测向方法,其特征在于:所述步骤2包括:
步骤2、对阵列接收数据预处理获得的信号子空间Es和噪声子空间En,计算扫描方向矢量a(θ)向信号子空间Es的投影矢量bs(θ):
bs(θ)=EsEs Ha(θ),-90°≤θ≤90°;
上标H是矩阵和矢量的共轭转置算子;对于信号入射方向角度θ1、θ2、...、θL,信号方向矢量a(θ1)、a(θ2)、…、a(θL)在信号子空间Es中,并与噪声子空间En正交,信号方向矢量a(θ1)、a(θ2)、…、a(θL)向信号子空间Es投影,投影矢量bs(θ1)、bs(θ2)、…、bs(θL)的元素绝对值较大;对于其它入射方向角度θ,扫描方向矢量a(θ)不完全在信号子空间Es中,甚至在噪声子空间En中,扫描方向矢量a(θ)向信号子空间Es的投影矢量bs(θ)的元素绝对值变小甚至为0。
4.根据权利要求1所述的比较扫描方向矢量投影极大值的阵列超分辨测向方法,其特征在于:所述步骤3包括:
步骤3、对阵列接收数据预处理获得的信号子空间Es和噪声子空间En,计算扫描方向矢量a(θ)向噪声子空间En的投影矢量bn(θ):
bn(θ)=EnEn Ha(θ),-90°≤θ≤90°;
对于信号入射方向角度θ1、θ2、...、θL,信号方向矢量a(θ1)、a(θ2)、…、a(θL)向噪声子空间En投影,投影矢量bn(θ1)、bn(θ2)、…、bn(θL)的元素绝对值为0;对于其它入射方向角度θ,扫描方向矢量a(θ)向噪声子空间En的投影矢量bn(θ)的元素绝对值变大。
5.根据权利要求1所述的比较扫描方向矢量投影极大值的阵列超分辨测向方法,其特征在于:所述步骤4包括:
步骤4、选取扫描方向矢量a(θ)对信号子空间Es的投影矢量bs(θ)的极大值:
vs(θ)=max|bs(θ)|=max|EsEs Ha(θ)|,-90°≤θ≤90°,
其中,||为将矢量元素取绝对值算子,max为取矢量元素极大值算子。
6.根据权利要求1所述的比较扫描方向矢量投影极大值的阵列超分辨测向方法,其特征在于:所述步骤5包括:
步骤5、选取扫描方向矢量a(θ)对噪声子空间En的投影矢量bn(θ)的极大值:
vn(θ)=max|bn(θ)|=max|EnEn Ha(θ)|,-90°≤θ≤90°。
7.根据权利要求1所述的比较扫描方向矢量投影极大值的阵列超分辨测向方法,其特征在于:所述步骤6包括:
步骤61、根据步骤4选取的扫描方向矢量a(θ)对信号子空间Es的投影矢量bs(θ)的极大值vs(θ)和步骤5选取的扫描方向矢量a(θ)对噪声子空间En的投影矢量bn(θ)的极大值vn(θ),构建如下的新的测向目标函数:
或
步骤62、当对θ在[-90°,90°]范围内进行扫描,对于θ的每一个值,计算F(θ),画出空间谱图F(θ),-90°≤θ≤90°;空间谱图的L个峰值对应的角度是信号入射方向角度θ1、θ2、...、θL的估计值,记为是通过谱峰搜索获得的。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210226222.XA CN114609581B (zh) | 2022-03-07 | 2022-03-07 | 一种比较扫描方向矢量投影极大值的阵列超分辨测向方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210226222.XA CN114609581B (zh) | 2022-03-07 | 2022-03-07 | 一种比较扫描方向矢量投影极大值的阵列超分辨测向方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114609581A CN114609581A (zh) | 2022-06-10 |
CN114609581B true CN114609581B (zh) | 2024-05-10 |
Family
ID=81860241
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210226222.XA Active CN114609581B (zh) | 2022-03-07 | 2022-03-07 | 一种比较扫描方向矢量投影极大值的阵列超分辨测向方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114609581B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2008106386A (ru) * | 2008-02-21 | 2009-08-27 | Государственное образовательное учреждение высшего профессионального образования "Московский государственный технический университ | Способ пеленгации источников радиоизлучения на одной частоте |
CN109633525A (zh) * | 2018-11-29 | 2019-04-16 | 上海无线电设备研究所 | 一种快速的均匀圆阵空间谱测向实现方法 |
CN110673086A (zh) * | 2019-10-31 | 2020-01-10 | 上海无线电设备研究所 | 一种基于数字阵列雷达的二维角度超分辨方法 |
CN114047481A (zh) * | 2021-11-09 | 2022-02-15 | 中国科学技术大学 | 一种基于子空间正交性的稳健自适应波束形成方法 |
CN114325565A (zh) * | 2021-12-27 | 2022-04-12 | 中国科学技术大学 | 一种基于子空间关系的阵列超分辨测向方法 |
-
2022
- 2022-03-07 CN CN202210226222.XA patent/CN114609581B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2008106386A (ru) * | 2008-02-21 | 2009-08-27 | Государственное образовательное учреждение высшего профессионального образования "Московский государственный технический университ | Способ пеленгации источников радиоизлучения на одной частоте |
CN109633525A (zh) * | 2018-11-29 | 2019-04-16 | 上海无线电设备研究所 | 一种快速的均匀圆阵空间谱测向实现方法 |
CN110673086A (zh) * | 2019-10-31 | 2020-01-10 | 上海无线电设备研究所 | 一种基于数字阵列雷达的二维角度超分辨方法 |
CN114047481A (zh) * | 2021-11-09 | 2022-02-15 | 中国科学技术大学 | 一种基于子空间正交性的稳健自适应波束形成方法 |
CN114325565A (zh) * | 2021-12-27 | 2022-04-12 | 中国科学技术大学 | 一种基于子空间关系的阵列超分辨测向方法 |
Non-Patent Citations (2)
Title |
---|
利用局域子空间投影提高子空间类DOA估计算法的谱分辨力;王布宏, 王永良, 陈辉;电子学报;20030325(第03期);全文 * |
提高空间谱算法分辨率的DOA估计;司伟建;蓝晓宇;刘学;;应用科学学报;20130730(第04期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN114609581A (zh) | 2022-06-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109633522B (zh) | 基于改进的music算法的波达方向估计方法 | |
US6653973B2 (en) | Adaptive digital beamforming radar method and system for maintaining multiple source angle super-resolution capability in jamming | |
US8669901B2 (en) | Method for determining azimuth and elevation angles of arrival of coherent sources | |
CN107092007A (zh) | 一种虚拟二阶阵列扩展的波达方向估计方法 | |
CN111580039A (zh) | 基于非圆信号的单基地展开互质阵列mimo雷达doa估计方法 | |
CN111693947A (zh) | 基于互质阵列doa估计的改进music方法 | |
US9444558B1 (en) | Synthetic robust adaptive beamforming | |
CN114325565B (zh) | 一种基于子空间关系的阵列超分辨测向方法 | |
Nickel | Radar target parameter estimation with array antennas | |
CN112363108A (zh) | 信号子空间加权超分辨的波达方向检测方法及系统 | |
CN114609581B (zh) | 一种比较扫描方向矢量投影极大值的阵列超分辨测向方法 | |
CN109061597B (zh) | 基于盲源分离与时频脊波域滤波的电离层杂波抑制方法 | |
CN113075613B (zh) | 一种卫星移动通信系统中的干扰测向方法 | |
CN106199550B (zh) | 一种波离角和二维波达角自动配对联合估计方法 | |
CN106844886B (zh) | 基于主分量分析的目标波达方向获取方法 | |
CN109861770B (zh) | 一种基于波束形成输出功率组合的宽带信号检测方法 | |
CN112731279B (zh) | 一种基于混合天线子阵列的到达角估计方法 | |
CN114265058A (zh) | Mimo雷达目标测角方法、装置、电子设备及存储介质 | |
CN112698263A (zh) | 一种基于正交传播算子的单基地互质mimo阵列doa估计算法 | |
CN115453454A (zh) | 基于扫描方向矢量与特征矢量内积的阵列超分辨测向方法 | |
Raia et al. | Performance of a hybrid non-contiguous planar array: effect of phase shifter quantization | |
Kazarinov et al. | The Experimental Research of DOA Estimation Based on Difference Co-array Method | |
CN111157968B (zh) | 一种基于稀疏mimo阵列的高精度测角方法 | |
CN109752688B (zh) | 一种针对传感器阵列系统的临近信源角度差值计算方法 | |
Tayem et al. | A unitary MUSIC-like algorithm for coherent sources |
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 |