CN104007429B - 基于极化分解的噪声稳健宽带全极化目标识别方法 - Google Patents
基于极化分解的噪声稳健宽带全极化目标识别方法 Download PDFInfo
- Publication number
- CN104007429B CN104007429B CN201410216267.4A CN201410216267A CN104007429B CN 104007429 B CN104007429 B CN 104007429B CN 201410216267 A CN201410216267 A CN 201410216267A CN 104007429 B CN104007429 B CN 104007429B
- Authority
- CN
- China
- Prior art keywords
- prime
- distance image
- training sample
- lambda
- image signal
- 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
- 230000010287 polarization Effects 0.000 title claims abstract description 210
- 238000000034 method Methods 0.000 title claims abstract description 47
- 238000000354 decomposition reaction Methods 0.000 title claims abstract description 18
- 238000012549 training Methods 0.000 claims abstract description 206
- 238000012360 testing method Methods 0.000 claims abstract description 203
- 230000004927 fusion Effects 0.000 claims abstract description 14
- 239000011159 matrix material Substances 0.000 claims description 119
- 239000013598 vector Substances 0.000 claims description 95
- 230000001427 coherent effect Effects 0.000 claims description 25
- 238000010606 normalization Methods 0.000 claims description 25
- 238000012545 processing Methods 0.000 claims description 24
- 238000012937 correction Methods 0.000 claims description 21
- 230000008569 process Effects 0.000 claims description 10
- 238000000605 extraction Methods 0.000 claims description 6
- 238000012935 Averaging Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000012512 characterization method Methods 0.000 claims description 4
- 125000004122 cyclic group Chemical group 0.000 claims description 2
- 230000000875 corresponding effect Effects 0.000 description 140
- 238000010586 diagram Methods 0.000 description 8
- 238000004088 simulation Methods 0.000 description 6
- 238000007796 conventional method Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000013075 data extraction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 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
- G01S7/414—Discriminating targets with respect to background clutter
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
本发明属于雷达目标识别技术领域,特别涉及基于极化分解的噪声稳健宽带全极化目标识别方法。该基于极化分解的噪声稳健宽带全极化目标识别方法包括以下步骤:将每组训练用全极化通道距离像信号作为对应的训练样本,针对每个训练样本,获取对应训练样本的平均距离像和分类特征;获取对应测试样本的平均距离像和分类特征;将对应测试样本的平均距离像分别与每个训练样本的平均距离像进行对齐;得出每个训练样本的对应测试样本融合特征;根据每个训练样本的对应测试样本融合特征、以及对应训练样本的分类特征,对对应测试样本进行分类识别。
Description
技术领域
本发明属于雷达目标识别技术领域,特别涉及基于极化分解的噪声稳健宽带全极化目标识别方法,可用于对不同类型飞机或车辆目标进行分类。
背景技术
宽带雷达目标识别可以利用目标的高分辨雷达回波信号,实现对目标类型的判定。高分辨雷达目标回波能够反应出目标更加详细的信息,具有低分辨雷达所不具有的独特优势。高分辨距离像提供了目标沿距离方向的几何结构信息,获取容易、处理方便,在目标识别方面有很好的应用。在单极化体制下,雷达发射和接收电磁波的极化方式是不变的。而目标在不同的极化电磁波照射下,电磁散射特性是不同的,通过不同极化方式的发射和接收,可以的到目标更加全面的极化信息。目标的极化信息和宽带信息从不同方面刻画了目标的散射特性,将两者结合起来也能更好的用于目标识别领域。
上述基于高分辨距离像、目标全极化信息以及两者相结合的方法都对噪声具有很强的敏感性。但是实际中由于接收机热噪声一定,目标距离雷达越远,回波信噪比越低。通常应用中,用于训练阶段的距离像数据可由合作的实验获取或者是直接由电磁仿真程序生成,其信噪比较高;而用于测试阶段的距离像数据经常在非合作情况下取得,无法保证高信噪比,目标数据提取的特征会发生变化。此时,训练和测试之间产生模型失配,导致了识别性能随着测试样本信噪比的降低而急剧恶化,造成识别率明显下降。
发明内容
本发明的目的在于提出基于极化分解的噪声稳健宽带全极化目标识别方法。本发明能够在低信噪比条件下,有效改善目标识别性能,从而提升识别距离。
为实现上述技术目的,本发明采用如下技术方案予以实现。
基于极化分解的噪声稳健宽带全极化目标识别方法包括以下步骤:
S1:利用雷达接收多组全极化通道距离像信号;获取全极化通道距离像信号的多个训练样本;
S2:针对每个训练样本,采用目标极化特征提取方法,获取每个训练样本的平均距离像和分类特征;
S3:将雷达接收的每组全极化通道距离像信号作为对应的测试样本,针对每个测试样本采用基于噪声修正的目标极化特征提取方法,获取每个测试样本的平均距离像和分类特征;针对每个测试样本,执行步骤S4至步骤S6;
S4:将对应测试样本的平均距离像分别与每个训练样本的平均距离像进行对齐;得出对应测试样本的平均距离像在与每个训练样本的平均距离像进行对齐时,向右移动的距离单元数;
S5:根据对应测试样本的平均距离像在与每个训练样本的平均距离像进行对齐时,向右移动的距离单元数,得出每个训练样本的对应测试样本融合特征;
S6:根据每个训练样本的对应测试样本融合特征、以及每个训练样本的分类特征,得出对应测试样本对应的目标类别。
本发明的有益效果为:本发明在信噪比变化时目标正确识别率没有大的变化,具有识别性能噪声稳健的优点,可用于在低信噪比宽带全极化目标识别。
附图说明
图1为本发明的基于极化分解的噪声稳健宽带全极化目标识别方法的流程图;
图2a,为现有方法和本发明得出的测试样本的极化相干矩阵的最大特征值所占比例随信噪比变化对比示意图;
图2b为现有方法和本发明得出的测试样本的极化相干矩阵对应的平均散射角随信噪比变化对比示意图;
图2c为现有方法和本发明得出的测试样本的极化相干矩阵对应的散射熵随信噪比变化对比示意图;
图3为现有方法和本发明得出的测试样本的极化相干矩阵的重构误差随信噪比变化的对比图;
图4为现有方法和本发明得出的测试样本对应的目标正确识别率随信噪比变化的对比示意图。
具体实施方式
下面结合附图对本发明作进一步说明:
参照图1,为本发明的基于极化分解的噪声稳健宽带全极化目标识别方法的流程图。该基于极化分解的噪声稳健宽带全极化目标识别方法包括以下步骤:
S1:利用雷达接收多组全极化通道距离像信号;获取全极化通道距离像信号的多个训练样本。具体说明如下:
上述全极化通道距离像信号包括HH极化通道距离像信号、HV极化通道距离像信号、VH极化通道距离像信号、以及VV极化通道距离像信号;在接收的全极化通道距离像信号中,选择信噪比最高的N1组全极化通道距离像信号作为对应的N1组训练用全极化通道距离像信号,N1为大于1的自然数。
全极化通道距离像信号的多个训练样本可以通过雷达获取,也可以通过仿真获取。在全极化通道距离像信号的每个训练样本中,目标方位角(目标相对于雷达的方位角)和目标类别(例如为直升机、喷气式飞机或螺旋桨飞机)均已知。
S2:针对每个训练样本,采用目标极化特征提取方法,获取每个训练样本的平均距离像和分类特征。具体说明如下:
在步骤S2中,获取每个训练样本的平均距离像的过程相同,获取每个训练样本的分类特征的过程相同。因此,下面以任一个训练样本为例,说明获取对应训练样本的平均距离像和分类特征的过程,获取对应训练样本的平均距离像和分类特征包括以下步骤:
S21:分别对对应训练样本的HH极化通道距离像信号、对应训练样本的HV极化通道距离像信号、对应训练样本的VH极化通道距离像信号、以及对应训练样本的VV极化通道距离像信号进行归一化处理,得到对应训练样本的HH极化通道归一化距离像信号、对应训练样本的HV极化通道归一化距离像信号、对应训练样本的VH极化通道归一化距离像信号、以及对应训练样本的VV极化通道归一化距离像信号;估计出对应训练样本的每个极化通道距离像信号的噪声平均功率。
具体地,在步骤S21中,对应训练样本的HH极化通道距离像信号表示为sHH,对应训练样本的HV极化通道距离像信号表示为sHV,对应训练样本的VH极化通道距离像信号表示为sVH,对应训练样本的VV极化通道距离像信号表示为sVV;其中,
sHH=[sHH1,...,sHHk,...,sHHM],sHV=[sHV1,...,sHVk,...,sHVM]
sVH=[sVH1,...,sVHk,...,sVHM],sVV=[sVV1,...,sVVk,...,sVVM]
其中,sHHk表示对应训练样本的HH极化通道第k个距离单元的距离像信号,sHVk表示对应训练样本的HV极化通道第k个距离单元的距离像信号,sVHk表示对应训练样本的VH极化通道第k个距离单元的距离像信号,sVVk表示对应训练样本的VV极化通道第k个距离单元的距离像信号,k取1至M,M为雷达接收的每个极化通道距离像信号(HH极化通道距离像信号sHH、HV极化通道距离像信号sHV、VH极化通道距离像信号sVH或VV极化通道距离像信号sVV)的长度;
分别对对应训练样本的HH极化通道距离像信号sHH、对应训练样本的HV极化通道距离像信号sHV、对应训练样本的VH极化通道距离像信号sVH、以及对应训练样本的VV极化通道距离像信号sVV进行归一化处理,得到sHH经归一化处理后的信号sHV经归一化处理后的信号sVH经归一化处理后的信号以及sVV经归一化处理后的信号
其中,表示对应训练样本的HH极化通道第k个距离单元经归一化处理之后的距离像信号,表示对应训练样本的HV极化通道第k个距离单元经归一化处理之后的距离像信号,表示对应训练样本的VH极化通道第k个距离单元经归一化处理之后的距离像信号,表示对应训练样本的VV极化通道第k个距离单元经归一化处理之后的距离像信号;以及分别为:
其中,‖·‖2表示求2-范数;
然后,估计出对应训练样本的每个极化通道距离像信号的噪声平均功率,对应训练样本的每个极化通道距离像信号的噪声平均功率均为σ2。
S22:针对每个训练样本,在所有训练样本选取与对应训练样本对应的多个训练样本,选取的每个训练样本的目标方位角与对应训练样本的目标方位角的角度差小于设定角度阈值;每个训练样本的目标方位角指:每个训练样本中目标相对于雷达的方位角;例如,设定角度范围为5度,则选取的任一个训练样本的目标方位角与对应训练样本的目标方位角的角度差小于5度。
根据选取的每个训练样本的HH极化通道归一化距离像信号、选取的每个训练样本的HV极化通道归一化距离像信号、选取的每个训练样本的VH极化通道归一化距离像信号、以及选取的每个训练样本的VV极化通道归一化距离像信号,得出对应训练样本的极化相干矩阵T、以及对应训练样本的平均距离像
具体地,在步骤S22中,根据以下公式得出对应训练样本的极化散射矩阵S:
S=[S1,…Sk,…,SM]
其中,Sk为:
根据以下公式得出对应训练样本的极化相干矩阵T:
T=[T1,…Tk,…,TM]
其中,|·|表示取模,上标*表示矩阵的共轭,在步骤S22中,对于选取的每个训练样本,均有对应的A值、B值和C值,<·>表示:根据步骤S22选取的多个训练样本求集合平均;例如,选取的训练样本的个数为5个,根据选取的这5个训练样本,得出对应的5个A值、5个B值和5个C值,此时AB*的值就对应有5个,<AB*>表示对这5个AB*的值作集合平均计算。
根据步骤S22选取的多个训练样本,对或进行集合平均计算,得出对应训练样本的平均距离像 或例如,选取的训练样本的个数为5个,根据选取的这5个训练样本,得出对应的5个的向量值,表示对这5个的向量值作集合平均计算(将5个的对应位置的元素作算术平均运算,得出新的元素,则得出的新的元素组合成)。
S23:根据对应训练样本的极化相干矩阵T、以及对应训练样本的每个极化通道距离像信号的噪声平均功率,得出对应训练样本的修正后极化相干矩阵根据对应训练样本的修正后极化相干矩阵得出对应训练样本的分类特征。
具体地,在步骤S23中,在雷达接收的4个极化通道距离像信号中,设定每个极化通道的噪声均为服从相同分布(高斯分布)的高斯白噪声,且设定噪声与信号互不相关。即<SZ*>=<ZS*>,其中S表示雷达接收的任一个极化通道中不含噪声的距离像信号,Z表示雷达接收的任一个极化通道的噪声,上标*表示矩阵的共轭,<·>表示求集合平均。
则雷达接收的不含噪声的距离像信号的极化相干矩阵TS可以表示为:
雷达接收的噪声的极化相干矩阵ZS可以表示为:
其中,SHH表示雷达接收的HH极化通道的不含噪声的距离像信号,SHV表示示雷达接收的HV极化通道的不含噪声的距离像信号,SVH表示雷达接收的VH极化通道的不含噪声的距离像信号,SVV表示雷达接收的VV极化通道的不含噪声的距离像信号;ZHH表示雷达接收的HH极化通道的噪声,ZHV表示雷达接收的HV极化通道的噪声,ZVH表示雷达接收的VH极化通道的噪声,ZVV表示雷达接收的VV极化通道的噪声。噪声的平均功率为为:
本发明实施例中,设定雷达接收的各个极化通道的噪声声不相关,即:<ZiZj *>=<ZjZi *>,其中,i∈{HH,HV,VH,VV},j∈{HH,HV,VH,VV},i≠j。
因此TZ为:
根据矩阵分解理论可得,
其中,T表示雷达接收的含有噪声的距离像信号的极化相干矩阵,λT1、λT2和λT3分别表示不Ts的三个特征值。
由上可知,对雷达接收的含有噪声的距离像信号的极化相干矩阵T做特征值分解之后,并根据噪声的平均功率,就可以对噪声进行修正。
然后,针对对应训练样本,分别对T1至TM进行修正,得出矩阵至然后利用矩阵至组合成对应训练样本的修正后极化相干矩阵对Tk进行修正包括以下步骤:
首先对Tk进行特征值分解,得出Tk=UkΛkUk -1,其中,上标-1表示求矩阵的逆,Λk为3×3维的矩阵,Λk为:
其中,λ1k为Tk的第1个特征值,λ2k为Tk的第2个特征值,λ3k为Tk的第3个特征值,λ1k≥λ2k≥λ3k≥0;Uk为3×3维的矩阵,Uk=[U1k,U2k,U3k],U1k为与λ1k对应的Tk的特征向量,U2k为与λ2k对应的Tk的特征向量,U3k为与λ3k对应的Tk的特征向量;U1k、U2k和U3k分别表示为如下形式:
其中,α1k、α2k、α3k、β1k、β2k、β3k、δ1k、δ2k、δ3k、γ1k、γ2k、以及γ3k均为角度参数。
根据对应训练样本的每个极化通道距离像信号的噪声平均功率σ2,对Λk进行修正,得到修正后特征值矩阵为:
则得出矩阵 此时有λ1k-σ2≥λ2k-σ2≥λ3k-σ2≥0。
由于k取1至M,因此通过遍历可得出矩阵至在矩阵至中,针对每个矩阵提取归一化后比例特征、归一化后散射熵特征和归一化后平均散射角特征;在矩阵中提取归一化后比例特征(用于表示对应距离单元最大特征值与所有特征值之和的比例)、归一化后散射熵特征和归一化后平均散射角特征包括以下步骤:
根据以下公式矩阵的比例特征Pk、散射熵特征Hk、以及平均散射角特征
其中g取1至3;由于k取1至M,因此通过遍历便得出P1至PM、H1至HM、以及至然后,计算出P1至PM的均值μP、P1至PM的标准差σP、H1至HM的均值μH、H1至HM的标准差σH、至的均值以及至的标准差
根据以下公式得出矩阵的归一化后比例特征归一化后散射熵特征和归一化后平均散射角特征
然后将矩阵至的归一化后比例特征组合成矢量,将矩阵至的归一化后散射熵特征组合成矢量将矩阵至的归一化后平均散射角特征组合成矢量矢量矢量和矢量分别为:
将矢量矢量和矢量组合成对应训练样本的分类特征C,C为:
S3:将雷达接收的每组全极化通道距离像信号作为对应的测试样本,针对每个测试样本采用基于噪声修正的目标极化特征提取方法,获取每个测试样本的平均距离像和分类特征;针对每个测试样本,执行步骤S4至步骤S6。具体说明如下:
在步骤S3中,获取每个测试训练样本的平均距离像和分类特征包括以下步骤:
S31:分别对每个测试样本的HH极化通道距离像信号、每个测试样本的HV极化通道距离像信号、每个测试样本的VH极化通道距离像信号、以及每个测试样本的VV极化通道距离像信号进行归一化处理,得到每个测试样本的HH极化通道归一化距离像信号、每个测试样本的HV极化通道归一化距离像信号、每个测试样本的VH极化通道归一化距离像信号、以及每个测试样本的VV极化通道归一化距离像信号;估计出每个测试样本的每个极化通道距离像信号的噪声平均功率。
具体地,步骤S31包括以下步骤:
在步骤S31中,对应测试样本的HH极化通道距离像信号表示为xHH,对应测试样本的HV极化通道距离像信号表示为xHV,对应测试样本的VH极化通道距离像信号表示为xVH,对应测试样本的VV极化通道距离像信号xVV;其中,
xHH=[xHH1,...,xHHk,...,xHHM],xHV=[xHV1,...,xHVk,...,xHVM]
xVH=[xVH1,...,xVHk,...,xVHM],xVV=[xVV1,...,xVVk,...,xVVM]
其中,xHHk表示对应测试样本的HH极化通道第k个距离单元的距离像信号,xHVk表示对应测试样本的HV极化通道第k个距离单元的距离像信号,xVHk表示对应测试样本的VH极化通道第k个距离单元的距离像信号,xVVk表示对应测试样本的VV极化通道第k个距离单元的距离像信号,k取1至M,M为雷达接收的每个极化通道距离像信号的长度;
分别对对应测试样本的HH极化通道距离像信号xHH、对应测试样本的HV极化通道距离像信号xHV、对应测试样本的VH极化通道距离像信号xVH、以及测试样本的VV极化通道距离像信号xVV进行归一化处理,得到xHH经归一化处理后的信号xHV经归一化处理后的信号xVH经归一化处理后的信号以及xVV经归一化处理后的信号
其中,表示对应测试样本的HH极化通道第k个距离单元经归一化处理之后的距离像信号,表示对应测试样本的HV极化通道第k个距离单元经归一化处理之后的距离像信号,表示对应测试样本的VH极化通道第k个距离单元经归一化处理之后的距离像信号,表示对应测试样本的VV极化通道第k个距离单元经归一化处理之后的距离像信号;以及分别为:
其中,‖·‖2表示求2-范数;
然后,估计出对应测试样本的每个极化通道距离像信号的噪声平均功率,对应测试样本的每个极化通道距离像信号的噪声平均功率均为σ′2;
S32:针对每个测试样本,在所有测试样本选取与对应测试样本对应的多个测试样本,选取的每个测试样本的雷达接收时刻与对应测试样本的雷达接收时刻的差小于设定时间长度;每个测试样本的雷达接收时刻指:雷达接收对应一组全极化通道距离像信号的时刻。例如,设定时间长短为5微秒,则选取的任一个测试样本的雷达接收时刻与对应测试样本的雷达接收时刻的差小于5微秒。
根据选取的每个测试样本的HH极化通道归一化距离像信号、选取的每个测试样本的HV极化通道归一化距离像信号、选取的每个测试样本的VH极化通道归一化距离像信号、以及选取的每个测试样本的VV极化通道归一化距离像信号,得出对应测试样本的极化相干矩阵T′、以及对应测试样本的平均距离像
具体地,在步骤S32中,根据以下公式得出对应训练样本的极化散射矩阵S′:
S′=[S′1,…,S′k,…,S′M]
其中,S′k为:
根据以下公式得出测试样本的极化相干矩阵T′:
T=[T′1,…,T′k,…,T′M]
其中,|·|表示取模,上标*表示矩阵的共轭,在步骤S32中,对于选取的每个测试样本,均有对应的A′值、B′值和C′值,<·>表示:根据步骤S32选取的多个测试样本求集合平均;例如,选取的测试样本的个数为5个,根据选取的这5个测试样本,得出对应的5个A′值、5个B′值和5个C′值,此时A′B′*的值就对应有5个,<A′B′*>表示对这5个A′B′*的值作集合平均计算。
根据步骤S32选取的多个测试样本,对或进行集合平均计算,得出对应测试样本的平均距离像 或例如,选取的测试样本的个数为5个,根据选取的这5个测试样本,得出对应的5个的向量值,表示对这5个的向量值作集合平均计算(将5个的对应位置的元素作算术平均运算,得出新的元素,则得出的新的元素组合成
S33:根据对应测试样本的极化相干矩阵T′、以及对应测试样本的每个极化通道距离像信号的噪声平均功率,得出对应测试样本的修正后极化相干矩阵根据对应测试样本的修正后极化相干矩阵得出对应测试样本的分类特征。
具体地,在步骤S33中,针对对应测试样本,分别对T′1至T′M进行修正,得出矩阵至然后利用矩阵至组合成对应测试样本的修正后极化相干矩阵对进行修正包括以下步骤:
首先对进行特征值分解,得出其中,上标-1表示求矩阵的逆,Λ′k为3×3维的矩阵,Λ′k为:
其中,λ′1k为的第1个特征值,λ′2k为的第2个特征值,λ′3k为的第3个特征值,λ′1k≥λ′2k≥λ′3k≥0;U′k为3×3维的矩阵,U′k=[U′1k,U′2k,U′3k],U′1k为与λ′1k对应的的特征向量,U′2k为与λ′2k对应的的特征向量,U′3k为与λ′3k对应的的特征向量;U′1k、U′2k和U′3k分别表示为如下形式:
其中,α′1k、α′2k、α′3k、β′1k、β′2k、β′3k、δ′1k、δ′2k、δ′3k、γ′1k、γ′2k、以及γ′3k均为角度参数。
根据对应测试样本的每个极化通道距离像信号的噪声平均功率σ′2,对Λ′k进行修正,得到修正后特征值矩阵为:
其中,λ′1k-σ′2≥λ′2k-σ′2≥λ′3k-σ′2≥0。然后得出矩阵
在矩阵至中,针对每个矩阵提取归一化后比例特征、归一化后散射熵特征和归一化后平均散射角特征;在矩阵中提取归一化后比例特征(表示对应距离单元最大特征值与特征值之和的比例)、归一化后散射熵特征和归一化后平均散射角特征包括以下步骤:
根据以下公式矩阵的比例特征P′k散射熵特征H′k、以及平均散射角特征
其中g取1至3;由于k取1至M,因此通过遍历便得出P′1至P′M、H′1至H′M以及至然后计算出P′1至P′M的均值μ'P、P′1至P′M的标准差σ'P、H′1至H′M的均值μ'H、H′1至H′M的标准差σ'H、至的均值以及至的标准差
根据以下公式得出矩阵的归一化后比例特征归一化后散射熵特征和归一化后平均散射角特征
然后将矩阵至的归一化后比例特征组合成矢量将矩阵至的归一化后散射熵特征组合成矢量将矩阵至的归一化后平均散射角特征组合成矢量矢量矢量和矢量分别为:
将矢量矢量和矢量组合成对应测试样本的分类特征C',C'为:
S4:将对应测试样本的平均距离像分别与每个训练样本的平均距离像进行对齐;得出对应测试样本的平均距离像在与每个训练样本的平均距离像进行对齐时,向右移动的距离单元数。具体说明如下:
在将对应测试样本的平均距离像分别与每个训练样本的平均距离像进行对齐时,使用相关系数最大准则;在上述对齐过程中,对应测试样本的平均距离像向右依次循环移位。
任一个训练样本的平均距离像表示为向量对应测试样本的平均距离像表示为向量则将向量按照相关系数最大准则与向量进行对齐;在向量与向量进行对齐的过程中,将向量向右依次循环移位,向量在每次移动时向右循环移动一个距离单元(移动一位)。例如在步骤S4中,在将对应测试样本的平均距离像分别与每个训练样本的平均距离像进行对齐,使用相关系数最大准则;在所述对齐过程中,对应测试样本的平均距离像向右依次循环移位。
任一个训练样本的平均距离像表示为向量测试样本的平均距离像表示为向量则将向量按照相关系数最大准则与向量进行对齐;在向量与向量进行对齐的过程中,将向量向右依次循环移位,向量在每次移动时向右循环移动一个距离单元(移动一位)。下面举例说明向右循环移动的含义,当时,在第一次向右循环移动后,变为向量[5,1,2,3,4],在第二次向右循环移动后,变为向量[4,5,1,2,3]。
向量每次向右循环移动m个距离单元(移动m位)后,转变为向量m取1至M;计算向量与向量的相关系数rm:
其中,k取1至M,sk为对应训练样本的平均距离像中第k个距离单元的取值(中第k个元素的值),s'(m)k为向量中第k个距离单元的取值(向量中第k个元素的值),为对应训练样本的平均距离像中所有元素的均值,为向量中所有元素的均值。
得出对应测试样本的分类特征所需移动距离单元数N,N为:r1至rM中最大值对应的m的取值;即
S5:根据对应测试样本的平均距离像在与每个训练样本的平均距离像进行对齐时,向右移动的距离单元数,得出每个训练样本的对应测试样本融合特征。具体说明如下:
在步骤S5中,将矢量向右循环移动N个距离单元(移动N位),得到新矢量将向右循环移动N个距离单元(移动N位),得到新矢量将向右循环移动N个距离单元(移动N位),得到新矢量则对应测试样本融合特征表示为C″,
S6:根据每个训练样本的对应测试样本融合特征、以及每个训练样本的分类特征,得出对应测试样本对应的目标类别。具体说明如下:
在步骤S6中,针对每个训练样本的对应测试样本融合特征、以及每个训练样本的分类特征,得出每个训练样本的欧式距离;得出每个训练样本的欧式距离包括以下步骤:根据以下公式得出任一个训练样本的欧式距离:
d=||C″||2
其中,C″表示与对应训练样本的对应的测试样本融合特征,C表示对应训练样本的分类特征;||·||2表示2-范数。
由于在本发明实施例中共有N1组训练用全极化通道距离像信号,则训练样本的个数为N1。在步骤S6中,通过遍历得出每个训练样本对应的欧式距离,最终得出N1个欧式距离。将N1个欧式距离中最小值对应的训练样本作为识别用训练样本;将识别用训练样本的目标类别作为测试样本的目标类别。例如,识别用训练样本的目标类别为直升机,则测试样本的目标类别为直升机。
本发明的效果可以通过以下仿真实验验证:
1)实验场景:
仿真实验所用的数据包括xpacth仿真的3类地面目标数据,分别为2类坦克目标数据和1类消防车目标数据。目标的姿态为相对雷达方位角从0°到180°。按照相对雷达方位角1.5°的间隔,测试样本和训练样本进行交替抽取仿真数据。训练样本为未经处理的仿真数据(信噪比大约为30dB)。测试样本是在仿真数据的基础上加上0-30db信噪比的噪声。
2)实验内容:
用现有方法提取训练样本和测试样本的极化相干矩阵,并用本发明的方法进行目标分类识别。
3)实验结果分析:
参照图2a,为现有方法和本发明得出的测试样本的极化相干矩阵的最大特征值所占比例随信噪比变化对比示意图。图2a中,现有方法得出的测试样本的极化相干矩阵的最大特征值所占比例指:未进行极化相干矩阵修正时得出的测试样本的极化相干矩阵的最大特征值所占比例(图2a中表示为“未修正”)。而本发明在图2a中表示为“修正后”。在图2a中,横坐标表示信噪比,单位为dB;纵坐标表示测试样本的极化相干矩阵的最大特征值所占比例。
参照图2b,为现有方法和本发明得出的测试样本的极化相干矩阵对应的平均散射角随信噪比变化对比示意图。图2b中,现有方法得出的测试样本对应的平均散射角指:未进行极化相干矩阵修正时得出的测试样本对应的平均散射角(图2b中表示为“未修正”)。而本发明在图2b中表示为“修正后”。在图2b中,横坐标表示信噪比,单位为dB;纵坐标表示测试样本对应的平均散射角,单位为度。
参照图2c,为现有方法和本发明得出的测试样本的极化相干矩阵对应的散射熵随信噪比变化对比示意图。图2c中,现有方法得出的测试样本对应的散射熵指:未进行极化相干矩阵修正时得出的测试样本对应的散射熵(图2c中表示为“未修正”)。而本发明在图2c中表示为“修正后”。在图2c中,横坐标表示信噪比,单位为dB;纵坐标表示测试样本对应的散射熵。
从图2a至图2c可以看出,经过修正后的测试样本的极化相干矩阵的三个特征(最大特征值所占比例、平均散射角以及散射熵),相对未修正前随信噪比变化的幅度减小,证明本发明提高了上述三种特征的噪声稳健性。例如,图2a所示,最大特征值所占比例这一特征未修正时随信噪比变化曲线的方差为0.0202,修正后曲线方差降为0.0025;同样的,如图2b所示,经过修正,对应特征的方差由0.044降为0.011;如图2(c)所示,对应特征的方差由0.0247降为0.0041。说明经过噪声修正,上述三个特征随信噪比的变化明显减小。
参照图3,为现有方法和本发明得出的测试样本的极化相干矩阵的重构误差随信噪比变化的对比图。在图3中,现有方法对应于“未修正”,本发明对应于“修正后”。在图3中,横坐标为信噪比,单位为为dB;纵坐标表示测试样本的极化相干矩阵的重构误差。从图3可以看出,本发明得出的测试样本的极化相干矩阵的重构误差明显小于现有方法得出的测试样本的极化相干矩阵的重构误差,本发明得出的测试样本的极化相干矩阵的重构误差更加接近于测试样本无噪声时极化相干矩阵的重构误差,从而证明了本发明的噪声修正方法的有效性。
参照图4,为现有方法和本发明得出的测试样本对应的目标正确识别率随信噪比变化的对比示意图。图4中,现有方法对应于“未修正”,本发明对应于“修正后”。在图4中,横坐标为信噪比,单位为为dB;纵坐标表示测试样本对应的目标正确识别率。从图4中可以看出,对于0-20db的信噪比,相比现有方法,本发明得出的测试样本对应的目标正确识别率有显著的提高,识别率平均提高20%以上,尤其在15db信噪比时最为明显,识别率提高了超过30%。若以75%作为识别率门限,本发明对信噪比的要求相对现有方法可以降低约20dB,对应的识别距离可以推远约2.2倍。因此,本发明提出的基于极化分解的噪声稳健宽带全极化目标识别方法对于推远识别距离具有重要意义。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (4)
1.基于极化分解的噪声稳健宽带全极化目标识别方法,其特征在于,包括以下步骤:
S1:利用雷达接收多组全极化通道距离像信号;获取全极化通道距离像信号的多个训练样本;
S2:针对每个训练样本,采用目标极化特征提取方法,获取每个训练样本的平均距离像和分类特征;
所述步骤S2具体包括以下步骤:
S21:分别对每个训练样本的HH极化通道距离像信号、每个训练样本的HV极化通道距离像信号、每个训练样本的VH极化通道距离像信号、以及每个训练样本的VV极化通道距离像信号进行归一化处理,得到每个训练样本的HH极化通道归一化距离像信号、每个训练样本的HV极化通道归一化距离像信号、每个训练样本的VH极化通道归一化距离像信号、以及每个训练样本的VV极化通道归一化距离像信号;估计出每个训练样本的每个极化通道距离像信号的噪声平均功率;
在步骤S21中,对应训练样本的HH极化通道距离像信号表示为sHH,对应训练样本的HV极化通道距离像信号表示为sHV,对应训练样本的VH极化通道距离像信号表示为sVH,对应训练样本的VV极化通道距离像信号表示为sVV;其中,
sHH=[sHH1,...,sHHk,...,sHHM],sHV=[sHV1,...,sHVk,...,sHVM]
sVH=[sVH1,...,sVHk,...,sVHM],sVV=[sVV1,...,sVVk,...,sVVM]
其中,sHHk表示对应训练样本的HH极化通道第k个距离单元的距离像信号,sHVk表示对应训练样本的HV极化通道第k个距离单元的距离像信号,sVHk表示对应训练样本的VH极化通道第k个距离单元的距离像信号,sVVk表示对应训练样本的VV极化通道第k个距离单元的距离像信号,k取1至M,M为雷达接收的每个极化通道距离像信号的长度;
分别对对应训练样本的HH极化通道距离像信号sHH、对应训练样本的HV极化通道距离像信号sHV、对应训练样本的VH极化通道距离像信号sVH、以及对应训练样本的VV极化通道距离像信号sVV进行归一化处理,得到sHH经归一化处理后的信号sHV经归一化处理后的信号sVH经归一化处理后的信号以及sVV经归一化处理后的信号
其中,表示对应训练样本的HH极化通道第k个距离单元经归一化处理之后的距离像信号,表示对应训练样本的HV极化通道第k个距离单元经归一化处理之后的距离像信号,表示对应训练样本的VH极化通道第k个距离单元经归一化处理之后的距离像信号,表示对应训练样本的VV极化通道第k个距离单元经归一化处理之后的距离像信号;以及分别为:
其中,||·||2表示求2-范数;
然后,估计出对应训练样本的每个极化通道距离像信号的噪声平均功率,对应训练样本的每个极化通道距离像信号的噪声平均功率均为σ2;
S22:针对每个训练样本,在所有训练样本中选取与对应训练样本对应的多个训练样本,选取的每个训练样本的目标方位角与对应训练样本的目标方位角的角度差小于设定角度阈值;每个训练样本的目标方位角指:每个训练样本中目标相对于雷达的方位角;
在步骤S22中,根据以下公式得出对应训练样本的极化散射矩阵S:
S=[S1,…Sk,…,SM]
其中,Sk为:
根据以下公式得出对应训练样本的极化相干矩阵T:
T=[T1,…Tk,…,TM]
其中,|·|表示取模,上标*表示矩阵的共轭,在步骤S22中,对于选取的每个训练样本,均有对应的A值、B值和C值,<·>表示:根据步骤S22选取的多个训练样本求集合平均;根据步骤S22选取的多个训练样本,对或进行集合平均计算,得出对应训练样本的平均距离像 或
根据选取的每个训练样本的HH极化通道归一化距离像信号、选取的每个训练样本的HV极化通道归一化距离像信号、选取的每个训练样本的VH极化通道归一化距离像信号、以及选取的每个训练样本的VV极化通道归一化距离像信号,得出对应训练样本的极化相干矩阵T、以及对应训练样本的平均距离像
S23:根据对应训练样本的极化相干矩阵T、以及对应训练样本的每个极化通道距离像信号的噪声平均功率,得出对应训练样本的修正后极化相干矩阵根据对应训练样本的修正后极化相干矩阵得出对应训练样本的分类特征;
在步骤S23中,针对对应训练样本,分别对T1至TM进行修正,得出矩阵至然后利用矩阵至组合成对应训练样本的修正后极化相干矩阵对Tk进行修正包括以下步骤:
首先对Tk进行特征值分解,得出Tk=UkΛkUk -1,其中,上标-1表示求矩阵的逆,Λk为3×3维的矩阵,Λk为:
其中,λ1k为Tk的第1个特征值,λ2k为Tk的第2个特征值,λ3k为Tk的第3个特征值;Uk为3×3维的矩阵,Uk=[U1k,U2k,U3k],U1k为与λ1k对应的Tk的特征向量,U2k为与λ2k对应的Tk的特征向量,U3k为与λ3k对应的Tk的特征向量;U1k、U2k和U3k分别表示为如下形式:
其中,α1k、α2k、α3k、β1k、β2k、β3k、δ1k、δ2k、δ3k、γ1k、γ2k、以及γ3k均为角度参数;
根据对应训练样本的每个极化通道距离像信号的噪声平均功率σ2,对Λk进行修正,得到修正后特征值矩阵为:
则得出矩阵
在矩阵至中,针对每个矩阵提取归一化后比例特征、归一化后散射熵特征和归一化后平均散射角特征;在矩阵中提取归一化后比例特征、归一化后散射熵特征和归一化后平均散射角特征包括以下步骤:
根据以下公式矩阵的比例特征Pk、散射熵特征Hk、以及平均散射角特征
其中,g取1至3;由于k取1至M,此时便得出P1至PM、H1至HM、以及至计算出P1至PM的均值μP、P1至PM的标准差σP、H1至HM的均值μH、H1至HM的标准差σH、至的均值以及至的标准差
根据以下公式得出矩阵的归一化后比例特征归一化后散射熵特征和归一化后平均散射角特征
然后将矩阵至的归一化后比例特征组合成矢量将矩阵至的归一化后散射熵特征组合成矢量将矩阵至的归一化后平均散射角特征组合成矢量矢量矢量和矢量分别为:
将矢量矢量和矢量组合成对应训练样本的分类特征C,C为:
S3:将雷达接收的每组全极化通道距离像信号作为对应的测试样本,针对每个测试样本采用基于噪声修正的目标极化特征提取方法,获取每个测试样本的平均距离像和分类特征;针对每个测试样本,执行步骤S4至步骤S6;
在步骤S3中,获取每个测试训练样本的平均距离像和分类特征包括以下步骤:
S31:分别对每个测试样本的HH极化通道距离像信号、每个测试样本的HV极化通道距离像信号、每个测试样本的VH极化通道距离像信号、以及每个测试样本的VV极化通道距离像信号进行归一化处理,得到每个测试样本的HH极化通道归一化距离像信号、每个测试样本的HV极化通道归一化距离像信号、每个测试样本的VH极化通道归一化距离像信号、以及每个测试样本的VV极化通道归一化距离像信号;估计出每个测试样本的每个极化通道距离像信号的噪声平均功率;
S32:针对每个测试样本,在所有测试样本选取与对应测试样本对应的多个测试样本,选取的每个测试样本的雷达接收时刻与对应测试样本的雷达接收时刻的差小于设定时间长度;每个测试样本的雷达接收时刻指:雷达接收对应一组全极化通道距离像信号的时刻;
根据选取的每个测试样本的HH极化通道归一化距离像信号、选取的每个测试样本的HV极化通道归一化距离像信号、选取的每个测试样本的VH极化通道归一化距离像信号、以及选取的每个测试样本的VV极化通道归一化距离像信号,得出对应测试样本的极化相干矩阵T′、以及对应测试样本的平均距离像
S33:根据对应测试样本的极化相干矩阵T′、以及对应测试样本的每个极化通道距离像信号的噪声平均功率,得出对应测试样本的修正后极化相干矩阵根据对应测试样本的修正后极化相干矩阵得出对应测试样本的分类特征;
S4:将对应测试样本的平均距离像分别与每个训练样本的平均距离像进行对齐;得出对应测试样本的平均距离像在与每个训练样本的平均距离像进行对齐时,向右移动的距离单元数;
S5:根据对应测试样本的平均距离像在与每个训练样本的平均距离像进行对齐时,向右移动的距离单元数,得出每个训练样本的对应测试样本融合特征;
S6:根据每个训练样本的对应测试样本融合特征、以及每个训练样本的分类特征,得出对应测试样本的目标类别。
2.如权利要求1所述的基于极化分解的噪声稳健宽带全极化目标识别方法,其特征在于,在步骤S1中,所述全极化通道距离像信号包括HH极化通道距离像信号、HV极化通道距离像信号、VH极化通道距离像信号、以及VV极化通道距离像信号,
在获取全极化通道距离像信号的多个训练样本的同时,获取每个训练样本对应的目标类别。
3.如权利要求1所述的基于极化分解的噪声稳健宽带全极化目标识别方法,其特征在于,在步骤S31中,对应测试样本的HH极化通道距离像信号表示为xHH,对应测试样本的HV极化通道距离像信号表示为xHV,对应测试样本的VH极化通道距离像信号表示为xVH,对应测试样本的VV极化通道距离像信号xVV;其中,
xHH=[xHH1,...,xHHk,...,xHHM],xHV=[xHV1,...,xHVk,...,xHVM]
xVH=[xVH1,...,xVHk,...,xVHM],xVV=[xVV1,...,xVVk,...,xVVM]
其中,xHHk表示对应测试样本的HH极化通道第k个距离单元的距离像信号,xHVk表示对应测试样本的HV极化通道第k个距离单元的距离像信号,xVHk表示对应测试样本的VH极化通道第k个距离单元的距离像信号,xVVk表示对应测试样本的VV极化通道第k个距离单元的距离像信号,k取1至M,M为雷达接收的每个极化通道距离像信号的长度;
分别对对应测试样本的HH极化通道距离像信号xHH、对应测试样本的HV极化通道距离像信号xHV、对应测试样本的VH极化通道距离像信号xVH、以及测试样本的VV极化通道距离像信号xVV进行归一化处理,得到xHH经归一化处理后的信号xHV经归一化处理后的信号xVH经归一化处理后的信号以及xVV经归一化处理后的信号
其中,表示对应测试样本的HH极化通道第k个距离单元经归一化处理之后的距离像信号,表示对应测试样本的HV极化通道第k个距离单元经归一化处理之后的距离像信号,表示对应测试样本的VH极化通道第k个距离单元经归一化处理之后的距离像信号,表示对应测试样本的VV极化通道第k个距离单元经归一化处理之后的距离像信号;以及分别为:
其中,||·||2表示求2-范数;
然后,估计出对应测试样本的每个极化通道距离像信号的噪声平均功率,对应测试样本的每个极化通道距离像信号的噪声平均功率均为σ′2;
在步骤S32中,根据以下公式得出对应训练样本的极化散射矩阵S′:
S′=[S′1,…,S′k,…,S′M]
其中,S′k为:
根据以下公式得出测试样本的极化相干矩阵T′:
T′=[T′1,…,T′k,…,T′M]
其中,|·|表示取模,上标*表示矩阵的共轭,在步骤S32中,对于选取的每个测试样本,均有对应的A′值、B′值和C′值,<·>表示:根据步骤S32选取的多个测试样本求集合平均;根据步骤S32选取的多个测试样本,对或进行集合平均计算,得出对应测试样本的平均距离像 或
在步骤S33中,针对对应测试样本,分别对T′1至T′M进行修正,得出矩阵至然后利用矩阵至组合成对应测试样本的修正后极化相干矩阵对进行修正包括以下步骤:
首先对进行特征值分解,得出其中,上标-1表示求矩阵的逆,Λ′k为3×3维的矩阵,Λ′k为:
其中,λ′1k为的第1个特征值,λ′2k为的第2个特征值,λ′3k为的第3个特征值;Uk′为3×3维的矩阵,U′k=[U′1k,U′2k,U′3k],U′1k为与λ′1k对应的的特征向量,U′2k为与λ′2k对应的的特征向量,U′3k为与λ′3k对应的的特征向量;U′1k、U′2k和U′3k分别表示为如下形式:
其中,α′1k、α′2k、α′3k、β′1k、β′2k、β′3k、δ′1k、δ′2k、δ′3k、γ′1k、γ′2k、以及γ′3k均为角度参数;
根据对应测试样本的每个极化通道距离像信号的噪声平均功率σ′2,对Λ′k进行修正,得到修正后特征值矩阵为:
则得出矩阵
在矩阵至中,针对每个矩阵提取归一化后比例特征、归一化后散射熵特征和归一化后平均散射角特征;在矩阵中提取归一化后比例特征、归一化后散射熵特征和归一化后平均散射角特征包括以下步骤:
根据以下公式矩阵的比例特征P′k、散射熵特征H′k、以及平均散射角特征
其中g取1至3;由于k取1至M,此时便得出P′1至P′M、H′1至H′M、以及至计算出P′1至P′M的均值μ'P、P′1至P′M的标准差σ'P、H′1至H′M的均值μ'H、H′1至H′M的标准差σ'H、至的均值μ'α、以及至的标准差
根据以下公式得出矩阵的归一化后比例特征归一化后散射熵特征和归一化后平均散射角特征
然后将矩阵至的归一化后比例特征组合成矢量将矩阵至的归一化后散射熵特征组合成矢量将矩阵至的归一化后平均散射角特征组合成矢量矢量矢量和矢量分别为:
将矢量矢量和矢量组合成对应测试样本的分类特征C',C'为:
4.如权利要求3所述的基于极化分解的噪声稳健宽带全极化目标识别方法,其特征在于,在步骤S4中,在将对应测试样本的平均距离像分别与每个训练样本的平均距离像进行对齐时,使用相关系数最大准则;在所述对齐过程中,对应测试样本的平均距离像向右依次循环移位;
任一个训练样本的平均距离像表示为向量对应测试样本的平均距离像表示为向量则将向量按照相关系数最大准则与向量进行对齐;在向量与向量进行对齐的过程中,将向量向右依次循环移位,向量在每次移动时向右循环移动一个距离单元;
向量每次向右循环移动m个距离单元后,转变为向量m取1至M;计算向量与向量的相关系数rm:
其中,k取1至M,sk为对应训练样本的平均距离像中第k个距离单元的取值,s'(m)k为向量中第k个距离单元的取值,为对应训练样本的平均距离像中所有元素的均值,为向量中所有元素的均值;
得出对应测试样本的分类特征所需移动距离单元数N,N为:r1至rM中最大值对应的m的取值;
在步骤S5中,将矢量向右循环移动N个距离单元,得到新矢量将向右循环移动N个距离单元,得到新矢量将向右循环移动N个距离单元,得到新矢量对应测试样本融合特征表示为C”,
在步骤S6中,针对每个训练样本的对应测试样本融合特征、以及每个训练样本的分类特征,得出每个训练样本的欧式距离;得出每个训练样本的欧式距离包括以下步骤:根据以下公式得出任一个训练样本的欧式距离:
d=||C”-C||2
其中,C”表示与对应训练样本的对应的测试样本融合特征,C表示对应训练样本的分类特征;||·||2表示2-范数;
选取欧式距离最小的训练样本作为识别用训练样本,将识别用训练样本的目标类别作为测试样本的目标类别。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410216267.4A CN104007429B (zh) | 2014-05-21 | 2014-05-21 | 基于极化分解的噪声稳健宽带全极化目标识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410216267.4A CN104007429B (zh) | 2014-05-21 | 2014-05-21 | 基于极化分解的噪声稳健宽带全极化目标识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104007429A CN104007429A (zh) | 2014-08-27 |
CN104007429B true CN104007429B (zh) | 2017-05-24 |
Family
ID=51368173
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410216267.4A Active CN104007429B (zh) | 2014-05-21 | 2014-05-21 | 基于极化分解的噪声稳健宽带全极化目标识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104007429B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104280724B (zh) * | 2014-10-23 | 2016-08-17 | 西安电子科技大学 | 基于散射点匹配的雷达高分辨距离像识别方法 |
CN105044697B (zh) * | 2015-06-26 | 2017-07-18 | 西安电子科技大学 | 基于广义特征分解的全极化高分辨距离像目标检测方法 |
CN109298402B (zh) * | 2018-09-14 | 2022-12-06 | 西安电子工程研究所 | 基于通道融合的极化特征提取方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SG191420A1 (en) * | 2010-12-29 | 2013-07-31 | Nestec Sa | Fiber and probiotics for reducing intestinal symptoms related to chronic stress |
CN103413146B (zh) * | 2013-08-23 | 2017-03-29 | 西安电子科技大学 | 基于Freeman熵和自学习的极化SAR图像精细分类方法 |
-
2014
- 2014-05-21 CN CN201410216267.4A patent/CN104007429B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN104007429A (zh) | 2014-08-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105445708B (zh) | 一种极化合成孔径雷达的定标方法 | |
Gu et al. | Information-theoretic compressive sensing kernel optimization and Bayesian Cramér–Rao bound for time delay estimation | |
CN103149561B (zh) | 一种基于场景块稀疏的稀疏微波成像方法 | |
CN103969640B (zh) | 双基地mimo雷达目标稀疏成像方法 | |
CN105425223B (zh) | 广义帕累托杂波下稀疏距离扩展雷达目标的检测方法 | |
Gu et al. | Radar target profiling and recognition based on TSI-optimized compressive sensing kernel | |
CN104199029B (zh) | 一种提高压缩感知雷达目标成像性能的测量矩阵设计方法 | |
CN105807267A (zh) | 一种mimo雷达扩展目标的检测方法 | |
CN105044693A (zh) | 基于辅助阵元的微波关联成像雷达幅相误差校正方法 | |
CN104931931A (zh) | 互耦条件下基于张量实值子空间的双基地mimo雷达角度估计方法 | |
CN102495393B (zh) | 基于子空间追踪的压缩感知雷达成像算法 | |
CN103792521B (zh) | 一种受Faraday旋转变化影响的雷达目标后向散射模拟方法 | |
CN104749564A (zh) | 海杂波Weibull幅度分布参数的多分位点估计方法 | |
CN108717189A (zh) | 基于压缩感知理论的双基地mimo雷达成像方法 | |
Gu et al. | Single snapshot imaging method in multiple‐input multiple‐output radar with sparse antenna array | |
CN104007429B (zh) | 基于极化分解的噪声稳健宽带全极化目标识别方法 | |
CN107942322A (zh) | 多径利用雷达时间反演目标距离估计方法 | |
CN104155651A (zh) | 一种用于极化雷达目标跟踪的概率数据关联方法 | |
CN104215947A (zh) | 一种双基地mimo雷达角度的估计方法 | |
CN104020456A (zh) | 一种基于多特显点线阵成像雷达系统幅相误差校正方法 | |
CN105403871A (zh) | 一种双基地mimo雷达阵列目标角度估计和互耦误差校准方法 | |
CN103267963B (zh) | 基于相位调制的多准则雷达发射信号设计方法 | |
Gu et al. | Information-theoretic waveform design for Gaussian mixture radar target profiling | |
CN109799484A (zh) | 一种外辐射源雷达系统多径抑制方法、系统及存储介质 | |
JP5699404B2 (ja) | レーダ受信信号処理装置とその方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |