CN107957574B - 基于ifft和混合匹配追踪的时分地基mimo滑坡雷达成像方法 - Google Patents
基于ifft和混合匹配追踪的时分地基mimo滑坡雷达成像方法 Download PDFInfo
- Publication number
- CN107957574B CN107957574B CN201711453838.6A CN201711453838A CN107957574B CN 107957574 B CN107957574 B CN 107957574B CN 201711453838 A CN201711453838 A CN 201711453838A CN 107957574 B CN107957574 B CN 107957574B
- Authority
- CN
- China
- Prior art keywords
- radar
- array
- receiving
- sparse
- formula
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 76
- 238000007906 compression Methods 0.000 claims abstract description 34
- 230000006835 compression Effects 0.000 claims abstract description 34
- 238000012544 monitoring process Methods 0.000 claims abstract description 32
- 238000000034 method Methods 0.000 claims abstract description 31
- 238000003491 array Methods 0.000 claims abstract description 21
- 238000012545 processing Methods 0.000 claims abstract description 19
- 239000011159 matrix material Substances 0.000 claims description 36
- 238000007781 pre-processing Methods 0.000 claims description 17
- 238000012937 correction Methods 0.000 claims description 9
- 238000013144 data compression Methods 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 6
- 230000009466 transformation Effects 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 5
- 230000001131 transforming effect Effects 0.000 claims description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 230000003321 amplification Effects 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000013508 migration Methods 0.000 claims description 3
- 230000005012 migration Effects 0.000 claims description 3
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 3
- 230000009191 jumping Effects 0.000 claims 1
- 230000000694 effects Effects 0.000 description 10
- 238000010586 diagram Methods 0.000 description 5
- 238000004088 simulation Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 239000000654 additive Substances 0.000 description 2
- 230000000996 additive effect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000012271 agricultural production Methods 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
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
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了基于IFFT和混合匹配追踪的时分地基MIMO滑坡雷达成像方法,用于解决现有山体滑坡监测领域成像方法所用天线数目多、数据处理复杂度高的问题,并针对多目标成像的伪影点问题引入混合匹配追踪算法,大大改善成像质量。实现步骤为:构造等效随机稀疏收发阵列;建立回波信号采集模型;距离向压缩处理;方位向压缩处理;完成二维目标高分辨成像。本发明具有以下优点:减少阵列数目,降低硬件成本;降低了数据处理复杂度;提高了距离向和方位向分辨率。
Description
技术领域
本发明属于山体滑坡监测雷达技术领域,涉及基于IFFT和混合匹配追踪的时分地基MIMO滑坡雷达成像方法。
背景技术
山体滑坡是常见的地质灾害之一,对铁路、公路建设及维护和工农业生产造成了巨大的危害。据不完全统计,每年自然灾害造成的人员死亡和失踪人数中约有六分之一来自于滑坡灾害,给人民群众生命财产安全造成巨大损失。因此山体滑坡实时监测研究是近年来灾害预防领域重要的研究方向,也具有现实意义。
地基雷达由于其稳定可控且可长期重复观测的特点,近年来已成为探测露天矿边坡形变、城区沉降、山体滑坡等领域的重要应用方向。本发明主要针对山体滑坡监测的关键技术——地基雷达成像进行研究。现有两大典型成像***应用在山体滑坡监测领域:意大利的Joint Research Centre研制的LISA(Linear SAR)成像***、佛罗伦萨大学和意大利IDS公司合作研制的GBInSAR(Ground Based InSAR)***IBIS(Image By InterferometrySurvey)。这些成像***通过雷达天线在水平直轨道运动获得较高的方位向分辨率,但是天线运动的同时会影响测量的精度,也限制了雷达数据的采集速度。
MIMO雷达利用其多发多收体制形成的虚拟天线阵元替代实际天线阵元,可以很好的解决合成孔径的问题。然而传统的成像算法,如距离多普勒算法和CS(Chirp Scaling)算法由于雷达结构的特殊性不再适用。后向投影算法,又称BP成像算法,是一种时域成像算法,具有很强的适应性,不会束缚MIMO雷达天线阵列的自由布阵。BP算法大量的运用在MIMO雷达成像,然而标准BP算法的成像效率较低。把波束形成算法引入BP成像中,可以进一步提高成像算法效率,但需要的天线数目仍然较多。从现有的文献来看,在山体滑坡监测雷达领域,采用的仍是经典的均匀线性布阵方式,稀疏阵列下的地基MIMO雷达成像技术仍待进一步研究。
发明内容
本发明的目的是提供基于IFFT和混合匹配追踪的时分地基MIMO滑坡雷达成像方法,用于解决现有山体滑坡监测领域成像方法所用天线数目多、数据处理复杂度高的问题,并针对多目标成像的伪影点问题引入混合匹配追踪算法,大大改善成像质量。
其原理是构造等效随机稀疏收发阵列,通过分时发射和接收步进频连续波信号建立时分地基MIMO雷达回波信号采集模型,对采集到的回波信号通过逆傅里叶变换完成距离向压缩,构造合适的测量矩阵和基于时延补偿因子的稀疏基,利用混合匹配追踪算法对距离压缩完的数据进行方位向压缩实现目标的二维高分辨成像。
为实现上述目的,本发明采用的以下技术方案。
基于IFFT和混合匹配追踪的时分地基MIMO滑坡雷达成像方法,包括以下步骤:
步骤1)构造等效随机稀疏收发阵列,所述构造等效随机稀疏收发阵列的方法为MIMO雷达是M发N收的均匀阵列,原始接收天线间隔为λ/2,发射天线间隔为N(λ/2),λ是发射波波长,并指定接收阵列稀疏布置,从原始N个接收阵列中随机选取N'个,为保证阵列孔径长度足够大,第一个和最后一个接收天线位置不变,根据得到的N'个稀疏接收阵元和M个均匀的发射阵元,利用多项式理论得到MN'个虚拟阵元;
步骤2)构造时分地基MIMO雷达回波信号采集模型,所述的信号为分时发射和接收步进频连续波信号,所述分时发射和接收步进频连续波信号的方法包括以下步骤:
步骤2-1,由信号产生单元产生步进频连续波信号,设发射信号中单个脉冲表示为S(t),
其中fi=f0+iΔf,f0为脉冲起始频率,Δf为频率步进量,rect(t)为单位矩形函数,Tr为发射信号脉冲宽度,Q为子脉冲个数;
步骤2-2,数字基带步进频信号经过数模转换后送入信号调制单元,雷达信号被调制到雷达工作频段后通过射频放大单元将信号功率放大,再由Tx分时选择器控制MIMO发射天线阵列分时将信号发射出去;
步骤2-3,分时选择器Rx控制MIMO接收天线阵列分时接收雷达回波信号,接收到的雷达回波信号非常微弱需放大后再进行正交解调,假设第n'个接收阵元接收到第m个发射阵元发射信号的回波,并将该通道称为第mn'个观测通道,将该观测通道回波信号幅度放大并经过解调采样后为,
其中,δ为目标散射系数,τmn'为该观测通道的时延,(x0,y0)、(xn',0)、(xm,0)分别为目标和收、发天线的坐标,c表示光速;
步骤2-4,对经过解调采样后的回波信号进行数据采集,将采集好的雷达数据送入信号处理单元进行后续的成像处理;
步骤3)通过采集到的回波信号完成距离向压缩,并对得到的压缩数据进行合理的相位校正,消除相位的平方项,完成方位向预处理,所述距离向压缩的方法为逆傅里叶变换,并对得到的压缩数据进行合理的相位校正,消除相位的平方项,完成方位向预处理,所述方位向预处理的方法包括以下步骤:
步骤3-1,雷达回波信号经过解调采样后进行数据采集得到一个原始雷达数据矩阵I[MN',Q],在小间距布阵条件下,MIMO雷达阵列长度远小于监测区域目标到阵列的距离,距离徙动不明显,因此,距离向只需要进行逆傅里叶变换(IFFT)即可得到较为准确的压缩数据,对I[MN',Q]矩阵按行分别进行IFFT得到第mn'行数据为,
步骤3-2,将式(3)所示的时延曲线按幂级数展开,省略高方次项,为方便表达,把完成距离向压缩后的雷达数据变换到极坐标下表示,监测区域目标(x0,y0)的波达角设为θ,根据真实滑坡监测区域的几何特点推出θ很小,sin2θ值非常小,sin2θ/R0的值更小,可以忽略不计,时延曲线将合理的近似为,
步骤3-3,对比传统的线性阵列,MIMO天线阵列稀疏布置的情况下会存在相位不连续现象,为提高雷达数据方位向压缩的质量,在方位向压缩之前进行预处理,即采用相位相乘的方法消除式(4)的第二项,因为R0是个变量,会造成算法的运算量巨大,所以对式(4)的第二项进行近似的校正,把式(4)第二项中的R0替换为R1,R1表示监测区域中心到雷达天线阵列坐标中心的距离,为一常量,因此对式(3)乘以校正因子完成方位向压缩预处理,那么得到,
步骤4)构造合适的观测矩阵和基于时延补偿因子的稀疏基,将经过距离向压缩的数据进行方位向压缩,对所有距离单元的数据进行重构处理,完成方位向数据压缩,所述方位向压缩的方法为混合匹配追踪算法,所述基于时延补偿因子的稀疏基的方法包括以下步骤:
步骤4-1,由步骤3处理得到的所有观测通道数据集合可看作一个大小为Q×MN'的矩阵,其中Q为距离单元数,MN'为等效稀疏收发虚拟天线个数,假设等效均匀收发天线得到的回波数据矩阵记为Y,那么X的每一行均可看作Y每一行的低维观测,其中观测矩阵Φ={φu,v}为广义单位阵,且
Φ矩阵由MN'行MN列元素构成,矩阵的每一行除了第δu个元素为1外,其余的全为0,δu由等效稀疏收发天线位置决定;
步骤4-2,根据滑坡监测区域几何形状,得出监测区域波达角的覆盖范围,雷达监测区域距离向从ymin到ymax,方位向从-x到x,则波达角的范围为(-arctan(x/ymin),arctan(x/ymin)),将其均分为MN份,得θl如式,
令Wl,mn=exp(-j2πf0sinθl(xm+xn)/c),可由Wl,mn组建MN行MN列时延补偿因子矩阵如式,
时延补偿因子与信号本身在某种程度上存在一一对应的关系,很全面的反应了散射点的特征信息,故可以把其作为稀疏变换的稀疏基进行后续的数据压缩处理;
步骤4-3,本文在构造稀疏接收阵列时,接收天线是随机选取的,等效的虚拟收发阵元的位置也是随机的,满足广义单位阵Φ与Ψ不相干的要求,因此Φ·Ψ满足RIP(Restricted Isometry Property)条件,由以上构造的稀疏变换矩阵和观测矩阵,利用混合正交匹配追踪算法对X的第q行求解以下问题:
然后对所有的q=1,2,...,Q求解式(9)问题,即可得到二维像,所述混合匹配追踪算法求解的方法包括以下步骤:
步骤4-4-1,采用标准正交匹配追踪算法得到稀疏解的初值,
δomp=omp(xq,T,K) (10),
步骤4-4-2,定义函数表示max_ind(y,P):返回P索引,找到向量y中幅值最大的元素所在的位置,由式(10)得到的稀疏解初值可定义初始支撑集为,
Λold=max_ind(|δomp|,K) (11),
步骤4-4-3,由式(11)得到的初始支撑集可以将残差初始化为,
其中xq为回波采样向量,恢复矩阵T=ΦΨ,K为稀疏度值;
步骤4-4-4,依然采用标准正交匹配追踪算法对式(12)所求的残差进行求值,可以得到,
步骤4-4-6,将原始测量信号向这2K个支撑集构成的子空间投影,可以得到更新的支撑集为,
步骤4-4-7,利用式(14)得到的新的支撑集对残差进行更新,
步骤5)获得目标的二维高分辨成像。
结合仿真实验,本发明的技术效果体现为:
根据山体滑坡的实际成像场景设置雷达仿真参数,雷达发射步进频信号,载频15GHz,带宽300MHz,子脉冲数4096个,这里按照步骤1稀疏布阵方法布置6发25收的阵列,均匀接收阵元的间隔为λ/2,发射阵元的间隔为50×(λ/2),阵列稀疏度为50%,雷达成像区域为距离向从1000m到2000m,方位向为-150m到150m,点目标的散射系数都设为1,附加噪声为加性高斯白噪声,信噪比为10dB;
首先分析单点目标成像,目标方位向为0m,距离向为1500m,传统波束成形算法、标准正交匹配追踪算法以及本发明提出的算法成像结果,三种方法均对单点聚焦正确成像,表明对于单点目标成像,标准正交匹配追踪算法和本发明成像效果一致,明显优于传统波束形成算法,标准正交匹配追踪算法对于单点目标成像没有伪影点出现;
其次分析多点目标成像,目标点位置信息表明了标准正交匹配追踪算法和本发明在稀疏阵列下多点目标聚焦成像的效果优于传统波束成形算法,但是存在明显的伪影点,影响多个目标之间的区分,这是由于标准正交匹配追踪算法在基信号选择时只能扩充而不能去除的策略造成的,本发明与标准正交匹配追踪算法不同之处在于迭代过程中结合了子空间追踪算法,去除了在前面迭代过程中被选择的病态索引,因此本发明的成像效果更好;
对比传统波束形成算法、标准正交匹配追踪算法、本发明的分辨率和目标杂波比(TCR)指标对比分析,目标杂波比对应用来衡量目标在背景杂波中的凸显程度,三种方法的距离向分辨率基本相同,这是因为距离向分辨率只与发射信号带宽有关,标准正交匹配追踪算法和本发明的方位向分辨率明显优于传统波束形成算法,这是因为阵列的稀疏先验信息用在压缩感知理论可以大大改善成像的质量,本发明的TCR值最高,点目标的散射特性保存最好,比标准正交匹配追踪算法的TCR值高4dB左右,极大的改善了多目标成像的伪影点问题。
与现有技术相比,本发明具有以下优点:
(1)在山体滑坡雷达监测领域,采用的仍是传统的均匀线性布阵方式,本发明采用稀疏布阵方式,大大减少了收发阵列的数目,降低了硬件成本;
(2)本发明采用的是较少的发射阵列均匀布阵,接收阵列在均匀布阵的基础上随机选取一半的位置数目作为接收阵列,从一定程度上降低了数据处理复杂度;
(3)在雷达天线数目低于传统均匀阵列的情况下,实现高分辨率成像。本发明通过发射步进频连续波信号获得高的距离向分辨率,方位向引入混合匹配追踪算法大大改善了多目标成像的伪影点问题,进一步提高了方位向分辨率。
附图说明:
图1为实施例MIMO稀疏布阵结构示意图;
图2为实施例稀疏阵列单点目标传统波束形成算法成像结果;
图3为实施例稀疏阵列单点目标标准正交匹配追踪算法成像结果;
图4为实施例稀疏阵列单点目标本发明算法成像结果;
图5为实施例稀疏阵列多目标传统波束形成算法成像结果;
图6为实施例稀疏阵列多目标标准正交匹配追踪算法成像结果;
图7为实施例稀疏阵列多目标本发明算法成像结果。
具体实施方式
本发明通过实施例,结合说明书附图对本发明内容作进一步详细说明,但不是对本发明的限定。
实施例
基于IFFT和混合匹配追踪的时分地基MIMO滑坡雷达成像方法:
步骤1,构造等效随机稀疏收发阵列,
根据山体滑坡监测场景设计满足要求的收发阵列,传统MIMO雷达是M发N收的均匀阵列,原始接收天线间隔为λ/2,发射天线间隔为N(λ/2),λ是发射波波长,并指定接收阵列稀疏布置,从原始N个接收阵列中随机选取N'个,为保证阵列孔径长度足够大,第一个和最后一个接收天线位置不变。根据得到的N'个稀疏接收阵元和M个均匀的发射阵元,利用多项式理论得到MN'个虚拟阵元,以2发6收的阵列为例,示意图如图1所示;
步骤2,构造时分地基MIMO雷达回波信号采集模型,
步骤2-1,由信号产生单元产生步进频连续波信号,设发射信号中单个脉冲表示为S(t),如(1)式,
其中fi=f0+iΔf,f0为脉冲起始频率,Δf为频率步进量,rect(t)为单位矩形函数,Tr为发射信号脉冲宽度,Q为子脉冲个数;
步骤2-2,数字基带步进频信号经过数模转换后送入信号调制单元,雷达信号被调制到雷达工作频段后通过射频放大单元将信号功率放大,再由Tx分时选择器控制MIMO发射天线阵列分时将信号发射出去,
步骤2-3,分时选择器Rx控制MIMO接收天线阵列分时接收雷达回波信号,接收到的雷达回波信号非常微弱需放大后再进行正交解调,假设第n'个接收阵元接收到第m个发射阵元发射信号的回波,并将该通道称为第mn'个观测通道,将该观测通道回波信号幅度放大并经过解调采样后为,
其中,δ为目标散射系数,τmn'为该观测通道的时延,(x0,y0)、(xn',0)、(xm,0)分别为目标和收、发天线的坐标,c表示光速;
步骤2-4,对经过解调采样后的回波信号进行数据采集,将采集好的雷达数据送入信号处理单元进行后续的成像处理;
步骤3)通过采集到的回波信号进行逆傅里叶变换完成距离向压缩,并对得到的压缩数据进行合理的相位校正,消除相位的平方项,完成方位向预处理,
传统MIMO雷达需要设计具有良好正交性的雷达信号,接收端通过设计相应的匹配滤波器进行匹配滤波,分离出MIMO雷达天线阵列所形成的多个观测通道信号,分选出各雷达各观测信号的同时,也实现了雷达数据距离向压缩,而本发明采用时分MIMO雷达体制,雷达通过时分方式发射和接收雷达信号,具体雷达工作方式在步骤2中进行了阐述,基于时分MIMO技术的地基雷达接收完数据后,各通道数据被很好的分离开来,但是距离向雷达数据还没有进行压缩处理,步进频连续波是一个频域信号,变换到时域可以形成一个sinc函数信号,即可以形成一个窄脉冲,斜距不同的目标窄脉冲出现的位置也不一样,从而在距离向分离出不同目标,所以本发明采用逆傅立叶变换算法对雷达数据进行距离向压缩,
步骤3-1,雷达回波信号经过解调采样后进行数据采集得到一个原始雷达数据矩阵I[MN',Q],在小间距布阵条件下,MIMO雷达阵列长度远小于监测区域目标到阵列的距离,距离徙动不明显,因此,距离向只需要进行逆傅里叶变换(IFFT)即可得到较为准确的压缩数据,对I[MN',Q]矩阵按行分别进行IFFT得到第mn'行数据为,
步骤3-2,将式(3)所示的时延曲线按幂级数展开,省略高方次项,为方便表达,把完成距离向压缩后的雷达数据变换到极坐标下表示,监测区域目标(x0,y0)的波达角设为θ,根据真实滑坡监测区域的几何特点推出θ很小,sin2θ值非常小,sin2θ/R0的值更小,可以忽略不计,时延曲线将合理的近似为,
步骤3-3,对比传统的线性阵列,MIMO天线阵列稀疏布置的情况下会存在相位不连续现象,为提高雷达数据方位向压缩的质量,在方位向压缩之前进行预处理,即采用相位相乘的方法消除式(4)的第二项,因为R0是个变量,会造成算法的运算量巨大,所以对式(4)的第二项进行近似的校正,把式(4)第二项中的R0替换为R1,R1表示监测区域中心到雷达天线阵列坐标中心的距离,为一常量,因此对式(3)乘以校正因子完成方位向压缩预处理,那么得到,
步骤4)构造合适的观测矩阵和基于时延补偿因子的稀疏基,将经过距离向压缩的数据进行方位向压缩,
若等效的虚拟收发天线是均匀布置的,在方位向用传统的波束形成算法进行数据压缩,即可得到高分辨率的图像,但本发明等效的虚拟收发天线是稀疏布置的,直接运用传统方法就会导致旁瓣很高,由于对山体滑坡监测区域目标成像时,聚焦性强的目标只占据监测区域的极小部分,满足压缩感知理论的稀疏性要求,稀疏阵列可认为是原始均匀阵列的低维观测,因此本发明结合压缩感知理论对方位向数据进行压缩处理,针对多目标的伪影点问题,引入混合匹配追踪算法,既保证了基信号选择的正交性又对支撑集选取过程进行回溯优化,大大改善了传统正交匹配追踪(Orthogonal Matching Pursuit,OMP)算法的成像质量,
步骤4-1,由步骤3处理得到的所有观测通道数据集合可看作一个大小为Q×MN'的矩阵,其中Q为距离单元数,MN'为等效稀疏收发虚拟天线个数,假设等效均匀收发天线得到的回波数据矩阵记为Y,那么X的每一行均可看作Y每一行的低维观测,其中观测矩阵Φ={φu,v}为广义单位阵,且
Φ矩阵由MN'行MN列元素构成,矩阵的每一行除了第δu个元素为1外,其余的全为0,δu由等效稀疏收发天线位置决定;
步骤4-2,根据滑坡监测区域几何形状,得出监测区域波达角的覆盖范围,雷达监测区域距离向从ymin到ymax,方位向从-x到x,则波达角的范围为(-arctan(x/ymin),arctan(x/ymin)),将其均分为MN份,得θl如式,
令Wl,mn=exp(-j2πf0sinθl(xm+xn)/c),可由Wl,mn组建MN行MN列时延补偿因子矩阵如式,
时延补偿因子与信号本身在某种程度上存在一一对应的关系,很全面的反应了散射点的特征信息,故可以把其作为稀疏变换的稀疏基进行后续的数据压缩处理;
步骤4-3,本文在构造稀疏接收阵列时,接收天线是随机选取的,等效的虚拟收发阵元的位置也是随机的,满足广义单位阵Φ与Ψ不相干的要求,因此Φ·Ψ满足RIP(Restricted Isometry Property)条件,由以上构造的稀疏变换矩阵和观测矩阵,利用混合正交匹配追踪算法对X的第q行求解以下问题,
然后对所有的q=1,2,...,Q求解式(9)问题,即可得到二维像;
步骤4-4,利用混合匹配追踪算法对式(9)进行求解的具体步骤如下,
步骤4-4-1,采用标准正交匹配追踪算法得到稀疏解的初值,
δomp=omp(xq,T,K) (10),
步骤4-4-2,定义函数表示max_ind(y,P):返回P索引,找到向量y中幅值最大的元素所在的位置,由式(10)得到的稀疏解初值可定义初始支撑集为,
Λold=max_ind(|δomp|,K) (11),
步骤4-4-3,由式(11)得到的初始支撑集可以将残差初始化为,
其中xq为回波采样向量,恢复矩阵T=ΦΨ,K为稀疏度值;
步骤4-4-4,依然采用标准正交匹配追踪算法对式(12)所求的残差进行求值,可以得到,
步骤4-4-6,将原始测量信号向这2K个支撑集构成的子空间投影,可以得到更新的支撑集为,
步骤4-4-7,利用式(14)得到的新的支撑集对残差进行更新,
步骤5)获得目标的二维高分辨成像。
以下结合仿真实验,对本发明的技术效果作进一步说明,
根据山体滑坡的实际成像场景设置雷达仿真参数,雷达发射步进频信号,载频15GHz,带宽300MHz,子脉冲数4096个,这里按照步骤1稀疏布阵方法布置6发25收的阵列,均匀接收阵元的间隔为λ/2,发射阵元的间隔为50×(λ/2),阵列稀疏度为50%,雷达成像区域为距离向从1000m到2000m,方位向为-150m到150m,点目标的散射系数都设为1,附加噪声为加性高斯白噪声,信噪比为10dB;
首先分析单点目标成像,目标方位向为0m,距离向为1500m,图2、图3、图4给出了传统波束形成算法、标准正交匹配追踪算法以及本发明提出的算法成像结果,从图2、图3、图4看出,三种方法均对单点聚焦正确成像,但是图2有低幅度值的旁瓣,图3和图4没有这种情况,且图2的分辨率低于图3和图4,表明对于单点目标成像,标准正交匹配追踪算法和本发明成像效果一致,明显优于传统波束形成算法,标准正交匹配追踪算法对于单点目标成像没有伪影点出现;
其次分析多点目标成像,目标点位置信息为:tg1-(0,1500),tg2-(-25,1505),tg3-(25,1505),tg4-(0,1510),tg5-(-50,1520),tg6-(50,1520),tg7-(0,1520),图5给出了传统波束形成算法效果图,图6给出了标准正交匹配追踪算法的成像效果图,图7给出了本发明成像效果图,可以看出,图5中由于阵元的缺失,不能很好的处理欠采样数据,呈现出的7个点目标图像模糊,只能大致辨别几个目标点的不同位置信息,且旁瓣水平比单点目标成像恶化明显,图6和图7表明了标准正交匹配追踪算法和本发明在稀疏阵列下多点目标聚焦成像的效果优于传统波束成形算法,但是图6存在明显的伪影点,影响多个目标之间的区分,这是由于标准正交匹配追踪算法在基信号选择时只能扩充而不能去除的策略造成的,本发明与标准正交匹配追踪算法不同之处在于迭代过程中结合了子空间追踪算法,去除了在前面迭代过程中被选择的病态索引,因此本发明的成像效果更好,为进一步分析本发明成像质量,结果见表1,
表1不同成像方法的分辨率和目标杂波比
表1分别给出了传统波束形成算法、标准正交匹配追踪算法、本发明的分辨率和目标杂波比(TCR)指标对比分析,目标杂波比对应用来衡量目标在背景杂波中的凸显程度,其值越大,聚焦度越高,以dB的形式定义为,
式中Γ为目标区域,Ω为杂波区域,分析表1可知,三种方法的距离向分辨率基本相同,这是因为距离向分辨率只与发射信号带宽有关,正交匹配追踪算法和本发明的方位向分辨率明显优于传统波束形成算法,这是因为阵列的稀疏先验信息用在压缩感知理论可以大大改善成像的质量,从表1可以看出,本发明的TCR值最高,点目标的散射特性保存最好,比标准正交匹配追踪算法的TCR值高4dB左右,极大的改善了多目标成像的伪影点问题。
Claims (5)
1.基于IFFT和混合匹配追踪的时分地基MIMO滑坡雷达成像方法,其特征在于包括以下步骤:
步骤1)构造等效随机稀疏收发阵列;所述步骤1)构造等效随机稀疏收发阵列的方法为MIMO雷达是M发N收的均匀阵列,原始接收天线间隔为λ/2,发射天线间隔为N(λ/2),λ是发射波波长,并指定接收阵列稀疏布置,从原始N个接收阵列中随机选取N'个,为保证阵列孔径长度足够大,第一个和最后一个接收天线位置不变,根据得到的N'个稀疏接收阵元和M个均匀的发射阵元,利用多项式理论得到MN'个虚拟阵元;
步骤2)构造时分地基MIMO雷达回波信号采集模型;所述步骤2)为分时发射和接收步进频连续波信号;
步骤3)通过采集到的回波信号完成距离向压缩,并对得到的压缩数据进行相位校正,消除相位的平方项,完成方位向预处理;所述步骤3)距离向压缩的方法为逆傅里叶变换,并对得到的压缩数据进行相位校正,消除相位的平方项,完成方位向预处理;
步骤4)构造观测矩阵和基于时延补偿因子的稀疏基,将经过方位向预处理的数据进行方位向压缩,对所有距离单元的数据进行重构处理,完成方位向数据压缩;所述步骤4)方位向压缩的方法为混合匹配追踪算法,采用的是基于时延补偿因子的稀疏基;
步骤5)获得目标的二维高分辨成像。
2.根据权利要求1所述的方法,其特征在于:所述分时发射和接收步进频连续波信号的方法包括以下步骤:
步骤2-1,由信号产生单元产生步进频连续波信号,设发射信号中单个脉冲表示为S(t),
其中fi=f0+iΔf,f0为脉冲起始频率,Δf为频率步进量,rect(t)为单位矩形函数,Tr为发射信号脉冲宽度,Q为子脉冲个数;
步骤2-2,数字基带步进频信号经过数模转换后送入信号调制单元,雷达信号被调制到雷达工作频段后通过射频放大单元将信号功率放大,再由Tx分时选择器控制MIMO发射天线阵列分时将信号发射出去;
步骤2-3,分时选择器Rx控制MIMO接收天线阵列分时接收雷达回波信号,接收到的雷达回波信号非常微弱需放大后再进行正交解调,假设第n'个接收阵元接收到第m个发射阵元发射信号的回波,并将第n'个接收阵元接收到第m个发射阵元发射信号的回波对应的通道称为第mn’个观测通道,将该观测通道回波信号幅度放大并经过解调采样后为,
i∈[0,1,2,...,Q-1],
其中,δ为目标散射系数,τmn′为该观测通道的时延,(x0,y0)、(xn′,0)、(xm,0)分别为目标和收、发天线的坐标,c表示光速;
步骤2-4,对经过解调采样后的回波信号进行数据采集,将采集好的雷达数据送入信号处理单元进行后续的成像处理。
3.根据权利要求1所述的方法,其特征在于:所述距离向压缩和方位向预处理的方法包括以下步骤:
步骤3-1,雷达回波信号经过解调采样后进行数据采集得到一个原始雷达数据矩阵I[MN′,Q],在小间距布阵条件下,MIMO雷达阵列长度远小于监测区域目标到阵列的距离,距离徙动不明显,因此,距离向只需要进行逆傅里叶变换(IFFT)即可得到较为准确的压缩数据,对I[MN′,Q]矩阵按行分别进行IFFT得到第mn'行数据为,
步骤3-2,将式(3)所示的时延曲线按幂级数展开,省略高方次项,为方便表达,把完成距离向压缩后的雷达数据变换到极坐标下表示,监测区域目标(x0,y0)的波达角设为θ,根据真实滑坡监测区域的几何特点推出θ很小,sin2θ值非常小,sin2θ/R0的值更小,可以忽略不计,时延曲线将合理的近似为,
步骤3-3,MIMO天线阵列稀疏布置的情况下会存在相位不连续现象,为提高雷达数据方位向压缩的质量,在方位向压缩之前进行预处理,即采用相位相乘的方法消除式(4)的第二项,因为R0是个变量,会造成算法的运算量巨大,所以对式(4)的第二项进行近似的校正,把式(4)第二项中的R0替换为R1,R1表示监测区域中心到雷达天线阵列坐标中心的距离,为一常量,因此对式(3)乘以校正因子完成方位向压缩预处理,那么得到,
Q表示子脉冲个数,m表示发射阵元序数,n'表示接收阵元的序数,f0表示脉冲起始频率,τmn'表示第mn'个通道产生的时延,λ代表波长,c代表光速。
4.根据权利要求1所述的方法,其特征在于:所述基于时延补偿因子的稀疏基的方法包括以下步骤:
步骤4-1,由步骤3)处理得到的所有观测通道数据集合X={Sr(tk,x1)H,...,Sr(tk,xMN')H}作为一个大小为Q×MN′的矩阵,其中Q为距离单元数,MN'为等效稀疏收发虚拟天线个数,假设等效均匀收发天线得到的回波数据矩阵记为Y,那么X的每一行均作为Y每一行的低维观测,其中观测矩阵Φ={φu,v}为广义单位阵,且
Φ矩阵由MN'行MN列元素构成,矩阵的每一行除了第δu个元素为1外,其余的全为0,δu由等效稀疏收发天线位置决定;
步骤4-2,根据滑坡监测区域几何形状,得出监测区域波达角的覆盖范围,雷达监测区域距离向从ymin到ymax,方位向从-x到x,则波达角的范围为(-arctan(x/ymin),arctan(x/ymin)),将其均分为MN份,得θl如式,
令Wl,mn=exp(-j2πf0sinθl(xm+xn)/c),由Wl,mn组建MN行MN列时延补偿因子矩阵如式,
时延补偿因子与信号本身在某种程度上存在一一对应的关系,很全面的反应了散射点的特征信息,故把其作为稀疏变换的稀疏基进行后续的数据压缩处理,
步骤4-3,在构造稀疏接收阵列时,接收天线是随机选取的,等效的虚拟收发阵元的位置也是随机的,满足广义单位阵Φ与Ψ不相干的要求,因此Φ·Ψ满足RIP条件,所述RIP为Restricted Isometry Property,由以上构造的稀疏变换矩阵和观测矩阵,利用混合正交匹配追踪算法对X的第q行求解以下问题:
然后对所有的q=1,2,...,Q求解式(9)问题,即可得到二维像;
步骤4-4,利用混合匹配追踪算法对式(9)进行求解,m表示发射阵元的序数,n表示接收阵元的序数,f0表示脉冲起始频率,c代表光速。
5.根据权利要求4所述的方法,其特征在于:所述混合匹配追踪算法求解的方法包括以下步骤:
步骤4-4-1,采用标准正交匹配追踪算法得到稀疏解的初值,
δomp=omp(xq,T,K) (10),
步骤4-4-2,定义函数表示max_ind(y,P):返回P索引,找到向量y中幅值最大的元素所在的位置,由式(10)得到的稀疏解初值可定义初始支撑集为,
Λold=max_ind(|δomp|,K) (11),
步骤4-4-3,由式(11)得到的初始支撑集将残差初始化为,
其中xq为回波采样向量,恢复矩阵T=ΦΨ,K为稀疏度值;
步骤4-4-4,依然采用标准正交匹配追踪算法对式(12)所求的残差进行求值,得到,
步骤4-4-6,将原始测量信号向这2K个支撑集构成的子空间投影得到更新的支撑集为,
步骤4-4-7,利用式(14)得到的新的支撑集对残差进行更新,
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711453838.6A CN107957574B (zh) | 2017-12-28 | 2017-12-28 | 基于ifft和混合匹配追踪的时分地基mimo滑坡雷达成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711453838.6A CN107957574B (zh) | 2017-12-28 | 2017-12-28 | 基于ifft和混合匹配追踪的时分地基mimo滑坡雷达成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107957574A CN107957574A (zh) | 2018-04-24 |
CN107957574B true CN107957574B (zh) | 2021-08-20 |
Family
ID=61957011
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711453838.6A Active CN107957574B (zh) | 2017-12-28 | 2017-12-28 | 基于ifft和混合匹配追踪的时分地基mimo滑坡雷达成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107957574B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108693511B (zh) * | 2018-05-25 | 2020-06-30 | 中国人民解放军国防科技大学 | 时分复用mimo雷达的运动目标角度计算方法 |
CN108710116B (zh) * | 2018-05-25 | 2021-01-12 | 中国人民解放军国防科技大学 | 一种mimo雷达的运动目标相位恢复方法 |
CN110011076B (zh) * | 2019-03-13 | 2020-11-13 | 成都聚利中宇科技有限公司 | 一种周期排列的稀疏阵列天线及排列方法 |
CN110609284B (zh) * | 2019-09-23 | 2023-08-18 | 成都盟升科技有限公司 | 一种sar中断数据恢复成像的方法 |
CN111239733B (zh) * | 2020-02-13 | 2023-03-28 | 电子科技大学 | 一种基于极化方式的短时移动正交信号宽幅成像方法 |
CN112014816B (zh) * | 2020-08-20 | 2022-09-23 | 桂林电子科技大学 | 基于改进层位追踪算法的双程走时计算方法 |
CN112986983B (zh) * | 2021-02-06 | 2023-12-19 | 中国人民解放军国防科技大学 | 一种全息穿透成像雷达杂波抑制方法 |
CN113311402B (zh) * | 2021-04-16 | 2022-03-01 | 苏州迈斯维通信技术有限公司 | 考虑发射波形非理想正交的mimo雷达目标测向方法及*** |
CN113835090B (zh) * | 2021-08-31 | 2024-04-12 | 电子科技大学 | 一种基于多通道sar***的高精度干涉相位获取方法 |
CN117480468A (zh) * | 2022-05-30 | 2024-01-30 | 京东方科技集团股份有限公司 | 阵列天线及人机交互装置 |
CN115480244B (zh) * | 2022-08-09 | 2024-05-31 | 中国人民解放军63921部队 | 一种同步轨道空间目标的二维成像方法和*** |
CN115166682B (zh) * | 2022-09-08 | 2022-11-22 | 南京安麦森电子科技有限公司 | 一种手持便携式雷达目标散射特性检测方法 |
CN116577782B (zh) * | 2023-04-23 | 2024-01-12 | 珠海微度芯创科技有限责任公司 | 一种毫米波一维线型阵列快速图像重建方法及*** |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104991249A (zh) * | 2015-08-12 | 2015-10-21 | 桂林电子科技大学 | 一种山体滑坡的mimo雷达监测***及监测方法 |
CN105652273A (zh) * | 2016-03-17 | 2016-06-08 | 哈尔滨工程大学 | 一种基于混合匹配追踪算法的mimo雷达稀疏成像算法 |
EP3094987A1 (en) * | 2014-01-13 | 2016-11-23 | Mitsubishi Electric Corporation | Method and system for reconstructing scene behind wall |
CN106249233A (zh) * | 2016-09-21 | 2016-12-21 | 清华大学 | 面向多极化穿墙雷达成像的前向混合匹配追踪方法和*** |
CN106338731A (zh) * | 2016-08-22 | 2017-01-18 | 西安电子科技大学 | 一种多子带并发的mimo‑sar雷达成像方法 |
CN106680815A (zh) * | 2016-12-15 | 2017-05-17 | 哈尔滨工程大学 | 基于张量稀疏表示的mimo雷达成像方法 |
-
2017
- 2017-12-28 CN CN201711453838.6A patent/CN107957574B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3094987A1 (en) * | 2014-01-13 | 2016-11-23 | Mitsubishi Electric Corporation | Method and system for reconstructing scene behind wall |
CN104991249A (zh) * | 2015-08-12 | 2015-10-21 | 桂林电子科技大学 | 一种山体滑坡的mimo雷达监测***及监测方法 |
CN105652273A (zh) * | 2016-03-17 | 2016-06-08 | 哈尔滨工程大学 | 一种基于混合匹配追踪算法的mimo雷达稀疏成像算法 |
CN106338731A (zh) * | 2016-08-22 | 2017-01-18 | 西安电子科技大学 | 一种多子带并发的mimo‑sar雷达成像方法 |
CN106249233A (zh) * | 2016-09-21 | 2016-12-21 | 清华大学 | 面向多极化穿墙雷达成像的前向混合匹配追踪方法和*** |
CN106680815A (zh) * | 2016-12-15 | 2017-05-17 | 哈尔滨工程大学 | 基于张量稀疏表示的mimo雷达成像方法 |
Non-Patent Citations (4)
Title |
---|
Hybrid matching pursuit for distributed through-wall radar imaging;Gang Li, et al.;《IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION》;20150430;p1701-1711 * |
基于压缩感知的稀疏阵列 MIMO 雷达成像方法;顾福飞等;《电子与信息学报》;20111031;第2452-2457页 * |
基于时分 MIMO 的地基雷达高分辨率成像研究;蒋留兵等;《电子与信息学报》;20160531;第1055-1063页 * |
基于混合匹配追踪算法的 MIMO 雷达稀疏成像方法;王伟等;《电子与信息学报》;20161031;第2415-2422页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107957574A (zh) | 2018-04-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107957574B (zh) | 基于ifft和混合匹配追踪的时分地基mimo滑坡雷达成像方法 | |
CN104714231B (zh) | 一种基于完全互补序列和相位补偿的mimo sar成像方法 | |
CN105699945B (zh) | 频控阵mimo雷达***中的波形优化设计方法 | |
CN104020469A (zh) | 一种mimo雷达距离-角度二维超分辨率成像算法 | |
CN107942295B (zh) | 一种前视阵列sar***的稀疏天线 | |
CN103336278A (zh) | 多视角观测下前视三维sar成像方法 | |
CN104133215B (zh) | 基于距离徙动微调和子带分割的同步轨道雷达成像方法 | |
CN103399315A (zh) | 实孔径相控阵雷达高分辨探测成像方法 | |
CN104898119A (zh) | 一种基于相关函数的动目标参数估计方法 | |
RU2546330C1 (ru) | Способ поляризационно-чувствительного радиоконтроля подвижных объектов | |
CN106199602B (zh) | 一种Staggered-SAR回波信号重建方法 | |
CN108693526B (zh) | 用于mimo雷达的运动目标相位补偿方法 | |
CN110596707B (zh) | 基于多快拍图像联合的mimo雷达三维成像方法 | |
RU2524401C1 (ru) | Способ обнаружения и пространственной локализации подвижных объектов | |
CN114545401A (zh) | 用于机载前视阵列雷达的空时自适应迭代超分辨成像方法 | |
CN104765033B (zh) | 利用互相关函数抑制步进频率成像中距离旁瓣的方法 | |
CN104280732A (zh) | 一种基于等效协同阵的穿墙雷达建筑布局成像方法 | |
RU2529483C1 (ru) | Способ скрытной радиолокации подвижных объектов | |
CN103969641A (zh) | 一种多波束发射三维成像方法 | |
RU2546329C1 (ru) | Способ поляризационно-чувствительного обнаружения подвижных объектов | |
CN113671485A (zh) | 基于admm的米波面阵雷达二维doa估计方法 | |
RU2524399C1 (ru) | Способ обнаружения малоразмерных подвижных объектов | |
CN111580040A (zh) | 双基地展开互质阵列mimo雷达dod和doa降维估计方法 | |
CN108710116B (zh) | 一种mimo雷达的运动目标相位恢复方法 | |
CN112068129B (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 |