CN109358312A - 确定入射信号来波方向的方法、装置、介质及设备 - Google Patents
确定入射信号来波方向的方法、装置、介质及设备 Download PDFInfo
- Publication number
- CN109358312A CN109358312A CN201811348931.5A CN201811348931A CN109358312A CN 109358312 A CN109358312 A CN 109358312A CN 201811348931 A CN201811348931 A CN 201811348931A CN 109358312 A CN109358312 A CN 109358312A
- Authority
- CN
- China
- Prior art keywords
- mutual coupling
- array
- aerial array
- spectral function
- covariance matrix
- 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
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)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明提供了一种确定入射信号来波方向的方法,基于等离子鞘环境,包括:建立等离子鞘环境下的天线阵列;采用空间平滑算法,对所述天线阵列进行解相关,获得所述天线阵列修正后的协方差矩阵;采用互耦校准算法,根据修正后的协方差矩阵,对所述天线阵列的空间谱函数进行处理,获得不含互耦系数的空间谱函数;搜索所述不含互耦系数的空间谱函数的极大值,确定入射信号的来波方向。通过结合空间平滑算法和互耦校准算法,能够得到不含互耦系数的空间谱函数,再通过不含互耦系数的空间谱函数确定入射信号的来波方向,能够减弱多径效应与天线耦合这两种误差对入射信号来波方向确定的影响,进而能估算出更加准确的来波方向。
Description
技术领域
本发明涉及卫星导航技术领域,具体涉及一种确定入射信号来波方向的方法、装置、介质及设备。
背景技术
飞行器重返大气层至落地的过程称为再入过程。再入过程中,高超声速飞行器与大气相互作用,在飞行器前形成弓形脱体激波,激波后气体温度急剧升高呈现“高温效应”,空气分子发生离解、电离反应;同时,飞行器表面的防热烧蚀材料在高温、高热流作用下烧蚀,将烧蚀产物释放到飞行器周围流场中,形成了成分多样的复杂的再入绕流流场。再入绕流流场是一团电离气体,当气体的电离度达到一定程度时,电离气体表现出等离子体所具有的集体行为而成为等离子体,此时再入绕流流场亦称为等离子体包覆流场、再入等离子体或等离子体鞘套。再入等离子体的形成,使得电磁波传播的功率衰减,且使电磁波产生反射、折射及散射。
信号穿过等离子体鞘后在飞行器表面和等离子体鞘之间产生多次反射,从而引起多径效应。这里的多径效应是指GPS接收机除直接接收到卫星信号外,尚可能接收到经等离子体鞘发射后,经过不同的传播路径到达天线信号。从而形成的定位误差,多径效应会降低接收机的定位精度,并使处理时间加长,在WAAS(Wide Area Augmentation System)和DGPS(Difference Global Positioning System)***中,多径效应误差将导致误差传递,严重降低***中用户的定位精度。
等离子鞘套的厚度和电子密度影响天线间的耦合特性,天线间的互耦影响着阵列天线的工作性能,因此,利用波束形成技术减弱等离子鞘的影响需考虑对耦合误差进行校正。在现有的等离子鞘模型中还没有可以校准阵列互耦的方法,如果能有效的校准等离子体鞘产生的多径效应和互耦影响,将为消除多径和耦合效应的影响起到重要的作用。
发明内容
针对现有技术中的缺陷,本发明提供一种确定入射信号来波方向的方法、装置、介质及设备,能够减弱多径效应与天线耦合这两种误差对入射信号来波方向确定的影响,进而能估算出更加准确的来波方向。
第一方面,本发明提供了一种确定入射信号来波方向的方法,基于等离子鞘环境,包括:
建立等离子鞘环境下的天线阵列;
采用空间平滑算法,对所述天线阵列进行解相关,获得所述天线阵列修正后的协方差矩阵;
采用互耦校准算法,根据修正后的协方差矩阵,对所述天线阵列的空间谱函数进行处理,获得不含互耦系数的空间谱函数;
搜索所述不含互耦系数的空间谱函数的极大值,确定入射信号的来波方向。
可选的,所述建立等离子鞘环境下的天线阵列,包括:
设等离子鞘环境下,M个天线阵元均匀排列,从远场θk(k=1,2,…N)处入射N个窄带信源(波长为λ),建立天线阵列,为:
X(t)=A(θ)S(t)+N(t)
其中,A(θ)是指方向矢量的阵列流形,S(t)是指信号,N(t)是指噪声。
可选的,所述采用空间平滑算法,对所述天线阵列进行解相关,获得所述天线阵列修正后的协方差矩阵,包括:
对所述阵列流形进行修正;
根据修正后的阵列流形,获得入射信号的协方差矩阵;
采用空间平滑算法,基于所述协方差矩阵,获得前后向平滑协方差矩阵;
根据所述前后向平滑协方差矩阵,修正所述协方差矩阵,获得所述天线阵列修正后的协方差矩阵。
可选的,在所述采用互耦校准算法,根据修正后的协方差矩阵,对所述天线阵列的空间谱函数进行处理,获得不含互耦系数的空间谱函数的步骤之前,还包括:
采用子空间理论,获得所述天线阵列的空间谱函数。
可选的,所述对所述阵列流形进行修正,包括:
阵列流形中的每个方向矢量可以修正为:
其中,C是均匀线阵的互耦矩阵,它是一个带状、对称Toeplitz矩阵;
设阵元间的间距为d=λ/2,互耦自由度为p(即两个阵元间间距大于d(p-1)时,互耦系数为0),其中,第一行互耦系数满足关系0<|cp-1|<…<c1<c0=1。
可选的,所述采用互耦校准算法,根据修正后的协方差矩阵,对所述天线阵列的空间谱函数进行处理,获得不含互耦系数的空间谱函数,包括:
根据所述互耦矩阵的Toeplitz性质,获得新互耦矩阵;
根据所述新互耦矩阵,获得新阵列流形;
根据所述修正后的协方差矩阵和所述新阵列流形,获得不含互耦系数的空间谱函数。
第二方面,本发明提供一种确定入射信号来波方向的装置,基于等离子鞘环境,包括:
阵列建立模块,用于建立等离子鞘环境下的天线阵列;
解相关模块,用于采用空间平滑算法,对所述天线阵列进行解相关,获得所述天线阵列修正后的协方差矩阵;
盲校准模块,用于采用互耦校准算法,根据修正后的协方差矩阵,对所述天线阵列的空间谱函数进行处理,获得不含互耦系数的空间谱函数;
方向确定模块,用于搜索所述不含互耦系数的空间谱函数的极大值,确定入射信号的来波方向。
可选的,所述阵列建立模块,具体用于:
设等离子鞘环境下,M个天线阵元均匀排列,从远场θk(k=1,2,…N)处入射N个窄带信源(波长为λ),建立天线阵列,为:
X(t)=A(θ)S(t)+N(t)
其中,A(θ)是指方向矢量的阵列流形,S(t)是指信号,N(t)是指噪声。
第三方面,本发明提供一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现上述一种确定入射信号来波方向的方法。
第四方面,本发明提供一种确定入射信号来波方向的处理设备,基于等离子鞘环境,包括:存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现上述一种确定入射信号来波方向的方法。
本发明通过结合空间平滑算法和互耦校准算法,能够得到不含互耦系数的空间谱函数,再通过不含互耦系数的空间谱函数确定入射信号的来波方向,能够减弱多径效应与天线耦合这两种误差对入射信号来波方向确定的影响,进而能估算出更加准确的来波方向。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案,下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍。在所有附图中,类似的元件或部分一般由类似的附图标记标识。附图中,各元件或部分并不一定按照实际的比例绘制。
图1为本发明提供的一种确定入射信号来波方向的方法的流程图;
图2为本发明提供的三种方法计算的空间谱曲线的对比图;
图3为本发明提供的未校准的阵列天线三维方向图;
图4为本发明提供的未校准的阵列天线二维方向图;
图5为本发明提供的经过校准的阵列天线三维方向图;
图6为本发明提供的经过校准的阵列天线二维方向图;
图7为本发明提供的一种确定入射信号来波方向的装置的示意图。
具体实施方式
下面将结合附图对本发明技术方案的实施例进行详细的描述。以下实施例仅用于更加清楚地说明本发明的技术方案,因此只是作为示例,而不能以此来限制本发明的保护范围。
需要注意的是,除非另有说明,本申请使用的技术术语或者科学术语应当为本发明所属领域技术人员所理解的通常意义。
本发明提供了一种确定入射信号来波方向的方法、装置、介质及设备。下面结合附图对本发明的实施例进行说明。
请参考图1,图1为本发明具体实施例提供的一种确定入射信号来波方向的方法的流程图,本实施例提供的一种确定入射信号来波方向的方法,基于等离子鞘环境,包括:
步骤S101:建立等离子鞘环境下的天线阵列。
步骤S102:采用空间平滑算法,对所述天线阵列进行解相关,获得所述天线阵列修正后的协方差矩阵。
步骤S103:采用互耦校准算法,根据修正后的协方差矩阵,对所述天线阵列的空间谱函数进行处理,获得不含互耦系数的空间谱函数。
步骤S104:搜索所述不含互耦系数的空间谱函数的极大值,确定入射信号的来波方向。
通过结合空间平滑算法和互耦校准算法,能够得到不含互耦系数的空间谱函数,再通过不含互耦系数的空间谱函数确定入射信号的来波方向,能够减弱多径效应与天线耦合这两种误差对入射信号来波方向确定的影响,进而能估算出更加准确的来波方向。
在本发明提供的一个具体实施例中,所述建立等离子鞘环境下的天线阵列,包括:设等离子鞘环境下,M个天线阵元均匀排列,从远场θk(k=1,2,…N)处入射N个窄带信源(波长为λ),建立天线阵列,为:
X(t)=A(θ)S(t)+N(t) (1)
其中,A(θ)是指方向矢量的阵列流形,S(t)是指信号,N(t)是指噪声。
在建立天线阵列时,可以通过HFSS仿真出等离子鞘环境下的天线阵列。
在本发明提供的一个具体实施例中,所述采用空间平滑算法,对所述天线阵列进行解相关,获得所述天线阵列修正后的协方差矩阵,包括:对所述阵列流形进行修正;根据修正后的阵列流形,获得入射信号的协方差矩阵;采用空间平滑算法,基于所述协方差矩阵,获得前后向平滑协方差矩阵;根据所述前后向平滑协方差矩阵,修正所述协方差矩阵,获得所述天线阵列修正后的协方差矩阵。
对阵列流形的修正:
当存在互耦误差时,阵列流形中的每个方向矢量可以修正为:
其中,C是均匀线阵的互耦矩阵,它是一个带状、对称Toeplitz矩阵。表示修正后的方向矢量,a(θ)表示修正前的方向矢量。
设阵元间的间距为d=λ/2。互耦自由度为p(即两个阵元间间距大于d(p-1)时,互耦系数为0)。其中,第一行互耦系数满足关系0<|cp-1|<…<c1<c0=1。
根据修正后的阵列流形,获得入射信号的协方差矩阵:
根据公式(2),可以得到阵列输出的协方差矩阵,具体为:
RXX=E[X(t)X(t)H]=CARSAHCH+σ2I (4)
其中,Rs表示入射信号的协方差矩阵;σ表示噪声功率;I表示单位矩阵。
采用空间平滑算法,基于所述协方差矩阵,获得前后向平滑协方差矩阵,具体过程为:
将M元的等距天线矩阵用滑动方式分成L个子阵,每个子阵有K个单元,其中K=M-L+1。那么,第l个向前子阵的输出为:
其中,
AM表示子阵的方向矢量阵列流形;s(t)表示子阵的输入信号;nl(t)表示子阵的噪声。
所以,第l个前向子阵的协方差矩阵为:
定义前向空间平滑协方差矩阵为:
同理,可得后向空间平滑协方差矩阵为:
因此,定义前后向平滑协方差矩阵为:
根据所述前后向平滑协方差矩阵,修正所述协方差矩阵,获得所述天线阵列修正后的协方差矩阵,具体过程为:
通常,Rfb只是Hermite矩阵,并不是Toeplitz矩阵。利用Toeplitz性质,对Rfb进行修正,RfbX是Rfb的无偏估计,因此,对公式(4)进行空间平滑后得:
其中,
在本发明提供的一个具体实施例中,在所述采用互耦校准算法,根据修正后的协方差矩阵,对所述天线阵列的空间谱函数进行处理,获得不含互耦系数的空间谱函数的步骤之前,还包括:采用子空间理论,获得所述天线阵列的空间谱函数。
根据子空间理论:
得到阵列空间谱函数PMUSIC(θ)为
其中,UN为噪声子空间。
对公式(13)中的角度进行变化,通过寻找波峰来确定入射信号的来波方向。
在本发明提供的一个具体实施例中,所述采用互耦校准算法,根据修正后的协方差矩阵,对所述天线阵列的空间谱函数进行处理,获得不含互耦系数的空间谱函数,包括:根据所述互耦矩阵的Toeplitz性质,获得新互耦矩阵;根据所述新互耦矩阵,获得新阵列流形;根据所述修正后的协方差矩阵和所述新阵列流形,获得不含互耦系数的空间谱函数。
根据所述互耦矩阵的Toeplitz性质,获得新互耦矩阵,具体过程为:
由于C是一个对称Toeplitz矩阵,从构成互耦矩阵的每个互耦系数入手,将C表示为:
得到新互耦矩阵,其中,cq为互耦系数,这时,阵列流形可以表示为新阵列流形,为:
其中,c=[c0,c1,c2,…,cp-1]T是互耦系数矩阵,W(θ)为M×p维由导向矢量组成的变换矩阵。
根据所述修正后的协方差矩阵和所述新阵列流形,获得不含互耦系数的空间谱函数,具体过程为:
对公式(10)进行特征值分解,得到噪声子空间UN,代入公式(13),得:
PMUSIC(θ)=1/||(Ca(θ))HUN||2 (16)
将公式(15)代入公式(16),得:
PMUSIC(θ)=1/||(W(θ)c)HUN||2
=1/(cHW(θ)HUNUN HW(θ)c)
=1/(cHΦ(θ)c) (17)
其中,Φ(θ)=W(θ)HUNUN HW(θ),与互耦系数不相关。
求PMUSIC(θ)的极大值问题可以转化为线性约束下的二次规划问题。即
mincHΦ(θ)c
s.t.εTc=1 (18)
其中,cε=1,ε=[1,0,…,0]T。
用拉格朗日求解法求得,
代入空间谱函数公式(17)得,
PMUSIC(θ)=εTΦ-1(θ)ε (20)
得到不含互耦系数的空间谱函数。利用构成的不含互耦系数的空间谱函数,搜索其极大值位置,确定入射信号的来波方向。
在使用本发明时,可以使用MATLAB仿真本发明的方法,验证该方法的准确性。具体过程为:使用HFSS仿真等离子鞘环境下的阵列模型,阵元数M=8,等离子体碰撞频率v=2.1e10Hz,电子密度Ne=1.6e12cm-3,厚度δ=0.2cm。所加源为GPS信号中心频率1.57542GHz,阵元间距为入射信号的半波长。在使用MATLAB仿真本发明时,互耦自由度选为3,互耦系数选取HFSS仿真的互耦参数为[1,-0.2232-0.3432i,-0.1143+0.1456i]T。噪声为零均值高斯白噪声,2个相干信号源,方位角分别为0和65,实验次数为100次,快拍数为200,信噪比为10dB。
结果表明,阵列互耦存在时,常规的空间平滑算法无法得到正确的谱峰而失效,本发明方法能够准确地分辨出来波方向为0和65。如图2所示。
通过本发明确定出阵列耦合的校准信息,即每个阵元需调节的幅度和相位。将校准信息加载到HFSS的阵列模型上,得到校准结果。
经过耦合校准之后的阵列天线得到了明显的改善,如图3-6所示。
本发明基于空间平滑的解相干处理技术,通过盲校准法校正阵列的耦合影响并估计出准确的波达方向。适合于卫星导航***中消除等离子鞘产生的多径效应和耦合的影响,提高卫星的定位精度。
以上,为本发明提供的一种确定入射信号来波方向的方法。
基于与上述一种确定入射信号来波方向的方法相同的发明构思,与之相对应的,本发明实施例还提供了一种确定入射信号来波方向的装置,如图2所示。由于装置实施例基本相似与方法实施例,所以描述得比较简单,相关之处参见方法实施例的部分说明即可。
本发明提供的一种确定入射信号来波方向的装置,基于等离子鞘环境,包括:
阵列建立模块101,用于建立等离子鞘环境下的天线阵列;
解相关模块102,用于采用空间平滑算法,对所述天线阵列进行解相关,获得所述天线阵列修正后的协方差矩阵;
盲校准模块103,用于采用互耦校准算法,根据修正后的协方差矩阵,对所述天线阵列的空间谱函数进行处理,获得不含互耦系数的空间谱函数;
方向确定模块104,用于搜索所述不含互耦系数的空间谱函数的极大值,确定入射信号的来波方向。
在本发明提供的一个具体实施例中,所述阵列建立模块101,具体用于:
设等离子鞘环境下,M个天线阵元均匀排列,从远场θk(k=1,2,…N)处入射N个窄带信源(波长为λ),建立天线阵列,为:
X(t)=A(θ)S(t)+N(t)
其中,A(θ)是指方向矢量的阵列流形,S(t)是指信号,N(t)是指噪声。
在本发明提供的一个具体实施例中,所述解相关模块102,包括:
阵列流形修正单元,用于对所述阵列流形进行修正;
协方差矩阵获得单元,用于根据修正后的阵列流形,获得入射信号的协方差矩阵;
平滑计算单元,用于采用空间平滑算法,基于所述协方差矩阵,获得前后向平滑协方差矩阵;
修正矩阵单元,用于根据所述前后向平滑协方差矩阵,修正所述协方差矩阵,获得所述天线阵列修正后的协方差矩阵。
在本发明提供的一个具体实施例中,所述装置,还包括:
谱函数计算模块,用于采用子空间理论,获得所述天线阵列的空间谱函数。
在本发明提供的一个具体实施例中,所述阵列流形修正单元,具体用于:
阵列流形中的每个方向矢量可以修正为:
其中,C是均匀线阵的互耦矩阵,它是一个带状、对称Toeplitz矩阵;
设阵元间的间距为d=λ/2,互耦自由度为p(即两个阵元间间距大于d(p-1)时,互耦系数为0),其中,第一行互耦系数满足关系0<|cp-1|<…<c1<c0=1。
在本发明提供的一个具体实施例中,所述盲校准模块103,包括:
性质分析单元,用于根据所述互耦矩阵的Toeplitz性质,获得新互耦矩阵;
新阵列流形获得单元,用于根据所述新互耦矩阵,获得新阵列流形;
谱函数计算单元,用于根据所述修正后的协方差矩阵和所述新阵列流形,获得不含互耦系数的空间谱函数。
以上,为本发明提供的一种确定入射信号来波方向的装置。
基于与上述一种确定入射信号来波方向的方法相同的发明构思,与之相对应的,本发明实施例还提供了一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现上述一种确定入射信号来波方向的方法。
由上述技术方案可知,本实施例提供的一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时,能够通过结合空间平滑算法和互耦校准算法,得到不含互耦系数的空间谱函数,再通过不含互耦系数的空间谱函数确定入射信号的来波方向,减弱多径效应与天线耦合这两种误差对入射信号来波方向确定的影响,进而能估算出更加准确的来波方向。
基于与上述一种确定入射信号来波方向的方法相同的发明构思,与之相对应的,本发明实施例还提供了一种确定入射信号来波方向的处理设备,包括:存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现上述一种确定入射信号来波方向的方法。
由上述技术方案可知,本实施例提供的一种确定入射信号来波方向的处理设备,通过结合空间平滑算法和互耦校准算法,能够得到不含互耦系数的空间谱函数,再通过不含互耦系数的空间谱函数确定入射信号的来波方向,能够减弱多径效应与天线耦合这两种误差对入射信号来波方向确定的影响,进而能估算出更加准确的来波方向。
本发明的说明书中,说明了大量具体细节。然而,能够理解,本发明的实施例可以在没有这些具体细节的情况下实践。在一些实例中,并未详细示出公知的方法、结构和技术,以便不模糊对本说明书的理解。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外,在不相互矛盾的情况下,本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
最后应说明的是:以上各实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述各实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的范围,其均应涵盖在本发明的权利要求和说明书的范围当中。
Claims (10)
1.一种确定入射信号来波方向的方法,基于等离子鞘环境,其特征在于,包括:
建立等离子鞘环境下的天线阵列;
采用空间平滑算法,对所述天线阵列进行解相关,获得所述天线阵列修正后的协方差矩阵;
采用互耦校准算法,根据修正后的协方差矩阵,对所述天线阵列的空间谱函数进行处理,获得不含互耦系数的空间谱函数;
搜索所述不含互耦系数的空间谱函数的极大值,确定入射信号的来波方向。
2.根据权利要求1所述的方法,其特征在于,所述建立等离子鞘环境下的天线阵列,包括:
设等离子鞘环境下,M个天线阵元均匀排列,从远场θk(k=1,2,…N)处入射N个窄带信源(波长为λ),建立天线阵列,为:
X(t)=A(θ)S(t)+N(t)
其中,A(θ)是指方向矢量的阵列流形,S(t)是指信号,N(t)是指噪声。
3.根据权利要求2所述的方法,其特征在于,所述采用空间平滑算法,对所述天线阵列进行解相关,获得所述天线阵列修正后的协方差矩阵,包括:
对所述阵列流形进行修正;
根据修正后的阵列流形,获得入射信号的协方差矩阵;
采用空间平滑算法,基于所述协方差矩阵,获得前后向平滑协方差矩阵;
根据所述前后向平滑协方差矩阵,修正所述协方差矩阵,获得所述天线阵列修正后的协方差矩阵。
4.根据权利要求1所述的方法,其特征在于,在所述采用互耦校准算法,根据修正后的协方差矩阵,对所述天线阵列的空间谱函数进行处理,获得不含互耦系数的空间谱函数的步骤之前,还包括:
采用子空间理论,获得所述天线阵列的空间谱函数。
5.根据权利要求3所述的方法,其特征在于,所述对所述阵列流形进行修正,包括:
阵列流形中的每个方向矢量可以修正为:
其中,C是均匀线阵的互耦矩阵,它是一个带状、对称Toeplitz矩阵;
设阵元间的间距为d=λ/2,互耦自由度为p(即两个阵元间间距大于d(p-1)时,互耦系数为0),其中,第一行互耦系数满足关系0<|cp-1|<…<c1<c0=1。
6.根据权利要求5所述的方法,其特征在于,所述采用互耦校准算法,根据修正后的协方差矩阵,对所述天线阵列的空间谱函数进行处理,获得不含互耦系数的空间谱函数,包括:
根据所述互耦矩阵的Toeplitz性质,获得新互耦矩阵;
根据所述新互耦矩阵,获得新阵列流形;
根据所述修正后的协方差矩阵和所述新阵列流形,获得不含互耦系数的空间谱函数。
7.一种确定入射信号来波方向的装置,基于等离子鞘环境,其特征在于,包括:
阵列建立模块,用于建立等离子鞘环境下的天线阵列;
解相关模块,用于采用空间平滑算法,对所述天线阵列进行解相关,获得所述天线阵列修正后的协方差矩阵;
盲校准模块,用于采用互耦校准算法,根据修正后的协方差矩阵,对所述天线阵列的空间谱函数进行处理,获得不含互耦系数的空间谱函数;
方向确定模块,用于搜索所述不含互耦系数的空间谱函数的极大值,确定入射信号的来波方向。
8.根据权利要求7所述的装置,其特征在于,所述阵列建立模块,具体用于:
设等离子鞘环境下,M个天线阵元均匀排列,从远场θk(k=1,2,…N)处入射N个窄带信源(波长为λ),建立天线阵列,为:
X(t)=A(θ)S(t)+N(t)
其中,A(θ)是指方向矢量的阵列流形,S(t)是指信号,N(t)是指噪声。
9.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,该程序被处理器执行时实现权利要求1-6之一所述的方法。
10.一种确定入射信号来波方向的处理设备,基于等离子鞘环境,包括:存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时实现权利要求1-6之一所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811348931.5A CN109358312B (zh) | 2018-11-13 | 2018-11-13 | 确定入射信号来波方向的方法、装置、介质及设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811348931.5A CN109358312B (zh) | 2018-11-13 | 2018-11-13 | 确定入射信号来波方向的方法、装置、介质及设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109358312A true CN109358312A (zh) | 2019-02-19 |
CN109358312B CN109358312B (zh) | 2023-05-26 |
Family
ID=65345003
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811348931.5A Active CN109358312B (zh) | 2018-11-13 | 2018-11-13 | 确定入射信号来波方向的方法、装置、介质及设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109358312B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110927659A (zh) * | 2019-11-25 | 2020-03-27 | 长江大学 | 互耦条件下任意阵列流形doa估计与互耦校准方法及*** |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101149429A (zh) * | 2006-09-18 | 2008-03-26 | 中国人民解放军空军工程大学导弹学院 | 适用于均匀圆阵的阵列互耦校正与信源测向方法 |
CN104931931A (zh) * | 2015-05-18 | 2015-09-23 | 哈尔滨工程大学 | 互耦条件下基于张量实值子空间的双基地mimo雷达角度估计方法 |
CN105974370A (zh) * | 2016-06-30 | 2016-09-28 | 哈尔滨工业大学 | 一种基于虚拟阵元空间平滑互耦校正的music空间谱估计方法 |
CN106154217A (zh) * | 2016-07-12 | 2016-11-23 | 南京邮电大学 | Ula和uca中互耦未知时基于空间谱伪峰消除的自校准方法 |
CN106324558A (zh) * | 2016-08-30 | 2017-01-11 | 东北大学秦皇岛分校 | 基于互质阵列的宽带信号doa估计方法 |
CN106443610A (zh) * | 2016-11-16 | 2017-02-22 | 西安电子科技大学 | Mimo雷达收发阵列互耦误差自校正方法 |
CN106788655A (zh) * | 2016-11-30 | 2017-05-31 | 电子科技大学 | 互耦条件下未知互耦信息的干扰相干稳健波束形成方法 |
CN107290730A (zh) * | 2017-05-10 | 2017-10-24 | 长江大学 | 互耦条件下双基地mimo雷达角度估算方法 |
-
2018
- 2018-11-13 CN CN201811348931.5A patent/CN109358312B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101149429A (zh) * | 2006-09-18 | 2008-03-26 | 中国人民解放军空军工程大学导弹学院 | 适用于均匀圆阵的阵列互耦校正与信源测向方法 |
CN104931931A (zh) * | 2015-05-18 | 2015-09-23 | 哈尔滨工程大学 | 互耦条件下基于张量实值子空间的双基地mimo雷达角度估计方法 |
CN105974370A (zh) * | 2016-06-30 | 2016-09-28 | 哈尔滨工业大学 | 一种基于虚拟阵元空间平滑互耦校正的music空间谱估计方法 |
CN106154217A (zh) * | 2016-07-12 | 2016-11-23 | 南京邮电大学 | Ula和uca中互耦未知时基于空间谱伪峰消除的自校准方法 |
CN106324558A (zh) * | 2016-08-30 | 2017-01-11 | 东北大学秦皇岛分校 | 基于互质阵列的宽带信号doa估计方法 |
CN106443610A (zh) * | 2016-11-16 | 2017-02-22 | 西安电子科技大学 | Mimo雷达收发阵列互耦误差自校正方法 |
CN106788655A (zh) * | 2016-11-30 | 2017-05-31 | 电子科技大学 | 互耦条件下未知互耦信息的干扰相干稳健波束形成方法 |
CN107290730A (zh) * | 2017-05-10 | 2017-10-24 | 长江大学 | 互耦条件下双基地mimo雷达角度估算方法 |
Non-Patent Citations (3)
Title |
---|
B. LIAO ETC.,: "DOA Estimation of Coherent Signals for Uniform Linear Arrays with Mutual Coupling", 《IEEE INTERNATIONAL SYMPOSIUM ON CIRCUITS AND SYSTEMS》 * |
司伟建 等: "一种阵列互耦矩阵与波达方向的级联估计方法", 《电子与信息学报》 * |
杨洋: "多径条件下波达方向估计与互耦自校正算法研究", 《中国优秀博硕士学位论文全文数据库(硕士)信息科技辑》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110927659A (zh) * | 2019-11-25 | 2020-03-27 | 长江大学 | 互耦条件下任意阵列流形doa估计与互耦校准方法及*** |
CN110927659B (zh) * | 2019-11-25 | 2022-01-14 | 长江大学 | 互耦条件下任意阵列流形doa估计与互耦校准方法及*** |
Also Published As
Publication number | Publication date |
---|---|
CN109358312B (zh) | 2023-05-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7855681B2 (en) | Systems and methods for determining element phase center locations for an array of antenna elements | |
US7969358B2 (en) | Compensation of beamforming errors in a communications system having widely spaced antenna elements | |
US10571544B2 (en) | Direction finding using signal power | |
US7970365B2 (en) | Systems and methods for compensating for transmission phasing errors in a communications system using a receive signal | |
Dai et al. | Spatial smoothing for direction of arrival estimation of coherent signals in the presence of unknown mutual coupling | |
CN110113085B (zh) | 一种基于协方差矩阵重构的波束形成方法及*** | |
CN106501770B (zh) | 基于幅相误差阵列的远近场宽带混合源中近场源定位方法 | |
CN109507698B (zh) | 卫星导航抗干扰导向矢量自动校正*** | |
CN107015191A (zh) | 一种在多径干扰环境下单偶极子极化敏感阵列降维doa估计方法 | |
CN113835063A (zh) | 一种无人机阵列幅相误差与信号doa联合估计方法 | |
CN112255629A (zh) | 基于联合uca阵列的序贯esprit二维不相干分布源参数估计方法 | |
CN103593509B (zh) | 一种飞行器表面由等离子鞘产生的多径干扰分析方法 | |
CN109358312A (zh) | 确定入射信号来波方向的方法、装置、介质及设备 | |
Al-Sadoon et al. | A more efficient AOA method for 2D and 3D direction estimation with arbitrary antenna array geometry | |
Liao et al. | A subspace-based method for DOA estimation of uniform linear array in the presence of mutual coupling | |
Chen et al. | Successive DSPE-based coherently distributed sources parameters estimation for unmanned aerial vehicle equipped with antennas array | |
Mills et al. | Fast iterative interpolated beamforming for interference DOA estimation in GNSS receivers using fully augmentable arrays | |
Hua et al. | Efficient two dimensional direction finding via auxiliary-variable manifold separation technique for arbitrary array structure | |
Yeom et al. | Performance analysis of beamspace MUSIC with beamforming angle | |
Song et al. | Robust adaptive beamforming under quadratic constraint with recursive method implementation | |
XU et al. | A maximum-likelihood based mutual coupling calibration algorithm in the presence of multipath for GPS antenna array | |
Brossard et al. | Calibration of radio interferometers using a sparse DOA estimation framework | |
Elbir et al. | Direction finding and localization for far-field sources with near-field multipath reflections | |
Al Jabr et al. | Modified UCA-ESPRIT for estimating DOA of coherent signals using one snapshot | |
Letrou | A Gaussian beam shooting scheme for fast multidimensional physical simulation of propagation channels in wireless communication systems |
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 |