CN101598784B - 基于fa模型的强噪声背景下雷达距离像统计识别方法 - Google Patents
基于fa模型的强噪声背景下雷达距离像统计识别方法 Download PDFInfo
- Publication number
- CN101598784B CN101598784B CN2009100232474A CN200910023247A CN101598784B CN 101598784 B CN101598784 B CN 101598784B CN 2009100232474 A CN2009100232474 A CN 2009100232474A CN 200910023247 A CN200910023247 A CN 200910023247A CN 101598784 B CN101598784 B CN 101598784B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msubsup
- mover
- overbar
- msub
- 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 74
- 238000012360 testing method Methods 0.000 claims abstract description 112
- 238000012549 training Methods 0.000 claims abstract description 31
- 239000011159 matrix material Substances 0.000 claims description 40
- 238000010606 normalization Methods 0.000 claims description 19
- 238000004364 calculation method Methods 0.000 claims description 15
- 230000008569 process Effects 0.000 claims description 13
- 238000013519 translation Methods 0.000 claims description 11
- 238000004422 calculation algorithm Methods 0.000 claims description 6
- 238000010998 test method Methods 0.000 claims description 3
- 239000013078 crystal Substances 0.000 claims 1
- 230000008901 benefit Effects 0.000 abstract description 4
- 238000009432 framing Methods 0.000 abstract description 2
- 101001048529 Plectranthus scutellarioides Hydroxyphenylpyruvate reductase Proteins 0.000 abstract 1
- 239000000523 sample Substances 0.000 description 57
- 230000009466 transformation Effects 0.000 description 12
- 230000035945 sensitivity Effects 0.000 description 11
- 230000006870 function Effects 0.000 description 6
- 238000002474 experimental method Methods 0.000 description 5
- 238000013179 statistical model Methods 0.000 description 5
- 238000000513 principal component analysis Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 230000007423 decrease Effects 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 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 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000000556 factor analysis Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000013074 reference sample Substances 0.000 description 1
- 238000010845 search algorithm Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000011426 transformation method Methods 0.000 description 1
- 238000005303 weighing Methods 0.000 description 1
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于FA模型的强噪声背景下雷达距离像统计识别方法,它涉及雷达自动目标识别技术领域,主要解决现有FA模型统计识别方法对噪声不稳健的问题。其训练阶段步骤为:对雷达连续HRRP分帧、平移对齐和强度归一化,并利用处理后的HRRP学习FA模型各方位帧的参数,并保存模板。其测试阶段步骤为:先对待测试样本强度归一化、平移对齐,然后估计它的信噪比范围,若该信噪比大于30dB,则计算各目标各帧的距离值,并判定类别属性;若该信噪比小于30dB,则改写现有距离值,通过最小化它,求解低信噪比条件下的噪声能量,最后计算各目标各帧的距离值,并判定类别属性。本发明具有对噪声稳健,计算量较小的优点,用于对雷达目标识别。
Description
技术领域
本发明属于雷达自动目标识别技术领域,具体的说是一种基于FA模型的强噪声背景下雷达距离像统计识别方法。
背景技术
一维高分辨距离像的自动目标技术可追溯到20世纪80年代。由于一维高分辨距离像HRRP能够提供目标沿距离方向几何结构信息,且具有易于获得和处理的独特优势;同时,雷达具有全天候、全天时的特点,雷达一维高分辨距离像自动目标识别受到广泛重视。
基于统计模型的识别方法是一种重要的雷达HRRP自动目标识别方法。一维高分辨距离像的统计识别是指根据测试样本在各类别下的类后验概率的大小确定该测试样本的类别归属的识别方法。相关文献中提出了利用多种统计建模进行目标识别的方法,这些模型包括独立高斯模型、独立Gamma模型,基于Gamma和混合高斯Gaussian Mixture的独立双分布模型、基于主分量分析PCA联合高斯模型、概率主分量分析PPCA联合高斯模型和因子分析FA联合高斯模型等。其中FA模型假设距离像回波数据x服从联合高斯分布,距离像x与隐变量y的关系为x=Ay+m+ε,其中x和y的维数分别为D和d,d<D,m是均值向量,A是加载矩阵。A的各列为正交向量,且ε服从零均值、协方差矩阵ψ为对角阵且对角线元素不相同的高斯分布:ε~N(0,ψ),其中ID为D×D单位矩阵;y服从零均值、协方差矩阵Id为对角阵且对角线元素相同的高斯分布:y~N(0,Id),其中Id为d×d单位矩阵。从而,类条件概率密度函数可表示为:
利用FA模型进行目标识别的步骤如下:
(一)训练过程
(1)按照目标所在的方位将目标的所有HRRP回波数据划分成多个数据段,每段称为一帧;
(2)将各帧内的HRRP回波数据平移对齐;
(3)将各帧内所有平移对齐后的HRRP数据进行强度归一化;
(4)分别对各个帧内平移对齐以及强度归一化后的数据建立一个FA模型,求取模型参数mjk,ψjk,Ajk并将其保存为模板 j=1,2,…,C,C是总的目标类别数,k=1,2,…,Kj,Kj是第j类目标的总的方位帧数目。
(二)测试过程
(1)对需要测试的HRRP回波数据,即测试样本,强度归一化,得到样本xtest;
(2)将xtest分别与各类目标模板中的均值向量mjk平移对齐,得到样本xtest jk*,j=1,2,…,C,k=1,2,…,Kj;
(3)利用训练过程得到的模板 计算样本xtest jk*对应于各类目标各帧的类条件概率密度:
(4)找出步骤(3)中计算出的类条件概率密度中最大的一个,若该类条件概率密度对应的模板属于第j类目标,j=1,2,…,C,则判定测试样本属于第j类目标,测试过程结束;
虽然传统FA模型考虑了噪声分量的统计建模,但是,实际战场环境复杂,雷达回波中噪声强度受目标距离雷达远近、目标特定方位反射特性、大气条件等影响,通常难于保持测试样本的信噪比与训练样本完全一致,尤其对在战场环境下的远距离的非合作目标,HRRP信噪比更低。因此,实际测试样本与训练样本的噪声分量是失配的;另外,由于距离像样本经过能量归一化预处理来克服强度敏感性,测试样本与训练样本的信号分量也是失配的,因而直接用传统的FA模型来对低信噪比条件下的样本进行识别,将会造成识别率下降,尤其当噪声较大时,识别率更低。
发明内容
本发明目的是克服上述已有技术的不足,提供一种基于FA模型的强噪声背景下雷达距离像统计识别方法,根据测试样本调节FA模型的模板参数mjk,ψjk,Ajk,使得FA模型的模板参数与测试样本相匹配,以实现在强噪声背景下能保持较高的识别率。
实现本发明目的的技术思路是:在训练阶段使统计模型的信号部分与噪声部分开;在测试阶段通过最小化距离值实现自适应学习,并将距离值改写成多项式的和,通过最小化各项,近似求出最小化的距离值,利用它判定测试样本所属的目标类别。具体实现过程如下:
A.训练步骤:
(A1)按照目标所在的方位将目标的所有高信噪比环境下获得的一维高分距离像HRRP回波数据划分成多个数据段,每段称为一帧;
(A2)将各帧内的HRRP回波数据平移对齐;
(A3)将各帧内所有平移对齐后的HRRP数据进行强度归一化;
(A4)分别对各个帧内强度归一化后的数据建立一个FA模型,求取模型参数均值mjk +,噪声协方差阵ψjk +,加载矩阵Ajk +,利用Jacobi算法联合对角化ψjk+和Ajk +TAjk +分别为QjkΛψ,jk +Qjk T和QjkΛA,jk +Qjk T,其中Qjk是列正交矩阵,Λψ,jk +和ΛA,jk +是对角矩阵,并保存为模板 j=1,2,…,C,C是总的目标类别数,k=1,2,…,Kj,Kj是第j类目标的总的方位帧数目;
B.测试步骤
(B1)对需要测试的样本进行强度归一化,得到归一化后的测试样本xtest;
(B2)将归一化后的测试样本xtest分别与各类目标模板中的均值向量mjk +平移对齐,得到对齐后的测试样本xtest jk*,=1,2,…,C,k=1,2,…,Kj;
(B3)估计待测试样本的信噪比范围,对于信噪比大于30dB的测试样本执行步骤(B4)至(B5),对于信噪比小于30dB的测试样本执行步骤(B6)至(B11);
(B4)利用训练过程得到的FA模板,计算测试样本对应于各类目标所有帧的距离值:
其中,mjk +,噪声协方差阵ψjk +,加载矩阵Ajk +为FA模型训练阶段存储模板参数;
(B5)找出步骤(B4)计算出的距离值中最小的一个,若该距离值对应的模板属于第j类目标,j=1,2,…,C,则判定测试样本属于第j类目标,测试过程结束;
(B6)将步骤(B4)中各类目标所有帧的距离值改写为:
其中,trace(ψjk +)表示求取矩阵ψjk +的迹,qp,jk是正交矩阵Qjk的列向量,Ps,jk是j类目标的第k帧的信号能量,Pw,jk +,Pw,jk -分别是j类目标的第k帧高、低信噪比条件下的噪声能量,D是距离像的维数,d是隐变量的维数,λA,p,jk +分别是训练阶段存储的模板参数Λψ,jk +和ΛA,jk +的第p个元素;
(B7)对步骤(B6)改写后的各帧距离值关于Pw,jk -求导数,并令之为零,最终等价为求解D个一元三次方程:
p=1,…,D
(B8)利用卡丹公式求解步骤(B7)中的一元三次方程,并通过判断方程的根与常用信噪比范围的关系,得到步骤(B6)中的距离值Ffa jk(xtest jk*)的近似最小值对应的参数z,把它记为zp,jk *,然后利用关系式 求得第j类目标第k帧的第p个距离单元对应的噪声能量: p=1,…,D;
(B9)求距离值近似最小值对应的低信噪比环境下的噪声能量:
(B10)计算测试样本到各类目标所有帧的距离值:
(B11)从步骤(B10)计算出的各类目标所有帧的距离值中找到最小的一个,若该距离值对应的模板属于第j类目标,j=1,2,…,C,则判定测试样本属于第j类目标,测试过程结束。
本发明与现有技术相比具有以下优点:
1.现有的FA统计模型没有考虑噪声环境对识别性能影响,如图2所示,当环境噪声变换时,识别率迅速下降;本发明对噪声稳健,可将识别的信噪比范围提高10dB以上。
2.本发明由于在测试过程中改写距离值时考虑了能量归一化预处理的影响,从而使模型对数据的描述准确度提高。
3.本发明由于在测试过程中通过改写距离值,并利用卡丹公式计算测试样本的噪声能量,所以与利用简单的一维迭代搜索算法来求测试样本噪声能量相比,计算量小,复杂度低,进而有更强的可实现性。
附图说明
图1是本发明的流程图;
图2本发明与现有基于FA统计模型的识别率随噪声变化仿真曲线比较;
图3在最优幂次变换下本发明与现有基于FA统计模型的识别方法的识别率随噪声变化仿真曲线比较。
具体实施方式
参照图1,本发明的统计识别方法包括训练和测试两个阶段,具体步骤如下:
步骤1,对雷达连续HRRP分帧。
将高信噪比环境下雷达录取到的一维高分辨距离像作为训练数据。按照目标所在的方位将目标的所有训练数据划分成等间隔的多个数据段,把每一段称为一帧;并按顺序标号存储帧样本。
步骤2,对各帧内的HRRP平移对齐。
由于目标沿雷达射线方向的平动会引起HRRP平移量的变化,所以会使发生平动的同一目标的HRRP变成两个相似性很小的样本,这对目标识别是不利的,称之为平移敏感性。对于分帧后的数据,各帧内样本具有平移敏感性。采用包络对齐方法来克服平移敏感性,具体方法为选取帧内某一样本作为基准,用xsd表示该样本,然后滑动帧内其它样本使它们与xsd的相关系数分别最大,得到对齐后的其它样本。存储基准样本和对齐后的其它样本作为平移对齐后的样本集合。
克服所述平移敏感性的方法除了采用包络对齐方法,也可以采用零相位绝对对齐方法、平移强度联合匹配方法、最大后验概率方法,但不限于这些方法。
步骤3,对各帧内的HRRP强度归一化。
由于目标距离雷达远近、系统损耗、电磁散射损耗等原因,雷达录取到的HRRP在强度上存在差异,这会给识别带来困难,称之为强度敏感性。平移对齐后的HRRP具有强度敏感性。强度敏感性可利用能量归一化准则来克服,具体而言,第j个目标第k个方位帧的第i个HRRP样本xi,jk可以归一化为:
从而使所有样本的能量都为1。
克服所述强度敏感性的方法除了采用能量归一化方法,也可以采用能量相等方法、模一归一化方法、平移强度联合匹配方法,但不限于这些方法。
步骤4,学习FA模型各方位帧的参数,并保存模板。
(4.1)建立FA模型。
设强度归一化后的样本xi,jk对应的隐变量和噪声变量分别是yi,jk,εi,jk。根据现有FA模型有:
其中,j=1,2…,C,C是目标总类别数目,k=1,2…,Njk,Njk是第j个目标第k个方位帧的样本数目,xi,jk,和yi,jk的维数分别为D和d,且d<D;上标+号表示高信噪比条件,mjk +是高信噪比条件下第j个目标第k个方位帧的均值向量,Ajk +是高信噪比条件下第j个目标第k个方位帧的加载矩阵;yi,jk服从均值为零向量协方差矩阵为Id的高斯分布,用符号表示为yi,jk~N(0,Id),其中Id是d×d单位矩阵,N(·,·)表示高斯分布;εi,jk服从均值为零向量协方差矩阵为ψjk +的高斯分布,用符号表示为εi,jk~N(0,ψjk +),其中ID是D×D单位矩阵。
(4.2)得到FA模型概率密度函数。
由于yi,jk~N(0,Id)和εi,jk~N(0,ψjk +),所以样本xi,jk服从高斯分布,它的均值和协方差矩阵分别为mjk +和ψjk ++Ajk +Ajk +T,用符号表示为xi,jk~N(mjk +,ψjk ++Ajk +Ajk +T),则类条件概率密度函数表示为:
(4.3)利用最大似然方法和EM算法求取模型参数。
由(3)式得到似然函数 最大化似然函数,
得到高信噪比条件下FA的模型参数mjk +:
利用经典的EM算法,求得FA的模型参数ψjk +,Ajk +。
(4.4)联合对角化ψjk +和Ajk +TAjk +。
利用Jacobi算法联合对角化ψjk +和Ajk +TAjk +分别为QjkΛψ,jk +Qjk T和QjkΛA,jk +Qjk T,即
其中Qjk是列正交矩阵,Λψ,jk +和ΛA,jk +是对角矩阵。
(4.5)存储模板。
将(4.4)得到的模型参数mjk +,Λψ,jk +,,ΛA,jk +,Qjk,Ajk +和ψjk +按照目标类别和帧序号存储为模板,模板用符号记为 j=1,2,…,C,C是总的目标类别数,
k=1,2,…,Kj,Kj是第j类目标的总的方位帧数目;
通过以上步骤1至步骤4完成本发明的训练阶段。
步骤5,对测试样本强度归一化。
由于测试样本与训练样本存在强度差异,需要对测试样本进行归一化,该归一化采用与训练阶段一致的准则,即:
其中xtest是雷达录取到得原始测试样本,xtest是归一化后的样本。
克服测试样本强度敏感性,除了能量归一化方法,也可以采用模一归一化方法、平移强度联合匹配方法,但不限于这些方法。需要注意的是,克服测试样本强度敏感性的方法要与训练阶段一致。
步骤6,测试样本与均值向量平移对齐。
由于测试样本与训练样本存在平移量的差异,需要对测试样本进行平移。该平移采用最大相关系数方法,具体为:将归一化后的测试样本xtest分别与各类目标模板中的均值向量mjk +滑动相关,相关系数最大时称为对齐,对齐后测试样本xtest jk*的平移量为:
xtest τ表示将xtest平移τ位,τ=1,2,...,D,克服平移敏感性,除了最大相关系数方法,也可以采用零相位绝对对齐方法、平移强度联合匹配方法、最大后验概率方法,但不限于这些方法。
步骤7,采用常规的雷达信噪比估计方法估计测试样本的信噪比范围,对于信噪比大于30dB的测试样本执行步骤8至步骤9,对于信噪比小于30dB的测试样本执行步骤10至步骤13。
步骤8,计算各目标各帧的距离值。
利用训练过程存储的FA模板 j=1,2,…,C,k=1,2,…,Kj,根据(3)式,得到测试样本对各目标各帧的类条件概率密度如下:
由于类条件概率密度Pjk(xtest jk*)计算复杂,采用计算其相对应的距离值:
利用(10)式计算出的距离值代替其相应的类条件概率密度进行识别。
步骤9,判定类别属性。
找出步骤8计算出的距离值中最小的一个,若该距离值对应的模板属于第j类目标,j=1,2,…,C,则判定测试样本属于第j类目标,测试过程结束。
步骤10,改写距离值。
(10.1)求FA模型参数在高、低信噪比条件下的关系
设同一无噪信号,在高信噪比环境下利用FA模型学习得到的均值向量为mjk +、加载矩阵为Ajk +,噪声协方差矩阵为ψjk +,低信噪比环境下学习得到的均值向量为mjk -、加载矩阵为Ajk -,噪声协方差矩阵为ψjk -,且mjk +与mjk -、Ajk +与Ajk -、Ψjk +与ψjk -形状一样,但存在尺度上的差异,令 根据(4)式,高、低信噪比条件下的mjk +与mjk -的关系式为:
其中,Ps,jk是j类目标的第k帧的信号能量,Pw,jk +,Pw,jk -分别是j类目标的第k帧高、低信噪比条件下的噪声能量。
根据(11)式,得到
由于没有强度归一化之前信号能量是恒定的,所以有如下等式:
将(11)式和 式代入(13)式,得到比例系数αjk 2为:
把(14)式代入 得到Ajk +与Ajk -的关系式如下:
由于强度归一化,所以有如下等式:
代入 以及(14)式、(16)式,得到
把(17)式代入 得到Ψjk +与Ψjk -的关系式如下:
(10.2)利用(10.1)中得到的高、低信噪比条件下模型参数的关系改写(10)式中的距离值为:
步骤11,最小化(19)式,求解低信噪比条件下的噪声能量。
(11.1)计算信号能量Ps,jk和高信噪比条件下的噪声能量Pw,jk +
信号能量Ps,jk和高信噪比条件下的噪声能量Pw,jk +的计算可通过下面两式计算得到:
其中,trace(Ajk +Ajk +T+mjk +mjk +T)表示对Ajk +Ajk +T+mjk +mjk +T求迹运算,trace(ψjk +)表示对ψjk +求迹运算。
(11.2)沿距离单元展开(19)式的距离值Ffa jk(xtest jk*):
式进一步可展开为
其中,第p个距离单元对应的FA距离定义为
由于雷达HRRP数据的各个距离单元的FA距离ffa p,jk增减形状大致相同,它们的和的最小值与ffa p,jk最小值对应的低信噪比条件下的噪声能量Pw,jk -是很接近的,所以最小化(19)式近似为最小化(24)式,而最小化(24)式可通过对其求导数,令导数为零,求得一元三次方程并求该方程的根来实现。
(11.3)计算与求解(24)式最小值等价的一元三次方程
对(24)式求导:
令 则有 把它代入(25)式,
得
令(26)式分子为零,得到一组三元一次方程:
求解(27)式的根,就能最小化(24)式。
(11.4)根据卡丹公式求(27)式的根
(27)式是一元三次方程,可用卡丹公式计算方程的根。首先给出求解三次方程的卡丹公式以及对应的根与系数的关系如下:
(a)当Δ>0时,方程有一个实根x1和两个互为共轭的复数根x2和x3;
(b)当Δ=0时,方程有一个实根x1和两个相等实数根x2和x3,x2=x3;
(c)当Δ<0时,方程有三个实根x1、x2和x3,且x1≠x2≠x3。
由(27)式知,一元三次方程系数为,
将上面的(27)式方程的系数代入卡丹公式,按照(28)式、(29)式和(30)式求出(27)式的三个跟z1,z2,z3。
(11.5)求(27)式的ffa p,jk在常用信噪比范围对应的区间 内最小值。
ffa p,jk在 上是连续函数,其一阶导数为令f(z)=0,根据(11.4)求得此方程的三个跟为z1,z2,z3,则fppcap,jk的最小值对应的点zp,jk *必为z1,z2z3, 之一,具体步骤如下:
(11.5a)当时,根据根与系数的关系得到其中z2和z3是两个互为共轭的复数或两个相等的实数。因为所以(z-z2)(z-z3)≥0,从而,当z<z1时,f(z)≤0,fppca p,jk随z是单调非增变化的;当z>z1时,f(z)≥0,fppca p,jk随z是单调非减变化的。不考虑取值范围,z1是一个极小值点。故此时ffa p,jk的最小值对应的点必在z1, 中。
当时,ffa p,jk最小值对应的点的计算公式为:
(11.5b)当时,根据根与系数的关系,有三个实根z1、z2和z3,且z1≠z2≠z3,由三次方根性质容易知道,z1>z2>z3,在不考虑取值范围时,z2是极大值点,z1和z3是极小值点。故此时,ffa p,jk的最小值对应的点必在z1,z3, 中。
(11.5c)当时,根据根与系数的关系,z2和z3是两个互为共轭的复数或两个相等的实数,则有(z-z2)(z-z3)≥0。当z<z1时,f(z)≥0,单调非减;当z>z1时,f(z)≤0,单调非增,不考虑取值范围,z1是一个极大值点。故此时ffa p,jk的最小值对应的点必为 之一。
(11.5d)当时,根据根与系数的关系,有三个实跟z1、z2和z3,且z1≠z2≠z3,由三次方根性质容易知道,z1>z2>z3,且不考虑取值范围时,z2是极小值点,z1和z3是极大值点。故此时,ffa p,jk的最小值对应的点必在z2, 中。
(11.6)求ffa p,jk的最小值对应的噪声能量
根据 则有
pw,(p),jk -*是第j类目标第k帧的第p维对应的噪声能量,p=d+1,…,D。
(11.7)求低信噪比环境下的噪声能量:
利用第d+1到D维对应的噪声能量的均值,近似求得低信噪比环境下的噪声能量:
低信噪比环境下的噪声能量进一步可采用 或 或 来近似求得。
步骤12,计算距离值。
根据(26)和(40)式,得到计算各个目标各帧的距离值的公式:
步骤13,判定类别属性。
从步骤14计算出的各类目标所有帧的距离值中找到最小的一个,若该距离值对应的模板属于第j类目标,j=1,2,…,C,则判定测试样本属于第j类目标,测试过程结束;
通过步骤5至步骤13完成本发明的测试阶段。
本发明的效果可以通过以下仿真实验说明:
试验一:所用实测数据包含三类目标:雅克42,安26和奖状飞机。训练样本大致包含全方位数据,近似认为是完备数据。由于本发明主要克服噪声对识别性能的影响,将雅克42训练数据分为35帧,安26数据分为50帧,奖状数据分为50帧。每帧训练样本数为1024。HRRP含有256个距离单元。
定义ISAR实测数据的各帧的平均信噪比为40dB,按照5dB间隔分别仿真了-10dB到36dB平均信噪比的训练数据、测试数据。利用现有FA模型,计算了不同噪声强度下的识别率并用虚线画在图2中。由于训练样本和噪声样本的噪声水平是一致的,故称这种方法为“匹配噪声”识别法。这种方法在实际操作是难以实现的,我们将它作为一种较理想的情况画出,用来比较衡量本发明。由图2可见,“匹配噪声”识别精度随信噪比下降而下降,到-10dB时,已几乎失效,因此,我们取噪声搜索范围为[Pw,jk +,10Ps,jk],这里Pw,jk +,Ps,jk表示利用训练数据根据(27)式计算得到的高信噪比条件下的噪声能量和根据(28)式计算得到的信号能量。[Pw,jk +,10Ps,jk]对应信噪比范围是[-10dB,r+],其中r+表示训练数据的信噪比。
用高信噪比条件下的模板按照现有方法测试低信噪比样本的识别方法被成为“失配噪声”识别方法,由图2可见,“失配噪声”方法随着噪声加大,识别率急剧下降;本发明方法相比于“失配噪声”方法能大幅提高识别率,若以正确识别率70%为可用前提,本发明能使可用噪声范围扩大约10dB多。
表1给出现有方法和本的样本匹配单个模板的计算时间,隐空间变量分别取50,100和150维,结果表示为平均时间±标准差,其实验所用软件为MATLABR2007b,所用电脑硬件配置为3.0-GHz Pentium-4处理器、2-GB内存和Windows XP操作系统。
表1各种方法的测试样本与单个模板匹配平均时间和标准差
常规算法 | 本发明 | |
10维 | 0.0005±8×10-5秒 | 0.0037±5×10-5秒 |
20维 | 0.0005±5×10-5秒 | 0.0037±3×10-5秒 |
40维 | 0.0004±8×10-5秒 | 0.0037±2×10-5秒 |
由表1可见,本发明有效提高了计算效率。
试验二:对试验一的数据做最优幂变换,重复试验一。
从识别率方面比较:对于现有FA方法,幂次变换能有效地抬高识别率,我们在图3中画出最优幂变换指数为0.2的下识别率随信噪比的变化关系。常见雷达HRRP信噪比范围大约为10dB~25dB,由图3知道,本发明方法获得了略微优于现有方法的识别性能。表2和表3分别给出信噪比10dB、15dB、20dB和25dB下幂变换前后本发明方法的识别率:
表2在不同信噪比下,本发明的平均识别率和混淆矩阵
表3不同信噪比下,最优幂变换指数为0.2时,本发明的平均识别率和混淆矩阵
可见,幂变换之后识别性能反而下降。
从计算速度方面比较:表4给出在最优幂变换指数为0.2时,现有方法和本发明的测试样本与单个模板匹配的平均时间和标准差,隐空间变量分别取10,20和40维,表格中形式为平均时间±标准差,其实验所用软件为MATLABR2007b,所用电脑硬件配置为3.0-GHz Pentium-4处理器、2-GB内存和Windows XP操作系统。
表4最优幂变换后,现有方法与本发明的计算时间
现有方法 | 本发明 | |
10维 | 0.0005±5×10-5秒 | 0.0037±4×10-5秒 |
20维 | 0.0005±8×10-5秒 | 0.0037±5×10-5秒 |
40维 | 0.0005±7×10-5秒 | 0.0037±5×10-5秒 |
比较表4与表1可见,幂变换占用更多时间资源。
试验二说明从识别率和计算速度来看,传统的幂变换方法对低信噪比统计识别是不利的,而本发明中由于没有采用幂变换,不仅能够节约计算时间,而且在低信噪比条件下的识别性能优于采用幂变换的方法。
Claims (5)
1.一种基于FA模型的强噪声背景下雷达距离像统计识别方法,包括
A.训练步骤:
(A1)按照目标所在的方位将目标的所有高信噪比环境下获得的一维高分距离像HRRP回波数据划分成多个数据段,每段称为一帧;
(A2)将各帧内的HRRP回波数据平移对齐;
(A3)将各帧内所有平移对齐后的HRRP数据进行强度归一化;
(A4)分别对各个帧内强度归一化后的数据建立一个FA模型,求取模型参数均值噪声协方差阵加载矩阵利用Jacobi算法联合对角化和分别为和其中Qjk是列正交矩阵,和是对角矩阵,并保存为模板 j=1,2,…,C,C是总的目标类别数,k=1,2,…,Kj,Kj是第j类目标的总的方位帧数目;
B.测试步骤
(B1)对需要测试的样本进行强度归一化,得到归一化后的测试样本
(B3)估计待测试样本的信噪比范围,对于信噪比大于30dB的测试样本执行步骤(B4)至(B5),对于信噪比小于30dB的测试样本执行步骤(B6)至(B11);
(B4)利用训练过程得到的FA模板,计算测试样本对应于各类目标所有帧的距离值:
(B5)找出步骤(B4)计算出的距离值中最小的一个,若该距离值对应的模板属于第j类目标,j=1,2,…,C,则判定测试样本属于第j类目标,测试过程结束;
(B6)将步骤(B4)中各类目标所有帧的距离值改写为:
其中,表示求取矩阵的迹,qp,jk是正交矩阵Qjk的列向量,Ps,jk是j类目标的第k帧的信号能量,分别是j类目标的第k帧高、低信噪比条件下的噪声能量,D是距离像的维数,d是隐变量的维数,分别是训练阶段存储的模板参数和的第p个元素;
p=1,…,D
(B8)利用卡丹公式求解步骤(B7)中的一元三次方程,并通过判断方程的根与常用信噪比范围的关系,得到步骤(B6)中的距离值的近似最小值对应的参数z,把它记为然后利用关系式求得第j类目标第k帧的第p个距离单元对应的噪声能量:p=1,…,D;
(B9)求距离值近似最小值对应的低信噪比环境下的噪声能量:
(B10)计算测试样本到各类目标所有帧的距离值:
(B11)从步骤(B10)计算出的各类目标所有帧的距离值中找到最小的一个,若该距离值对应的模板属于第j类目标,j=1,2,…,C,则判定测试样本属于第j类目标,测试过程结束。
2.根据权利要求1所述的基于FA模型的强噪声背景下雷达距离像统计识别方法,其中步骤(B9)求距离值近似最小值对应的低信噪比环境下的噪声能量,采用如下公式计算:
3.根据权利要求1所述的基于FA模型的强噪声背景下雷达距离像统计识别方法,其中步骤(B9)求距离值近似最小值对应的低信噪比环境下的噪声能量,采用如下公式计算:
4.根据权利要求1所述的基于FA模型的强噪声背景下雷达距离像统计识别方法,其中步骤(B9)求距离值近似最小值对应的低信噪比环境下的噪声能量,采用如下公式计算:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100232474A CN101598784B (zh) | 2009-07-08 | 2009-07-08 | 基于fa模型的强噪声背景下雷达距离像统计识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100232474A CN101598784B (zh) | 2009-07-08 | 2009-07-08 | 基于fa模型的强噪声背景下雷达距离像统计识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101598784A CN101598784A (zh) | 2009-12-09 |
CN101598784B true CN101598784B (zh) | 2012-06-27 |
Family
ID=41420288
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009100232474A Active CN101598784B (zh) | 2009-07-08 | 2009-07-08 | 基于fa模型的强噪声背景下雷达距离像统计识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101598784B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107817490A (zh) * | 2017-10-09 | 2018-03-20 | 南京航空航天大学 | 一种帧内最大相似像的雷达一维距离像帧中心提取方法 |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102184408B (zh) * | 2011-04-11 | 2012-12-26 | 西安电子科技大学 | 基于自回归模型的雷达高分辨距离像目标识别方法 |
CN102636777B (zh) * | 2012-04-29 | 2013-08-14 | 西安电子科技大学 | 基于噪声先验的概率主成分分析雷达目标稳健识别方法 |
CN103235296B (zh) * | 2013-05-05 | 2014-09-17 | 西安电子科技大学 | 基于功率谱特征修正的噪声背景下雷达目标识别方法 |
CN105116400B (zh) * | 2015-08-31 | 2018-04-10 | 西安电子科技大学 | 基于iMMFA模型的雷达高分辨距离像的目标识别方法 |
CN109490840B (zh) * | 2018-11-22 | 2022-05-31 | 中国人民解放军海军航空大学 | 基于改进稀疏自编码模型的雷达目标hrrp的降噪和重构方法 |
CN111308440A (zh) * | 2020-03-17 | 2020-06-19 | 上海无线电设备研究所 | 超宽带雷达目标检测方法、电子设备和存储介质 |
CN111693139B (zh) * | 2020-06-19 | 2022-04-22 | 浙江讯飞智能科技有限公司 | 声音强度测量方法、装置、设备及存储介质 |
CN112882010B (zh) * | 2021-01-12 | 2022-04-05 | 西安电子科技大学 | 基于信噪比领域知识网络的高分辨距离像目标识别方法 |
-
2009
- 2009-07-08 CN CN2009100232474A patent/CN101598784B/zh active Active
Non-Patent Citations (10)
Title |
---|
保铮 |
保铮.雷达高分辨距离像目标识别研究进展.《电子与信息学报》.2005,第27卷(第8期),1328-1334. |
候庆禹.一种利用强度信息的雷达HRRP自动目标识别方法.《电子学报》.2009,第27卷(第3期),459-463. |
刘宏伟 |
刘宏伟;杜兰;袁莉;保铮.雷达高分辨距离像目标识别研究进展.《电子与信息学报》.2005,第27卷(第8期),1328-1334. * |
杜兰 |
杜兰.雷达高分辨距离像目标识别方法研究.《中国博士学位论文全文数据库》.2008,(第4期),I136-41. * |
袁莉 |
陈凤 |
陈凤;杜兰;刘宏伟;保铮;候庆禹.一种利用强度信息的雷达HRRP自动目标识别方法.《电子学报》.2009,第27卷(第3期),459-463. * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107817490A (zh) * | 2017-10-09 | 2018-03-20 | 南京航空航天大学 | 一种帧内最大相似像的雷达一维距离像帧中心提取方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101598784A (zh) | 2009-12-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101598784B (zh) | 基于fa模型的强噪声背景下雷达距离像统计识别方法 | |
CN108133232B (zh) | 一种基于统计字典学习的雷达高分辨距离像目标识别方法 | |
CN110084221B (zh) | 一种基于深度学习的带中继监督的序列化人脸关键点检测方法 | |
CN102175999B (zh) | 基于谱图特征的雷达目标高分辨距离像识别方法 | |
CN101598783B (zh) | 基于ppca模型的强噪声背景下雷达距离像统计识别方法 | |
CN102169177B (zh) | 基于时域特征的雷达目标高分辨距离像识别方法 | |
CN102628938B (zh) | 基于噪声先验的联合高斯模型雷达目标稳健识别方法 | |
CN110045015B (zh) | 一种基于深度学习的混凝土结构内部缺陷检测方法 | |
CN111079847B (zh) | 一种基于深度学习的遥感影像自动标注方法 | |
CN106054155B (zh) | 基于卷积因子分析模型的雷达高分辨距离像目标识别方法 | |
CN104280724B (zh) | 基于散射点匹配的雷达高分辨距离像识别方法 | |
CN101853509B (zh) | 基于Treelets和模糊C-均值聚类的SAR图像分割方法 | |
CN111507047B (zh) | 一种基于SP-CUnet的逆散射成像方法 | |
CN114726692B (zh) | 一种基于SEResNet-LSTM的辐射源调制方式识别方法 | |
CN107085206A (zh) | 一种基于自适应稀疏保持投影的一维距离像识别方法 | |
US20230243744A1 (en) | Method and system for automatically detecting and reconstructing spectrum peaks in near infrared spectrum analysis of tea | |
CN111581909A (zh) | 基于改进的自适应重要性采样算法的sram良率评估方法 | |
CN103454623B (zh) | 基于复数agc模型的噪声背景下雷达目标识别方法 | |
CN106951822B (zh) | 一种基于多尺度稀疏保持投影一维距离像融合识别方法 | |
CN102636777B (zh) | 基于噪声先验的概率主成分分析雷达目标稳健识别方法 | |
Chung et al. | Recursive K-distribution parameter estimation | |
CN113240081B (zh) | 针对雷达载频变换的高分辨距离像目标稳健识别方法 | |
CN115932770A (zh) | 雷达辐射源个体精准智能识别方法、系统、设备及终端 | |
CN115294398A (zh) | 一种基于多姿态角联合学习的sar图像目标识别方法 | |
CN110941908B (zh) | 一种基于核密度估计的海杂波分布建模方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C41 | Transfer of patent application or patent right or utility model | ||
TR01 | Transfer of patent right |
Effective date of registration: 20151023 Address after: 710075, B, building 38, Tengfei innovation center, six hi tech Zone, Xi'an hi tech Zone, Shaanxi, 606 Patentee after: XI'AN CETC XIDIAN UNIVERSITY RADAR TECHNOLOGY COLLABORATIVE INNOVATION RESEARCH INSTITUTE CO., LTD. Address before: Xi'an City, Shaanxi province Taibai Road 710071 No. 2 Patentee before: Xidian University |