一种目标微动特征相干激光探测效果量化评价方法 |
|||||||
申请号 | CN201510997458.3 | 申请日 | 2015-12-24 | 公开(公告)号 | CN105629254A | 公开(公告)日 | 2016-06-01 |
申请人 | 中国人民解放军电子工程学院; | 发明人 | 胡以华; 李政; 郭力仁; 赵楠翔; 石亮; 王金诚; 王阳阳; 徐世龙; 董晓; | ||||
摘要 | 本 发明 提供一种目标微动特征相干激光探测效果量化评价方法,包括:采用平滑伪Winger-Ville分布 算法 对采集的目标微动特征相干激光探测回波 信号 进行时频分析处理,得到所述回波信号的时频分布矩阵;计算所述回波信号的频域聚集度;计算所述回波信号的分量分辨 力 ;计算所述回波信号的 信噪比 ;计算总指标。本发明给出了定量评价目标微 多普勒效应 探测及特征提取效果的方法,得到对目标微动特征探测效果规范、客观、量化的评价,解决了现行定性判断方法局限性大、规范性差的不足,为后续微多普勒效应探测影响因素、影响规律的研究及最优系统参数的选择奠定了 基础 。 | ||||||
权利要求 | 1.一种目标微动特征相干激光探测效果量化评价方法,其特征在于,包括以下步骤: |
||||||
说明书全文 | 一种目标微动特征相干激光探测效果量化评价方法技术领域[0001] 本发明涉及目标微多普勒特征探测技术领域,具体是一种目标微动特征相干激光探测效果量化评价方法。 背景技术[0002] 目标的精确探测和识别向来是各国国防研究的重点方向,然而,随着近年来隐身技术、诱饵技术等的发展和应用,传统的探测识别手段受到严重干扰,工作效能大打折扣,难以提供有效判断目标种类和真假的信息。21世纪初,美国海军研究实验室的V.C.Chen教授率先提出了微多普勒效应的概念,开辟了目标探测识别的新途径。 [0003] 目标与雷达存在相对径向运动会产生多普勒频移,而微多普勒效应就是指在此同时,由于目标自身或某些组成部分存在的微动,而对回波信号产生的附加频率调制,最终导致的以多普勒频移为中心的频率展宽效应。微多普勒效应被认为是反映目标微动的唯一特征,包含了目标的特有信息。使用激光雷达进行探测,其波长短的特点能够测得更为显著的微多普勒效应,从而获得更高的探测分辨率和灵敏度。 [0004] 实现目标微多普勒效应的探测,首先需要研究探测系统参数、结构、传输介质等因素对目标回波影响的规律,找到最佳的系统参数和结构的设计方案以及对外界干扰的补偿措施,这就要求在不同实验条件下对回波中目标微多普勒特征的探测效果进行对比分析。目前,在微多普勒特征提取中最常用的方法就是对信号进行时频分析,然后从时频分布中提取目标信息。系统探测效果的好坏会直观地反映到时频分布结果上,因此时频分布的可读性直接决定了微多普勒特征提取的准确性,这对于目标的探测和识别至关重要。 [0005] 目标微多普勒效应的回波相干探测信号为非平稳信号,所以在对其处理中,研究者们普遍采用了时频分析方法来提取目标微动特征(微多普勒特征),然而目前开展的对回波特征影响因素的研究中(比如频率源相位噪声、激光谱线宽度、探测距离、探测系统结构和参数等因素),在分析不同因素导致的不同目标微动特征探测效果时,研究人员主要还是通过定性的观测并根据经验来对比评价不同条件下得到的目标微动特征时频分布结果。 [0006] 这种定性的对比方法个体差异大,规范性差,有较大的局限性,比如:对于目标微动特征探测效果差异明显的两个时频分布结果,通过观察能够判断出哪种情况下得到的时频分布在提取目标微动特征时效果更好,但是难以具体给出效果好的比差的强多少,无法与之前定量变化的探测条件相对应;而对于探测效果差异不明显的两个时频分布结果,只通过观测则很难准确判断出哪种情况下的探测效果更好,甚至出现误判。为解决这一问题,我们需要科学、定量的评价方法来判断目标微多普勒特征的探测及提取效果。 发明内容[0007] 本发明的目的在于提供一种目标微动特征相干激光探测效果量化评价方法,基于平滑伪Winger-Ville分布(SPWVD)这种已被证明具有较好性能的时频分析方法,提出一种测量标准及其计算方法来定量评价目标微动特征的探测及提取效果,以解决定性分析的不足。 [0008] 本发明的技术方案为: [0009] 一种目标微动特征相干激光探测效果量化评价方法,包括以下步骤: [0011] (2)将所述回波信号的频域聚集度作为目标微动特征相干激光探测效果量化评价的第一个指标,并根据所述时频分布矩阵P(t,f)对所述频域聚集度进行计算; [0012] (3)将所述回波信号的分量分辨力作为目标微动特征相干激光探测效果量化评价的第二个指标,并根据所述时频分布矩阵P(t,f)对所述分量分辨力进行计算; [0013] (4)将所述回波信号的信噪比作为目标微动特征相干激光探测效果量化评价的第三个指标,并根据所述时频分布矩阵P(t,f)对所述信噪比进行计算; [0014] (5)根据所述频域聚集度、分量分辨力和信噪比,采用以下公式对目标微动特征相干激光探测效果量化评价的总指标进行计算: [0015] 若所述回波信号为单分量信号,则所述总指标的计算公式为: [0016] [0017] 若所述回波信号为多分量信号,则所述总指标的计算公式为: [0018] [0019] 其中,Q表示目标微动特征相干激光探测效果量化评价的总指标,a表示频域聚集度,r表示分量分辨力,SNR表示信噪比。 [0020] 所述的目标微动特征相干激光探测效果量化评价方法,步骤(2)中,根据所述时频分布矩阵P(t,f)对所述频域聚集度进行计算,包括以下步骤: [0021] (21)计算瞬时频域聚集度: [0022] 对于单分量信号,包括: [0023] a1、在时频分布矩阵P(t,f)中提取tn时刻的频率-能量分布一维矩阵P(tn,f); [0024] b1、找到P(tn,f)中的能量峰值ASmax及其对应的频率fm; [0025] c1、找到能量峰值ASmax两侧与fm距离最近的对应能量为 的两频率fp和fq; [0026] d1、求出tn时刻信号的3dB带宽BW(tn),即为tn时刻信号的瞬时频域聚集度a(tn): [0027] a(tn)=BW(tn)=fp-fq,fp>fq [0028] 对于多分量信号,包括: [0029] a2、在时频分布矩阵P(t,f)中提取tn时刻的频率-能量分布一维矩阵P(tn,f); [0030] b2、找到P(tn,f)中最大的m个能量峰值及其对应的频率,m为分量数目,m≥2; [0031] c2、依次求出m个能量峰值对应的m个3dB带宽BW1(tn)、BW2(tn)、…、BWm(tn),再求出m个3dB带宽的平均值,即为tn时刻信号的瞬时频域聚集度a(tn): [0032] [0033] (22)计算整体频域聚集度: [0034] 求出整个时间轴上所有时刻信号的瞬时频域聚集度a(tn),再对其求平均值,即得到信号的整体频域聚集度a: [0035] [0036] 其中,N表示整个时频分布中的离散时间点总数。 [0037] 所述的目标微动特征相干激光探测效果量化评价方法,步骤(3)中,根据所述时频分布矩阵P(t,f)对所述分量分辨力进行计算,包括以下步骤: [0038] (31)计算瞬时分量分辨力: [0039] 对于含有m个分量的多分量信号,m≥2,包括: [0040] a、针对其中两个分量,在时频分布矩阵P(t,f)中提取tn时刻的频率-能量分布一维矩阵P(tn,f); [0041] b、在P(tn,f)中找到最大的两个能量峰值ASmax1和ASmax2及其对应的频率f1、f2; [0042] c、求出能量峰值ASmax2对应的3dB带宽BW2(tn); [0043] d、根据以下公式求出两个分量在tn时刻的分辨力: [0044] [0045] e、按照步骤a~d,求出所有分量两两之间的瞬时分辨力rd(tn),再对所有rd(tn)求平均值,即得到信号的瞬时分量分辨力r(tn): [0046] [0047] 其中, 表示从m个分量中选2个得到的组合数; [0048] (32)计算整体分量分辨力: [0049] 求出整个时间轴上所有时刻信号的瞬时分量分辨力r(tn),再对其求平均值,即得到信号的整体分量分辨力r: [0050] [0051] 其中,N表示整个时频分布中的离散时间点总数。 [0052] 所述的目标微动特征相干激光探测效果量化评价方法,步骤(4)中,根据所述时频分布矩阵P(t,f)对所述信噪比进行计算,包括以下步骤: [0053] (41)计算瞬时信噪比: [0054] 对于单分量信号,包括: [0055] a1、在时频分布矩阵P(t,f)中提取tn时刻的频率-能量分布一维矩阵P(tn,f); [0056] b1、找到P(tn,f)中的能量峰值ASmax及其对应的频率; [0057] c1、求出tn时刻信号的3dB带宽BW(tn),找到3dB带宽内信号对应的频率坐标; [0058] d1、再求出P(tn,f)中3dB带宽以外所有频率坐标对应的能量平均值E(Anoise),作为噪声能量; [0059] e1、采用以下公式计算tn时刻信号的瞬时信噪比SNR(tn): [0060] [0061] 对于多分量信号,包括: [0062] a2、在时频分布矩阵P(t,f)中提取tn时刻的频率-能量分布一维矩阵P(tn,f); [0063] b2、找到P(tn,f)中最大的m个能量峰值ASmax1、ASmax2、…、ASmaxm及其对应的频率,m为分量数目,m≥2; [0064] c2、依次求出m个能量峰值对应的m个3dB带宽BW1(tn)、BW2(tn)、…、BWm(tn),找到m个3dB带宽内信号对应的频率坐标; [0065] d2、再求出P(tn,f)中m个3dB带宽以外所有频率坐标对应的能量平均值E(Anoise),作为噪声能量; [0066] e2、采用以下公式计算tn时刻信号的瞬时信噪比SNR(tn): [0067] [0068] (42)计算整体信噪比: [0069] 求出整个时间轴上所有时刻信号的瞬时信噪比SNR(tn),再对其求平均值,即得到信号的整体信噪比SNR: [0070] [0071] 其中,N表示整个时频分布中的离散时间点总数。 [0072] 本发明的有益效果为: [0073] 由上述技术方案可知,本发明给出了定量评价目标微多普勒效应探测及特征提取效果的方法,得到对目标微动特征探测效果规范、客观、量化的评价,解决了现行定性判断方法局限性大、规范性差的不足,为后续微多普勒效应探测影响因素、影响规律的研究及最优系统参数的选择奠定了基础。附图说明 [0074] 图1是本发明的评价标准定量计算流程图; [0075] 图2是本发明定义对单分量信号评价指标的求解方法,其中,图2(a)为时频分布,图2(b)为各指标计算参数说明; [0076] 图3是本发明定义对多分量信号评价指标的求解方法,其中,图3(a)为时频分布,图3(b)为各指标计算参数说明。 具体实施方式[0077] 下面结合附图和具体实施例进一步说明本发明。 [0078] 如图1所示,一种目标微动特征相干激光探测效果量化评价方法,包括以下步骤: [0079] S1、利用平滑伪Winger-Ville分布算法对采集的包含目标微多普勒效应的回波信号进行时频分析处理,得到含有目标微动特征信息的信号时频分布图和时频分布矩阵P(t,f)。 [0080] SPWVD的原理是在Winger-Ville分布基础上对频域和时域分别加窗进行平滑,来抑制交叉项的影响。由于SPWVD对交叉项的抑制效果较好,在时频分布结果中,交叉项基本不会对信号特征的提取带来影响,所以在考虑性能评价指标时就忽略了交叉项这个指标,但加窗平滑使能量分布展宽,时频聚集性有所下降,所以聚集度和分辨力这两个指标需要重点分析。 [0081] SPWVD的计算公式为: [0082] [0083] 其中,g(μ)是对频率加的窗函数,h(τ)是对时间加的窗函数,它们都是实的偶窗函* *数,且有g(0)=h(0)=1;s(·)s(·)表示回波信号及其共轭,所以s(·)s(·)表示某一过去时间信号与未来时间信号的自相关函数。 [0084] 计算得到的时频分布矩阵P(t,f)=SPWVDs(t,f),其横坐标表示时间,纵坐标表示频率,矩阵中的数值表示信号在时频面上的能量分布。 [0085] S2、计算目标微动特征相干激光探测效果量化评价的第一个指标——频域聚集度: [0086] 频域聚集度,其定义是信号能量在时频分布中的展宽程度,反映瞬时频率的不确定度。 [0087] S21、计算瞬时频域聚集度: [0088] 对于单分量信号,即只含有一个分量的信号,图2(a)给出了其SPWVD分布,其对应时频分布矩阵P(t,f),在P(t,f)中提取tn时刻时间剖面上的频率-能量分布一维矩阵P(tn,f),即图中垂直于时间轴的直线,具体分布如图2(b)所示,找到P(tn,f)中的能量峰值ASmax及其对应的频率fm;找到峰值两侧与fm距离最近的对应能量为 的两频率fp和fq,二者间距即为tn时刻信号的瞬时频域聚集度a(tn),相当于该时刻信号的3dB带宽BW(tn): [0089] a(tn)=BW(tn)=fp-fq,fp>fq [0090] 对于多分量信号,即含有多个分量的信号,其瞬时频域聚集度的求法在单分量信号的基础上进行扩展,在求峰值时由原来的求一个峰值改为求最大的前m个峰值(m为分量数目,m≥2),然后依次求m个峰值对应的3dB带宽BW1(tn)、BW2(tn)、…、BWm(tn),tn时刻信号的瞬时频域聚集度a(tn)为各个分量3dB带宽的平均值: [0091] [0092] 图3(a)给出了两分量信号的SPWVD分布,按照单分量信号的计算方法,求出图3(b)中的BW1(tn)、BW2(tn),此时信号的瞬时频域聚集度为: [0093] [0094] S22、计算整体频域聚集度: [0095] 按照步骤S21中描述的方法求出整个时间轴上所有时刻信号的瞬时频域聚集度a(tn),再对其求平均值,即得到信号的整体频域聚集度a: [0096] [0097] 其中,N表示整个时频分布中的离散时间点总数。 [0098] S3、计算目标微动特征相干激光探测效果量化评价的第二个指标——分量分辨力: [0099] 分量分辨力,是指多分量信号自项在频域的分辨能力,单分量信号不存在分量分辨的情况。 [0100] S31、计算瞬时分量分辨力: [0101] 以两分量信号为例,在图3(a)中提取tn时刻时间剖面上的频率-能量分布一维矩阵P(tn,f),即图中直线位置,剖面的能量分布如图3(b)所示,在P(tn,f)中找到最大的两个峰值ASmax1和ASmax2及峰值对应的频率值f1、f2,根据步骤S21所描述的方法,求出f2对应的分量在tn时刻的3dB带宽BW2(tn),再根据以下公式求出两个分量在tn时刻的分辨力: [0102] [0103] 当分量1的峰值与分量2的峰值距离小于分量2的3dB带宽,即rd(tn)<0时,认为两个分量不可分辨。 [0104] 对于多分量信号,先求出各个分量两两之间的瞬时分辨力rd(tn),再对所有rd(tn)求平均值,即得到信号的瞬时分量分辨力r(tn): [0105] [0106] 其中,m为分量数目,m≥2, 表示从m个分量中选2个得到的组合数。 [0107] S32、计算整体分量分辨力: [0108] 按照步骤S31中描述的方法求出整个时间轴上所有时刻信号的瞬时分量分辨力r(tn),再对其求平均值,即得到信号的整体分量分辨力r: [0109] [0110] 其中,N表示整个时频分布中的离散时间点总数。 [0111] S4、计算目标微动特征相干激光探测效果量化评价的第三个指标——信噪比: [0112] 信噪比,其定义为时频分布图中信号成分峰值能量与噪声成分平均能量之比,反映时频分布中信号成分的纯净度和清晰度,直接决定了微多普勒特征提取的精确程度。 [0113] S41、计算瞬时信噪比: [0114] 对于单分量信号,如图2(b)所示,在其时频分布矩阵P(t,f)中提取tn时刻时间剖面上的频率-能量分布一维矩阵P(tn,f),找到P(tn,f)中的能量峰值ASmax,根据步骤S21所描述的方法,求出tn时刻信号的3dB带宽BW(tn),找到3dB带宽内信号对应的频率坐标;再求出P(tn,f)中3dB带宽以外所有频率坐标对应的能量平均值E(Anoise),认为是噪声能量。采用以下公式计算tn时刻信号的瞬时信噪比SNR(tn): [0115] [0116] 对于含有m个分量(m≥2)的多分量信号,其瞬时信噪比定义为m个信号成分峰值能量的均值与噪声成分平均能量之比: [0117] [0118] 以两分量信号为例,如图3(b)所示,找到P(tn,f)中的能量峰值ASmax1和ASmax2,此时信号的瞬时信噪比为: [0119] [0120] S42、计算整体信噪比: [0121] 按照步骤S41中描述的方法求出整个时间轴上所有时刻信号的瞬时信噪比SNR(tn),再对其求平均值,即得到信号的整体信噪比SNR: [0122] [0123] 其中,N表示整个时频分布中的离散时间点总数。 [0124] S5、计算目标微动特征相干激光探测效果量化评价的总指标: [0125] 对于单分量信号,由于不存在对分量的分辨,所以总指标Q的计算公式为: [0126] [0127] 对于多分量信号,总指标Q的计算公式为: [0128] [0129] Q值越大,说明从时频分布结果中能够提取出目标微动特征的效果越好。 [0130] 对于研究某一参数α对目标微多普勒效应探测的影响,在α的可变化范围内依次求出其对应的Q值,最大Q值对应的就是目标微动特征探测中该参数的最佳取值。 |