CN106772224A - 一种采用时频分析的l型阵列二维波达方向估计算法 - Google Patents
一种采用时频分析的l型阵列二维波达方向估计算法 Download PDFInfo
- Publication number
- CN106772224A CN106772224A CN201611238085.2A CN201611238085A CN106772224A CN 106772224 A CN106772224 A CN 106772224A CN 201611238085 A CN201611238085 A CN 201611238085A CN 106772224 A CN106772224 A CN 106772224A
- Authority
- CN
- China
- Prior art keywords
- time
- frequency
- information source
- matrix
- array
- 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.)
- Withdrawn
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)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种采用时频分析的L型阵列二维波达方向估计算法,主要解决时频谱混叠和空间邻近信源的二维DOA估计问题,并且适用于低信噪比和欠定情况,其实现步骤是:1)利用两个均匀线阵形成xz平面的L型阵列,构造总体接收数据;2)采用STFT将总体接收数据变换到时频域;3)选取各信源在二维方向上的单信源时频点,建立单信源的时频接收数据矩阵;4)构造时频互相关矩阵;5)根据共轭对称性及旋转不变性,计算孔径扩大后的新接收矩阵;6)基于传播算子原则,定义孔径扩大后的传播算子;7)构造角度选择矩阵,依据子阵间的旋转不变性,计算二维DOA。本发明提高了对时频谱混叠和空间邻近信源的估计精度和成功率,对噪声鲁棒性强,并且能够减少所需阵元数。
Description
技术领域
本发明属于阵列信号处理领域,特别涉及一种基于时频分析的二维波达方向估计仿真方法,可解决时频谱混叠和空间邻近信源的方位估计。
背景技术
二维波达方向(DOA)估计在阵列信号处理领域非常重要,它是通过处理空间中按一定方式排列的阵元上的接收数据,估计多个来波信号的方位信息,如方位角和仰角。二维DOA估计技术可应用于无线通信,雷达,声呐定位等应用,受到国内外学者广泛研究。研究过程中,发展出许多阵列模型,如均匀矩形阵列,双平行均匀线阵,均匀圆阵及L型阵,其中L型阵凭借较为简单的结构和参数设置,易于传统算法的移植,更因其能够提供更大的阵列孔径,使得该类型阵列在二维DOA估计效果上更佳,从而受到大量关注。
目前,已经有大量针对二维DOA估计的算法提出。子空间类算法主要将接收数据协方差矩阵分解为信号子空间和噪声子空间,其中多重信号分类MUSIC算法思路简单易于操作,但需要在二维方向上谱峰搜索;旋转不变ESPRIT算法无需谱峰搜索,但要求额外的角度配对过程。虽然子空间类算法具有超分辨率,但当信源角度邻近或信噪比较低时,估计性能会急速下降。另一类基于传播算子PM思想的算法仅需线性操作即可获得信号子空间,运算复杂度降低。经过不断发展,已出现许多具有高分辨率、抗噪性强、无需谱峰搜索和能自动配对的变形算法,然而这些算法要求信源个数需小于阵元个数,不适用于欠定情况,在工程应用中受限。除此之外,当多个信源DOA在二维空间十分邻近时,估计效果衰退。
实际应用中,多个待估计信源的时频谱混叠且空间邻近,即单个时间采样点上可能对应多个空间邻近信源的频率。在一维DOA估计中,通过引入时频分析,滤除时频域中多信号项和噪声项时频点,仅采用单信源时频点参与角度计算,已能达到优良的估计效果。虽然随着国内外学者对二维DOA估计研究的不断深入,各方面性能都得以大幅提升,但目前为止,仍未提出有效针对该类信源的二维DOA估计算法。此外,现有的大多数二维DOA估计算法在低信噪比和欠定情况下效果不佳。
发明内容
本发明主要针对现有技术存在的问题,提出了一种经短时傅里叶变换(Short-time Fourier transform,STFT)将时域信号变换到时频域,采用二维方向上各信源对应的单信源时频点STFT值形成接收数据,再结合PM法则和阵列的旋转不变性进行二维角度自动配对的仿真方法。本发明着重解决时频谱混叠和空间邻近信源,对噪声具有高鲁棒性,且能够在阵元数小于待估计信源数情况下保持良好的估计效果。
本发明提出的仿真方法主要包括以下步骤:
1)两个均匀线阵正交摆放,形成xz平面上的L型阵列,并假设K个时频谱混叠信源以方位角θ和仰角φ入射到此阵列上。其中每个方向上有M个阵元,阵元间距为d(0<d≤λ/2,λ为入射窄带信号波长),θ代表入射信源与x轴正方向的夹角,φ代表入射信源与z轴正方向夹角;
2)对阵列二维方向上的接收数据分别采样,得到输出数据X(t),Z(t),获得总体接收数据
3)采用时频分析工具STFT,将Y(t)变换到时频域Wy(t,f),Wx(t,f),Wz(t,f)分别是总体输出数据、x方向输出数据和z方向输出数据的STFT值;
4)利用合理的阈值和k-means聚类法选取时频域中每个时间采样点上有且只有一个频率的时频点。所形成的单信源时频点集用Φ*表示。已知信源数K的情况下,再次利用k-means聚类法将点集Φ*划分为K类,每类代表每个信源对应的单信源时频点集Φi,i=1,…,K;
5)依据Φi,i=1,…,K,建立单个信源的时频接收数据矩阵Xtf=Ytf(1:M,:),Ztf=Ytf(M+1:2M,:)代表x,z方向上各信源的单信源时频点STFT值。D(i:j,:)代表矩阵D的第i到j行数据组成的子阵;
6)计算时频域互相关矩阵:
其中E{}代表期望运算符,H代表共轭转置运算符;
7)根据线性阵列的共轭对称特性,以及均匀线阵的旋转不变性,计算扩大阵列孔径后的新接收矩阵:
其中Y1=Rxz(:,1:M-1),Y2=Rxz(:,2:M),Y4=JMY1 *JM-1,JM为副对角线全1的M×M维置换矩阵,*代表共轭运算符,T代表转置运算符;
8)基于步骤7中所获Y,通过传播算子法原则,求出传播算子P,并定义一个扩展的时频传播算子:
Pe=[IK,P]H;
其中IK为K×K维单位矩阵;
9)定义两个方位角选择矩阵Gx1,Gx2,两个仰角选择矩阵Gz1,Gz2,结合ESPRIT思想建立旋转不变等式,推导出方位角和仰角计算公式。
若信源数目K=1,顺序执行步骤1)至9);若K>1,执行至第4)步后步骤5)-9)重复K次,共获得K对计算结果。
步骤4)所述的求取各信源在二维方向上的单信源时频点,按如下步骤进行:
步骤4a)、选择高能量信号项时频点:
式中,WY(t,v)表示总体输出数据在时间点t时,能获得的最大STFT值,且该值出现于频率v处。δ0定义为高能量时频点选定阈值,取值范围为0.3~0.5,在此范围内取值大小仅影响选取的时频点个数,对最终估计效果无影响。满足条件的点集为Φ,计算Φ中每个点的空间矢量:
分别代表矩阵Y(t)各行,即各阵元接收数据的STFT值。采用k-means聚类算法将计算出的空间矢量分为K0类,K0>K,得到K0个方向矢量:
其中γk是第k类时频点集合,是集合γk中时频点的个数;方向矢量构成矩阵
步骤4b)、计算Φ中每个时频点上最可能出现的两个信源所对应的导向矢量an1,an2:
其中代表指向的噪声子空间的正交投影矩阵,其元素是矩阵A0中的随机两列,I为单位矩阵,根据导向矢量an1,an2可以获得对应的两个STFT值:
其中A2=[an1,an2],代表Moore-Penrose伪逆运算符;
步骤4c)、选取单信源时频点并分类,步骤如下:
单信源时频点上仅有信号项和噪声项,因此W1,W2的数量级有明显差异:
阈值0反映信号与噪声STFT值相差程度,设为5,τ0反映噪声门限值,其值设定与噪声大小有关;将满足上述条件的时频点集用Φ*表示,再次使用k-means聚类法将Φ*中的点分为K类,每一类用集合Φi,i=1,…,K表示;至此,各个信源的单信源时频点被选取出来。
步骤5)所述的依据Φi,i=1,…,K,建立单个信源的时频接收数据矩阵按如下步骤进行:
步骤5a)、建立一个三维矩阵Dy{m}(t,f),存放2M个阵元的时频分布:
定义索引矢量其各元素代表点集Φi中每个单信源时频点在Φ*上对应的坐标,坐标定义为:例如有一个3×3矩阵,坐标5代表第二列的第二个元素;坐标7代表第三列的第一个元素。Ni为Φi中时频点数,时频接收数据矩阵可用如下计算式获得:
Ytf(m,:)=Dy{m}(qi),m=1,…,2M
步骤5b)、分离两个方向的时频数据:
Xtf=Ytf(1:M,:),Ztf=Ytf(M+1:2M,:);
x和z方向上的接收信号都来自于相同的信号,因此x和z方向上的信号具有相同的时频分布,对总体时频接收数据WY(t,f)进行一次单信源分离,即可获得x,z方向上时频坐标相同的匹配时频点。
与现有算法相比,本发明具有下述优点:
1)能够解决实际应用中时频谱混叠和空间邻近信源的二维DOA估计问题。本发明将时域信源转换到时频域中,放大信源特征,利用合理的阈值设定和k-means聚类法,选取信源在二维方向上具有相同坐标的单信源时频点进行计算,消除了多信号时频点和噪声时频点带来的影响。
2)在低信噪比情况下仍具有良好的估计效果。本发明一方面利用噪声在时频域上均匀分布,而信号在时频处有聚集性的特点集中信号强度,另一方面过滤能量低的噪声时频点,提高对噪声的鲁棒性。
3)能够在欠定情况下保持高估计精度。本发明分别对单个信源进行二维DOA估计,则可以放宽各方向阵元数M必须大于信源数K的要求,仅需少数阵元就可以获得较高精度的估计结果。
附图说明
图1是本发明的实现流程图;
图2是本发明的阵列及参数示意图;
图3是单脉冲(AM),线性调频(LFM),正弦调频(SFM)以及频移键控(FSK)信号在进行单信源时频点选取前后的时频分布图;
图4是本发明与现有两种二维DOA估计算法在不同信噪比条件下均方根误差对比图;
图5是本发明与现有两种二维DOA估计算法在不同仰角角度间隔下估计成功率对比图;
图6是本发明与现有两种二维DOA估计算法在非欠定(a)和欠定(b)情况下估计效果图,其中,图6(a)为M=10,K=4时二维DOA估计图,图6(b)为M=3,K=4时二维DOA估计图。
具体实施方式
以下参照附图,对本发明的技术方案和效果作进一步的详细说明。
参照图1,本发明的实现步骤如下:
步骤1:利用两个均匀线阵形成xz平面的L型阵列。
在x和z方向上,分别以相同间隔d放置阵元,各方向上有M个阵元,并且规定原点阵元处于z方向上。两均匀线阵正交摆放形成xz平面上的L型阵列。假设有K个频率混叠远场信号(N代表快拍采样个数)以方位角θk和仰角φk入射到此阵列上,且信号在传播过程中加入均值为0的复高斯白噪声;
步骤2:对x,z方向的均匀线阵接收数据,分别进行时域采样,获得输出数据X(n),Z(n),n=1,…,N,构造总体接收数据Y(n):
其中Y(n)维度为2M×N,X(n)=Y(1:M,:),Z(n)=Y(M+1:2M,:);
步骤3:采用STFT将时域总体接收数据变换到时频域。
总体接收数据Y(n)的STFT值可以用X(n),Z(n)的STFT值表示:
其中WY(t,f),WX(t,f),WZ(t,f),WS(t,f)分别表示Y(n),X(n),Z(n),S(n)的短时傅里叶值,Ax,Az代表x和z阵上形成接收数据的导向矢量阵。至此,接收信号的信息可以通过时频点的分布来表示;
步骤4:选取各信源在二维方向上的单信源时频点。
4a)选择高能量信号项时频点:
式中WY(t,v)表示总体输出数据在时间点t时,能获得的最大STFT值,且该值出现于频率v处。δ0定义为高能量时频点选定阈值,取值范围为0.3~0.5,在此范围内取值大小仅影响选取的时频点个数,对最终估计效果无影响。满足条件的点集为Φ,计算Φ中每个点的空间矢量:
分别代表矩阵Y(t)各行,即各阵元接收数据的STFT值。采用k-means聚类算法将计算出的空间矢量分为K0类,K0>K,得到K0个方向矢量:
其中γk是第k类时频点集合,Nγk是集合γk中时频点的个数;方向矢量构成矩阵
4b)计算Φ中每个时频点上最可能出现的两个信源所对应的导向矢量an1,an2:
其中代表指向的噪声子空间的正交投影矩阵,其元素是矩阵A0中的随机两列,I为单位矩阵,根据导向矢量an1,an2可以获得对应的两个STFT值:
其中A2=[an1,an2],代表Moore-Penrose伪逆运算符;
4c)选取单信源时频点并分类,步骤如下:
单信源时频点上仅有信号项和噪声项,因此W1,W2的数量级有明显差异:
阈值λ0反映信号与噪声STFT值相差程度,设为5,τ0反映噪声门限值,其值设定与噪声大小有关;将满足上述条件的时频点集用Φ*表示,再次使用k-means聚类法将Φ*中的点分为K类,每一类用集合Φi,i=1,…,K表示;至此,各个信源的单信源时频点被选取出来。
步骤5、根据上述选择结果,建立单信源的时频接收数据矩阵。
5a)建立一个三维矩阵Dy{m}(t,f),存放2M个阵元的时频分布:
定义索引矢量其各元素代表点集Φi中每个单信源时频点在Φ*上对应的坐标,坐标定义为:例如有一个3×3矩阵,坐标5代表第二列的第二个元素;坐标7代表第三列的第一个元素。Ni为Φi中时频点数,时频接收数据矩阵可用如下计算式获得:
Ytf(m,:)=Dy{m}(qi),m=1,…,2M
5b)分离两个方向的时频数据:
Xtf=Ytf(1:M,:),Ztf=Ytf(M+1:2M,:);
x和z方向上的接收信号都来自于相同的信号,因此x和z方向上的信号具有相同的时频分布,对总体时频接收数据WY(t,f)进行一次单信源分离,即可获得x,z方向上时频坐标相同的匹配时频点。
步骤6:求解时频互相关矩阵:
步骤7:根据线性阵列的共轭对称特性,以及均匀线阵的旋转不变性,计算扩大阵列孔径后的新接收矩阵。
7a)划分时频互相关子阵:
Y1=Rxz(:,1:M-1),Y2=Rxz(:,2:M),
由于接收数据的导向矢量阵是范德蒙特矩阵,且线性阵列具有共轭对称的特性,因此建立两个新的互相关子阵:
其中JM为副对角线全1的M×M维置换矩阵。
7b)构造类似孔径扩大后的新接收矩阵:
新接收矩阵维度扩大为4M×(M-1)。
步骤8、根据传播算子原理,构成Y的新导向矢量阵可以被传播算子P分为两部分,传播算子的计算公式如下:
其中Ryy1,Ryy2为自相关矩阵Ryy=E{YYH}前K行和后4M-K行,为伪逆运算符。为了后面利用旋转不变性建立子阵关系,这里定义一个扩大的传播算子:
Pe=[IK,P]H
步骤9、定义两个方位角选择矩阵和两个仰角选择矩阵,用于建立孔径扩大后子阵间旋转不变关系:
其中0(M-1)×1,0M×M代表(M-1)×1维和M×M的全0矩阵。根据旋转不变等式,可以推导出θk,φk的计算式:
其中Φx为对角阵,对角元素由矩阵奇异值分解后的特征值按降序排列构成;Γ由矩阵奇异值分解后的特征矢量阵构成。
若信源数目K=1,顺序执行步骤1)至9);若K>1,执行至第4)步后步骤5)-9)重复K次,共获得K对计算结果。
本发明的效果可通过以下仿真说明:
仿真条件:可能采用的信号为AM,LFM,SFM以及FSK,各信号时频混叠情况如图3上层所示,经过单信源选取后时频分布如图3下层所示。采样点数为512,观测空域角度范围为[-90°,90°]。均方根误差的计算式为:
其中J表示实验次数,φk,θk分别表示二维角度估计值与实际值。估计成功率表示J次独立实验中检测成功的次数,其中一次实验检测成功的定义为该次所有的DOA估计值均在真实值±1°范围内。
仿真1:信噪比对算法的影响
假设二维各方向M=6,K=3,信源为LFM,SFM和FSK信号,入射角度θ=[98.49°,101.42°,103.08°],φ=[91.72°,65.84°,76.16°]。信噪比由-5-15dB变化,本发明与JSVD,Dong三种算法在实验次数J为1000时,均方根误差结果如图4。
图4中横坐标表示信噪比值,纵坐标表示均方根误差。
从图4中可以看出本发明在各信噪比下,甚至是低信噪比时,均方根误差都明显小于其他两种算法。
仿真2:不同仰角角度间隔下的估计成功率
假设二维各方向M=5,信噪比为10dB。K=2,分别是SFM,FSK,入射角度为θ=[55.49°,65.42°],φ=[140.72°,140.72°+Δφ],Δφ由1°-20°变化。本发明与JSVD,Dong三种算法在实验次数J为100时,估计成功率结果如图5。
图5中横坐标表示仰角角度间隔,纵坐标表示估计成功率。
从图5中可以看出,本发明与其他两种算法相比,大幅提高了估计成功率,即便两个信源邻近时,依然能够分辨二维DOA。
仿真3:阵元数对算法的影响
假设信噪比为10dB,K=4,分别为AM,LFM,SFM和FSK,入射角度为θ=[30°,35°,40°,41°],φ=[104°,70°,75°,80°]。本发明与其他两种算法在M=10,M=3(欠定情况,即阵元数小于信源数情况)时,仿真结果如图6(a),6(b)所示。
图6(a)(b)中横坐标表示方位角,纵坐标表示仰角,圆圈代表信源的真实DOA。
从6(a)中可以看出,阵元数大于信源数时,本发明比其他两种算法估计精度更高;
从6(b)中可以看出,本发明同样适用于欠定情况,且保持较高的估计精度,其他两种算法已失去估计功能。
Claims (4)
1.一种采用时频分析的L型阵列二维波达方向估计算法,其特征在于,包括以下步骤:
步骤1)、两个均匀线阵正交摆放,形成xz平面上的L型阵列,并假设K个时频谱混叠信源以方位角θ和仰角φ入射到此阵列上;其中每个方向上有M个阵元,阵元间距为d,0<d≤λ/2,λ为入射窄带信号波长,θ代表入射信源与x轴正方向的夹角,φ代表入射信源与z轴正方向夹角;
步骤2)、对阵列二维方向上的接收数据分别采样,得到输出数据X(t),Z(t),获得总体接收数据
步骤3)、采用时频分析工具STFT,将Y(t)变换到时频域Wy(t,f),Wx(t,f),Wz(t,f)分别表示总体输出数据、x方向输出数据和z方向输出数据的STFT值;
步骤4)、根据信号和噪声的差异程度及噪声功率,设定阈值,并结合和k-means聚类选取时频域中每个时间采样点上有且只有一个频率的时频点;所形成的单信源时频点集用Φ*表示;已知信源数K的情况下,再次利用k-means聚类法将点集Φ*划分为K类,每类代表每个信源对应的单信源时频点集Φi,i=1,…,K;
步骤5)、依据Φi,i=1,…,K,建立单个信源的时频接收数据矩阵Xtf=Ytf(1:M,:),Ztf=Ytf(M+1:2M,:)代表x,z方向上各信源的单信源时频点STFT值;D(i:j,:)代表矩阵D的第i到j行数据组成的子阵;
步骤6)、计算时频域互相关矩阵:
其中,E{}代表期望运算符,H代表共轭转置运算符;
步骤7)、根据线性阵列的共轭对称特性,以及均匀线阵的旋转不变性,计算扩大阵列孔径后的新接收矩阵:
其中JM为副对角线全1的M×M维置换矩阵,*代表共轭运算符,T代表转置运算符;
步骤8)、基于步骤7中所获Y,通过传播算子法原则,求出传播算子P,并定义一个扩展的时频传播算子:
Pe=[IK,P]H,
其中,IK为K×K维单位矩阵;
步骤9)、定义两个方位角选择矩阵Gx1,Gx2,两个仰角选择矩阵Gz1,Gz2,结合ESPRIT思想建立旋转不变等式,推导出方位角和仰角计算公式。
2.根据权利要求1所述的一种采用时频分析的L型阵列二维波达方向估计算法,其特征在于:
若信源数目K=1,顺序执行步骤1)至9);若K>1,执行至第4)步后步骤5)-9)重复K次,共获得K对计算结果。
3.根据权利要求1所述的一种采用时频分析的L型阵列二维波达方向估计算法,其特征在于:步骤4)所述的求取各信源在二维方向上的单信源时频点,按如下步骤进行:
步骤4a)、选择高能量信号项时频点:
式中,WY(t,v)表示总体输出数据在时间点t时,能获得的最大STFT值,且该值出现于频率v处;δ0定义为高能量时频点选定阈值,取值范围为0.3~0.5;满足条件的点集为Φ,计算Φ中每个点的空间矢量:
分别代表矩阵Y(t)各行,即各阵元接收数据的STFT值;采用k-means聚类算法将计算出的空间矢量分为K0类,K0>K,得到K0个方向矢量:
其中γk是第k类时频点集合,是集合γk中时频点的个数;方向矢量构成矩阵
步骤4b)、计算Φ中每个时频点上最可能出现的两个信源所对应的导向矢量an1,an2:
其中代表指向的噪声子空间的正交投影矩阵,其元素是矩阵A0中的随机两列,I为单位矩阵,根据导向矢量an1,an2可以获得对应的两个STFT值:
其中A2=[an1,an2],◇代表Moore-Penrose伪逆运算符;
步骤4c)、选取单信源时频点并分类,步骤如下:
单信源时频点上仅有信号项和噪声项,因此W1,W2的数量级有明显差异:
阈值λ0反映信号与噪声STFT值相差程度,设为5,τ0反映噪声门限值,其值设定与噪声大小有关;将满足上述条件的时频点集用Φ*表示,再次使用k-means聚类法将Φ*中的点分为K类,每一类用集合Φi,i=1,…,K表示;至此,各个信源的单信源时频点被选取出来。
4.根据权利要求1所述的一种采用时频分析的L型阵列二维波达方向估计算法,其特征在于:步骤5)所述的依据Φi,i=1,…,K,建立单个信源的时频接收数据矩阵按如下步骤进行:
步骤5a)、建立一个三维矩阵Dy{m}(t,f),存放2M个阵元的时频分布:
定义索引矢量其各元素代表点集Φi中每个单信源时频点在Φ*上对应的坐标,Ni为Φi中时频点数,时频接收数据矩阵可用如下计算式获得:
Ytf(m,:)=Dy{m}(qi),m=1,…,2M
步骤5b)、分离两个方向的时频数据:
Xtf=Ytf(1:M,:),Ztf=Ytf(M+1:2M,:);
x和z方向上的接收信号都来自于相同的信号,因此x和z方向上的信号具有相同的时频分布,对总体时频接收数据WY(t,f)进行一次单信源分离,即可获得x,z方向上时频坐标相同的匹配时频点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611238085.2A CN106772224A (zh) | 2016-12-28 | 2016-12-28 | 一种采用时频分析的l型阵列二维波达方向估计算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611238085.2A CN106772224A (zh) | 2016-12-28 | 2016-12-28 | 一种采用时频分析的l型阵列二维波达方向估计算法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106772224A true CN106772224A (zh) | 2017-05-31 |
Family
ID=58925037
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611238085.2A Withdrawn CN106772224A (zh) | 2016-12-28 | 2016-12-28 | 一种采用时频分析的l型阵列二维波达方向估计算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106772224A (zh) |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109507634A (zh) * | 2018-11-08 | 2019-03-22 | 中国电子科技集团公司第二十八研究所 | 一种任意传感器阵列下的基于传感算子的盲远场信号波达方向估计方法 |
CN109507675A (zh) * | 2019-01-07 | 2019-03-22 | 中国科学院声学研究所东海研究站 | 基于频分系统实现水下多目标时延估计处理的方法 |
CN109521392A (zh) * | 2018-10-24 | 2019-03-26 | 华南理工大学 | 基于非圆信号和l型线阵的水下一维doa估计方法和装置 |
CN109541524A (zh) * | 2018-12-14 | 2019-03-29 | 哈尔滨工业大学 | 一种基于欠采样的信号载频与二维doa联合估计方法 |
CN109814064A (zh) * | 2019-02-28 | 2019-05-28 | 中国电子科技集团公司第三十六研究所 | 一种基于三阵元l型直角阵干涉仪测向方法和装置 |
CN110045322A (zh) * | 2019-03-21 | 2019-07-23 | 中国人民解放军战略支援部队信息工程大学 | 一种基于高分辨测向语图智能识别的短波自动测向方法 |
CN110133577A (zh) * | 2019-05-23 | 2019-08-16 | 成都理工大学 | 基于时频残差网络的单基互质mimo阵列doa估计算法 |
CN110208738A (zh) * | 2019-06-13 | 2019-09-06 | 哈尔滨工业大学 | 基于阵列调制宽带转换器的信号频率与二维doa联合估计方法 |
CN110749857A (zh) * | 2019-09-12 | 2020-02-04 | 宁波大学 | 基于秩损法的均匀矩形阵列二维非圆信号doa估计方法 |
CN111521968A (zh) * | 2020-05-22 | 2020-08-11 | 南京理工大学 | 基于目标空间分集的欠定doa估计方法 |
CN114036458A (zh) * | 2021-10-22 | 2022-02-11 | 哈尔滨工程大学 | 一种非合作水声信号时频信息获取方法 |
US20220120847A1 (en) * | 2020-10-20 | 2022-04-21 | Aptiv Technologies Limited | Radar System with Modified Orthogonal Linear Antenna Subarrays |
US11644565B2 (en) | 2020-10-13 | 2023-05-09 | Aptiv Technologies Limited | Radar system with sparse primary array and dense auxiliary array |
US11714180B2 (en) | 2021-01-29 | 2023-08-01 | Aptiv Technologies Limited | Radar system to detect angles in bistatic and monostatic scenarios |
US11774570B2 (en) | 2020-02-04 | 2023-10-03 | Aptiv Technologies Limited | Radar device |
US11808846B2 (en) | 2021-02-12 | 2023-11-07 | Aptiv Technologies Limited | Angle-finding process for sparse uniform arrays |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102253363A (zh) * | 2011-03-29 | 2011-11-23 | 西安交通大学 | 基于l型阵列的相干信号二维波达方向估计装置及其方法 |
CN104678368A (zh) * | 2015-02-16 | 2015-06-03 | 零八一电子集团有限公司 | 一维相扫三坐标雷达空域自适应干扰抑制方法 |
CN104730491A (zh) * | 2015-03-06 | 2015-06-24 | 中国计量学院 | 一种基于l型阵的虚拟阵列doa估计方法 |
CN105158735A (zh) * | 2015-06-03 | 2015-12-16 | 西安电子科技大学 | 基于压缩采样阵列的空频二维谱估计方法 |
CN105203989A (zh) * | 2015-09-28 | 2015-12-30 | 天津伽利联科技有限公司 | 一种l型阵列相干信号波达方向的估计方法 |
CN105577587A (zh) * | 2016-02-03 | 2016-05-11 | 清华大学 | 一种60GHz室内场景下MIMO系统信道估计方法及装置 |
-
2016
- 2016-12-28 CN CN201611238085.2A patent/CN106772224A/zh not_active Withdrawn
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102253363A (zh) * | 2011-03-29 | 2011-11-23 | 西安交通大学 | 基于l型阵列的相干信号二维波达方向估计装置及其方法 |
CN104678368A (zh) * | 2015-02-16 | 2015-06-03 | 零八一电子集团有限公司 | 一维相扫三坐标雷达空域自适应干扰抑制方法 |
CN104730491A (zh) * | 2015-03-06 | 2015-06-24 | 中国计量学院 | 一种基于l型阵的虚拟阵列doa估计方法 |
CN105158735A (zh) * | 2015-06-03 | 2015-12-16 | 西安电子科技大学 | 基于压缩采样阵列的空频二维谱估计方法 |
CN105203989A (zh) * | 2015-09-28 | 2015-12-30 | 天津伽利联科技有限公司 | 一种l型阵列相干信号波达方向的估计方法 |
CN105577587A (zh) * | 2016-02-03 | 2016-05-11 | 清华大学 | 一种60GHz室内场景下MIMO系统信道估计方法及装置 |
Cited By (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109521392A (zh) * | 2018-10-24 | 2019-03-26 | 华南理工大学 | 基于非圆信号和l型线阵的水下一维doa估计方法和装置 |
CN109521392B (zh) * | 2018-10-24 | 2022-12-16 | 华南理工大学 | 基于非圆信号和l型线阵的水下一维doa估计方法和装置 |
CN109507634B (zh) * | 2018-11-08 | 2020-08-11 | 中国电子科技集团公司第二十八研究所 | 一种任意传感器阵列下的基于传播算子的盲远场信号波达方向估计方法 |
CN109507634A (zh) * | 2018-11-08 | 2019-03-22 | 中国电子科技集团公司第二十八研究所 | 一种任意传感器阵列下的基于传感算子的盲远场信号波达方向估计方法 |
CN109541524A (zh) * | 2018-12-14 | 2019-03-29 | 哈尔滨工业大学 | 一种基于欠采样的信号载频与二维doa联合估计方法 |
CN109541524B (zh) * | 2018-12-14 | 2020-07-24 | 哈尔滨工业大学 | 一种基于欠采样的信号载频与二维doa联合估计方法 |
CN109507675A (zh) * | 2019-01-07 | 2019-03-22 | 中国科学院声学研究所东海研究站 | 基于频分系统实现水下多目标时延估计处理的方法 |
CN109507675B (zh) * | 2019-01-07 | 2020-10-16 | 中国科学院声学研究所东海研究站 | 基于频分系统实现水下多目标时延估计处理的方法 |
CN109814064A (zh) * | 2019-02-28 | 2019-05-28 | 中国电子科技集团公司第三十六研究所 | 一种基于三阵元l型直角阵干涉仪测向方法和装置 |
CN110045322A (zh) * | 2019-03-21 | 2019-07-23 | 中国人民解放军战略支援部队信息工程大学 | 一种基于高分辨测向语图智能识别的短波自动测向方法 |
CN110045322B (zh) * | 2019-03-21 | 2021-01-22 | 中国人民解放军战略支援部队信息工程大学 | 一种基于高分辨测向语图智能识别的短波自动测向方法 |
CN110133577A (zh) * | 2019-05-23 | 2019-08-16 | 成都理工大学 | 基于时频残差网络的单基互质mimo阵列doa估计算法 |
CN110208738B (zh) * | 2019-06-13 | 2022-12-06 | 哈尔滨工业大学 | 基于阵列调制宽带转换器的信号频率与二维doa联合估计方法 |
CN110208738A (zh) * | 2019-06-13 | 2019-09-06 | 哈尔滨工业大学 | 基于阵列调制宽带转换器的信号频率与二维doa联合估计方法 |
CN110749857A (zh) * | 2019-09-12 | 2020-02-04 | 宁波大学 | 基于秩损法的均匀矩形阵列二维非圆信号doa估计方法 |
CN110749857B (zh) * | 2019-09-12 | 2024-03-08 | 深圳万知达科技有限公司 | 基于秩损法的均匀矩形阵列二维非圆信号doa估计方法 |
US11774570B2 (en) | 2020-02-04 | 2023-10-03 | Aptiv Technologies Limited | Radar device |
CN111521968A (zh) * | 2020-05-22 | 2020-08-11 | 南京理工大学 | 基于目标空间分集的欠定doa估计方法 |
US11644565B2 (en) | 2020-10-13 | 2023-05-09 | Aptiv Technologies Limited | Radar system with sparse primary array and dense auxiliary array |
US20220120847A1 (en) * | 2020-10-20 | 2022-04-21 | Aptiv Technologies Limited | Radar System with Modified Orthogonal Linear Antenna Subarrays |
US11619705B2 (en) * | 2020-10-20 | 2023-04-04 | Aptiv Technologies Limited | Radar system with modified orthogonal linear antenna subarrays |
US11921228B2 (en) | 2020-10-20 | 2024-03-05 | Aptiv Technologies Limited | Radar system with modified orthogonal linear antenna subarrays |
US11714180B2 (en) | 2021-01-29 | 2023-08-01 | Aptiv Technologies Limited | Radar system to detect angles in bistatic and monostatic scenarios |
US11808846B2 (en) | 2021-02-12 | 2023-11-07 | Aptiv Technologies Limited | Angle-finding process for sparse uniform arrays |
CN114036458A (zh) * | 2021-10-22 | 2022-02-11 | 哈尔滨工程大学 | 一种非合作水声信号时频信息获取方法 |
CN114036458B (zh) * | 2021-10-22 | 2024-04-30 | 哈尔滨工程大学 | 一种非合作水声信号时频信息获取方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106772224A (zh) | 一种采用时频分析的l型阵列二维波达方向估计算法 | |
CN103018730B (zh) | 分布式子阵波达方向估计方法 | |
CN106526530B (zh) | 基于传播算子的2-l型阵列二维doa估计算法 | |
CN101592721B (zh) | 基于特征值重构的相干信号到达角估计方法 | |
CN106019234B (zh) | L型天线阵的低计算复杂度二维波达方向估计方法 | |
CN106646376A (zh) | 基于加权修正参数的p范数噪声源定位识别方法 | |
CN104678372B (zh) | 正交频分复用雷达超分辨距离与角度值联合估计方法 | |
CN103344940B (zh) | 低复杂度的doa估计方法及系统 | |
CN104515969B (zh) | 一种基于六角形阵列的相干信号二维doa估计方法 | |
CN101351058B (zh) | 一种麦克风阵列及实现语音信号增强的方法 | |
CN103323827B (zh) | 基于快速傅里叶变换的mimo雷达系统角度估计方法 | |
CN109655799A (zh) | 基于iaa的协方差矩阵向量化的非均匀稀疏阵列测向方法 | |
CN103353595A (zh) | 基于阵列内插压缩感知的米波雷达测高方法 | |
CN103744061A (zh) | 基于迭代最小二乘方法的mimo雷达doa估计方法 | |
CN110515038B (zh) | 一种基于无人机-阵列的自适应无源定位装置及实现方法 | |
CN102694588B (zh) | 一种基于共轭扩展的波达方向估计方法 | |
CN106019214A (zh) | 宽带相干信号源doa估计方法 | |
CN106019213A (zh) | 一种部分稀疏l阵及其二维doa估计方法 | |
CN106021637A (zh) | 互质阵列中基于迭代稀疏重构的doa估计方法 | |
CN103116162B (zh) | 基于目标空间稀疏性的高分辨声呐定位方法 | |
CN101644760B (zh) | 一种适用于高分辨阵列的快速鲁棒的信源个数检测方法 | |
CN102662158B (zh) | 一种对传感器天线阵列接收信号的快速处理方法 | |
CN112130111A (zh) | 一种大规模均匀十字阵列中单快拍二维doa估计方法 | |
CN103760527A (zh) | 单基地mimo雷达相干源波达方向估计方法 | |
CN104933290A (zh) | 双l型拉伸正交电偶对阵列的多参数联合估计四元数方法 |
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 | ||
WW01 | Invention patent application withdrawn after publication | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20170531 |