一种ISAR回波信号的稳健时频分析方法技术

本发明专利技术涉及一种ISAR回波信号的稳健时频分析方法。该方法包括：获取ISAR回波信号；根据所述ISAR回波信号确定所述ISAR回波信号的时间片信号集；确定所述时间片信号集中每一时间片信号的频域表示基；根据所述频域表示基构造时频过完备卷积框架；根据所述时频过完备卷积框架以及所述ISAR回波信号建立联合优化模型；根据所述联合优化模型，利用迭代阈值算法进行求解，确定所述ISAR回波信号的时频谱。本发明专利技术有效提升了ISAR回波信号的时频分辨率。

【技术实现步骤摘要】
一种ISAR回波信号的稳健时频分析方法
本专利技术涉及信号分析处理领域，特别是涉及一种ISAR回波信号的稳健时频分析方法。
技术介绍
时间与频率是描述ISAR回波信号最重要的两个物理量，时频分析通过研究信号频率随时间的变化规律，提供了ISAR回波信号时间与频率的联合分布。经过多年的研究，传统时频分析已经在ISAR回波信号分析中获得广泛的应用。常用的时频分析方法主要分为三类：线性时频分析、双线性时频分析方法以及数据驱动的时频分析方法。线性时频分析主要包含短时傅立叶变换、连续小波变换、Gabor变换、S变换等，此类时频分析算法简单易于实现，但时频分辨率受到限制，对于ISAR回波信号中的强非线性、突变等自适应能力不足。双线性时频分析主要包含魏格纳分布、平滑伪魏格纳分布、科恩类时频分布等，此类时频分析显著提升信号的时频分辨率，但对于ISAR回波信号中的多分量成分存在交叉项影响，使得频谱图上存在假频、模糊等现象。数据驱动的时频分析方法主要包含经验模态分解、变分模态分解等方法，其能够很好的处理非平稳信号，克服交叉项干扰等不利因素，但是此类方法理论基础尚不完备，在ISAR回波信号时频分析中的稳定性不足。此外，一些ISAR回波信号的实际应用条件，如信号部分缺失、信噪比过低等，都会对其时频分析都会产生严重影响。稳健时频分析通过挖掘ISAR回波信号中时频谱的先验信息提升时频分析的性能与稳健性，其主要涉及两个关键因素：时频表示基与先验约束。传统时频表示基不论调和分析类还是过完备字典，都难以兼顾表示性能与表示效率；而分析模型又往往针对多个时间片信号分别施加先验约束，稳健性有待提高。
技术实现思路
本专利技术的目的是提供一种ISAR回波信号的稳健时频分析方法，能够有效提升ISAR回波信号的时频分辨率。为实现上述目的，本专利技术提供了如下方案：一种ISAR回波信号的稳健时频分析方法，包括：获取ISAR回波信号；根据所述ISAR回波信号确定所述ISAR回波信号的时间片信号集；确定所述时间片信号集中每一时间片信号的频域表示基；根据所述频域表示基构造时频过完备卷积框架；根据所述时频过完备卷积框架以及所述ISAR回波信号建立联合优化模型；根据所述联合优化模型，利用迭代阈值算法进行求解，确定所述ISAR回波信号的时频谱。可选的，所述根据所述频域表示基构造时频过完备卷积框架，具体包括：利用公式x＝Σs构造时频过完备卷积框架；其中，x为ISAR回波信号，Σ为时频过完备卷积框架，为卷积算子，CM×M为M×M维的复矩阵空间，s为矩阵S＝[s1,s2,…,sN]∈CM×N逐列首尾相连矢量化得到的向量，N为ISAR回波信号的长度。可选的，所述根据所述时频过完备卷积框架以及所述ISAR回波信号建立联合优化模型，具体包括：利用公式确定所述联合优化模型的逼近项；利用公式确定s的先验约束项；根据所述逼近项和所述先验约束项确定联合优化模型。可选的，所述根据所述逼近项和所述先验约束项确定联合优化模型，具体包括：利用公式确定联合优化模型，其中，λ>0为联合优化模型的正则化参数。根据本专利技术提供的具体实施例，本专利技术公开了以下技术效果：本专利技术所提供的一种ISAR回波信号的稳健时频分析方法，首先，根据ISAR回波信号构建时频过完备卷积框架；其次，建立联合优化模型，进而对联合优化模型进行求解，得到ISAR回波信号的时频谱，实现ISAR回波信号的稳健时频分析。其中，时频过完备卷积框架能够有效提升ISAR回波信号的时频分辨率，并且能避免ISAR回波信号中的假频、模糊等交叉项干扰。联合优化建模与数值算法设计能够克服ISAR回波信号中存在的信号部分缺失、信噪比过低等问题。附图说明为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。图1为本专利技术所提供的一种ISAR回波信号的稳健时频分析方法流程示意图；图2为本专利技术所提供的联合优化模型数值算法流程图；图3为本专利技术实施例中的时域信号示意图；图4为本专利技术实施例中构建的时频过完备卷积框架示意图；图5为本专利技术实施例中多种时频分析方法的处理结果示意图。具体实施方式下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。本专利技术的目的是提供一种ISAR回波信号的稳健时频分析方法，能够有效提升ISAR回波信号的时频分辨率。为使本专利技术的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本专利技术作进一步详细的说明。图1为本专利技术所提供的一种ISAR回波信号的稳健时频分析方法流程示意图，如图1所示，本专利技术所提供的一种ISAR回波信号的稳健时频分析方法，包括：S101，获取ISAR回波信号x＝[x[1],x[2],…,x[N]]T∈RN；其中，N为回波信号长度，RN为N维实空间。S102，根据所述ISAR回波信号确定所述ISAR回波信号的时间片信号集。ISAR回波信号的时间片信号集为：其中，M为时间片长度。S103，确定所述时间片信号集中每一时间片信号的频域表示基Φ。其中，Φ＝[φ1,φ2,…,φM]∈CM×M，时间片信号的频域表示基Φ可以设为傅立叶基、小波基等多种形式。其中，CM×M为M×M维的复矩阵空间，则经过时频表示，时间片信号xn在Φ上的线性表示为：xn＝Φsn,n＝1,2,…,N。其中，sn为表示系数，即xn的频谱，而S＝[s1,s2,…,sN]∈CM×N即为时频分析的时频谱；相应的，若在上述表示模型下建立约束模型进而重构表示系数，则可获得传统稳健时频分析。S104，根据所述频域表示基构造时频过完备卷积框架。S104具体包括：利用公式x＝Σs构造时频过完备卷积框架。其中，x为ISAR回波信号，Σ为时频过完备卷积框架，为卷积算子，CM×M为M×M维的复矩阵空间，s为矩阵S＝[s1,s2,…,sN]∈CM×N逐列首尾相连矢量化得到的向量，N为ISAR回波信号的长度。构造时频过完备卷积框架的具体过程为：联合上述N个线性表示模型，将矢量形式改写为矩阵形式X＝ΦS。其中，X＝[x1,x2,…,xN]为ISAR回波信号x的汉克矩阵，S为x的时频图。进一步上式可表示为X＝ΦSE。其中，E＝[e1,e2,…,eN]∈RN×N为一单位矩阵，RN×N为N×N维的实矩阵空间。依据卷积框架基本理论可知，IS本文档来自技高网...

