CN107194329A - 一种基于自适应局部稀疏保持投影的一维距离像识别方法 - Google Patents
一种基于自适应局部稀疏保持投影的一维距离像识别方法 Download PDFInfo
- Publication number
- CN107194329A CN107194329A CN201710311830.XA CN201710311830A CN107194329A CN 107194329 A CN107194329 A CN 107194329A CN 201710311830 A CN201710311830 A CN 201710311830A CN 107194329 A CN107194329 A CN 107194329A
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- sparse
- dimensional
- matrix
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 59
- 230000003044 adaptive effect Effects 0.000 title claims abstract description 18
- 238000012549 training Methods 0.000 claims abstract description 72
- 239000011159 matrix material Substances 0.000 claims abstract description 51
- 238000012360 testing method Methods 0.000 claims abstract description 22
- 238000012706 support-vector machine Methods 0.000 claims abstract description 11
- 238000005259 measurement Methods 0.000 claims abstract description 5
- 239000013598 vector Substances 0.000 claims description 36
- 238000013519 translation Methods 0.000 claims description 22
- 230000008569 process Effects 0.000 claims description 10
- 238000004321 preservation Methods 0.000 claims description 2
- 230000014616 translation Effects 0.000 description 13
- 238000000605 extraction Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- 230000035945 sensitivity Effects 0.000 description 4
- 238000000513 principal component analysis Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 238000002592 echocardiography Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/08—Feature extraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/12—Classification; Matching
- G06F2218/16—Classification; Matching by matching signal segments
- G06F2218/20—Classification; Matching by matching signal segments by applying autoregressive analysis
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Image Analysis (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于自适应局部稀疏保持投影的一维距离像识别方法,该方法首先对实测一维距离像信号样本进行预处理;然后通过稀疏保持投影法得到稀疏系数矩阵,由局部保持投影法得到局部相似矩阵;接着将稀疏保持投影方程和局部保持投影方程及自适应最大间距准则进行融合,建立联合约束方程组,得到自适应局部稀疏保持投影矩阵;最后通过投影矩阵将训练样本和测试样本投影到低维空间,用支持向量机对其进行训练和分类。本发明基于稀疏保持投影、局部保持投影和自适应最大间距准则,充分利用样本稀疏重构、近邻关系中包含的识别信息结合自适应最大间距准则提取样本低维特征,提高了一维距离像信号的识别精度、降低了特征维度、增强了抗干扰性。
Description
技术领域
本发明涉及一种基于自适应局部稀疏保持投影的一维距离像识别方法,具体涉及一种在干扰环境下快速准确识别雷达目标一维距离像的技术,属于雷达一维信号识别技术领域。
背景技术
雷达自动目标识别是雷达信号处理领域的一个重要研究方向。随着雷达自动目标识别技术的广泛运用,人们对于雷达识别精度、实时性和抗干扰性的要求越来越高,而一般的雷达目标二维像如SAR、ISAR由于其自身成像机理的原因,在获取、储存和应用时数据量较大,处理耗时较长成为其实用过程中的一大障碍。雷达高分辨率一维距离像(HRRP)作为一种一维信号,由目标散射中心向雷达反射的回波所构成,其反映了目标的结构特征中的散射点中心强度、散射体中心位置、目标大小等信息,对目标识别具有重要价值,而由于其一维信号的特性,对储存和处理的要求相对较低,现已成为最具实际应用价值的雷达自动目标识别方案。
虽然HRRP相对雷达二维图像数据量较小,然而由于目标尺寸大小和雷达分辨率的影响,如今的雷达HRRP信号越来越精细,维度也越来越高;此外,由于其对方位角敏感,每个雷达目标需要获取许多不同方位的HRRP信号对其进行识别,因此研究HRRP信号低维且易于识别的特征至关重要。一直以来,国内外学者提出了多种特征提取方法,如传统的主成分分析法(PCA)、局部保持投影法(LPP)以及伴随着稀疏表示发展而提出的稀疏保持投影法(SPP)等。这些方法都利用了信号本身或者其空间结构特性,将其作为识别信息融入提取的低维特征中,虽然能在一定程度上提高识别率、降低特征维度,但其都未深入挖掘信号间的联系关系,对其包含的识别信息也利用有限,导致识别率提升有限且抗噪声能力不强,难以在实际环境中达到令人满意的识别效果。
发明内容
本发明所要解决的技术问题是:提供一种基于自适应局部稀疏保持投影的一维距离像识别方法,能够以低维特征获得满意的识别精度,该方法识别阶段耗时短且鲁棒性强,具有一定的实际应用价值。
本发明为解决上述技术问题采用以下技术方案:
一种基于自适应局部稀疏保持投影的一维距离像识别方法,包括如下步骤:
步骤1,获取雷达目标一维距离像并分成两部分,一部分作为训练样本集,另一部分作为测试样本集;
步骤2,对训练样本集X=[x1,x2,…,xN],提取其归一化幅度特征集并进行平移对齐,得到平移对齐后的训练样本归一化幅度特征集H=[h1,h2,…,hN],N表示训练样本总数;
步骤3,采用稀疏保持投影法对步骤2得到的平移对齐后的训练样本归一化幅度特征集H中各特征样本进行稀疏重构,建立稀疏重构方程,得到稀疏系数矩阵R;
步骤4,采用局部保持投影法对步骤2得到的平移对齐后的训练样本归一化幅度特征集H中各特征样本进行局部权重度量,得到局部相似矩阵S;
步骤5,结合稀疏系数矩阵R、局部相似矩阵S以及自适应最大间距准则建立联合约束方程组,并对方程组进行化简,利用拉格朗日乘数法求解得到自适应局部稀疏保持投影矩阵W;
步骤6,利用投影矩阵W对步骤2得到的平移对齐后的训练样本归一化幅度特征集H进行低维空间投影,得到训练样本低维特征向量集P=[p1,p2,…,pN];
步骤7,由训练样本低维特征向量集P对支持向量机进行训练,得到训练好的支持向量机;
步骤8,对测试样本y提取其归一化幅度特征并与步骤2得到的平移对齐后的训练样本归一化幅度特征集H进行平移对齐,得到平移对齐后的测试样本归一化幅度特征hy;
步骤9,利用步骤5求得的投影矩阵W对测试样本归一化幅度特征hy进行低维投影,得到测试样本低维特征向量py;
步骤10,利用训练好的支持向量机对测试样本低维特征向量进行分类,输出测试样本分类结果。
作为本发明的一种优选方案,所述步骤3的具体过程为:
根据稀疏保持投影法,对步骤2得到的平移对齐后的训练样本归一化幅度特征集H中任一特征样本hi(i=1,2,…,N),用除其自身以外的剩余特征样本对其进行稀疏重构,建立稀疏重构方程:
s.t.||Hri-hi||≤ε
1=eTri
其中,为满足约束条件的稀疏表示系数向量,N表示训练样本总数,ri=[ri,1,…,ri,i-1,0,ri,i+1,…,ri,N]T表示稀疏表示系数向量,ri,j表示训练样本特征向量hj对重构hi的贡献量(j=1,2,…,N,j≠i),||·||1表示取1范数,ε为噪声松弛量,e表示所有元素均为1的列向量;计算所有特征样本满足约束条件的稀疏表示系数向量得到稀疏系数矩阵
作为本发明的一种优选方案,所述步骤4的具体过程为:
采用局部保持投影法对步骤2得到的平移对齐后的训练样本归一化幅度特征集H中各特征样本进行局部权重度量,得到局部相似矩阵S:
其中,sij表示矩阵S的第i行第j列元素,||hi-hj||表示hi和hj的欧几里得距离,α为阈值且α>0,t为热核参数,i=1,2,…,N,j=1,2,…,N,N表示训练样本总数。
作为本发明的一种优选方案,所述步骤5的具体过程为:
结合稀疏系数矩阵R、局部相似矩阵S以及自适应最大间距准则建立联合约束方程组,则自适应局部稀疏保持投影矩阵W满足:
s.t.WTW=I
其中,Sβ=I-R-RT+RTR,H为平移对齐后的训练样本归一化幅度特征集,I表示单位矩阵,tr表示矩阵的迹,Sb和Sw分别表示样本类间散布矩阵和类内散布矩阵,Sb和Sw的表达式分别为:
其中,C表示训练样本目标类别总数,分别表示训练样本中属于类别t1、类别t2的样本归一化幅度特征均值向量,表示训练样本中属于类别t1的样本数目,表示属于类别t1的第k个样本归一化幅度特征,f(t1,t2),分别定义为:
其中,为间的欧氏距离,表示间的欧氏距离;由以上公式求得投影矩阵W,通过方程组化简得:
max tr[WT(Sb-Sw)W-WTH(Sβ-S)HTW]
s.t.WTW=I
利用拉格朗日乘数法,上述方程可转换成如下特征值问题的求解:
[(Sb-Sw)-H(Sβ-S)HT]w=λHHTw
则λ和w即为该问题的特征值和特征向量,求解以上方程所得到的最大的d个特征值所对应的特征向量即构成最终投影矩阵W=[w1,w2,…,wd]。
作为本发明的一种优选方案,所述步骤6的具体过程为:
利用投影矩阵W对步骤2得到的平移对齐后的训练样本归一化幅度特征集H进行低维空间投影:P=WTH,得到训练样本低维特征向量集P=[p1,p2,…,pN],N表示训练样本总数。
本发明采用以上技术方案与现有技术相比,具有以下技术效果:
1、本发明提出的识别方法利用稀疏保持投影和局部保持投影,保留信号样本邻域空间固有联系作为识别信息,结合自适应最大间距准则达到类间可分类内聚集,可以自适应调整各样本类别贡献权重,相比传统方案识别精度更高。
2、本发明提出的方法通过稀疏保持投影、局部保持投影及最大间距准则联合约束降维,在取得很好识别精度的同时避免了特征的大量冗余,达到降维的目的,在一定程度上提高了识别效率。
3、本发明提出的方法通过结合样本间稀疏表示关系且自适应调整贡献权重,在信噪比较低的情况下依然拥有足够的鉴别信息,对噪声抗性强,特征鲁棒性好。
4、本发明提出的一维距离像识别方法可以通过适当变换应用到多种一维或二维信号处理问题,如人脸识别、雷达SAR、ISAR图像的识别等,应用范围广泛。
附图说明
图1是本发明种基于自适应局部稀疏保持投影的一维距离像识别方法的流程图。
图2是本发明实施例处理的一维距离像原始信号的示意图。
图3是本发明实施例提取一维距离像归一化幅度特征的示意图。
图4是本发明由稀疏保持投影获得的实施例样本稀疏系数示意图。
图5是本发明方法在不同信噪比条件下取得的分类精度和其他方法的对比结果图。
具体实施方式
下面详细描述本发明的实施方式,所述实施方式的示例在附图中示出。下面通过参考附图描述的实施方式是示例性的,仅用于解释本发明,而不能解释为对本发明的限制。
本发明提出一种基于自适应局部稀疏保持投影的一维距离像识别方法实现低维特征提取达到雷达在干扰环境下的稳健识别。由于结合了稀疏保持投影、局域保持投影及自适应最大间距准则,充分挖掘信号内部结构信息,将其融入到低维特征提取过程中,本发明可在维数较低的情况下提取信号易于识别的特征,控制了计算量的同时提高了识别精度。后期通过线性支持向量机来达到分类目的,相比传统的低维特征提取方法达到了更好的识别精度,对噪声的鲁棒性也得到提高。具体流程如图1所示。
本发明结合实施例对具体流程进行阐述,现有一飞机一维距离像回波数据如图2所示,实际情况中不同型号的飞机回波各不相同,同一型号飞机不同角度的回波信号也有所差异,本发明主要目的就是完成这类回波信号的特征提取和类别判断。为了解决该类一维回波信号的识别分类问题,本发明方法包括以下几个步骤:
训练阶段:
第1步:对训练样本集X=[x1,x2,…,xN](xi∈Rl),l表示原始向量维数,这里所用一维距离像信号维度为256,提取其归一化幅度特征集如图3所示;
其中,|·|表示取模,||·||2表示取2范数,N表示训练样本总数;由于幅度特征的平移敏感性,需要对其进行平移对齐操作,这里采用较为常用的最大相关平移对齐法,信号和的互相关系数为:
式中,<·>表示内积运算;则每个信号以前一个幅度特征作为标准,后面依次平移直到与前一个幅度特征信号的互相关系数达到最大,设前一个信号幅度特征为则需平移t个距离单位,t满足:
由此得到平移对齐后的训练样本归一化幅度特征向量集H=[h1,h2,…,hN]。
第2步:基于第1步获得的训练样本幅度特征向量集H,对于其任一特征hi,用除其自身以外的剩余训练样本幅度特征对其进行稀疏重构,建立稀疏重构方程:
得到满足约束条件的稀疏表示系数向量式中e表示所有元素均为1的列向量,ri=[ri,1,…,ri,i-1,0,ri,i+1,…,ri,N]T表示稀疏表示系数向量,ri,j表示训练样本特征向量hj对重构hi的贡献量,ε为噪声松弛量,在本实施例中,ε取0.05,||·||1表示取1范数;计算所有训练样本稀疏表示系数向量得到稀疏系数矩阵图4给出了本实验中稀疏保持投影获得的示例样本稀疏系数示意图,可以发现对于一个训练样本特征,用剩余样本中很少的一部分便能对其进行精确稀疏重构,其中,稀疏系数作为该样本的重要识别信息将带入后续步骤融入最终低维特征提取当中。
第3步:基于第1步获得的训练样本幅度特征向量集H,由局部保持投影法,得到局部相似矩阵S:
这里sij代表矩阵S的第i行第j列元素,||hi-hj||代表hi和hj的欧几里得距离,α为一个极小的阈值且α>0,t是热核参数,由此可以得到局部相似矩阵S。
第4步:结合稀疏保持投影、局部保持投影方程和自适应最大间距准则(AMMC)建立联合约束方程组,则投影矩阵W应满足:
s.t.WTW=I
其中第一个约束条件来自SPP算法,且Sβ=I-R-RT+RTR,第二个约束条件来自LPP算法,第三个约束条件来自AMMC算法,I表示单位矩阵,tr表示矩阵的迹,Sb和Sw分别表示样本类间散布矩阵和类内散布矩阵,Sb和Sw的表达式分别为:
其中,C表示训练样本目标类别总数,表示训练样本中属于类别t1和类别t2的样本均值向量,表示训练样本中属于类别t1的样本数目,表示属于类别t1的第k个样本,这里,f(t1,t2),分别定义为:
其中,为间的欧氏距离,表示间的欧氏距离;由以上(6)(7)(8)可求得投影矩阵W,通过方程组合化简得:
利用拉格朗日乘数法,上述方程(9)可转换成如下特征值问题的求解:
[(Sb-Sw)-H(Sβ-S)HT]w=λHHTw (10)
则λ和w即为该问题的特征值和特征向量,求解以上方程所得到的最大的d个特征值所对应的特征向量即构成最终投影矩阵W=[w1,w2,…,wd],本实施例中取d=20。
第5步:由第4步求得的投影矩阵W对训练样本集进行低维空间投影:
P=WTH (11)
由此得到训练样本低维空间投影特征向量集P,将其作为最终特征进行分类器的训练。
第6步:由训练样本投影特征向量集P进行线性支持向量机的参数训练。
测试阶段:
第1步:对测试样本y提取其归一化幅度特征并与平移对齐后训练样本归一化幅度特征集H进行平移对齐,得到平移对齐后的幅度特征hy,方法如上训练阶段第1步所示。
第2步:利用训练阶段求得的投影矩阵W对测试样本幅度特征hy进行低维投影,得到测试样本低维特征向量py:
py=WThy (12)
第3步:利用训练完成的分类器对测试样本低维特征py进行分类识别,输出测试样本分类结果。
训练阶段第6步和测试阶段第3步中,选用线性支持向量机作为特征分类器;由于本发明所提取的低维特征已具有一定的线性可分性且线性支持向量机分类精度较高,参数训练时间较短,满足一维距离像分类过程中的实时性要求。
图5给出了本发明方法和其他主要特征降维算法在不同信噪比下的识别精度对比示意图(d为最终特征维数),由以上数据可以看出,本发明所提出的一维距离像识别方法相比传统的PCA、LPP、SPP等方法平均识别精度有4-12个百分点的提高,因步骤中利用了样本局部稀疏重构信息、类间类内信息在一定程度上丰富了特征的信息量,克服了单一LPP法对噪声和稀疏点的敏感性以及SPP对姿态的大变化的敏感性,结合两者优势,减小了噪声和姿态影响,使其在复杂环境中识别率更加稳定,其相比于直接利用原始信号256维幅度特征识别大幅降低了特征维度,缩短了识别时间,在实际工程应用中有很好的发展前景。
以上实施例仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明保护范围之内。
Claims (5)
1.一种基于自适应局部稀疏保持投影的一维距离像识别方法,其特征在于,包括如下步骤:
步骤1,获取雷达目标一维距离像并分成两部分,一部分作为训练样本集,另一部分作为测试样本集;
步骤2,对训练样本集X=[x1,x2,…,xN],提取其归一化幅度特征集并进行平移对齐,得到平移对齐后的训练样本归一化幅度特征集H=[h1,h2,…,hN],N表示训练样本总数;
步骤3,采用稀疏保持投影法对步骤2得到的平移对齐后的训练样本归一化幅度特征集H中各特征样本进行稀疏重构,建立稀疏重构方程,得到稀疏系数矩阵R;
步骤4,采用局部保持投影法对步骤2得到的平移对齐后的训练样本归一化幅度特征集H中各特征样本进行局部权重度量,得到局部相似矩阵S;
步骤5,结合稀疏系数矩阵R、局部相似矩阵S以及自适应最大间距准则建立联合约束方程组,并对方程组进行化简,利用拉格朗日乘数法求解得到自适应局部稀疏保持投影矩阵W;
步骤6,利用投影矩阵W对步骤2得到的平移对齐后的训练样本归一化幅度特征集H进行低维空间投影,得到训练样本低维特征向量集P=[p1,p2,…,pN];
步骤7,由训练样本低维特征向量集P对支持向量机进行训练,得到训练好的支持向量机;
步骤8,对测试样本y提取其归一化幅度特征并与步骤2得到的平移对齐后的训练样本归一化幅度特征集H进行平移对齐,得到平移对齐后的测试样本归一化幅度特征hy;
步骤9,利用步骤5求得的投影矩阵W对测试样本归一化幅度特征hy进行低维投影,得到测试样本低维特征向量py;
步骤10,利用训练好的支持向量机对测试样本低维特征向量进行分类,输出测试样本分类结果。
2.根据权利要求1所述基于自适应局部稀疏保持投影的一维距离像识别方法,其特征在于,所述步骤3的具体过程为:
根据稀疏保持投影法,对步骤2得到的平移对齐后的训练样本归一化幅度特征集H中任一特征样本hi(i=1,2,…,N),用除其自身以外的剩余特征样本对其进行稀疏重构,建立稀疏重构方程:
<mrow>
<msub>
<mover>
<mi>r</mi>
<mo>&OverBar;</mo>
</mover>
<mi>i</mi>
</msub>
<mo>=</mo>
<munder>
<mi>argmin</mi>
<msub>
<mi>r</mi>
<mi>i</mi>
</msub>
</munder>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>r</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<msub>
<mo>|</mo>
<mn>1</mn>
</msub>
</mrow>
s.t.||Hri-hi||≤ε
1=eTri
其中,为满足约束条件的稀疏表示系数向量,N表示训练样本总数,ri=[ri,1,…,ri,i-1,0,ri,i+1,…,ri,N]T表示稀疏表示系数向量,ri,j表示训练样本特征向量hj对重构hi的贡献量(j=1,2,…,N,j≠i),||·||1表示取1范数,ε为噪声松弛量,e表示所有元素均为1的列向量;计算所有特征样本满足约束条件的稀疏表示系数向量得到稀疏系数矩阵
3.根据权利要求1所述基于自适应局部稀疏保持投影的一维距离像识别方法,其特征在于,所述步骤4的具体过程为:
采用局部保持投影法对步骤2得到的平移对齐后的训练样本归一化幅度特征集H中各特征样本进行局部权重度量,得到局部相似矩阵S:
<mrow>
<msub>
<mi>s</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>h</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>h</mi>
<mi>j</mi>
</msub>
<mo>|</mo>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>/</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>h</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>h</mi>
<mi>j</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
<mo><</mo>
<mi>&alpha;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>h</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>h</mi>
<mi>j</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
<mo>&GreaterEqual;</mo>
<mi>&alpha;</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,sij表示矩阵S的第i行第j列元素,||hi-hj||表示hi和hj的欧几里得距离,α为阈值且α>0,t为热核参数,i=1,2,…,N,j=1,2,…,N,N表示训练样本总数。
4.根据权利要求1所述基于自适应局部稀疏保持投影的一维距离像识别方法,其特征在于,所述步骤5的具体过程为:
结合稀疏系数矩阵R、局部相似矩阵S以及自适应最大间距准则建立联合约束方程组,则自适应局部稀疏保持投影矩阵W满足:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>min</mi>
<mi> </mi>
<mi>t</mi>
<mi>r</mi>
<mo>&lsqb;</mo>
<msup>
<mi>W</mi>
<mi>T</mi>
</msup>
<msub>
<mi>HS</mi>
<mi>&beta;</mi>
</msub>
<msup>
<mi>H</mi>
<mi>T</mi>
</msup>
<mi>W</mi>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>max</mi>
<mi> </mi>
<mi>t</mi>
<mi>r</mi>
<mo>&lsqb;</mo>
<msup>
<mi>W</mi>
<mi>T</mi>
</msup>
<msup>
<mi>HSH</mi>
<mi>T</mi>
</msup>
<mi>W</mi>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>max</mi>
<mi> </mi>
<mi>t</mi>
<mi>r</mi>
<mo>&lsqb;</mo>
<msup>
<mi>W</mi>
<mi>T</mi>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>S</mi>
<mi>b</mi>
</msub>
<mo>-</mo>
<msub>
<mi>S</mi>
<mi>w</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>W</mi>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
s.t.WTW=I
其中,Sβ=I-R-RT+RTR,H为平移对齐后的训练样本归一化幅度特征集,I表示单位矩阵,tr表示矩阵的迹,Sb和Sw分别表示样本类间散布矩阵和类内散布矩阵,Sb和Sw的表达式分别为:
<mrow>
<msub>
<mi>S</mi>
<mi>b</mi>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>C</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<msub>
<mi>t</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>C</mi>
</munderover>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>t</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>m</mi>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
</msub>
<mo>-</mo>
<msub>
<mi>m</mi>
<msub>
<mi>t</mi>
<mn>2</mn>
</msub>
</msub>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>m</mi>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
</msub>
<mo>-</mo>
<msub>
<mi>m</mi>
<msub>
<mi>t</mi>
<mn>2</mn>
</msub>
</msub>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
</mrow>
<mrow>
<msub>
<mi>S</mi>
<mi>w</mi>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>C</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>c</mi>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
</msub>
</munderover>
<mi>g</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>h</mi>
<mrow>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>h</mi>
<mrow>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>m</mi>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
</msub>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>h</mi>
<mrow>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>m</mi>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
</msub>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
</mrow>
其中,C表示训练样本目标类别总数,分别表示训练样本中属于类别t1、类别t2的样本归一化幅度特征均值向量,表示训练样本中属于类别t1的样本数目,表示属于类别t1的第k个样本归一化幅度特征,f(t1,t2),分别定义为:
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>t</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>1</mn>
<mo>/</mo>
<mi>d</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>m</mi>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
</msub>
<mo>,</mo>
<msub>
<mi>m</mi>
<msub>
<mi>t</mi>
<mn>2</mn>
</msub>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>g</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>h</mi>
<mrow>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>d</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>h</mi>
<mrow>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>m</mi>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
</msub>
<mo>)</mo>
</mrow>
</mrow>
其中,为间的欧氏距离,表示间的欧氏距离;由以上公式求得投影矩阵W,通过方程组化简得:
max tr[WT(Sb-Sw)W-WTH(Sβ-S)HTW]
s.t.WTW=I
利用拉格朗日乘数法,上述方程可转换成如下特征值问题的求解:
[(Sb-Sw)-H(Sβ-S)HT]w=λHHTw
则λ和w即为该问题的特征值和特征向量,求解以上方程所得到的最大的d个特征值所对应的特征向量即构成最终投影矩阵W=[w1,w2,…,wd]。
5.根据权利要求1所述基于自适应局部稀疏保持投影的一维距离像识别方法,其特征在于,所述步骤6的具体过程为:
利用投影矩阵W对步骤2得到的平移对齐后的训练样本归一化幅度特征集H进行低维空间投影:P=WTH,得到训练样本低维特征向量集P=[p1,p2,…,pN],N表示训练样本总数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710311830.XA CN107194329B (zh) | 2017-05-05 | 2017-05-05 | 一种基于自适应局部稀疏保持投影的一维距离像识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710311830.XA CN107194329B (zh) | 2017-05-05 | 2017-05-05 | 一种基于自适应局部稀疏保持投影的一维距离像识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107194329A true CN107194329A (zh) | 2017-09-22 |
CN107194329B CN107194329B (zh) | 2020-12-08 |
Family
ID=59874034
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710311830.XA Active CN107194329B (zh) | 2017-05-05 | 2017-05-05 | 一种基于自适应局部稀疏保持投影的一维距离像识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107194329B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107909077A (zh) * | 2017-10-10 | 2018-04-13 | 安徽信息工程学院 | 半监督情况下基于稀疏化理论的特征选择方法 |
CN108122006A (zh) * | 2017-12-20 | 2018-06-05 | 南通大学 | 基于不等权局部保持嵌入的故障诊断方法 |
CN108647690A (zh) * | 2017-10-17 | 2018-10-12 | 南京工程学院 | 用于非约束人脸识别的判别稀疏保持投影方法 |
CN109212501A (zh) * | 2018-08-29 | 2019-01-15 | 杭州电子科技大学 | 基于局部嵌入的雷达高分辨距离像目标识别方法 |
CN109242010A (zh) * | 2018-08-27 | 2019-01-18 | 电子科技大学 | 一种稀疏学习rcs序列特征提取方法 |
CN110221267A (zh) * | 2019-06-17 | 2019-09-10 | 电子科技大学 | 一种目标库属性判别邻域保持学习子空间特征提取方法 |
CN113255745A (zh) * | 2021-05-13 | 2021-08-13 | 南京邮电大学 | 一种基于局部保留投影法的信道状态信息定位方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102902979A (zh) * | 2012-09-13 | 2013-01-30 | 电子科技大学 | 一种合成孔径雷达自动目标识别的方法 |
US8483492B2 (en) * | 2005-10-25 | 2013-07-09 | William Marsh Rice University | Method and apparatus for signal detection, classification and estimation from compressive measurements |
CN103941244A (zh) * | 2014-04-23 | 2014-07-23 | 电子科技大学 | 一种雷达目标一维距离像局部最优子空间识别方法 |
CN103729652B (zh) * | 2014-01-17 | 2016-08-31 | 重庆大学 | 基于稀疏保持流形嵌入的高光谱遥感影像分类方法 |
CN105975940A (zh) * | 2016-05-09 | 2016-09-28 | 大连大学 | 一种基于稀疏双向二维局部判别投影的掌纹图像识别方法 |
CN106257488A (zh) * | 2016-07-07 | 2016-12-28 | 电子科技大学 | 一种基于邻域特征空间鉴别分析的雷达目标识别方法 |
-
2017
- 2017-05-05 CN CN201710311830.XA patent/CN107194329B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8483492B2 (en) * | 2005-10-25 | 2013-07-09 | William Marsh Rice University | Method and apparatus for signal detection, classification and estimation from compressive measurements |
CN102902979A (zh) * | 2012-09-13 | 2013-01-30 | 电子科技大学 | 一种合成孔径雷达自动目标识别的方法 |
CN103729652B (zh) * | 2014-01-17 | 2016-08-31 | 重庆大学 | 基于稀疏保持流形嵌入的高光谱遥感影像分类方法 |
CN103941244A (zh) * | 2014-04-23 | 2014-07-23 | 电子科技大学 | 一种雷达目标一维距离像局部最优子空间识别方法 |
CN105975940A (zh) * | 2016-05-09 | 2016-09-28 | 大连大学 | 一种基于稀疏双向二维局部判别投影的掌纹图像识别方法 |
CN106257488A (zh) * | 2016-07-07 | 2016-12-28 | 电子科技大学 | 一种基于邻域特征空间鉴别分析的雷达目标识别方法 |
Non-Patent Citations (8)
Title |
---|
HUANHUAN ZHANG等: ""Adaptive neighborhood-preserving discriminant projection method for HRRP-based radar target recognition"", 《IEEE ANTENNAS WIRELESS PROPAGATION LETTERS》 * |
JIE GUI等: ""Locality preserving discriminant projection for face and palmprint recognition"", 《NEUROCOMPUTING》 * |
JIWEN LU等: ""Adaptive maximum margin criterion for image classification"", 《2011 IEEE INTERNATIONAL CONFERENCE ON MULTIMEDIA AND EXPO》 * |
JUN KONG等: ""A Novel Approach for Face Recognition Based on Supervised Locality Preserving Projection and Maximum Margin Criterion"", 《2009 INTERNATIONAL CONFERENCE ON COMPUTER ENGINEERING AND TECHNOLOGY》 * |
KUN CHEN等: ""The target HRRP recognition with supervised discrimination sparse neighborhood preserving embedding"", 《INTERNATIONAL JOURNAL OF ADVANCED RESEARCH IN SCIENCE,ENGINEERING AND TECHNOLOGY》 * |
WANG YI等: ""A Novel Feature Extraction Algorithm of Acoustic Targets Based on Locality Preserving Discriminant Projections"", 《2010 INTERNATIONAL CONFERENCE ON COMPUTER APPLICATION AND SYSTEM MODELING (ICCASM 2010)》 * |
YUE JIANG等: ""Target recognition of radar HRRP using manifold learning with feature weighting"", 《2016 IEEE INTERNATIONAL WORKSHOP ON ELECTROMAGNETICS: APPLICATIONS AND STUDENT INNOVATION COMPETITION (IWEM)》 * |
张欢欢: ""目标宽带电磁散射分析与识别技术研究"", 《中国博士学位论文全文数据库 信息科技辑》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107909077A (zh) * | 2017-10-10 | 2018-04-13 | 安徽信息工程学院 | 半监督情况下基于稀疏化理论的特征选择方法 |
CN108647690A (zh) * | 2017-10-17 | 2018-10-12 | 南京工程学院 | 用于非约束人脸识别的判别稀疏保持投影方法 |
CN108647690B (zh) * | 2017-10-17 | 2021-12-31 | 南京工程学院 | 基于判别稀疏保持投影的非约束人脸图像降维方法 |
CN108122006A (zh) * | 2017-12-20 | 2018-06-05 | 南通大学 | 基于不等权局部保持嵌入的故障诊断方法 |
CN109242010A (zh) * | 2018-08-27 | 2019-01-18 | 电子科技大学 | 一种稀疏学习rcs序列特征提取方法 |
CN109212501A (zh) * | 2018-08-29 | 2019-01-15 | 杭州电子科技大学 | 基于局部嵌入的雷达高分辨距离像目标识别方法 |
CN109212501B (zh) * | 2018-08-29 | 2020-11-06 | 杭州电子科技大学 | 基于局部嵌入的雷达高分辨距离像目标识别方法 |
CN110221267A (zh) * | 2019-06-17 | 2019-09-10 | 电子科技大学 | 一种目标库属性判别邻域保持学习子空间特征提取方法 |
CN110221267B (zh) * | 2019-06-17 | 2022-08-05 | 电子科技大学 | 一种目标库属性判别邻域保持学习子空间特征提取方法 |
CN113255745A (zh) * | 2021-05-13 | 2021-08-13 | 南京邮电大学 | 一种基于局部保留投影法的信道状态信息定位方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107194329B (zh) | 2020-12-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107194329B (zh) | 一种基于自适应局部稀疏保持投影的一维距离像识别方法 | |
CN108133232B (zh) | 一种基于统计字典学习的雷达高分辨距离像目标识别方法 | |
CN107085206B (zh) | 一种基于自适应稀疏保持投影的一维距离像识别方法 | |
CN101908138B (zh) | 基于噪声独立成分分析的合成孔径雷达图像目标识别方法 | |
CN107133648B (zh) | 基于自适应多尺度融合稀疏保持投影的一维距离像识别方法 | |
CN107992891B (zh) | 基于光谱矢量分析多光谱遥感图像变化检测方法 | |
CN107132535A (zh) | 基于变分贝叶斯学习算法的isar稀疏频带成像方法 | |
CN106951915B (zh) | 一种基于类别置信度的一维距离像多分类器融合识别法 | |
CN103955701B (zh) | 多层次结合的多视合成孔径雷达图像目标识别方法 | |
CN107729926B (zh) | 一种基于高维空间变换的数据扩增方法、机器识别系统 | |
CN105809198A (zh) | 基于深度置信网络的sar图像目标识别方法 | |
CN104122554B (zh) | 高分辨率sar图像目标的属性散射中心特征提取方法 | |
CN102262736B (zh) | 空间目标图像分类与识别方法 | |
CN105741333B (zh) | 一种Video-SAR图像实时压缩重构方法 | |
CN106096517A (zh) | 一种基于低秩矩阵与特征脸的人脸识别方法 | |
CN106203256A (zh) | 一种基于稀疏保持典型相关分析的低分辨率人脸识别方法 | |
CN107831473B (zh) | 基于高斯过程回归的距离-瞬时多普勒图像序列降噪方法 | |
CN107219510B (zh) | 基于无限最大间隔线性判别投影模型的雷达目标识别方法 | |
CN103218623B (zh) | 基于自适应近邻保持鉴别投影的雷达目标特征提取方法 | |
CN109001702B (zh) | 无载波超宽带雷达人体动作识别方法 | |
CN104680182B (zh) | 基于nsct和判别字典学习的极化sar分类方法 | |
CN108734115A (zh) | 一种基于标签一致字典学习的雷达目标识别方法 | |
CN105116400A (zh) | 基于iMMFA模型的雷达高分辨距离像的目标识别方法 | |
CN109343043A (zh) | 一种基于选择主成分分析的雷达hrrp目标识别方法 | |
CN107045120B (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 | ||
TR01 | Transfer of patent right |
Effective date of registration: 20220107 Address after: 210000 Room 215, block B, Zhongyue Industrial Park, No. 129, Wutang village, Gulou District, Nanjing, Jiangsu Province Patentee after: NANJING YUDADIAN INFORMATION TECHNOLOGY CO.,LTD. Address before: No. 29, Qinhuai District, Qinhuai District, Nanjing, Jiangsu Patentee before: Nanjing University of Aeronautics and Astronautics |
|
TR01 | Transfer of patent right |