CN109212526B - 用于高频地波雷达的分布式阵列目标角度测量方法 - Google Patents
用于高频地波雷达的分布式阵列目标角度测量方法 Download PDFInfo
- Publication number
- CN109212526B CN109212526B CN201811207040.8A CN201811207040A CN109212526B CN 109212526 B CN109212526 B CN 109212526B CN 201811207040 A CN201811207040 A CN 201811207040A CN 109212526 B CN109212526 B CN 109212526B
- Authority
- CN
- China
- Prior art keywords
- array
- target
- signal
- range
- doppler
- 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
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/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/9004—SAR image acquisition techniques
- G01S13/9011—SAR image acquisition techniques with frequency domain processing of the SAR signals in azimuth
-
- 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/9094—Theoretical aspects
-
- 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
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/418—Theoretical aspects
Abstract
本发明提出了一种用于高频地波雷达的分布式阵列目标角度测量方法,其特征在于包括:步骤一、雷达各子阵阵元接收回波距离多普勒处理;步骤二、数据协方差矩阵稀疏表示建模;步骤三、待求解优化问题转化;步骤四、内循环更新信号功率估计值;步骤五、外循环更新代价函数并最终获得目标角度估计;相比传统雷达阵列,本发明采用的分布式阵列布阵灵活,选址方便,且仅需少量阵元即可达到较大孔径传统阵列的测向性能;而本发明提出的角度测量方法最少仅需距离多普勒域的单次快拍,无需已知目标个数和任何人工参数调节,测向性能也明显优于传统的雷达测角方法。
Description
技术领域
本发明涉及雷达目标测角技术领域,具体的说是一种用于高频地波雷达的分布式阵列目标角度测量方法。
背景技术
雷达目标测角技术是雷达信号处理领域的一个重要研究课题。传统的高频雷达测角一般仍采用波束形成技术,但受限于阵列孔径,其测角精度和分辨力较低。为提高测向性能,基于现代超分辨谱估计技术的测角方法被提出,然而,由于高频雷达系统存在目标个数众多且数量未知、回波信号相干性强等特殊性,此类方法的应用仍受到诸多限制,例如需要已知目标个数、处理相干信号时的解相干预处理步骤要求特定的阵列形式等。尽管近年来涌现的一部分基于稀疏表示的测向方法能够解决以上问题,但其局限性在于,其性能在很大程度上依赖于一个或多个可调的超参数,在未知噪声方差的实际情况下很难将该参数调整到最优以获得满意的性能。另外,为进一步提高测角性能,传统雷达阵列需要较多的接收阵元以加大阵列孔径,但这将导致苛刻的场地需求并显著提高硬件成本。
发明内容
本发明的目的在于改善现有传统雷达阵列布置方案以及测向算法的局限性,提供了一种用于高频地波雷达的分布式阵列目标角度测量方法,该方法最少仅需距离多普勒域的单次快拍,无需已知目标个数和任何人工的参数调节,能够估计同一距离多普勒单元上的多个相干目标角度,能够利用更少的接收阵元获得接近较大孔径传统雷达阵列的测向性能。该方法适用于目标对于整个分布式阵列满足远场假设的条件下,在处理中仅需传输子站中距离多普勒处理后的少量数据。
一种用于高频地波雷达的分布式阵列目标角度测量方法,所述目标角度测量方法包括如下步骤:
步骤一:针对平面内由K(K>1)个子阵组成的分布式雷达阵列(子站内阵元间距不大于工作波长的一半),对各接收通道的雷达回波进行距离多普勒处理,针对感兴趣的目标提取相应距离多普勒单元上的L个数据快拍(L≥1),收集各子站信息用于后续的目标角度测量;
步骤二:利用所有接收通道获得的目标距离多普勒域数据,建立协方差矩阵稀疏表示模型,将目标角度估计问题建模为稀疏功率矢量重构问题,以负对数随机最大似然函数为初始代价函数,优化求解各角度网格上的信号功率;
步骤三:利用Majorization-Minimization方法,通过迭代更新并最小化一个替代的凸代价函数达到最小化原始非凸代价函数的目的,为便于求解,通过引入中间变量进一步将该问题转化为一个交替迭代优化问题;
步骤四:计算中间变量的最优解并将其代入代价函数,进一步对其求导并令导数为零,得到待求解的功率矢量的更新表达式,迭代更新中间变量和功率矢量,直至得到当前代价函数下的功率矢量最优解,再利用获得的功率矢量估计值更新代价函数;
步骤五:重复步骤四所述过程,直至算法收敛,获得功率矢量的最终估计值,根据该矢量中的峰值位置确定感兴趣的目标入射角度。
进一步地,步骤一所述阵列接收数据的距离多普勒处理包括:
第一步:由K(K>1)个子阵构成分布式雷达接收阵列,其中第k(k=1,…,K)个子站中包含Mk个阵元,整个接收阵列包含的总阵元数为由于高频地波雷达针对的是海面目标,无需估计其俯仰角,因此阵列排布采用简单的线阵,要求子阵内部阵元间距不超过半波长,且目标对于整个接收阵列满足远场入射条件,即可认为目标到达各子阵的入射角度相同。雷达发射信号表示为s(t),则第m个阵元接收到的第n个目标回波可表示为:
smn(t)=μns(t-Δtn-τmn)
其中,m=1,2,…,M为阵元编号索引,n=1,2,…,N0为目标编号索引,N0为总目标数;μn为第n个目标信号的传播衰减因子,Δtn为该目标回波到达参考阵元时相对于发射信号的延迟,其大小与目标距离Rn和径向速度vn有关:Δtn=2(Rn-vnt)/c,c为光速;为回波到达第m个阵元时相对于参考阵元的延迟,xm为阵元位置,为带估计的目标方位角。
对于高频雷达而言,其发射信号载频f0一般远大于带宽,可视为窄带信号,因此第m个阵元的接收回波可进一步写为:
一般而言,对于高频雷达系统常用的脉冲截断线性调频信号和相位编码信号而言,其距离变换可分别通过快速傅里叶变换和相关处理实现,而多普勒变换往往通过快速傅里叶变换实现。鉴于以上距离多普勒处理过程是一个线性操作,对sm(t)进行距离多普勒变换后可得:
对L(L≥1)个批次的回波数据分别进行上述变换,可得(r,d)单元上L次快拍的距离多普勒域数据:
X=AS+N
其中,X=[x(1) … x(L)],x(l)=[x1(l)|r,d … xM(l)|r,d]T,S=[s(1) … s(L)],s(l)=[s1′(l)|r,d … s′N(l)|r,d]T,l为数据批次索引;A=[a1 … aN]为M×N维阵列流型矩阵,N为M×L维观测噪声。
进一步地,步骤二所述协方差矩阵稀疏表示建模包括:
首先,要想有效地利用完整的阵列孔径,需对所有子阵数据进行联合处理,而接收阵列为由多个子阵构成的分布式阵列,即稀疏非均匀阵列,为解决该阵列形式下的强相关信号角度估计问题,采用基于稀疏信号表示的角度估计方法:将感兴趣的空间范围按方位角划分为Z(Z>M且Z>>N)个网格,使目标信号的真实方位角包含在网格中,利用相应网格点对应的信号导向矢量构成字典A′=[a1 … aZ]。
其次,假设信号S的协方差矩阵为对角阵,且由于基于傅里叶变换的距离多普勒处理过程存在对偶性,可认为观测噪声N也服从高斯分布,则步骤一所得距离多普勒域数据的协方差矩阵可表示为:
其中,IM为M阶单位矩阵,P=diag(p),p=[p1 … pZ+M],pz(z=1,2,…,Z)为各角度网格上的信号功率,pZ+m=σm(m=1,2,…,M)为各接收通道噪声方差。显然,仅在信号真实方位角对应的网格上pz非零,否则pz为零,因此功率矢量p是一个稀疏向量,通过重构p即可确定真实信号的方位角。
最后,考虑最小化如下负对数随机最大似然函数以求解功率矢量p:
进一步地,步骤三所述的优化问题转化过程包括:
因此,令
进一步地,步骤四所述的功率矢量求解过程包括:
其中,aZ+m=em为M阶单位矩阵IM的第m列,cq为Copt的第q(q=1,2,…,Z+M)行。
第三,由于pq的独立性,可将上述优化问题分解为Z+M个独立问题分别求解,分别对pq求导并令导数为零,可得第i次外循环的第j+1次内循环迭代中,第z(z=1,2,…,Z)个角度网格上的信号功率更新表达式:
以及第m(m=1,2,…,M)个接收通道的噪声方差更新表达式:
而当各通道噪声方差相等时,有:
进一步地,步骤五所述的最终功率矢量估计和目标角度测量过程包括:
首先,以上述步骤四所述的循环迭代作为所提方法的内循环求解过程,令i→i+1,重复步骤四作为所提方法的外循环求解过程,通过嵌套循环直至算法收敛,得到最终的功率矢量估计值。
最后,以信号功率矢量估计值p=[p1 … pZ]作为空间谱,利用其中的峰值位置确定感兴趣的(r,d)距离多普勒单元内的目标信号角度。
本发明有益效果:
本发明提出的分布式阵列测角方法基于距离多普勒域数据的协方差矩阵稀疏表示模型,算法实施无需已知目标个数、无需人工超参数调节,既充分利用了高频地波雷达距离多普勒处理后信噪比提高、目标在距离多普勒域分离后单个单元内目标个数少、子阵之间数据传输量小等优点,算法本身又能够较好地适用于距离多普勒处理后可利用快拍数少甚至为1、目标信号相关性强等条件,相比于传统的雷达测角方法,其性能有明显的提升。另外,与相同孔径的传统雷达阵列相比,本发明中采用的分布式接收阵列需要的阵元数更少,在节约系统硬件成本的同时又能达到接近的角度估计性能,且布阵灵活,为雷达阵地选址提供了很大的便利。
附图说明
图1为本发明所提方法流程图。
图2为分布式雷达阵列角度测量系统示意图。
图3为角度测量性能仿真结果。
具体实施方式
下面结合具体实施例对本发明做进一步说明,但本发明不受实施例的限制。
一种用于高频地波雷达的分布式阵列目标角度测量方法,如图1所示,所述目标角度测量方法包括如下步骤:
步骤一:针对平面内由K(K>1)个子阵组成的分布式雷达阵列(子站内阵元间距不大于工作波长的一半),对各接收通道的雷达回波进行距离多普勒处理,针对感兴趣的目标提取相应距离多普勒单元上的L个数据快拍(L≥1),收集各子站信息用于后续的目标角度测量;
步骤二:利用所有接收通道获得的目标距离多普勒域数据,建立协方差矩阵稀疏表示模型,将目标角度估计问题建模为稀疏功率矢量重构问题,以负对数随机最大似然函数为初始代价函数,优化求解各角度网格上的信号功率;
步骤三:利用Majorization-Minimization方法,通过迭代更新并最小化一个替代的凸代价函数达到最小化原始非凸代价函数的目的,为便于求解,通过引入中间变量进一步将该问题转化为一个交替迭代优化问题;
步骤四:计算中间变量的最优解并将其代入代价函数,进一步对其求导并令导数为零,得到待求解的功率矢量的更新表达式,迭代更新中间变量和功率矢量,直至得到当前代价函数下的功率矢量最优解,再利用获得的功率矢量估计值更新代价函数;
步骤五:重复步骤四所述过程,直至算法收敛,获得功率矢量的最终估计值,根据该矢量中的峰值位置确定感兴趣的目标入射角度。
为方便表述与理解,此处对下文所用符号进行统一说明:矢量和矩阵以粗斜体表示;操作符(·)T、(·)H和(·)-1分别表示矩阵的转置、共轭转置和取逆;符号|·|、||·||2、E{·}、tr(·)和diag(·)分别表示取行列式、取2范数、取期望、取迹和对角化操作。
具体而言,步骤一所述阵列接收数据的距离多普勒处理包括:
第一步:如图2所示,由K(K>1)个子阵构成分布式雷达接收阵列,其中第k(k=1,…,K)个子站中包含Mk个阵元,整个接收阵列包含的总阵元数为由于高频地波雷达针对的是海面目标,无需估计其俯仰角,因此阵列排布采用线阵,要求子阵内部阵元间距不超过半波长,且目标对于整个接收阵列满足远场入射条件,即可认为目标达到各子阵的入射角度相同。雷达发射信号表示为s(t),则第m个阵元接收到的第n个目标回波可表示为:
smn(t)=μns(t-Δtn-τmn)
其中,m=1,2,…,M为阵元编号索引,n=1,2,…,N0为目标编号索引,N0为总目标数;μn为第n个目标信号的传播衰减因子,Δtn为该目标回波到达参考阵元时相对于发射信号的延迟,其大小与目标距离Rn和径向速度vn有关:Δtn=2(Rn-vnt)/c,c为光速;为回波到达第m个阵元时相对于参考阵元的延迟,xm为阵元位置,为带估计的目标方位角。
对于高频雷达而言,其发射信号载频f0一般远大于带宽,可视为窄带信号,因此第m个阵元的接收回波可进一步写为:
一般而言,对于高频雷达系统常用的脉冲截断线性调频信号和相位编码信号而言,其距离变换可分别通过快速傅里叶变换和相关处理实现,而多普勒变换往往通过快速傅里叶变换实现。鉴于以上距离多普勒处理过程是一个线性操作,对sm(t)进行距离多普勒变换后可得:
对L(L≥1)批次的回波数据分别进行上述变换,可得(r,d)单元上L次快拍的距离多普勒域数据:
X=AS+N
其中,X=[x(1) … x(L)],x(l)=[x1(l)|r,d … xM(l)|r,d]T,S=[s(1) … s(L)],s(l)=[s1′(l)|r,d … s′N(l)|r,d]T,l为数据批次索引;A=[a1 … aN]为M×N维阵列流型矩阵,N为M×L维观测噪声。
进一步地,步骤二所述协方差矩阵稀疏表示建模包括:
首先,要想有效地利用完整的阵列孔径,需对所有子阵数据进行联合处理,而接收阵列为由多个子阵构成的分布式阵列,即稀疏非均匀阵列,为解决该阵列形式下的强相关信号角度估计问题,采用基于稀疏信号表示的角度估计方法:将感兴趣的空间范围按方位角划分为Z(Z>M且Z>>N)个网格,使目标信号的真实方位角包含在网格中,利用相应网格点对应的信号导向矢量构成字典A′=[a1…aZ]。
其次,假设信号S的协方差矩阵为对角阵,且由于基于傅里叶变换的距离多普勒处理过程存在对偶性,可认为观测噪声N也服从高斯分布,则步骤一所得距离多普勒域数据的协方差矩阵可表示为:
其中,IM为M阶单位矩阵,P=diag(p),p=[p1 … pZ+M],pz(z=1,2,…,Z)为各角度网格上的信号功率,pZ+m=σm(m=1,2,…,M)为各接收通道噪声方差。显然,仅在信号真实方位角对应的网格上pz非零,否则pz为零,因此功率矢量p是一个稀疏向量,通过重构p即可确定真实信号的方位角。
最后,考虑最小化如下负对数随机最大似然函数以求解功率矢量p:
进一步地,步骤三所述的优化问题转化过程包括:
因此,令
进一步地,步骤四所述的功率矢量求解过程包括:
其中,aZ+m=em为M阶单位矩阵IM的第m列,cq为Copt的第q(q=1,2,…,Z+M)行。
第三,由于pq的独立性,可将上述优化问题分解为Z+M个独立问题分别求解,分别对pq求导并令导数为零,可得第i次外循环的第j+1次内循环迭代中,第z(z=1,2,…,Z)个角度网格上的信号功率更新表达式:
以及第m(m=1,2,…,M)个接收通道的噪声方差更新表达式:
而当各通道噪声方差相等时,有:
进一步地,步骤五所述的最终功率矢量估计和目标角度测量过程包括:
首先,以上述步骤四所述的循环迭代作为所提方法的内循环求解过程,令i→i+1,重复步骤四作为所提方法的外循环求解过程,通过嵌套循环直至算法收敛,得到最终的功率矢量估计值。
最后,以信号功率矢量估计值p=[p1 … pZ]作为空间谱,利用其中的峰值位置确定感兴趣的(r,d)距离多普勒单元内的目标信号角度。
本发明所提的一种用于高频地波雷达的分布式阵列目标角度测量方法,无需已知目标个数、无需人工参数调节,能够估计同一距离多普勒单元上的多个相干目标角度,采用的分布式阵列结构能够利用更少的接收阵元获得接近较大孔径传统雷达阵列的测向性能。该方法适用于目标对于整个分布式阵列满足远场假设的条件下,最少仅需距离多普勒域的单次快拍,在处理中仅需传输子站中距离多普勒处理后的少量数据。
实施例:
利用所提的用于高频地波雷达的分布式阵列目标角度测量方法得到的测角均方误差曲线如图3所示,其仿真条件如下:雷达发射信号为线性调频中断连续波,载波频率为3MHz,带宽为75KHz,调频周期为200ms,脉冲重复周期为4ms,脉冲宽度为0.4ms,采样频率为37.5KHz,每一批次回波数据相干积累采用的总调频周期数为500,积累时间为100s;接收阵列包括阵元间距为半波长的20阵元均匀线阵、8阵元均匀线阵、以及如图1所示的本发明所提出的分布式阵列,该分布式阵列与20阵元均匀线阵孔径相同,但其由两个4阵元的子阵构成,总阵元数仅为8;在距离为100km、速度为25m/s的距离多普勒单元上存在两个目标,其方位角分别为88°、92°,取该单元上距离多普勒域的单个数据快拍进行目标测向,仿真噪声背景为高斯白噪声,蒙特卡洛仿真次数为500。从图3给出的角度测量均方误差随相干积累后信噪比变化曲线可以看出:本发明采用的分布式阵列尽管阵元数仅为10,但其测向性能接近由20个阵元构成的传统雷达均匀线阵,在性能接近的条件下大大减少了接收阵元数从而节约了系统硬件成本;进一步地,在相同的阵列配置下,本发明所提方法的角度测量精度和分辨力明显高于传统的MUSIC算法。
虽然本发明已以较佳的实施例公开如上,但其并非用以限定本发明,任何熟悉此技术的人,在不脱离本发明的精神和范围内,都可以做各种改动和修饰,因此本发明的保护范围应该以权利要求书所界定的为准。
Claims (3)
1.一种用于高频地波雷达的分布式阵列目标角度测量方法,其特征在于,所述目标角度测量方法包括如下步骤:
步骤一:针对平面内由K个子阵组成的分布式雷达阵列,其中,K>1,子站内阵元间距不大于工作波长的一半,对各接收通道的雷达回波进行距离多普勒处理,针对感兴趣的目标提取相应距离多普勒单元上的L个数据快拍,L≥1,收集各子站信息用于后续的目标角度测量;
步骤二:利用所有接收通道获得的目标距离多普勒域数据,建立协方差矩阵稀疏表示模型,将目标角度估计问题建模为稀疏功率矢量重构问题,以负对数随机最大似然函数为初始代价函数,优化求解各角度网格上的信号功率;
步骤三:利用Majorization-Minimization方法,通过迭代更新并最小化一个替代的凸代价函数达到最小化原始非凸代价函数的目的,为便于求解,通过引入中间变量进一步将该问题转化为一个交替迭代优化问题;
步骤四:计算中间变量的最优解并将其代入代价函数,进一步对其求导并令导数为零,得到待求解的功率矢量的更新表达式,迭代更新中间变量和功率矢量,直至得到当前代价函数下的功率矢量最优解,再利用获得的功率矢量估计值更新代价函数;
步骤五:重复步骤四过程,直至算法收敛,获得功率矢量的最终估计值,根据该矢量中的峰值位置确定感兴趣的目标入射角度;
步骤二中协方差矩阵稀疏表示建模包括:
首先,要想有效地利用完整的阵列孔径,需对所有子阵数据进行联合处理,而接收阵列为由多个子阵构成的分布式阵列,即稀疏非均匀阵列,为解决该阵列形式下的强相关信号角度估计问题,采用基于稀疏信号表示的角度估计方法:将感兴趣的空间范围按方位角划分为Z个网格,Z>M且Z>>N,使目标信号的真实方位角包含在网格中,利用相应网格点对应的信号导向矢量构成字典A′=[a1…aZ],
其次,假设信号S的协方差矩阵为对角阵,且由于基于傅里叶变换的距离多普勒处理过程存在对偶性,可认为观测噪声N也服从高斯分布,则步骤一所得距离多普勒域数据的协方差矩阵可表示为:
其中,IM为M阶单位矩阵,P=diag(p),p=[p1…pZ+M],pz(z=1,2,…,Z)为各角度网格上的信号功率,pZ+m=σm(m=1,2,…,M)为各接收通道噪声方差,显然,仅在信号真实方位角对应的网格上pz非零,否则pz为零,因此功率矢量p是一个稀疏向量,通过重构p即可确定真实信号的方位角;
最后,考虑最小化如下负对数随机最大似然函数以求解功率矢量p:
步骤三优化问题转化过程包括:
因此,令
其中,aZ+m=em为M阶单位矩阵IM的第m列,cq为Copt的第q(q=1,2,…,Z+M)行;第三,由于pq的独立性,可将上述优化问题分解为Z+M个独立问题分别求解,分别对pq求导并令导数为零,可得第i次外循环的第j+1次内循环迭代中,第z(z=1,2,…,Z)个角度网格上的信号功率更新表达式:
以及第m(m=1,2,…,M)个接收通道的噪声方差更新表达式:
2.根据权利要求1所述一种用于高频地波雷达的分布式阵列目标角度测量方法,其特征在于,步骤一中阵列接收数据的距离多普勒处理包括:
第一步:由K个子阵构成分布式雷达接收阵列,其中K>1,第k(k=1,…,K)个子站中包含Mk个阵元,整个接收阵列包含的总阵元数为由于高频地波雷达针对的是海面目标,无需估计其俯仰角,因此阵列排布采用简单的线阵,要求子阵内部阵元间距不超过半波长,且目标对于整个接收阵列满足远场入射条件,即可认为目标到达各子阵的入射角度相同,雷达发射信号表示为s(t),则第m个阵元接收到的第n个目标回波可表示为:
smn(t)=μns(t-Δtn-τmn)
其中,m=1,2,…,M为阵元编号索引,n=1,2,…,N0为目标编号索引,N0为总目标数;μn为第n个目标信号的传播衰减因子,Δtn为该目标回波到达参考阵元时相对于发射信号的延迟,其大小与目标距离Rn和径向速度vn有关:Δtn=2(Rn-vnt)/c,c为光速;
对于高频雷达而言,其发射信号载频f0一般远大于带宽,可视为窄带信号,因此第m个阵元的接收回波可进一步写为:
对于高频雷达系统常用的脉冲截断线性调频信号和相位编码信号而言,其距离变换可分别通过快速傅里叶变换和相关处理实现,而多普勒变换往往通过快速傅里叶变换实现,鉴于以上距离多普勒处理过程是一个线性操作,对sm(t)进行距离多普勒变换后可得:
对L(L≥1)个批次的回波数据分别进行上述变换,可得(r,d)单元上L次快拍的距离多普勒域数据:
X=AS+N
3.根据权利要求1所述一种用于高频地波雷达的分布式阵列目标角度测量方法,其特征在于,步骤五最终功率矢量估计和目标角度测量过程包括:
首先,以上述步骤四的循环迭代作为所提方法的内循环求解过程,令i→i+1,重复步骤四作为所提方法的外循环求解过程,通过嵌套循环直至算法收敛,得到最终的功率矢量估计值,最后,以信号功率矢量估计值p=[p1…pZ]作为空间谱,利用其中的峰值位置确定感兴趣的(r,d)距离多普勒单元内的目标信号角度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811207040.8A CN109212526B (zh) | 2018-10-17 | 2018-10-17 | 用于高频地波雷达的分布式阵列目标角度测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811207040.8A CN109212526B (zh) | 2018-10-17 | 2018-10-17 | 用于高频地波雷达的分布式阵列目标角度测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109212526A CN109212526A (zh) | 2019-01-15 |
CN109212526B true CN109212526B (zh) | 2022-07-01 |
Family
ID=64980685
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811207040.8A Active CN109212526B (zh) | 2018-10-17 | 2018-10-17 | 用于高频地波雷达的分布式阵列目标角度测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109212526B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110231590B (zh) * | 2019-06-11 | 2023-02-28 | 东南大学 | 一种基于dft变换的阵列目标角度测量方法 |
CN112147598A (zh) * | 2019-06-27 | 2020-12-29 | 武汉爱速达机器人科技有限公司 | 一种基于直角墙面的激光标定方法 |
CN110208768A (zh) * | 2019-06-27 | 2019-09-06 | 哈尔滨工业大学 | 超并行自主产生式的相控阵雷达实时多目标回波模拟方法 |
CN111427022A (zh) * | 2020-05-08 | 2020-07-17 | 北京理工大学重庆创新中心 | 一种基于极大似然估计的阵列雷达测角方法 |
CN111693958B (zh) * | 2020-06-09 | 2024-04-19 | 周口师范学院 | 一种被动雷达三维空间测向方法及系统 |
CN112255625B (zh) * | 2020-09-01 | 2023-09-22 | 杭州电子科技大学 | 基于深度学习的二维角度依赖误差下的一维线阵测向方法 |
CN112198489B (zh) * | 2020-09-10 | 2022-02-08 | 北京理工大学 | 基于改进的最大似然算法的机扫雷达角度超分辨测角方法 |
CN112731326A (zh) * | 2020-12-24 | 2021-04-30 | 深圳大学 | 非高斯噪声下的雷达信号波形与目标角度联合估计方法 |
CN112882007B (zh) * | 2021-01-12 | 2023-06-23 | 北京交通大学 | 基于稀布阵雷达的单脉冲多目标超分辨测角实现方法 |
CN112946599B (zh) * | 2021-02-04 | 2022-05-06 | 哈尔滨工业大学(威海) | 一种基于稀布阵列的雷达空间谱估计方法 |
CN113406620B (zh) * | 2021-04-28 | 2022-08-19 | 北京理工大学 | 一种阵列分解的分布式阵列测角方法 |
CN115047456B (zh) * | 2022-08-16 | 2022-10-25 | 中国人民解放军国防科技大学 | 无人驾驶车三维全景感知方法和雷达系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103954950A (zh) * | 2014-04-25 | 2014-07-30 | 西安电子科技大学 | 一种基于样本协方差矩阵稀疏性的波达方向估计方法 |
CN104155648A (zh) * | 2014-08-26 | 2014-11-19 | 国家海洋局第一海洋研究所 | 基于阵列数据重排的高频地波雷达单次快拍music测向方法 |
-
2018
- 2018-10-17 CN CN201811207040.8A patent/CN109212526B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103954950A (zh) * | 2014-04-25 | 2014-07-30 | 西安电子科技大学 | 一种基于样本协方差矩阵稀疏性的波达方向估计方法 |
CN104155648A (zh) * | 2014-08-26 | 2014-11-19 | 国家海洋局第一海洋研究所 | 基于阵列数据重排的高频地波雷达单次快拍music测向方法 |
Non-Patent Citations (4)
Title |
---|
基于压缩感知的高频地波雷达二维DOA估计;赵春雷 等;《系统工程与电子技术》;20170430;第39卷(第4期);全文 * |
基于改进S变换的高频地波雷达阵列信号DOA估计;田潇潇等;《中国海洋大学学报(自然科学版)》;20170215(第02期);全文 * |
改进的高频天波雷达空间谱重构算法;侯煜冠等;《中国科学技术大学学报》;20080315(第03期);全文 * |
时―空欠采样下的频率和DOA联合估计算法;黄翔东等;《通信学报》;20160525(第05期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN109212526A (zh) | 2019-01-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109212526B (zh) | 用于高频地波雷达的分布式阵列目标角度测量方法 | |
CN109212527B (zh) | 用于高频地波雷达的大孔径分布式多站目标定位方法 | |
CN110208735B (zh) | 一种基于稀疏贝叶斯学习的相干信号doa估计方法 | |
CN108459307B (zh) | 基于杂波的mimo雷达收发阵列幅相误差校正方法 | |
CN108983168B (zh) | 基于孔径补全的压缩感知mimo雷达目标探测方法 | |
CN109507666B (zh) | 基于离网变分贝叶斯算法的isar稀疏频带成像方法 | |
CN112415476B (zh) | 一种基于稀疏贝叶斯学习的字典失配杂波空时谱估计方法 | |
CN109212500A (zh) | 一种基于稀疏重构的ka-stap杂噪协方差矩阵高精度估计方法 | |
CN109613532B (zh) | 一种机载雷达实时多普勒波束锐化超分辨成像方法 | |
CN113126087B (zh) | 一种星载干涉成像高度计天线 | |
CN113189592B (zh) | 考虑幅相互耦误差的车载毫米波mimo雷达测角方法 | |
CN111736131A (zh) | 一种剔除一比特信号谐波虚假目标的方法及相关组件 | |
CN112305537B (zh) | 单比特随机频控阵雷达目标距离-角度联合估计方法 | |
CN109085556B (zh) | 一种基于一二阶峰比值的高频地波雷达浪场形成方法 | |
CN115453528A (zh) | 基于快速sbl算法实现分段观测isar高分辨成像方法及装置 | |
CN112147608A (zh) | 一种快速高斯网格化非均匀fft穿墙成像雷达bp方法 | |
CN110579737A (zh) | 一种杂波环境中基于稀疏阵列的mimo雷达宽带doa计算方法 | |
CN104698448B (zh) | 运动平台下基于流形分离的共形阵列稳健估角方法 | |
CN113608218A (zh) | 一种基于后向投影原理的频域干涉相位稀疏重构方法 | |
CN111812612A (zh) | 基于子阵正交lfm信号的mimo雷达部分相关波形设计方法 | |
CN109884621B (zh) | 雷达高度表回波相参积累方法 | |
CN116577749A (zh) | 一种天线方向图未知展宽下的扫描雷达超分辨方法 | |
CN113359196B (zh) | 基于子空间法和dbf的多目标生命体征探测方法 | |
CN113671485B (zh) | 基于admm的米波面阵雷达二维doa估计方法 | |
CN110231590B (zh) | 一种基于dft变换的阵列目标角度测量方法 |
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 |