CN108445465B - 一种基于Arnoldi的特征向量相位补偿方法 - Google Patents

一种基于Arnoldi的特征向量相位补偿方法 Download PDF

Info

Publication number
CN108445465B
CN108445465B CN201810225591.0A CN201810225591A CN108445465B CN 108445465 B CN108445465 B CN 108445465B CN 201810225591 A CN201810225591 A CN 201810225591A CN 108445465 B CN108445465 B CN 108445465B
Authority
CN
China
Prior art keywords
arnoldi
vector
phase compensation
matrix
equal
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
Application number
CN201810225591.0A
Other languages
English (en)
Other versions
CN108445465A (zh
Inventor
刘秋晨
王勇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN201810225591.0A priority Critical patent/CN108445465B/zh
Publication of CN108445465A publication Critical patent/CN108445465A/zh
Application granted granted Critical
Publication of CN108445465B publication Critical patent/CN108445465B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41Details 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

一种基于Arnoldi的特征向量相位补偿方法,本发明涉及基于Arnoldi的特征向量相位补偿方法。本发明的目的是为了解决现有基于特征向量的相位补偿方法计算量大的问题。一、取M个目标回波,对复距离像取模得到目标回波的一维距离像,并利用积累互相关方法进行包络对齐,得到包络对齐后的距离维压缩数据;二、求得采样协方差矩阵;三、利用Arnoldi迭代求解采样协方差矩阵的最大特征值对应的特征向量,提取最大特征值对应的特征向量中各元素的相位信息,对一中的包络对齐后的距离维压缩数据进行相位补偿,得到相位补偿后的数据;四、利用距离多普勒方法对相位补偿后的数据进行ISAR成像。本发明用于逆合成孔径雷达领域。

Description

一种基于Arnoldi的特征向量相位补偿方法
技术领域
本发明属于逆合成孔径雷达(ISAR)领域,涉及基于Arnoldi的特征向量相位补偿方法。
背景技术
逆合成孔径雷达(ISAR)利用雷达和目标的相对运动来实现对空间目标的高分辨成像,如今已广泛运用于军事和民用领域。运动补偿是ISAR成像的一个重要步骤,其又分为包络对齐和相位补偿两部分。包络对齐为相位补偿的基础。目前利用相邻一维像的相似性进行包络对齐的方法应用最为普遍,因此以积累互相关的方法来进行一维距离像的包络对齐。基于特征向量的相位补偿方法作为一种全局相位补偿方法既可应用于全孔径又可运用于稀疏孔径,且拥有很好的相位误差估计性能,但是该方法需要对数据的采样协方差矩阵进行特征分解来进一步获取相位误差信息。对于一个高分辨的ISAR图像,用于成像的回波个数较多,因此所对应的采样协方差矩阵维度较高,对其进行特征分解的计算量较大,从而导致其在实际工程中的应用受到限制。在基于特征向量的相位补偿方法中包含相位误差信息的仅为对应于采样协方差矩阵最大特征值的特征向量,若能仅对该向量进行求解便能大大降低计算量。因此,基于Arnoldi的特征向量相位补偿方法的研究尤为重要。
发明内容
本发明的目的是为了解决现有基于特征向量的相位补偿方法计算量大的问题,从而提出了一种基于Arnoldi的特征向量相位补偿方法。
一种基于Arnoldi的特征向量相位补偿方法具体过程为:
步骤一、取M个目标回波,经过距离维压缩后的第m次目标回波沿距离n分布的复距离像为s(m,n),对复距离像取模得到目标回波的一维距离像,并利用积累互相关方法进行包络对齐,得到包络对齐后的距离维压缩数据;
其中,0≤m≤(M-1),0≤n≤(N-1),M为方位向点数(即目标回波个数),N为距离向点数;
步骤二、根据包络对齐后的距离维压缩数据求得采样协方差矩阵;
步骤三、利用Arnoldi迭代求解采样协方差矩阵的最大特征值对应的特征向量,提取最大特征值对应的特征向量中各元素的相位信息,利用各元素的相位信息对步骤一中的包络对齐后的距离维压缩数据进行相位补偿,得到相位补偿后的数据;
步骤四、利用距离多普勒(RD)方法对相位补偿后的数据进行ISAR成像;
所述ISAR为逆合成孔径雷达。
本发明的有益效果为:
对于一个高分辨的ISAR图像,用于成像的回波个数较多,所对应的采样协方差矩阵的维度较高,现有基于特征向量的相位补偿方法直接对高维度的采样协方差矩阵进行特征分解,求出所有的特征向量后再提取最大特征值所对应的特征向量,计算量较大。本发明利用Arnoldi迭代仅对包含相位误差信息的特征向量进行求解,即利用Arnoldi迭代仅对采样协方差矩阵最大特征值对应的特征向量进行求解,避免了直接对高维度的采样协方差矩阵进行特征分解,从而减少了计算量。本发明相对于现有基于特征向量的相位补偿方法在保证了相位误差估计性能的同时,降低了计算复杂度,减小了计算量,更适用于实时ISAR成像。如图2为经过现有基于特征向量的相位补偿方法进行相位补偿后的飞机ISAR成像结果,图3为经过本发明所提出的基于Arnoldi的特征向量相位补偿方法进行相位补偿后的飞机ISAR成像结果,从图2图3可以看出,本发明所提出的基于Arnoldi的特征向量相位补偿方法和现有基于特征向量的相位补偿方法拥有一致的相位补偿效果。如表1为本发明所提出的基于Arnoldi的特征向量相位补偿方法与现有基于特征向量的相位补偿方法的计算时长对比。表1中数据为100次蒙特卡洛实验的总计算时长,从表1可以看出,本发明所提出的方法的计算时长为5.9635s,远小于现有基于特征向量的相位补偿方法的计算时长243.1921s。因此,本发明所提出的基于Arnoldi的特征向量相位补偿方法更适合用于实时ISAR成像系统的相位补偿。
附图说明
图1为本发明流程图;
图2为现有基于特征向量的相位补偿方法的ISAR成像结果图;
图3为本发明基于Arnoldi的特征向量相位补偿方法的ISAR成像结果图。
具体实施方式
具体实施方式一:结合图1说明本实施方式,本实施方式的一种基于Arnoldi的特征向量相位补偿方法具体过程为:
步骤一、取M个目标回波,经过距离维压缩后的第m次目标回波沿距离n分布的复距离像为s(m,n),对复距离像取模得到目标回波的一维距离像,并利用积累互相关方法进行包络对齐,得到包络对齐后的距离维压缩数据;
其中,0≤m≤(M-1),0≤n≤(N-1),M为方位向点数(即目标回波个数),N为距离向点数;
步骤二、根据包络对齐后的距离维压缩数据求得采样协方差矩阵;
步骤三、利用Arnoldi迭代求解采样协方差矩阵的最大特征值对应的特征向量,提取最大特征值对应的特征向量中各元素的相位信息,利用各元素的相位信息对步骤一中的包络对齐后的距离维压缩数据进行相位补偿,得到相位补偿后的数据;
步骤四、利用距离多普勒(RD)方法对相位补偿后的数据进行ISAR成像;
所述ISAR为逆合成孔径雷达。
Arnoldi为Arnoldi迭代算法。
具体实施方式二:本实施方式与具体实施方式一不同的是:所述步骤一中取M个目标回波,经过距离维压缩后的第m次目标回波沿距离n分布的复距离像为s(m,n),对复距离像取模得到目标回波的一维距离像,并利用积累互相关方法进行包络对齐,得到包络对齐后的距离维压缩数据;具体过程为:
步骤一一、经过距离维压缩后的第m次目标回波沿距离n分布的复距离像表示为s(m,n);
其中,0≤m≤(M-1),0≤n≤(N-1),M为方位向点数(即目标回波个数),N为距离向点数,N、M取值为正整数;
步骤一二、对复距离像取模,得到目标回波的一维距离像:
p(m,n)=abs{s(m,n)} (1)
其中,p(m,n)为目标回波的一维距离像,abs{·}为取绝对值;
步骤一三、利用积累互相关方法将当前一维距离像与前l个一维距离像进行互相关:
Figure BDA0001601242760000031
其中,l为脉冲积累数,R(τ)为当前一维距离像与前l个一维距离像的互相关函数,τ为时延,
Figure BDA0001601242760000032
为前l个一维距离像的和,p*(m+a,n)为p(m+a,n)的共轭,p(m+l,n)为当前所要对齐的一维距离像,a为叠加变量,a=0,1,...,l-1;
对于积累互相关法,脉冲积累数有一个最佳脉冲积累数值范围,即当脉冲积累数取值恰当时补偿精度更高。脉冲积累数的取值范围和成像期间接收到的回波数相关,因此,在发明的实施例中,因用于成像的脉冲回波数取2048,为使成像效果较好,积累互相关法的脉冲积累数取100。
对τ进行搜索,(将不同的τ带入式(2))计算R(τ)的最大值(峰值)所对应的时延τ,将R(τ)的最大值(峰值)所对应的时延τ作为补偿,来实现包络的对齐,得到包络对齐后的距离维压缩数据。
其它步骤及参数与具体实施方式一相同。
具体实施方式三:本实施方式与具体实施方式一或二不同的是:所述步骤二中根据包络对齐后的距离维压缩数据求得采样协方差矩阵;具体过程为:
步骤二一、第d个距离单元的距离维压缩数据的复包络sd表示为:
Figure BDA0001601242760000041
其中,ad为第d个距离单元的强散射点的散射率,d=1,2,...,N-1,v为相位误差向量,nd为高斯白噪声,
Figure BDA0001601242760000042
为第e+1个方位向单元与第一个方位向单元之间的相位差,e=1,2,...,M-1,
Figure BDA0001601242760000043
为虚数单位;
s(m,n)为经过距离维压缩后的第m次回波沿距离n分布的复距离像表达式,0≤m≤(M-1),0≤n≤(N-1),因此经过距离维压缩后的回波数据是一个M×N的矩阵,该矩阵的每一行为一个复距离像序列(取模后为一维距离像序列),每一列为一个距离单元的复包络序列。
因此,M×N维的采样矩阵S写作:
S=[s0 s1 ... sN-1]M×N (4)
步骤二二、采样协方差矩阵C为:
Figure BDA0001601242760000044
其中,采样协方差矩阵C为Hermitian(厄尔米特)矩阵,E[·]为取期望,SH为S的共轭转置,
Figure BDA0001601242760000045
为sd的共轭转置。
其它步骤及参数与具体实施方式一或二相同。
具体实施方式四:本实施方式与具体实施方式一至三之一不同的是:所述步骤三中利用Arnoldi迭代求解采样协方差矩阵的最大特征值对应的特征向量,提取最大特征值对应的特征向量中各元素的相位信息,利用各元素的相位信息对步骤一中的包络对齐后的距离维压缩数据进行相位补偿,得到相位补偿后的数据;具体过程为:
原理:
Rayleigh-Ritz(RR)逼近引理表明,Hermitian矩阵C的大特征值λ1>λ2>...>λc和大特征值对应的特征向量u1,u2,...,uc可以分别由矩阵C的RR值和RR向量逼近。而C的RR值和RR向量可由b×b维矩阵
Figure BDA0001601242760000051
的大特征值
Figure BDA0001601242760000052
和大特征值对应的特征向量
Figure BDA0001601242760000053
求出。因此,矩阵C的大特征值可由式(10)和式(11)近似求得,其中b>c。
Figure BDA0001601242760000054
Figure BDA0001601242760000055
由于C为Hermitian矩阵,因此Arnoldi迭代可用来计算得到矩阵Qb以及矩阵
Figure BDA0001601242760000056
Figure BDA0001601242760000057
进行特征分解即可进一步求取采样协方差矩阵C的最大特征值对应的特征向量。由于
Figure BDA0001601242760000058
的维度远小于矩阵C,因此,对
Figure BDA0001601242760000059
特征分解的计算量远小于直接对矩阵C进行特征分解的计算量。
其中,Arnoldi迭代推导过程为:
设酉矩阵Q=[q1 q2 ... qM]M×M
Figure BDA00016012427600000510
为Hessenberg矩阵,形如
Figure BDA00016012427600000511
其中,C为Hermitian矩阵,q1,q2,...,qM为Arnoldi向量,hik为Hessenberg矩阵
Figure BDA00016012427600000512
的第i行第k列元素;
由于QQH=I,对
Figure BDA00016012427600000513
等式两边左乘Q,得
Figure BDA00016012427600000514
其中I为单位矩阵;
对比
Figure BDA0001601242760000061
的两边,有
Figure BDA0001601242760000062
其中,qk为第k个Arnoldi向量,由
Figure BDA0001601242760000063
表示,qi为第i个Arnoldi向量,1≤i≤k+1;
由式(13)可得:
Figure BDA0001601242760000064
其中,
Figure BDA0001601242760000065
i=1,2,...,k;
若rk≠0,那么
qk+1=rk/hk+1,k (15)
其中,hk+1,k=||rk||2
因此,Arnoldi迭代步骤可以概括为:
Figure BDA0001601242760000066
其中,q1为给定的单位二范数向量,b为Arnoldi迭代长度,qk,k=1,...,b为迭代求出的Arnoldi向量。
通过迭代即可求得矩阵Qb=[q1 q2 ... qb]M×b以及矩阵
Figure BDA0001601242760000067
若q1选取不当,会带来一定的估计误差,Arnoldi的估计误差可以通过增加迭代长度b来减小,即迭代出更多的Arnoldi向量。
具体过程:
步骤三一、给定初始Arnoldi向量q1
其中,q1为单位二范数向量,
Figure BDA0001601242760000071
Figure BDA0001601242760000072
为M×1维复数向量空间;
步骤三二、求出Arnoldi向量q2
步骤三三、将步骤三二求出的Arnoldi向量q2以及步骤三一给定的初始Arnoldi向量q1,根据步骤三二求得Arnoldi向量q3,并依次求出Arnoldi向量q4,...,qb
其中,b为Arnoldi迭代长度(即要求迭代出的Arnoldi向量个数);
步骤三四、根据求出的Arnoldi向量q2,...,qb以及步骤三一给定的初始Arnoldi向量q1,构造矩阵Qb=[q1 q2 … qb]M×b,根据Qb计算得到矩阵
Figure BDA0001601242760000073
其中,b<<M,Qb为由求得的Arnoldi向量所构造的矩阵,
Figure BDA0001601242760000074
为b×b维的Hessenberg矩阵,
Figure BDA0001601242760000075
为Qb的共轭转置;
步骤三五、对矩阵
Figure BDA0001601242760000076
进行特征分解,求解出
Figure BDA0001601242760000077
的最大特征值所对应的特征向量
Figure BDA0001601242760000078
步骤三六、通过
Figure BDA0001601242760000079
求解得到采样协方差矩阵C的最大特征值所对应的特征向量u1
步骤三七、提取采样协方差矩阵C的最大特征值所对应的特征向量u1的各元素的相位信息,根据各元素的相位信息对步骤一中的包络对齐后的距离维压缩数据进行相位补偿,得到相位补偿后的数据。
其它步骤及参数与具体实施方式一至三之一相同。
具体实施方式五:本实施方式与具体实施方式一至四之一不同的是:所述步骤三二中求出Arnoldi向量q2;具体过程为:
设酉矩阵Q=[q1 q2 ... qM]M×M,Hessenberg矩阵
Figure BDA00016012427600000710
Figure BDA00016012427600000711
具体形式为
Figure BDA0001601242760000081
其中,C为采样协方差矩阵(此处矩阵C为步骤二求得的采样协方差矩阵,该采样协方差矩阵具有Hermitian属性),q1,q2,...,qM为Arnoldi向量,hik为Hessenberg矩阵
Figure BDA00016012427600000810
的第i行第k列元素,1≤i≤M,1≤k≤M,QH为Q的共轭转置;
由于QQH=I,对
Figure BDA0001601242760000082
等式两边左乘Q,得
Figure BDA0001601242760000083
其中I为单位矩阵;
对比
Figure BDA0001601242760000084
的两边,有
Figure BDA0001601242760000085
其中,qk为第k个Arnoldi向量,由
Figure BDA0001601242760000086
表示,qi为第i个Arnoldi向量,i=1,2,...,k+1;
将式(7)中
Figure BDA0001601242760000087
的hk+1,kqk+1项提取到等式一边得:
Figure BDA0001601242760000088
其中,
Figure BDA0001601242760000089
i=1,2,...,k,rk=hk+1,kqk+1,rk为中间变量;
若rk≠0,则
qk+1=rk/hk+1,k (9)
其中,hk+1,k=||rk||2
||·||2为二范数;
由式(6)至式(9)知,第k+1个Arnoldi向量qk+1由前k个Arnoldi向量q1,...,qk求出。
其它步骤及参数与具体实施方式一至四之一相同。
采用以下实施例验证本发明的有益效果:
实施例一:
本实施例一种基于Arnoldi的特征向量相位补偿方法,具体是按照以下步骤制备的:
采用实测数据来验证本发明的有效性,本发明相对于现有基于特征向量的相位补偿方法在保证相位补偿性能的同时,降低了计算量,更适用于实时ISAR成像。
实测数据由载频为5520MHz、带宽为400MHz、脉冲重复频率为400Hz、脉宽为25.6us的雷达系统测得。用于成像的总回波数为2048,积累互相关法中脉冲积累数为100,基于Arnoldi的特征向量相位补偿方法中Arnoldi迭代长度为6。
图2为现有基于特征向量的相位补偿方法的ISAR成像结果,图3为本发明基于Arnoldi的特征向量相位补偿方法的ISAR成像结果。从图2图3可以看出,本发明所提出的基于Arnoldi的特征向量相位补偿方法的ISAR成像效果和现有基于特征向量的相位补偿方法几乎一致。
表1为现有基于特征向量的相位补偿方法和本发明所提出的基于Arnoldi的特征向量相位补偿方法的计算时长对比。由于本发明所提出的相位补偿方法仅在计算包含相位误差信息的特征向量步骤区别于现有方法,因此,表1中计算时长的对比仅涉及该步骤的计算复杂度,表中数据为100次蒙特卡洛实验的总计算时长。从表1可以看出,本发明所提出的基于Arnoldi的特征向量相位补偿方法的计算量远小于现有基于特征向量的相位补偿方法。因此,本发明所提出的基于Arnoldi的特征向量相位补偿方法更适合用于实时ISAR成像系统的相位补偿。
表1计算时长对比(s)
Figure BDA0001601242760000091
表2为不同Arnoldi迭代长度下的ISAR图像熵。经过现有基于特征向量的相位补偿方法所成图像的熵为8.8979。由表2可以看出Arnoldi迭代长度b越长,所对应的ISAR图像的熵越小,即本发明所提出的基于Arnoldi的特征向量相位补偿方法的相位补偿效果越好,但当熵降低到一定程度时便不再继续下降,此时经过本发明所提出的相位补偿方法所成图像的熵近似等于经过现有基于特征向量的相位补偿方法所成图像的熵。即当参数选择恰当时,本发明拥有和现有方法一致的相位补偿性能,但计算速度大大加快。
表2不同Arnoldi迭代长度下的ISAR图像熵
Figure BDA0001601242760000092
Figure BDA0001601242760000101
表3为不同Arnoldi迭代长度下的计算时长。由于本发明所提出的相位补偿方法仅在计算包含相位误差信息的特征向量步骤区别于现有方法,因此,图中计算时长的数据仅涉及该步骤的计算时长,图中数据均为100次蒙特卡洛实验的总计算时长。现有基于特征向量的相位补偿算法的计算时长为243.1921s。由表3可以看出,随着Arnoldi迭代长度b的增加,相应的计算时长增加。因此,恰当的参数选择对本发明的效率有着较为重要意义。由于Arnoldi迭代长度d远小于采样协方差矩阵维度2048,因此,结合表2可知,当参数选择合适时,本发明所提出的基于Arnoldi的特征向量相位补偿方法拥有和现有基于特征向量的相位补偿方法相似的相位误差补偿效果,但大大减少了计算时长。
表3不同Arnoldi迭代长度下的计算时长(s)
Figure BDA0001601242760000102
从以上结果可以看出,本发明所提出的基于Arnoldi的特征向量相位补偿方法在保证相位误差估计精度的同时大量减少了计算量,相比于现有基于特征向量的相位补偿方法更适合应用于实时系统。
本发明还可有其它多种实施例,在不背离本发明精神及其实质的情况下,本领域技术人员当可根据本发明作出各种相应的改变和变形,但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。

Claims (4)

1.一种基于Arnoldi的特征向量相位补偿方法,其特征在于:所述方法具体过程为:
步骤一、取M个目标回波,经过距离维压缩后的第m次目标回波沿距离n分布的复距离像为s(m,n),对复距离像取模得到目标回波的一维距离像,并利用积累互相关方法进行包络对齐,得到包络对齐后的距离维压缩数据;
其中,0≤m≤(M-1),0≤n≤(N-1),M为方位向点数,N为距离向点数;
步骤二、根据包络对齐后的距离维压缩数据求得采样协方差矩阵;
步骤三、利用Arnoldi迭代求解采样协方差矩阵的最大特征值对应的特征向量,提取最大特征值对应的特征向量中各元素的相位信息,利用各元素的相位信息对步骤一中的包络对齐后的距离维压缩数据进行相位补偿,得到相位补偿后的数据;
步骤四、利用距离多普勒方法对相位补偿后的数据进行ISAR成像;
所述ISAR为逆合成孔径雷达;
所述步骤一中取M个目标回波,经过距离维压缩后的第m次目标回波沿距离n分布的复距离像为s(m,n),对复距离像取模得到目标回波的一维距离像,并利用积累互相关方法进行包络对齐,得到包络对齐后的距离维压缩数据;具体过程为:
步骤一一、经过距离维压缩后的第m次目标回波沿距离n分布的复距离像表示为s(m,n);
其中,0≤m≤(M-1),0≤n≤(N-1),M为方位向点数,N为距离向点数,N、M取值为正整数;
步骤一二、对复距离像取模,得到目标回波的一维距离像:
p(m,n)=abs{s(m,n)} (1)
其中,p(m,n)为目标回波的一维距离像,abs{·}为取绝对值;
步骤一三、利用积累互相关方法将当前一维距离像与前l个一维距离像进行互相关:
Figure FDA0003153180210000011
其中,l为脉冲积累数,R(τ)为当前一维距离像与前l个一维距离像的互相关函数,τ为时延,
Figure FDA0003153180210000012
为前l个一维距离像的和,p*(m+a,n)为p(m+a,n)的共轭,p(m+l,n)为当前所要对齐的一维距离像,a为叠加变量,a=0,1,...,l-1;
对τ进行搜索,计算R(τ)的最大值所对应的时延τ,将R(τ)的最大值所对应的时延τ作为补偿,来实现包络的对齐,得到包络对齐后的距离维压缩数据。
2.根据权利要求1所述一种基于Arnoldi的特征向量相位补偿方法,其特征在于:所述步骤二中根据包络对齐后的距离维压缩数据求得采样协方差矩阵;具体过程为:
步骤二一、第d个距离单元的距离维压缩数据的复包络sd表示为:
Figure FDA0003153180210000021
其中,ad为第d个距离单元的强散射点的散射率,d=0,1,2,...,N-1,v为相位误差向量,nd为高斯白噪声,
Figure FDA0003153180210000022
为第e+1个方位向单元与第一个方位向单元之间的相位差,e=1,2,...,M-1,
Figure FDA0003153180210000023
为虚数单位;
因此,M×N维的采样矩阵S写作:
S=[s0 s1 ... sN-1]M×N (4)
步骤二二、采样协方差矩阵C为:
Figure FDA0003153180210000024
其中,采样协方差矩阵C为Hermitian矩阵,E[·]为取期望,SH为S的共轭转置,
Figure FDA0003153180210000025
为sd的共轭转置。
3.根据权利要求2所述一种基于Arnoldi的特征向量相位补偿方法,其特征在于:所述步骤三中利用Arnoldi迭代求解采样协方差矩阵的最大特征值对应的特征向量,提取最大特征值对应的特征向量中各元素的相位信息,利用各元素的相位信息对步骤一中的包络对齐后的距离维压缩数据进行相位补偿,得到相位补偿后的数据;具体过程为:
步骤三一、给定初始Arnoldi向量q1
其中,q1为单位二范数向量,
Figure FDA0003153180210000026
Figure FDA0003153180210000027
为M×1维复数向量空间;
步骤三二、求出Arnoldi向量q2
步骤三三、将步骤三二求出的Arnoldi向量q2以及步骤三一给定的初始Arnoldi向量q1,根据步骤三二求得Arnoldi向量q3,并依次求出Arnoldi向量q4,...,qb
其中,b为Arnoldi迭代长度;
步骤三四、根据求出的Arnoldi向量q2,...,qb以及步骤三一给定的初始Arnoldi向量q1,构造矩阵Qb=[q1 q2 … qb]M×b,根据Qb计算得到矩阵
Figure FDA0003153180210000031
其中,b<<M,Qb为由求得的Arnoldi向量所构造的矩阵,
Figure FDA0003153180210000032
为b×b维的Hessenberg矩阵,
Figure FDA0003153180210000033
为Qb的共轭转置;
步骤三五、对矩阵
Figure FDA0003153180210000034
进行特征分解,求解出
Figure FDA0003153180210000035
的最大特征值所对应的特征向量
Figure FDA0003153180210000036
步骤三六、通过
Figure FDA0003153180210000037
求解得到采样协方差矩阵C的最大特征值所对应的特征向量u1
步骤三七、提取采样协方差矩阵C的最大特征值所对应的特征向量u1的各元素的相位信息,根据各元素的相位信息对步骤一中的包络对齐后的距离维压缩数据进行相位补偿,得到相位补偿后的数据。
4.根据权利要求3所述一种基于Arnoldi的特征向量相位补偿方法,其特征在于:所述步骤三二中求出Arnoldi向量q2;具体过程为:
设酉矩阵Q=[q1 q2 … qM]M×M,Hessenberg矩阵
Figure FDA0003153180210000038
Figure FDA0003153180210000039
具体形式为
Figure FDA00031531802100000310
其中,C为采样协方差矩阵,q1,q2,...,qM为Arnoldi向量,hik为Hessenberg矩阵
Figure FDA00031531802100000311
的第i行第k列元素,1≤i≤M,1≤k≤M,QH为Q的共轭转置;
由于QQH=I,对
Figure FDA00031531802100000312
等式两边左乘Q,得
Figure FDA00031531802100000313
其中I为单位矩阵;
对比
Figure FDA0003153180210000041
的两边,有
Figure FDA0003153180210000042
其中,qk为第k个Arnoldi向量,qi为第i个Arnoldi向量,i=1,2,...,k+1,1≤k≤(M-1);
将式(7)中
Figure FDA0003153180210000043
的hk+1,kqk+1项提取到等式一边得:
Figure FDA0003153180210000044
其中,
Figure FDA0003153180210000045
rk=hk+1,kqk+1,rk为中间变量;
qk+1=rk/hk+1,k (9)
其中,hk+1,k=||rk||2
||·||2为二范数;
由式(6)至式(9)知,第k+1个Arnoldi向量qk+1由前k个Arnoldi向量q1,...,qk求出。
CN201810225591.0A 2018-03-19 2018-03-19 一种基于Arnoldi的特征向量相位补偿方法 Active CN108445465B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810225591.0A CN108445465B (zh) 2018-03-19 2018-03-19 一种基于Arnoldi的特征向量相位补偿方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810225591.0A CN108445465B (zh) 2018-03-19 2018-03-19 一种基于Arnoldi的特征向量相位补偿方法

Publications (2)

Publication Number Publication Date
CN108445465A CN108445465A (zh) 2018-08-24
CN108445465B true CN108445465B (zh) 2021-09-07

Family

ID=63195727

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810225591.0A Active CN108445465B (zh) 2018-03-19 2018-03-19 一种基于Arnoldi的特征向量相位补偿方法

Country Status (1)

Country Link
CN (1) CN108445465B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5627543A (en) * 1994-08-05 1997-05-06 Deutsche Forschungsanstalt Fur Luft-Und Raumfahrt E.V. Method of image generation by means of two-dimensional data processing in connection with a radar with synthetic aperture
JPH11248834A (ja) * 1998-03-06 1999-09-17 Mitsubishi Electric Corp レーダ装置
JP2004020380A (ja) * 2002-06-17 2004-01-22 Mitsubishi Electric Corp レーダ装置およびisar画像の視認性向上方法
CN103675816A (zh) * 2013-11-20 2014-03-26 中国科学院电子学研究所 一种基于压缩感知理论的干涉逆合成孔径雷达成像方法
CN103901429A (zh) * 2014-04-09 2014-07-02 西安电子科技大学 基于稀疏孔径的机动目标逆合成孔径雷达成像方法
CN105425238A (zh) * 2016-01-08 2016-03-23 中国人民解放军国防科学技术大学 一种逆合成孔径雷达成像并行处理方法
CN106291551A (zh) * 2016-07-28 2017-01-04 南京航空航天大学 一种基于gpu的并行结构isar距离对准方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5627543A (en) * 1994-08-05 1997-05-06 Deutsche Forschungsanstalt Fur Luft-Und Raumfahrt E.V. Method of image generation by means of two-dimensional data processing in connection with a radar with synthetic aperture
JPH11248834A (ja) * 1998-03-06 1999-09-17 Mitsubishi Electric Corp レーダ装置
JP2004020380A (ja) * 2002-06-17 2004-01-22 Mitsubishi Electric Corp レーダ装置およびisar画像の視認性向上方法
CN103675816A (zh) * 2013-11-20 2014-03-26 中国科学院电子学研究所 一种基于压缩感知理论的干涉逆合成孔径雷达成像方法
CN103901429A (zh) * 2014-04-09 2014-07-02 西安电子科技大学 基于稀疏孔径的机动目标逆合成孔径雷达成像方法
CN105425238A (zh) * 2016-01-08 2016-03-23 中国人民解放军国防科学技术大学 一种逆合成孔径雷达成像并行处理方法
CN106291551A (zh) * 2016-07-28 2017-01-04 南京航空航天大学 一种基于gpu的并行结构isar距离对准方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
逆合成孔径雷达成像包络对齐性能评估方法研究;黄璐 等;《无线电工程》;20170331;第47卷(第3期);全文 *

Also Published As

Publication number Publication date
CN108445465A (zh) 2018-08-24

Similar Documents

Publication Publication Date Title
CN109100718B (zh) 基于贝叶斯学习的稀疏孔径isar自聚焦与横向定标方法
CN104020469B (zh) 一种mimo雷达距离-角度二维超分辨率成像算法
CN103901429B (zh) 基于稀疏孔径的机动目标逆合成孔径雷达成像方法
Li et al. Doppler keystone transform: An approach suitable for parallel implementation of SAR moving target imaging
CN105652273B (zh) 一种基于混合匹配追踪算法的mimo雷达稀疏成像算法
CN108459307B (zh) 基于杂波的mimo雷达收发阵列幅相误差校正方法
CN106680815B (zh) 基于张量稀疏表示的mimo雷达成像方法
CN110244303B (zh) 基于sbl-admm的稀疏孔径isar成像方法
CN108008386B (zh) 一种基于单快拍music算法的距离向处理方法
CN105445704B (zh) 一种sar图像中的雷达动目标抑制方法
CN104698431B (zh) 基于模糊分量doa估计的多通道sar方位解模糊方法
CN104977582A (zh) 一种实现扫描雷达方位超分辨成像的解卷积方法
CN110113085A (zh) 一种基于协方差矩阵重构的波束形成方法及系统
CN105301589B (zh) 高分辨宽测绘带sar地面运动目标成像方法
CN110109098B (zh) 一种扫描雷达快速超分辨成像方法
CN114371478B (zh) 基于单天线系统的机载雷达对舰船目标三维成像方法
CN112305495A (zh) 一种基于原子范数最小的互质阵列协方差矩阵重构方法
CN104777479A (zh) 基于多核dsp的前侧视sar实时成像方法
CN105223554A (zh) 基于多普勒三通道联合处理的空时自适应单脉冲测角方法
CN105785326A (zh) 一种非正侧视阵雷达杂波谱配准的优化方法
CN107703478A (zh) 基于互相关矩阵的扩展孔径二维doa估计方法
CN112147608A (zh) 一种快速高斯网格化非均匀fft穿墙成像雷达bp方法
CN109946672B (zh) 基于被动孔径合成稀疏阵列的doa估计方法
CN106546964A (zh) 一种机载雷达杂波秩估计方法及其系统
CN109946663A (zh) 一种线性复杂度的Massive MIMO目标空间方位估计方法和装置

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