【技术保护点】
1.一种ISAR回波信号的稳健时频分析方法，其特征在于，包括：/n获取ISAR回波信号；/n根据所述ISAR回波信号确定所述ISAR回波信号的时间片信号集；/n确定所述时间片信号集中每一时间片信号的频域表示基；/n根据所述频域表示基构造时频过完备卷积框架；/n根据所述时频过完备卷积框架以及所述ISAR回波信号建立联合优化模型；/n根据所述联合优化模型，利用迭代阈值算法进行求解，确定所述ISAR回波信号的时频谱。/n

【技术特征摘要】
1.一种ISAR回波信号的稳健时频分析方法，其特征在于，包括：
获取ISAR回波信号；
根据所述ISAR回波信号确定所述ISAR回波信号的时间片信号集；
确定所述时间片信号集中每一时间片信号的频域表示基；
根据所述频域表示基构造时频过完备卷积框架；
根据所述时频过完备卷积框架以及所述ISAR回波信号建立联合优化模型；
根据所述联合优化模型，利用迭代阈值算法进行求解，确定所述ISAR回波信号的时频谱。


2.根据权利要求1所述的一种ISAR回波信号的稳健时频分析方法，其特征在于，所述根据所述频域表示基构造时频过完备卷积框架，具体包括：
利用公式x＝Σs构造时频过完备卷积框架；
其中，x为ISAR回波信号，Σ为时频过完备卷积框架，为卷...

【专利技术属性】
技术研发人员：王泽龙，刘吉英，谭欣桐，余奇，郭文博，舒小虎，
申请(专利权)人：中国人民解放军国防科技大学，
类型：发明
国别省市：湖南;43