CN107064896B - 基于截断修正sl0算法的mimo雷达参数估计方法 - Google Patents
基于截断修正sl0算法的mimo雷达参数估计方法 Download PDFInfo
- Publication number
- CN107064896B CN107064896B CN201710201492.4A CN201710201492A CN107064896B CN 107064896 B CN107064896 B CN 107064896B CN 201710201492 A CN201710201492 A CN 201710201492A CN 107064896 B CN107064896 B CN 107064896B
- Authority
- CN
- China
- Prior art keywords
- matrix
- mimo radar
- singular
- target
- algorithm
- 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
- 238000000034 method Methods 0.000 title claims abstract description 77
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 59
- 238000012937 correction Methods 0.000 title claims abstract description 23
- 239000011159 matrix material Substances 0.000 claims abstract description 109
- 239000013598 vector Substances 0.000 claims abstract description 29
- 238000001514 detection method Methods 0.000 claims abstract description 14
- 230000009466 transformation Effects 0.000 claims abstract description 5
- 230000008447 perception Effects 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 7
- 238000000354 decomposition reaction Methods 0.000 claims description 6
- 230000001575 pathological effect Effects 0.000 claims description 5
- 238000002945 steepest descent method Methods 0.000 claims description 4
- 239000000654 additive Substances 0.000 claims description 3
- 230000000996 additive effect Effects 0.000 claims description 3
- 230000014759 maintenance of location Effects 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 3
- 238000013461 design Methods 0.000 claims description 2
- 238000003491 array Methods 0.000 claims 1
- 230000000717 retained effect Effects 0.000 claims 1
- 238000004364 calculation method Methods 0.000 description 5
- 238000004088 simulation Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 230000001629 suppression Effects 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
- 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
本发明公开了一种基于截断修正SL0算法的MIMO雷达参数估计方法,首先改善MIMO雷达感知矩阵的病态性,利用SVD反变换从修正后的奇异值及其对应的左右奇异矩阵中获得非病态感知矩阵;其次利用SL0算法对MIMO雷达目标参数进行估计,并在计算初值和梯度投影值时,将病态感知矩阵的伪逆由获得的非病态感知矩阵的伪逆代替;最后根据SL0算法得到的目标场景向量估计值中非零元素的位置确定MIMO雷达目标的角度、距离和多普勒信息。本发明解决了病态感知矩阵下MIMO雷达目标参数估计问题。降低了MIMO雷达目标参数估计的复杂度,有利于工程实现,适合于在现代战争中要求实时性较高的目标探测场合。
Description
技术领域
本发明涉及一种基于截断修正SL0算法的MIMO雷达目标参数估计方法,属于MIMO雷达目标参数估计技术领域。
背景技术
多输入多输出雷达(Multiple Input and Multiple Output,MIMO)是一种新体制雷达系统。与相控阵雷达相比,MIMO雷达采用波形分集技术,提高了目标分辨率,增强了系统参数识别能力,在参数估计、噪声抑制和目标探测方面具有很大优势。压缩感知(Compressed Sensing,CS)是一种新兴的信号采样和重建理论,不同于传统的奈奎斯特(Nyquist)采样定理,它通过随机采样的少量观测值就能实现稀疏信号重构,是目前信号处理领域的研究热点。在实际雷达探测区域,目标呈稀疏分布,其回波信号呈稀疏性,故可以将CS理论应用于雷达目标探测。稀疏重构问题等价于l0范数最小化问题,此问题是NP-hard,其求解难度随着维度增加而增大,因此利用l0范数最小化难以对高维度稀疏重构问题进行求解。现有技术中,提出一种利用正则化迭代重加权最小化方法(RegularizedIterative reweighted Minimization Approach,RIRMA)实现MIMO雷达目标参数估计,该方法通过一系列迭代加权lq(0<q≤1)范数来逼近l0范数,只需少量观测值即可精确重构MIMO雷达的目标信号,但RIRMA方法在每次迭代时需对更新后的大维度矩阵进行求逆,导致该方法计算速度较慢。SL0(Smoothed l0norm)算法是一种匹配度高和计算效率高的稀疏重构算法,它由一系列高斯函数来趋近l0范数,从而将l0范数最小化问题转化成平滑函数最小化问题,并利用最速下降法和梯度投影法求解该最小化问题,从而快速实现稀疏信号的重构。还有将SL0算法应用于MIMO雷达的目标参数估计,采用收敛性更好的双曲正切函数来近似l0范数,然后利用修正牛顿法求解该极值问题,并采用正则化方法改善SL0算法中的病态问题,提高了MIMO雷达目标参数估计的实时性。
为保证稀疏重构算法对MIMO雷达目标参数估计的精度,目标场景一般被划分成精细的栅格,则MIMO雷达的感知矩阵不可避免地存在近似线性相关的列,从而导致该矩阵呈病态。将SL0算法应用于MIMO雷达能明显提高其目标参数的估计速度,但是SL0算法中初值和梯度投影步骤都需要对病态感知矩阵求伪逆,则当MIMO雷达的接收信号中存在微小的噪声扰动时,就会引起初值和梯度投影计算误差较大,进而导致SL0算法失效。现有技术中,将截断奇异值分解(Truncated Singular Value Decomposition,TSVD)技术引入SL0算法,对感知矩阵进行奇异值分解,通过剔除对噪声干扰较为敏感的较小奇异值及其对应的左右奇异矩阵,改善感知矩阵的病态性,从而提高了SL0算法的稳健性。但是该方法仅仅剔除感知矩阵中较小的奇异值而未对较大奇异值进行修正,导致感知矩阵病态性的改善效果有限。因此,为了能将SL0算法应用于MIMO雷达目标参数估计问题中,对MIMO雷达感知矩阵的奇异值进行截断修正处理以改善其病态性是非常有必要的。
发明内容
目的:为了克服SL0算法在MIMO雷达目标参数估计中因感知矩阵的病态性而导致其失效的问题,本发明提供基于截断修正SL0算法的MIMO雷达参数估计方法。
技术方案:为解决上述技术问题,本发明采用的技术方案为:
一种基于截断修正SL0算法的MIMO雷达参数估计方法,包括如下步骤:
步骤一:建立MIMO雷达接收信号向量模型;
步骤二:通过对MIMO雷达感知矩阵的奇异值进行二次截断和一次修正处理来改善感知矩阵的病态性,并利用SVD反变换获得非病态感知矩阵;
步骤三:利用SL0算法对MIMO雷达目标参数进行估计,并在计算初值和梯度投影值的步骤中,将病态感知矩阵的伪逆由上面获得的非病态感知矩阵的伪逆代替;
步骤四:最终输出求解后的目标场景向量α,根据最终输出求解后的目标场景向量α中的非零元素位置来确定各个目标的角度、距离和多普勒频移值。
作为优选方案,所述步骤一包括如下步骤:
1a:建立MIMO雷达的接收信号为式中,(·)*表示复共轭转置,E是加性噪声,αp,k,h(p=1,...,P,k=1,...,K,h=1,...,H)表示目标的复散射系数,JP是用来描述从不同距离单元返回的信号所用的转移矩阵;
1b:将上式简化成向量形式,定义
式中,vec(·)表示矩阵向量化运算,A是感知矩阵,α是目标场景向量;
1c:得到接收信号的向量模型为y=Aα+e式中,e=vec(E)。
作为优选方案,所述步骤二包括如下步骤:
2a:对A作奇异值分解:A=UΣV*,其中为由奇异值σi(i=1,2,…,Mr(N+P-1))构成的对角矩阵,U、V分别为奇异值对应的左奇异矩阵、右奇异矩阵。
2b:选取一次截断门限保留大于等于的奇异值,剔除小于的奇异值,即
2c:对保留的大于等于的奇异值再次截断并修正,选择二次截断门限将保留的奇异值分成大于等于和小于的两部分,利用Tikhonov正则化准则对大于等于的奇异值进行修正,将小于的奇异值修正为与接近的常数值,即q1和q2为修正参数;
2d:定义U1和V1分别为修正后的奇异值所对应的左奇异矩阵和右奇异矩阵,其中U1=U(:,1:t1)由矩阵U的前t1列组成的矩阵,V1=V(1:t1,:)由矩阵V的前t1行组成的矩阵,Σ1由经修正后的奇异值σ'm(1≤m≤t1)构成的对角矩阵;利用SVD反变换获得非病态感知矩阵A1=U1Σ1V1 *,矩阵A1的伪逆
作为优选方案,所述步骤三包括如下步骤:
3a:初始化:
(1)令初值内循环次数为L、步长μ,其中L,μ>0;
(2)选取一个序列[δ1,δ2,...,δJ],且δ1>δ2>...>δJ,其中,3b:算法迭代:
for j=1,2,...,J
(1)令δ=δj
(2)利用最速下降法求解Fδ(α)的最小值,并将其投影到可行集上;
初始化:
for l=1,...,L
(a)设φ=[α1exp(-α1 2/2δ2),...,αzexp(-αZ 2/2δ2)]T,其中Z=P·K·H,P、K和H分别是在目标探测场景中所划分的距离单元个数,角度单元个数,以及多普勒单元个数;
(3)令
3c:当或δj≤δJ时,转至步骤四,否则,令j=j+1,δj=ρδj-1,返回步骤二;其中η为一个正数,用于判断算法是否终止,ρ为衰减因子,0<ρ<1。
有益效果:本发明提供的基于截断修正SL0算法的MIMO雷达参数估计方法,该方法在截断奇异值分解(Truncated Singular Value Decomposition,TSVD)法基础上,通过设置截断门限将保留的感知矩阵奇异值分成较大和较小两部分,并分别采用不同的修正准则对它们进行修正,然后利用SVD反变换从修正后的奇异值及其对应的左右奇异矩阵中获得非病态感知矩阵,在SL0算法中以截断修正后的非病态感知矩阵来重构MIMO雷达的目标信号,从而显著提高了MIMO雷达目标参数估计的精度和速度。
与现有技术相比,本发明具有如下优点:
1、在利用SL0算法对MIMO雷达目标参数进行估计时,由于MIMO雷达的感知矩阵呈病态性,则雷达接收信号中微小的噪声干扰会导致SL0算法中初值和梯度投影计算误差较大,进而使得SL0算法失效;本发明方法通过对MIMO雷达感知矩阵的奇异值进行两次截断和一次修正处理,明显改善了感知矩阵的病态性,从而提高了SL0算法中初值和梯度投影的计算精度,使得SL0算法能够以较高的精度重构出MIMO雷达的稀疏目标信号,解决了病态感知矩阵下MIMO雷达目标参数估计问题。
2、本发明方法能利用SL0算法对MIMO雷达的目标参数进行快速估计,降低了MIMO雷达目标参数估计的复杂度,有利于工程实现,适合于在现代战争中要求实时性较高的目标探测场合。
附图说明
图1为不同方法在SNR=0dB时距离-角度二维估计图;
图2为不同方法在SNR=0dB时距离-多普勒二维估计图;
图3为不同方法的重构信噪比和回波信噪比的变化关系图;
图4为不同方法的运行时间和回波信噪比的变化关系图。
具体实施方式
下面结合附图对本发明作更进一步的说明。
一种基于截断修正SL0算法的MIMO雷达参数估计方法,包括如下步骤:
一、MIMO雷达接收信号模型
假设MIMO雷达的发射阵列和接收阵列分别由Mt个发射阵元和Mr个接收阵元组成,其中,发射阵元和接收阵元间隔分别为dt和dr。发射阵列的发射信号矩阵表示为
式中,sm=[sm(1),sm(2),...,sm(N)]T表示第m个发射阵元的发射信号,N为发射信号长度。
将雷达目标探测场景划分为Z(Z=P·K·H)个离散的距离-角度-多普勒单元,其中,P是距离单元个数,K是角度单元个数,H是多普勒单元个数。τp(1≤p≤P)表示第p个距离单元的延迟,且目标回波间最大延时单元(系统第一个距离单元的发射信号与最后一个距离单元反射信号的时延)为P-1;θk(1≤k≤K)表示第k个角度单元对应的目标角度;ωh(1≤h≤H)为第h个多普勒单元对应的目标多普勒频移。则对应ωh的多普勒频移信号矩阵为
式中,
和d(ωh)分别表示对应ωh的多普勒频移信号向量和多普勒频移向量。
令
式中,是补零发射信号矩阵,0(P-1)×Mt为维度(P-1)×Mt的零矩阵。发射阵列和接收阵列的导向向量分别表示为
式中,λ0表示载波波长。因此,MIMO雷达的接收信号为
式中,(·)*表示复共轭转置,E是加性噪声,αp,k,h(p=1,...,P,k=1,...,K,h=1,...,H)表示目标的复散射系数。JP是用来描述从不同距离单元返回的信号所用的转移矩阵,其表达式如下
将式(8)简化成向量形式,定义
式中,vec(·)表示矩阵向量化运算。A是感知矩阵,α是目标场景向量。则接收信号的向量模型为
y=Aα+e (12)
式中,e=vec(E)。由于α中非零元素只占少数,式(12)为稀疏表示的MIMO雷达接收信号模型,可以利用稀疏重构算法估计α,根据α中非零元素的位置从而获得目标的参数信息(距离、角度、多普勒)。
感知矩阵A存在近似线性相关的列会导致矩阵病态。A中列向量bp,k,h等价于划分栅格{(τp,θk,ωh)}上单位幅度假想目标的回波信号。为了使邻近的角度(距离或多普勒)单元所对应的感知矩阵A中列向量线性无关,则其角度(距离或多普勒)单元划分间隔需大于或等于角度(距离或多普勒)分辨率。以角度分辨率为例,MIMO雷达接收阵列在波束指向角为θk时的角分辨率为定义Δθ=θk+1-θk为角度划分间隔大小,则当满足Δθ≥θ0.5s时,可认为A中相邻角度θk和θk+1所对应的列向量近似无关。假设MIMO雷达接收阵元数Mr=8,阵元间隔dr=0.5λ0,θk=0°,则此时的角分辨率为θ0.5s=12.7°,即当角度划分间隔|Δθ|≥12.7°时,A中相邻角度θk和θk+1所对应的列向量近似无关。而在实际雷达探测场景中,为保证稀疏重构算法的精确性,其角度单元划分间隔一般远小于角分辨率,则感知矩阵A中不可避免地会存在线性相关的列。同样,当距离和多普勒单元精细划分时也会导致MIMO雷达感知矩阵A中存在近似线性相关的列,且A的条件数非常大,从而呈严重的病态性。
二、通过对MIMO雷达感知矩阵的奇异值进行二次截断和一次修正处理来改善感知矩阵的病态性,并利用SVD反变换获得非病态感知矩阵。
具体步骤如下:
(1)对A作奇异值分解:A=UΣV*,其中为由奇异值σi(i=1,2,…,Mr(N+P-1))构成的对角矩阵,U、V分别为奇异值对应的左奇异矩阵、右奇异矩阵。
(2)选取一次截断门限保留大于等于的奇异值,剔除小于的奇异值,即
(3)对保留的大于等于的奇异值再次截断并修正:选择二次截断门限将保留的奇异值分成大于等于和小于的两部分,利用Tikhonov正则化准则对大于等于的奇异值进行修正,将小于的奇异值修正为与接近的常数值,即q1和q2为修正参数。
(4)定义U1和V1分别为修正后的奇异值所对应的左奇异矩阵和右奇异矩阵,其中U1=U(:,1:t1)由矩阵U的前t1列组成的矩阵,V1=V(1:t1,:)
由矩阵V的前t1行组成的矩阵,Σ1由经修正后的奇异值σ'm(1≤m≤t1)构成的对角矩阵。利用SVD反变换获得非病态感知矩阵A1=U1Σ1V1 *,矩阵A1的伪逆为
三、利用SL0算法对MIMO雷达目标参数进行估计,并在计算初值和梯度投影值的步骤中,将病态感知矩阵的伪逆由上面获得的非病态感知矩阵的伪逆代替。
具体步骤如下:
(1)初始化:
(a)令初值内循环次数为L、步长μ,其中L,μ>0。
(b)选取一个序列[δ1,δ2,...,δJ],且δ1>δ2>...>δJ。其中,
(2)算法迭代:
for j=1,2,...,J
①令δ=δj
②利用最速下降法求解Fδ(α)的最小值,并将其投影到可行集上
初始化:
for l=1,...,L
(a)设φ=[α1exp(-α1 2/2δ2),...,αzexp(-αZ 2/2δ2)]T,其中Z=P·K·H,P、K和H分别是在目标探测场景中所划分的距离单元个数,角度单元个数,以及多普勒单元个数。
③令
(3)当或δj≤δJ时,转至步骤(4),否则,令j=j+1,δj=ρδj-1,返回步骤(2)。其中η为一个正数,用于判断算法是否终止,ρ为衰减因子,0<ρ<1。
(4)最终输出求解后的目标场景向量α。根据最终输出求解后的目标场景向量α中的非零元素位置来确定各个目标的角度、距离和多普勒频移值。
本发明的技术效果可以通过以下仿真结果进一步说明。设MIMO雷达的发射阵元数和接收阵元数分别为Mt=8和Mr=8,发射阵元间隔和接收阵元间隔分别为dt=2.5λ0、dr=0.5λ0;采用相互正交的噪声调频信号作为发射阵列的发射信号,发射波形采样个数N=32;回波噪声采用均值为零、方差为σn 2的加性高斯白噪声,定义回波信噪比为
SNR=10log10{tr(S*S)/Nσn 2} (13)
式中,tr(·)表示矩阵求迹运算。假设雷达目标探测场景有P=12个距离单元,雷达扫描角度范围[-30°,30°],以1°为间隔将其划分成61个角度单元,即K=61;目标的多普勒频移用角度表示,即Φh=ωhN(180°/π),目标多普勒范围[-25°,25°],以5°为间隔将其划分成11个多普勒单元,即H=11。
在迭代加权lq算法中,设置迭代次数l=8,迭代范数q=0.9,ε=0.01,η=0.01。SL0_TSVD算法表示采用TSVD方法解决SL0算法中的病态问题;SL0_Tikhonov正则化法表示采用Tikhonov正则化方法解决SL0算法中的病态问题。为了表述方便,本发明算法简称为TMSL0算法。在SL0算法、SL0_TSVD算法、SL0_Tikhonov正则化法及TMSL0算法中,设置σJ=0.01,ρ=0.8,η=0.01,内循环次数L=50,步长μ=2。在TMSL0算法中,设置C=25,满足条件0<C<100,以确定截断门限修正参数q1=115,q2=95。在SL0_TSVD方法和SL0_Tikhonov方法中截断参数和正则化参数的选取都是采用经典的L曲线法。
定义重构信噪比如下,
式中,α为真实目标场景向量;目标场景向量估计值;||·||2表示2-范数。
仿真内容1:MIMO雷达在回波信噪比SNR=0dB时的距离-角度-多普勒估计图
图1为MIMO雷达目标的距离-角度估计图,回波信噪比为0dB。图1(a)为真实目标在多普勒单元5°处的距离-角度单元估计图,图1(b)、图1(c)、图1(d)、图1(e)和图1(f)分别是SL0算法、迭代加权lq方法、SL0_Tikhonov方法、SL0_TSVD方法和TMSL0方法估计得到的距离-角度估计图。图2(a)为真实目标在角度单元-10°处的距离-多普勒单元估计图,回波信噪比为0dB。图2(b)、图2(c)、图2(d)、图2(e)和图2(f)分别是由上述五种方法获得的的距离-多普勒单元估计图。由图1和图2可知,当感知矩阵病态时,SL0算法失效,导致其无法有效重构MIMO雷达的目标信号;而迭代加权lq方法、SL0_Tikhonov方法、SL0_TSVD方法虽都对病态感知矩阵作了不同的修正处理,在一定程度上改善了MIMO雷达感知矩阵的病态性,但是这些方法在参数估计时存在旁瓣电平较高,并在真实目标附近存在许多虚假目标的问题;TMSL0方法对感知矩阵奇异值进行二次截断和一次修正处理,进一步改善了MIMO雷达感知矩阵的病态性,从而能够以较高的精度重构目标信号,该方法的目标参数估计性能最好。
仿真内容2:五种方法的重构信噪比与回波信噪比的变化关系
设回波信噪比在-10dB~20dB之间变化,重复100次单独实验。图3为五种算法重构性能与回波信噪比的变化关系。由图3可知,由于SL0_Tikhonov方法和SL0_TSVD方法都不同程度地改善了感知矩阵的病态性,它们的重构信噪比都明显高于SL0方法,而本发明提出的TMSL0方法重构性能优于SL0_Tikhonov方法和SL0_TSVD方法,相比之下,TMSL0方法更适合应用于MIMO雷达目标参数估计。
仿真内容3:四种方法运行时间与回波信噪比的变化关系
设回波信噪比在-10dB~20dB之间变化,重复100次单独实验。图4表示四种方法运行时间与回波信噪比的变化关系。迭代加权lq方法的每次迭代都需要对更新后的大维度矩阵进行求逆运算,导致利用该方法对MIMO雷达目标参数估计耗时较长。SL0_Tikhonov方法、SL0_TSVD方法及TMSL0方法的运行时间明显低于迭代加权lq方法,在较低的信噪比情况下,TSVD方法和Tikhonov方法由于采用L曲线法求解截断参数和正则化参数不稳定,对感知矩阵的病态性改善效果不理想,导致SL0算法中初值求解的精度较低,则需要更多的迭代次数才能达到退出条件,因此这两种方法运行时间较长,而当信噪比增加后,这三种方法的运行时间接近。
以上所述仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (5)
1.一种基于截断修正SL0算法的MIMO雷达参数估计方法,其特征在于:包括如下步骤:
步骤一:建立MIMO雷达接收信号向量模型;
步骤二:通过对MIMO雷达感知矩阵的奇异值进行二次截断和一次修正处理来改善感知矩阵的病态性,并利用SVD反变换获得非病态感知矩阵;
步骤三:利用SL0算法对MIMO雷达目标参数进行估计,并在计算初值和梯度投影值的步骤中,将病态感知矩阵的伪逆由上面获得的非病态感知矩阵的伪逆代替;
步骤四:最终输出求解后的目标场景向量α,根据最终输出求解后的目标场景向量α中的非零元素位置来确定各个目标的角度、距离和多普勒频移值;
所述步骤一包括如下步骤:
1a:建立MIMO雷达的接收信号为式中,(·)*和(·)T分别表示矩阵或矢量的复共轭转置和转置;P是雷达目标探测区域内所划分的距离单元个数,K是雷达目标探测区域内所划分的角度单元个数,H是雷达目标探测区域内所划分的多普勒单元个数;E是加性噪声;αp,k,h,p=1,...,P;k=1,...,K;h=1,...,H,表示目标的复散射系数,其中;JP是用来描述从不同距离单元返回的信号所用的转移矩阵;aT(θk)和aR(θk)分别是发射阵列和接收阵列的导向向量,其中,θk表示第k,1≤k≤K,个角度单元对应的目标角度;是补零发射信号矩阵,其中,ωh为第h个多普勒单元对应的目标多普勒频移,h=1,...,H,Sd(ωh)为对应ωh的多普勒频移信号矩阵,为维度(P-1)×Mt的零矩阵,Mt为发射阵元个数,N为发射信号长度;
1b:将上式简化成向量形式,定义
式中,vec(·)表示矩阵向量化运算,A是感知矩阵;α是目标场景向量;N为发射信号长度;Mr为接收阵元的数量;
1c:得到接收信号的向量模型为y=Aα+e,式中,e=vec(E);
所述步骤二包括如下步骤:
2a:对矩阵A作奇异值分解:A=UΣV*,其中为由奇异值σi,i=1,2,…,Mr(N+P-1),构成的对角矩阵,U、V分别为奇异值对应的左奇异矩阵、右奇异矩阵,A是感知矩阵;
2b:选取一次截断门限保留大于等于的奇异值,剔除小于的奇异值,即
2c:对保留的大于等于的奇异值再次截断并修正,选择二次截断门限将保留的奇异值分成大于等于和小于的两部分,利用Tikhonov正则化准则对大于等于的奇异值进行修正,将小于的奇异值修正为与接近的常数值,即q1和q2为修正参数;
2d:定义U1和V1分别为修正后的奇异值所对应的左奇异矩阵和右奇异矩阵,其中U1=U(:,1:t1)由矩阵U的前t1列组成的矩阵,V1=V(1:t1,:)由矩阵V的前t1行组成的矩阵,Σ1是修正后的奇异值σ'm,1≤m≤t1,构成的对角矩阵;利用SVD反变换获得非病态感知矩阵A1=U1Σ1V1 *,矩阵A1的伪逆
2.根据权利要求1所述的基于截断修正SL0算法的MIMO雷达参数估计方法,其特征在于:所述步骤三包括如下步骤:
3a:初始化:
(1)令初值内循环次数为L、步长为μ,其中L,μ>0,U1和V1分别为修正后的奇异值所对应的左奇异矩阵和右奇异矩阵;y为向量形式表示的MIMO雷达接收信号;Σ1是修正后的奇异值σ'm,1≤m≤t1,构成的对角矩阵;
(2)选取一个序列[δ1,δ2,...,δJ],且δ1>δ2>...>δJ,其中,
3b:算法迭代:
forj=1,2,...,J
(1)令δ=δj
(2)利用最速下降法求解Fδ(α)的最小值,并将其投影到可行集上;
初始化:
forl=1,...,L
(a)设φ=[α1exp(-α1 2/2δ2),...,αZexp(-αZ 2/2δ2)]T,其中Z=P·K·H,P、K和H分别是在雷达目标探测区域中所划分的距离单元个数,角度单元个数,以及多普勒单元个数;
(b)
(c)
(3)令
3c:当或δj≤δJ时,转至步骤四,否则,令j=j+1,δj=ρδj-1,返回步骤3b;其中η为一个正数,用于判断算法是否终止,ρ为衰减因子,0<ρ<1。
3.根据权利要求2所述的基于截断修正SL0算法的MIMO雷达参数估计方法,其特征在于:所述η=0.01,ρ=0.8。
4.根据权利要求2所述的基于截断修正SL0算法的MIMO雷达参数估计方法,其特征在于:所述内循环次数L=50,步长μ=2。
5.根据权利要求2所述的基于截断修正SL0算法的MIMO雷达参数估计方法,其特征在于:所述修正参数q1=115,q2=95,C=25。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710201492.4A CN107064896B (zh) | 2017-03-30 | 2017-03-30 | 基于截断修正sl0算法的mimo雷达参数估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710201492.4A CN107064896B (zh) | 2017-03-30 | 2017-03-30 | 基于截断修正sl0算法的mimo雷达参数估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107064896A CN107064896A (zh) | 2017-08-18 |
CN107064896B true CN107064896B (zh) | 2019-12-10 |
Family
ID=59601409
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710201492.4A Active CN107064896B (zh) | 2017-03-30 | 2017-03-30 | 基于截断修正sl0算法的mimo雷达参数估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107064896B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108197414A (zh) * | 2018-01-30 | 2018-06-22 | 江苏建筑职业技术学院 | 一种改进的结构方程模型参数估计算法 |
CN110780274B (zh) * | 2019-11-04 | 2022-04-01 | 电子科技大学 | 一种用于扫描雷达的改进l1正则化方位超分辨成像方法 |
CN112816971B (zh) * | 2020-12-31 | 2023-04-07 | 湖州师范学院 | 分离式mimo雷达波系统中目标位置的估计方法和装置 |
CN113376630A (zh) * | 2021-07-05 | 2021-09-10 | 中国人民解放军空军工程大学 | 雷达成像的方法、装置、电子设备和计算机可读存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101719118A (zh) * | 2009-11-24 | 2010-06-02 | 中国科学院对地观测与数字地球科学中心 | 一种用于克服有理函数模型病态性的改进的奇异值修正方法 |
CN103337087A (zh) * | 2013-07-04 | 2013-10-02 | 西北工业大学 | 一种基于伪逆自适应算法的压缩感知图像重构方法 |
CN105930310A (zh) * | 2016-05-06 | 2016-09-07 | 南京信息工程大学 | 一种自适应正则化平滑l0范数方法 |
CN106125041A (zh) * | 2016-07-26 | 2016-11-16 | 清华大学 | 基于子空间加权稀疏恢复的宽带信号源定位方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8912951B2 (en) * | 2012-10-09 | 2014-12-16 | Raytheon Company | Moving target detection using a two-dimensional folding approach |
CA2891839C (en) * | 2014-05-16 | 2023-02-14 | Mohamed Mabrouk | Remote sensing of human breathing at a distance |
-
2017
- 2017-03-30 CN CN201710201492.4A patent/CN107064896B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101719118A (zh) * | 2009-11-24 | 2010-06-02 | 中国科学院对地观测与数字地球科学中心 | 一种用于克服有理函数模型病态性的改进的奇异值修正方法 |
CN103337087A (zh) * | 2013-07-04 | 2013-10-02 | 西北工业大学 | 一种基于伪逆自适应算法的压缩感知图像重构方法 |
CN105930310A (zh) * | 2016-05-06 | 2016-09-07 | 南京信息工程大学 | 一种自适应正则化平滑l0范数方法 |
CN106125041A (zh) * | 2016-07-26 | 2016-11-16 | 清华大学 | 基于子空间加权稀疏恢复的宽带信号源定位方法 |
Non-Patent Citations (4)
Title |
---|
sparse learning via iterative minimization with application to MIMO radar imaging;Xing Tan等;《IEEE transaction on signal processing》;20110331;第59卷(第3期);1088-1101 * |
基于SL0算法的快速局部稀疏多径信道估计;刘婷等;《计算机工程与设计》;20140331;第35卷(第3期);785-790 * |
基于SL0范数的改进稀疏信号重构算法;冯俊杰等;《数据采集与处理》;20160131;178-183 * |
奇异值分解法在病态问题中的应用;徐文等;《测绘通报》;20160131(第1期);62-63,83 * |
Also Published As
Publication number | Publication date |
---|---|
CN107064896A (zh) | 2017-08-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110275166B (zh) | 基于admm的快速稀疏孔径isar自聚焦与成像方法 | |
CN110113085B (zh) | 一种基于协方差矩阵重构的波束形成方法及系统 | |
CN110208735B (zh) | 一种基于稀疏贝叶斯学习的相干信号doa估计方法 | |
CN107064896B (zh) | 基于截断修正sl0算法的mimo雷达参数估计方法 | |
CN105652273B (zh) | 一种基于混合匹配追踪算法的mimo雷达稀疏成像算法 | |
CN105929386B (zh) | 一种基于高阶累积量的波达估计方法 | |
CN107167783B (zh) | 一种共形阵列杂波协方差矩阵的稀疏重构方法 | |
CN106646414B (zh) | 基于Bi-CGSTAB和SL0算法的MIMO雷达目标参数估计方法 | |
CN107870315B (zh) | 一种利用迭代相位补偿技术估计任意阵列波达方向方法 | |
CN106226754A (zh) | 基于时间反演的低仰角波达方向估计方法 | |
CN109298383A (zh) | 一种基于变分贝叶斯推断的互质阵波达方向角估计方法 | |
CN107102298B (zh) | 基于迭代互耦校正的雷达协方差矩阵重构波束形成方法 | |
CN113189592B (zh) | 考虑幅相互耦误差的车载毫米波mimo雷达测角方法 | |
CN106872934B (zh) | L型电磁矢量传感器阵列解相干esprit参数估计方法 | |
CN110161489A (zh) | 一种基于伪框架的强弱信号测向方法 | |
CN108303683A (zh) | 单基地mimo雷达实值esprit非圆信号角度估计方法 | |
CN111337873A (zh) | 一种基于稀疏阵的doa估计方法 | |
CN109557504B (zh) | 一种近场窄带信号源的定位方法 | |
CN110196417B (zh) | 基于发射能量集中的双基地mimo雷达角度估计方法 | |
CN106908754B (zh) | L型声矢量传感器阵列esprit解相干参数估计方法 | |
Bingbing | DOA estimation of the coherent signals using beamspace matrix reconstruction | |
CN115236584A (zh) | 基于深度学习的米波雷达低仰角估计方法 | |
CN108828586B (zh) | 一种基于波束域的双基地mimo雷达测角优化方法 | |
CN107167782B (zh) | 基于信杂噪比最大的雷达三维异构阵稀疏重构方法 | |
CN117092585A (zh) | 单比特量化DoA估计方法、系统和智能终端 |
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 |