CN107622035A - 一种基于模拟退火的多项式相位信号自适应时频变换方法 - Google Patents
一种基于模拟退火的多项式相位信号自适应时频变换方法 Download PDFInfo
- Publication number
- CN107622035A CN107622035A CN201710914700.5A CN201710914700A CN107622035A CN 107622035 A CN107622035 A CN 107622035A CN 201710914700 A CN201710914700 A CN 201710914700A CN 107622035 A CN107622035 A CN 107622035A
- Authority
- CN
- China
- Prior art keywords
- mrow
- frequency
- signal
- msub
- polynomial phase
- 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
Landscapes
- Complex Calculations (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提供一种基于模拟退火的多项式相位信号自适应时频变换方法,能够完成多项式相位信号的时频分解,其中分解得到的每一个信号分量都为任一时刻都只对应一个频点的单分量,然后利用各信号分量、各时刻的瞬时频率取值,通过仅保留主瓣响应的Sinc函数直接计算生成相应时刻对应的信号频率分布,克服了传统的时频变换中一个时刻对应多个频点的非单分量存在交叉项的缺陷,最终输出无任何交叉项干扰且时频联合分辨率较优的时频分布;本发明原理简单,操作方便,可有效克服经典时频分析方法交叉项干扰的不利影响以及时频联合分辨率的损失,能够有效提升非平稳多项式相位信号时频分析的质量和效益。
Description
技术领域
本发明属于信号分析领域,尤其涉及一种基于模拟退火的多项式相位信号自适应时频变换方法。
背景技术
许多天然和人工的信号,譬如语音、生物医学信号、在色散媒质中传播的波、机械振动、动物叫声,音乐、雷达、声纳信号等,都是典型的非平稳信号,其特点是持续时间有限,并且频率是时变的,具有非平稳、非线性、非均匀、非结构、非确定、非可积、非可逆、非晶态、非规则、非连续、非光滑、非周期、非对称等特点。时频联合分析(joint time-frequencyanalysis,简称时频分析)正是着眼于真实信号组成成分的时变特征,将一个一维的时间信号以二维的时间-频率密度函数形式表示出来,旨在揭示信号中包含了多少频率分量,以及每一分量是怎样随时间变化的。
1948年,法国学者J.Ville将匈牙利布达佩斯出生的美籍物理学家E.P.Wigner在1932年提出的Wigner分布引入信号处理领域,得到了称为“Wigner-Ville分布”(Wigner-Ville distribution,WVD)。后续学者起而效仿,提出了一些新型的时频分布。整个时频分析的历史,几乎就是一部与WVD的不足作斗争的历史。按照各派的本质特征,可将形形色色的时频分布归入如下几类:(1)线性时频表示;(2)Cohen类双线性时频分布;(3)仿射类双线性时频分布;(4)重排类双线性时频分布;(5)自适应核函数类时频分布;(6)参数化时频分布。
其中,线性时频变换Gabor变换、STFT的时频分辨率受制于窗函数的形状和宽度。小波分析本质上是一种时间-尺度分析,更适于分析具有自相似结构的信号(如分形)和突变(瞬态)信号,而从刻画信号的时变结构角度看,小波变换的结果往往难于解释。
Cohen类双线性时频分布的实质,是将信号的能量(信号的某种平方形式)分布于时频平面内,其基础则为WVD。但WVD不是线性的,即两信号之和的WVD并非每一个信号的WVD之和,其中多出一个附加项。交叉项对WVD影响之烈,可见一斑。交叉项是实的,混迹于自项成分之间,且幅度较大;另外,交叉项是振荡型的,每两个信号分量就会产生一个交叉项。若信号有N个分量,则会产生个交叉项。交叉项的存在严重地干扰着人们对WVD的解释;当信号的组成成分变得复杂时,WVD给出的时频分布甚至变得毫无意义。
为了解决交叉项干扰的影响,人们相继提出了若干种不同的时频分析方法,其中Cohen类核函数时频分析方法通过设计二维核函数(二维滤波器),产生具有所需特性的时频分布。但是,该类时频分布用平滑的方法抑制交叉项,是以牺牲整个分布的时频分辨率为代价的。
前面介绍的各种分布(小波变换除外)对信号的时间和频率局部特性的刻画,是通过时移及频移变换实现的;与此对照的是,仿射类的分布则是通过时移和伸缩变换实现的。
仿射类中最著名的分布当推尺度图(Scale gram),即小波变换的平方。由于这一类分布的基础仍是WVD,因此,WVD自然地成为其成员之一。事实上,WVD正是连接Cohen类和仿射类的纽带。前者基于对WVD进行的时频平滑(time-frequency smoothing),而后者则基于仿射平滑(affine smoothing)。
Cohen类和仿射类时频分布,通过对WVD进行时频平滑及时间-尺度平滑处理,如图4所示,可以极大地抑制交叉项干扰,但仍有不少交叉项残留,而且有些分布还会引入一些新的交叉项。为进一步提高这两类分布的性能,K.Kdoera等首先提出对时频平面进行重排的思想,此后F.Auger和P.Flandrin拓展并完善了重排的方法。
在Cohen类、仿射类和重排类的双线性时频分布中,每一种分布均与一固定的核函数相对应,正是该核函数决定了相应分布的交叉项抑制特性。不言而喻,一种核函数只对一类信号有效,因而所有上述三类中的双线性时频分布均缺乏对信号的适应性。
前面介绍的各种时频分析方法,均为非参数化的方法,它们均没有先验地假定信号是由何种模型信号组成的。而参数化时频分析(parametric time frequency analysis)方法,则根据对信号层次结构的分析,构造出与信号层次结构最佳匹配的信号模型,因而能浓缩信号的信息,简化信号的表示,并由此得到信号的时频分布。
在线性时频表示(原子分解)方法中,若选用的原子与信号的主要成分相似,则仅需少数原子的线性组合,就能比较精确地表示信号,分解的结果将是稀疏的(sparse)。反之,如果原子的性状与信号的主要结构相去甚远或迥异,那么就需要用大量的,甚至无穷多的原子,才能足够精确地组装成原信号,信号的信息将弥散在太多的原子上,不利于有效地表示信号。所以在采用原子分解方法时,必须根据信号的局部结构特征,自适应地选择原子的组合,以期用尽可能少的原子来分解信号。1994年S.Qian和D.Chen开创了参数化时频分析之先河(其创造性思想可追溯到1988年),提出了“自适应展开’,(adaptive expansion)算法;1993年S.Mallat和Z.Zhang提出了自适应展开的姊妹算法—“matchingpursuit”算法。这两种算法尽管名称不同,但本质上没有不可调和的异点,实质上等价。
自适应匹配投影塔形分解算法的实质,是用原子的时频能量分布逼近原信号的时频能量分布。由于S.Qian和S.Mallat采用的是频率不变的Gabor原子,因此迭代算法对时频平面的划分是一种格型分割。这种算法对时不变的频率分量效果很好,但当待分析的信号是Chirp信号时,这种匹配相当于零阶曲线逼近,势必造成分解过程存在许多截断和分量之间的混合畸变。为克服这一缺陷,S.Mann和S.Haykin等人几乎同时提出采用经伸缩、时移、频移和频率倾斜的Gauss函数一一“chirplet”作为原子,以取代频率不变的Gabor原子(S.Mann另增加了时间倾斜操作),并用内积法得到所谓的“chirplet变换”。
基于自适应匹配投影塔形分解算法的Chirplet变换,在本质上是对时频平面上的任意一条能量曲线用一组任意倾斜的线段进行线性逼近。显然,一阶逼近比零阶逼近能更紧凑地表达chirp类信号。尽管学者们一直热爱线性,然而大自然通常是非线性的。当信号的频率成分随时间非线性地变化时,譬如一类天然存在的或人工产生的Doppler信号,用频率随时间作线性(零阶或一阶)变化的原子来表征,势必造成原子数目的增加,从而既影响对分解结果的理解和途释,又影响分解结果的数据压缩能力。
但上述自适应时频分解方法基于已知的信号层次结构或者信号模型,在信号参数上是自适应分解的,但在信号类型上是非盲的。为进一步改进上述参数化时频分析方法还需要进行新的尝试。
时频分析在语音识别,雷达信号处理和图像处理,地震信号处理,信号重构以及扩频通信中的干扰抑制等方面,已有不少成功应用。总体说来,时频分析的应用领域大致有如下四类:一是时变谱分析;二是由时频分布间接计算出某些物理量;三是利用时频分布作为信号所携带的信息的载体(而不关心它是否真能表示能量密度);四是信号的重构、压缩和编码等。
从考察信号的频率成分随时间的演化特性角度来说,小波变换的结果令人费解,尽管这一领域炙手可热。小波变换是以时间和尺度为参数,在时间-尺度平面的不同位置上具有不同的分辨率,因而是一种多分辨率分析方法。小波分析得益于小波原子的完备性、自相似性和多分辨性,它能获得成功的两个最重要的原因,是它拥有塔形快速算法和良好的时频局域特性;缺点则是一旦母小波选择不当,应用效果会大受影响。从信号压缩及消除交叉项干扰角度看,参数化的时频分析方法较好,但求其原子模型的参数,也非易事。WVD及所有其它的Cohen类时频分布,都可用于分析窄带信号,不过它们均不太适于分析多分量的宽带信号以及雷达和声纳信号,若采用仿射类的分布,则存在分辨率空变的问题。客观地说,各种时频分析技术难分优劣,关键是其适合何种类型的信号。一条屡试不爽的经验是,可以首先用运算速度很快的STFT的谱图尝试,如图3所示;若需要较高的时频分辨率,可以采用参数化时频分析方法。
但是,现有各种时频分析方法尚难完全适用于相位调制可表示为有限项多项式级数的多项式相位信号。
专利“一种模型驱动的多项式相位信号自适应时频变换方法”针对上述现实问题,提出了一种可行的新思路和新方法,但并未具体给出其中各种现代优化算法的具体的实用操作和运算过程。本发明针对这一问题,给出了一种用于非平稳多项式相位信号的自适应时频分解新算法,该算法是“模型驱动的多项式相位信号自适应时频分解新方法”的一种具体实现方式和实用运算过程及计算流程,可以自适应的完成多项式相位信号的时频分解,并输出无任何交叉项干扰且时频联合分辨率较优的时频分布。
发明内容
为解决上述问题,本发明提供一种基于模拟退火的多项式相位信号自适应时频变换方法,利用各信号分量、各时刻的瞬时频率取值,通过仅保留主瓣响应的Sinc函数直接计算生成相应时刻对应的信号频率分布,克服了传统的时频变换中一个时刻对应多个频点的非单分量存在交叉项的缺陷,最终输出无任何交叉项干扰且时频联合分辨率较优的时频分布。
一种基于模拟退火的多项式相位信号自适应时频变换方法,包括以下步骤:
步骤1:利用模拟退火算法对雷达接收的原始多项式相位信号s(t)进行多次分解,每次分解均确定一个信号分量的最优模型阶数和最优待定系数集,具体的:
步骤11:根据多项式相位模型,生成原始多项式相位信号s(t)对应的参考函数hp,并初始化参考函数hp的模型阶数N1=1,并计算原始多项式相位信号s(t)的能量E0,随机生成参考函数hp的两个待定系数集{an},n=0,1,2,...,N1,其中待定系数集{an}作为模拟退火算法中的初始解;
步骤12:将两个待定系数集{an}对应的参考函数hp的共轭分别与原始多项式相位信号s(t)相乘得到混合调制信号x(t),对混合调制信号x(t)实施傅利叶变换,得到变换结果X(f),计算变换结果X(f)其中一个频域特征的特征值,并将该特征值作为评价变换结果X(f)的适应度,根据各个变换结果X(f)的适应度优劣,其中适应度较优的变换结果X(f)对应的初始解为当前最优解Gbest,适应度较差的为次级最优解PreBest,同时将当前最优解Gbest赋值给Prelocal;
步骤13:根据Prelocal,按照设定的退火步长因子Tstep计算新的可能解Nextp,并按步骤12的方法计算可能解Nextp的适应度fitnext,如果可能解Nextp的适应度fitnext优于步骤12所述的当前最优解Gbest的适应度,则将其替换为当前最优解Gbest,并将原最优解替换为次级最优解PreBest,否则当前最优解Gbest和次级最优解PreBest保持不变;
步骤14:随机生成另一个新的可能解Nextq,并将其对应的待定系数集{an}按步骤12的方法,计算该新的可能解Nextq对应变换结果X(f)的适应度fitnew,如果适应度fitnew优于适应度fitnext,则将新的可能解Nextq替换为Prelocal,否则根据替换概率pω与随机概率rand的大小来决定是否将新的可能解Nextq替换为Prelocal,其中如果替换概率pω大于随机概率rand,则新的可能解Nextq替换为Prelocal,否则Prelocal保持不变;
步骤15:重复步骤13-步骤14,直至循环次数达到设定的马尔科夫链的数值Lmarkov,得到迭代完成的当前最优解Gbest与次优解PreBest;
步骤16:计算步骤15得到的当前最优解Gbest与次优解PreBest的适应度差值Δfit,并计算退火温度Tt,其中退火温度Tt随着步骤16的重复次数增多而衰减;如果适应度差值Δfit的绝对值大于设定的退火的容限ξ,且退火温度Tt大于设定的于终止温度Tlow,则重复步骤13-步骤15,直到适应度差值Δfit的绝对值不大于设定的退火的容限ξ,或退火温度Tt不大于设定的于终止温度Tlow,从而得到最终的最优解Gbest';
步骤17:令模型阶数N1依次从2取到最大可能阶数order_max,重复步骤11至步骤16,从而得到不同模型阶数下的最终的最优解Gbest′;
步骤18:从步骤17得到所有最优解Gbest'中选出适应度最优的解,从而确定原始多项式相位信号s(t)当前信号分量的最优模型阶数Np以及Np对应的最优待定系数集{an}max;
步骤2:利用最优模型阶数Np和最优待定系数集{an}max确定的信号分量的共轭与原始多项式相位信号s(t)相乘并实施傅利叶变换,得到频谱X'(f),将频谱X'(f)包络最大值处的强度复数取值置零后,实施逆傅立叶变换,从而得到时域信号y(t);
步骤3:利用模型阶数为Np和最优待定系数集{an}max确定的信号分量与时域信号y(t)相乘得到残差信号z(t),本次分解结束;
步骤4:计算残差信号z(t)的能量Ed,与原始多项式相位信号s(t)能量E0取比值R,如果比值R小于设定门限γ或分解次数达到设定的上限数量Nmax,则停止分解,得到各分量信号的最优模型阶数Np和最优待定系数集{an}max;否则,利用步骤3计算的残差信号z(t)替换步骤1中的原始多项式相位信号s(t)重新计算混合调制信号x(t),重复步骤1至步骤3,直至比值R小于设定门限γ或分解次数达到设定的上限数量Nmax;
步骤5:将最优模型阶数Np和最优待定系数集{an}max对应的各个信号分量按分解顺序编号为hc,其中c=1,2,...,C,其中C为自适应分解的次数;
步骤6:按单分量信号瞬时频率的物理定义,得到各信号分量hc对应的瞬时频率曲线fc(t),并根据所有信号分量hc的瞬时频率曲线fc(t)确定原始多项式相位信号s(t)的频率分布范围,最后对所述频率分布范围依次进行离散化、保留主瓣响应以及叠加操作后,得到原始多项式相位信号s(t)的时频联合分布f(t)。
进一步地,所述原始多项式相位信号s(t)对应的参考函数hp具体为:
其中,order_max为原始多项式相位信号s(t)的最大可能阶数。
进一步地,所述原始多项式相位信号s(t)的频率分布范围具体为:
所有信号分量hc的瞬时频率曲线fc(t)中的瞬时频率最小值为原始多项式相位信号s(t)的最小频率fimin,瞬时频率最大值为原始多项式相位信号s(t)的最大频率fimax。
进一步地,所述待定系数集{an}具体为矢量其中N1代表当前优化的模型阶数,a0为固定取值,设为原始多项式相位信号s(t)直接傅利叶变换所得峰值的所在采样点的位置序号;为待优化的待定系数,由随机函数Rand产生,取值范围为[-M,M],其中M为原始多项式相位信号s(t)的采样点数目,具体的:
其中rand(2,N1)表示随机产生2个维数为N1的分布于[0,1]的随机数,从而得到模型阶数为N1的初始解。
进一步地,步骤12所述的频域特征为频谱峰值、频谱的对比度或频谱的谱形熵中的任意一个,其中各频域特征的特征值频谱峰值Xp、频谱的对比度Xc以及频谱的谱形熵Xe分别为:
Xp=max{abs[X(f)]}
Xe=sum(-pS·logpS)
其中
sum(·)代表向量的求和操作,abs(·)为取绝对值;
当选择频谱峰值Xp或频谱对比度Xc作为适应度时,该值越大,评价结果越优,当选择频谱的谱形熵Xe作为适应度时,该值越小,评价结果越优。
进一步地,步骤13所述的按照设定的退火步长因子Tstep计算新的可能解Nextp,具体方法为:
Nextp=Prelocal+Tstep·[2M·rand(1,N1)-M]
其中rand(1,N1)表示随机产生1个维数为N1的分布于[0,1]的随机数,M为原始多项式相位信号s(t)的采样点数目。
进一步地,步骤14所述的替换概率pω具体计算方法为:
pω=exp[-(fitnext-fitnew)]。
进一步地,步骤16所述的退火温度Tt具体计算方法为:
Tt=KtT0
其中,K为设定的模拟退火的温度衰减参数,T0为设定的模拟退火初始温度,t为步骤16的重复次数,至少为1。
进一步地,步骤6所述的对所述频率分布范围依次进行离散化具体为:
fis=[fimin fimin+Δf fimin+2Δf fimin+3Δf ... fimax]T
其中,fis为离散频率向量的瞬时频率值,Δf为设定的频率分辨率,具体的:
其中,fimin为原始多项式相位信号s(t)的最小频率,fimax为原始多项式相位信号s(t)的最大频率,K为设定的离散频率维数,i为虚部单位。
进一步地,步骤6所述的对所述频率分布范围进行保留主瓣响应具体计算方法为:
其中,Ac为第c个信号分量hc的频谱包络最大值处的强度复数取值,i为虚部单位,Δf为设定的频率分辨率,fis为离散频率向量fs的瞬时频率值、fic(t)为各信号分量hc对应的瞬时频率曲线fc(t)的瞬时频率值。
有益效果:
本发明提供一种基于模拟退火的多项式相位信号自适应时频变换方法,能够完成多项式相位信号的时频分解,其中分解得到的每一个信号分量都为任一时刻都只对应一个频点的单分量,然后利用各信号分量、各时刻的瞬时频率取值,通过仅保留主瓣响应的Sinc函数直接计算生成相应时刻对应的信号频率分布,克服了传统的时频变换中一个时刻对应多个频点的非单分量存在交叉项的缺陷,最终输出无任何交叉项干扰且时频联合分辨率较优的时频分布;
本发明原理简单,操作方便,可有效克服经典时频分析方法交叉项干扰的不利影响以及时频联合分辨率的损失,能够有效提升非平稳多项式相位信号时频分析的质量和效益。
附图说明
图1为本发明基于模拟退火的自适应时频变换算法流程图;
图2为本发明其中一个多分量多项式相位信号的时域波形图;
图3为现有技术中经典STFT分析时频图;
图4为现有技术在中经典WVD时频图;
图5为本发明基于模拟退火的自适应时频变换算法所得的时频图。
具体实施方式
下面结合附图并举实施例,对本发明进行详细叙述。
本发明首先通过模拟退火估计出多项式相位信号各组成分量的模型阶数及对应的各阶相位参数,然后利用“清洁”的思想,提取出相应的多项式相位信号分量,并从原始信号中剔除该信号分量得到残差信号,接着以迭代“清洁”的方式,重复利用模拟退火针对残差信号实施自适应模型定阶及最优模型参数优化,并逐步提取出各信号分量,如此反复,直至残差信号能量低于预设的门限或提取出的信号分量数目达到预设的最大值。随后,利用模拟退火分解所得的各分量信号的相位参数构建相位-时间历史,并根据单分量信号瞬时频率的物理定义,对相位历史求导,以解析的方式直接得到各信号分量的频率-时间历史,即时频变化曲线,根据信号各分量的时频变化曲线,确定整个信号的最大及最小频率,利用其确定整个信号的频率变化范围,按应用所需的频率分辨率需求,对该频率范围进行离散化,得到所需维数的离散频率向量,最后,利用各分量、各时刻的瞬时频率取值,通过仅保留主瓣响应的Sinc函数直接计算生成相应时刻对应的信号频率分布;如此反复,直至所有信号分量所有时刻的频率分布均生成完毕,将其按时间先后顺序存储并逐信号分量累加,即可得到最终多项式相位信号的时频联合分布;本发明利用模拟退火过程实现,原理简单,操作方便,可有效克服经典时频分析方法交叉项干扰的不利影响以及时频联合分辨率的损失,提升非平稳多项式相位信号时频分析的质量和效益。
如图1所示,一种基于模拟退火的多项式相位信号自适应时频变换方法,包括以下步骤:
步骤1:利用模拟退火算法对雷达接收的原始多项式相位信号s(t)进行多次分解,每次分解均确定一个信号分量的最优模型阶数和最优待定系数集,具体的:
步骤11:根据多项式相位模型,生成原始多项式相位信号s(t)对应的参考函数hp,并初始化参考函数hp的模型阶数N1=1,并计算原始多项式相位信号s(t)的能量E0,随机生成参考函数hp的两个待定系数集{an},n=0,1,2,...,N1,其中待定系数集{an}作为模拟退火算法中的初始解;
步骤12:将两个待定系数集{an}对应的参考函数hp的共轭分别与原始多项式相位信号s(t)相乘得到混合调制信号x(t),对混合调制信号x(t)实施傅利叶变换,得到变换结果X(f),计算变换结果X(f)其中一个频域特征的特征值,并将该特征值作为评价变换结果X(f)的适应度,根据各个变换结果X(f)的适应度优劣,其中适应度较优的变换结果X(f)对应的初始解为当前最优解Gbest,适应度较差的变换结果X(f)对应的初始解为次级最优解PreBest,同时将当前最优解Gbest赋值给Prelocal;
步骤13:根据Prelocal,按照设定的退火步长因子Tstep计算新的可能解Nextp,并按步骤12的方法计算可能解Nextp的适应度fitnext,如果可能解Nextp的适应度指标fitnext优于步骤12所述的当前最优解Gbest的适应度,则将其替换为当前最优解Gbest,并将原最优解替换为次级最优解PreBest,否则当前最优解Gbest和次级最优解PreBest保持不变;
步骤14:实施Metropolis过程:随机生成另一个新的可能解Nextq,并将其对应的待定系数集{an}按步骤12的方法,计算该新的可能解Nextq对应变换结果X(f)的适应度fitnew,如果适应度fitnew优于适应度fitnext,则将新的可能解Nextq替换为Prelocal,否则根据替换概率pω与随机概率rand的大小来决定是否将新的可能解Nextq替换为Prelocal,其中如果替换概率pω大于随机概率rand,则新的可能解Nextq替换为Prelocal,否则Prelocal保持不变;
步骤15:重复步骤13-步骤14,直至循环次数达到设定的马尔科夫链的数值Lmarkov,得到迭代完成的当前最优解Gbest与次优解PreBest;
步骤16:计算步骤15得到的当前最优解Gbest与次优解PreBest的适应度差值Δfit,并计算退火温度Tt,其中退火温度Tt随着步骤16的重复次数增多而衰减;如果适应度差值Δfit的绝对值大于设定的退火的容限ξ,且退火温度Tt大于设定的于终止温度Tlow,则重复步骤13-步骤15,直到适应度差值Δfit的绝对值不大于设定的退火的容限ξ,或退火温度Tt不大于设定的于终止温度Tlow,从而得到最终的最优解Gbest';
步骤17:令模型阶数N1依次从2取到最大可能阶数order_max,重复步骤11至步骤16,从而得到不同模型阶数下的最终的最优解Gbest′;
步骤18:从步骤17得到所有最优解Gbest'中选出适应度最优的解,从而确定原始多项式相位信号s(t)当前信号分量的最优模型阶数Np以及Np对应的最优待定系数集{an}max;
步骤2:利用最优模型阶数Np和最优待定系数集{an}max确定的信号分量的共轭与原始多项式相位信号s(t)相乘并实施傅利叶变换,得到频谱X'(f),将频谱X'(f)包络最大值处的强度复数取值置零后,实施逆傅立叶变换,从而得到时域信号y(t),如图2所示;
步骤3:利用模型阶数为Np和最优待定系数集{an}max确定的信号分量与时域信号y(t)相乘得到残差信号z(t),本次分解结束;
步骤4:计算残差信号z(t)的能量Ed,与原始多项式相位信号s(t)能量E0取比值R,如果比值R小于设定门限γ或分解次数达到设定的上限数量Nmax,则停止分解,得到各分量信号的最优模型阶数Np和最优待定系数集{an}max;否则,利用步骤3计算的残差信号z(t)替换步骤1中的原始多项式相位信号s(t)重新计算混合调制信号x(t),重复步骤1至步骤3,直至比值R小于设定门限γ或分解次数达到设定的上限数量Nmax;
步骤5:将最优模型阶数Np和最优待定系数集{an}max对应的各个信号分量按分解顺序编号为hc,其中c=1,2,...,C,其中C为自适应分解的次数;
步骤6:按单分量信号瞬时频率的物理定义,得到各信号分量hc对应的瞬时频率曲线fc(t),并根据所有信号分量hc的瞬时频率曲线fc(t)确定原始多项式相位信号s(t)的频率分布范围,最后对所述频率分布范围依次进行离散化、保留主瓣响应以及叠加操作后,得到原始多项式相位信号s(t)的时频联合分布f(t),如图5所示。
进一步地,所述原始多项式相位信号s(t)对应的参考函数hp具体为:
其中,待定系数集{an},n=0,1,2,...N1,中的待定系数an初始值随机生成,order_max为原始多项式相位信号s(t)的最大可能阶数。
进一步地,所述原始多项式相位信号s(t)的频率分布范围具体为:
所有信号分量hc的瞬时频率曲线fc(t)中的瞬时频率最小值为原始多项式相位信号s(t)的最小频率fimin,瞬时频率最大值为原始多项式相位信号s(t)的最大频率fimax。
进一步地,所述待定系数集{an}具体为矢量其中N1代表当前优化的模型阶数,a0为固定取值,设为原始多项式相位信号s(t)直接傅利叶变换所得峰值的所在采样点的位置序号;为待优化的待定系数,由随机函数Rand产生,取值范围为[-M,M],其中M为原始多项式相位信号s(t)的采样点数目,具体的:
其中rand(2,N1)表示随机产生2个维数为N1的分布于[0,1]的随机数,从而得到模型阶数为N1的初始解。
进一步地,步骤12所述的频域特征为频谱峰值、频谱的对比度或频谱的谱形熵中的任意一个,其中各频域特征的特征值频谱峰值Xp、频谱的对比度Xc以及频谱的谱形熵Xe分别为:
Xp=max{abs[X(f)]}
Xe=sum(-pS·logpS)
其中
sum(·)代表向量的求和操作,abs(·)为取绝对值;
当选择频谱峰值Xp或频谱对比度Xc作为适应度时,该值越大,评价结果越优,当选择频谱的谱形熵Xe作为适应度时,该值越小,评价结果越优。
进一步地,步骤13所述的按照设定的退火步长因子Tstep计算新的可能解Nextp,具体方法为:
Nextp=Prelocal+Tstep·[2M·rand(1,N1)-M]
其中rand(1,N1)表示随机产生1个维数为N1的分布于[0,1]的随机数,M为原始多项式相位信号s(t)的采样点数目。
进一步地,步骤14所述的替换概率pω具体计算方法为:
pω=exp[-(fitnext-fitnew)]。
进一步地,步骤16所述的退火温度Tt具体计算方法为:
Tt=KtT0
其中,K为设定的模拟退火的温度衰减参数,T0为设定的模拟退火初始温度,t为步骤16的重复次数,至少为1。
进一步地,步骤6所述的对所述频率分布范围依次进行离散化具体为:
fis=[fimin fimin+Δf fimin+2Δf fimin+3Δf ... fimax]T
其中,fis为离散频率向量的瞬时频率值,Δf为设定的频率分辨率,具体的:
其中,fimin为原始多项式相位信号s(t)的最小频率,fimax为原始多项式相位信号s(t)的最大频率,K为设定的离散频率维数,i为虚部单位。
进一步地,步骤6所述的对所述频率分布范围进行保留主瓣响应具体计算方法为:
其中,前一个因式代表Sinc函数直接频率分布;后一个因式代表仅保留主瓣响应的操作,Ac为第c个信号分量hc的频谱包络最大值处的强度复数取值,i为虚部单位,Δf为设定的频率分辨率,fis为离散频率向量fs的瞬时频率值、fic(t)为各信号分量hc对应的瞬时频率曲线fc(t)的瞬时频率值。
当然,本发明还可有其他多种实施例,在不背离本发明精神及其实质的情况下,熟悉本领域的技术人员当可根据本发明作出各种相应的改变和变形,但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。
Claims (10)
1.一种基于模拟退火的多项式相位信号自适应时频变换方法,其特征在于,包括以下步骤:
步骤1:利用模拟退火算法对雷达接收的原始多项式相位信号s(t)进行多次分解,每次分解均确定一个信号分量的最优模型阶数和最优待定系数集,具体的:
步骤11:根据多项式相位模型,生成原始多项式相位信号s(t)对应的参考函数hp,并初始化参考函数hp的模型阶数N1=1,并计算原始多项式相位信号s(t)的能量E0,随机生成参考函数hp的两个待定系数集{an},n=0,1,2,...,N1,其中待定系数集{an}作为模拟退火算法中的初始解;
步骤12:将两个待定系数集{an}对应的参考函数hp的共轭分别与原始多项式相位信号s(t)相乘得到混合调制信号x(t),对混合调制信号x(t)实施傅利叶变换,得到变换结果X(f),计算变换结果X(f)其中一个频域特征的特征值,并将该特征值作为评价变换结果X(f)的适应度,根据各个变换结果X(f)的适应度优劣,其中适应度较优的变换结果X(f)对应的初始解为当前最优解Gbest,适应度较差的为次级最优解PreBest,同时将当前最优解Gbest赋值给Prelocal;
步骤13:根据Prelocal,按照设定的退火步长因子Tstep计算新的可能解Nextp,并按步骤12的方法计算可能解Nextp的适应度fitnext,如果可能解Nextp的适应度fitnext优于步骤12所述的当前最优解Gbest的适应度,则将其替换为当前最优解Gbest,并将原最优解替换为次级最优解PreBest,否则当前最优解Gbest和次级最优解PreBest保持不变;
步骤14:随机生成另一个新的可能解Nextq,并将其对应的待定系数集{an}按步骤12的方法,计算该新的可能解Nextq对应变换结果X(f)的适应度fitnew,如果适应度fitnew优于适应度fitnext,则将新的可能解Nextq替换为Prelocal,否则根据替换概率pω与随机概率rand的大小来决定是否将新的可能解Nextq替换为Prelocal,其中如果替换概率pω大于随机概率rand,则新的可能解Nextq替换为Prelocal,否则Prelocal保持不变;
步骤15:重复步骤13-步骤14,直至循环次数达到设定的马尔科夫链的数值Lmarkov,得到迭代完成的当前最优解Gbest与次优解PreBest;
步骤16:计算步骤15得到的当前最优解Gbest与次优解PreBest的适应度差值Δfit,并计算退火温度Tt,其中退火温度Tt随着步骤16的重复次数增多而衰减;如果适应度差值Δfit的绝对值大于设定的退火的容限ξ,且退火温度Tt大于设定的于终止温度Tlow,则重复步骤13-步骤15,直到适应度差值Δfit的绝对值不大于设定的退火的容限ξ,或退火温度Tt不大于设定的于终止温度Tlow,从而得到最终的最优解Gbest';
步骤17:令模型阶数N1依次从2取到最大可能阶数order_max,重复步骤11至步骤16,从而得到不同模型阶数下的最终的最优解Gbest′;
步骤18:从步骤17得到所有最优解Gbest'中选出适应度最优的解,从而确定原始多项式相位信号s(t)当前信号分量的最优模型阶数Np以及Np对应的最优待定系数集{an}max;
步骤2:利用最优模型阶数Np和最优待定系数集{an}max确定的信号分量的共轭与原始多项式相位信号s(t)相乘并实施傅利叶变换,得到频谱X'(f),将频谱X'(f)包络最大值处的强度复数取值置零后,实施逆傅立叶变换,从而得到时域信号y(t);
步骤3:利用模型阶数为Np和最优待定系数集{an}max确定的信号分量与时域信号y(t)相乘得到残差信号z(t),本次分解结束;
步骤4:计算残差信号z(t)的能量Ed,与原始多项式相位信号s(t)能量E0取比值R,如果比值R小于设定门限γ或分解次数达到设定的上限数量Nmax,则停止分解,得到各分量信号的最优模型阶数Np和最优待定系数集{an}max;否则,利用步骤3计算的残差信号z(t)替换步骤1中的原始多项式相位信号s(t)重新计算混合调制信号x(t),重复步骤1至步骤3,直至比值R小于设定门限γ或分解次数达到设定的上限数量Nmax;
步骤5:将最优模型阶数Np和最优待定系数集{an}max对应的各个信号分量按分解顺序编号为hc,其中c=1,2,...,C,其中C为自适应分解的次数;
步骤6:按单分量信号瞬时频率的物理定义,得到各信号分量hc对应的瞬时频率曲线fc(t),并根据所有信号分量hc的瞬时频率曲线fc(t)确定原始多项式相位信号s(t)的频率分布范围,最后对所述频率分布范围依次进行离散化、保留主瓣响应以及叠加操作后,得到原始多项式相位信号s(t)的时频联合分布f(t)。
2.如权利要求1所述的一种基于模拟退火的多项式相位信号自适应时频变换方法,其特征在于,所述原始多项式相位信号s(t)对应的参考函数hp具体为:
<mrow>
<msub>
<mi>h</mi>
<mi>p</mi>
</msub>
<mo>=</mo>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mi>j</mi>
<mo>&CenterDot;</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<msub>
<mi>N</mi>
<mn>1</mn>
</msub>
</munderover>
<msub>
<mi>a</mi>
<mi>n</mi>
</msub>
<msup>
<mi>t</mi>
<mi>n</mi>
</msup>
<mo>&rsqb;</mo>
<mo>,</mo>
<mn>0</mn>
<mo>&le;</mo>
<msub>
<mi>N</mi>
<mn>1</mn>
</msub>
<mo>&le;</mo>
<mi>o</mi>
<mi>r</mi>
<mi>d</mi>
<mi>e</mi>
<mi>r</mi>
<mo>_</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
其中,order_max为原始多项式相位信号s(t)的最大可能阶数。
3.如权利要求1所述的一种基于模拟退火的多项式相位信号自适应时频变换方法,其特征在于,所述原始多项式相位信号s(t)的频率分布范围具体为:
所有信号分量hc的瞬时频率曲线fc(t)中的瞬时频率最小值为原始多项式相位信号s(t)的最小频率fimin,瞬时频率最大值为原始多项式相位信号s(t)的最大频率fimax。
4.如权利要求1所述的一种基于模拟退火的多项式相位信号自适应时频变换方法,其特征在于,所述待定系数集{an}具体为矢量其中N1代表当前优化的模型阶数,a0为固定取值,设为原始多项式相位信号s(t)直接傅利叶变换所得峰值的所在采样点的位置序号;为待优化的待定系数,由随机函数Rand产生,取值范围为[-M,M],其中M为原始多项式相位信号s(t)的采样点数目,具体的:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>a</mi>
<msub>
<mi>N</mi>
<mn>1</mn>
</msub>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mn>2</mn>
<mi>M</mi>
<mo>&CenterDot;</mo>
<mi>r</mi>
<mi>a</mi>
<mi>n</mi>
<mi>d</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>,</mo>
<msub>
<mi>N</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>M</mi>
</mrow>
其中rand(2,N1)表示随机产生2个维数为N1的分布于[0,1]的随机数,从而得到模型阶数为N1的初始解。
5.如权利要求1所述的一种基于模拟退火的多项式相位信号自适应时频变换方法,其特征在于,步骤12所述的频域特征为频谱峰值、频谱的对比度或频谱的谱形熵中的任意一个,其中各频域特征的特征值频谱峰值Xp、频谱的对比度Xc以及频谱的谱形熵Xe分别为:
Xp=max{abs[X(f)]}
<mrow>
<msub>
<mi>X</mi>
<mi>c</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
<mo>{</mo>
<mi>a</mi>
<mi>b</mi>
<mi>s</mi>
<mo>&lsqb;</mo>
<mi>X</mi>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>}</mo>
<mo>-</mo>
<mi>min</mi>
<mo>{</mo>
<mi>a</mi>
<mi>b</mi>
<mi>s</mi>
<mo>&lsqb;</mo>
<mi>X</mi>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>}</mo>
</mrow>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
<mo>{</mo>
<mi>a</mi>
<mi>b</mi>
<mi>s</mi>
<mo>&lsqb;</mo>
<mi>X</mi>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>}</mo>
</mrow>
</mfrac>
</mrow>
Xe=sum(-pS·logpS)
其中
<mrow>
<msub>
<mi>p</mi>
<mi>S</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>a</mi>
<mi>b</mi>
<mi>s</mi>
<mo>&lsqb;</mo>
<mi>X</mi>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<munder>
<mo>&Sigma;</mo>
<mi>f</mi>
</munder>
<mi>a</mi>
<mi>b</mi>
<mi>s</mi>
<mo>&lsqb;</mo>
<mi>X</mi>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mfrac>
</mrow>
sum(·)代表向量的求和操作,abs(·)为取绝对值;
当选择频谱峰值Xp或频谱对比度Xc作为适应度时,该值越大,评价结果越优,当选择频谱的谱形熵Xe作为适应度时,该值越小,评价结果越优。
6.如权利要求1所述的一种基于模拟退火的多项式相位信号自适应时频变换方法,其特征在于,步骤13所述的按照设定的退火步长因子Tstep计算新的可能解Nextp,具体方法为:
Nextp=Prelocal+Tstep·[2M·rand(1,N1)-M]
其中rand(1,N1)表示随机产生1个维数为N1的分布于[0,1]的随机数,M为原始多项式相位信号s(t)的采样点数目。
7.如权利要求1所述的一种基于模拟退火的多项式相位信号自适应时频变换方法,其特征在于,步骤14所述的替换概率pω具体计算方法为:
pω=exp[-(fitnext-fitnew)]。
8.如权利要求1所述的一种基于模拟退火的多项式相位信号自适应时频变换方法,其特征在于,步骤16所述的退火温度Tt具体计算方法为:
Tt=KtT0
其中,K为设定的模拟退火的温度衰减参数,T0为设定的模拟退火初始温度,t为步骤16的重复次数,至少为1。
9.如权利要求1所述的一种基于模拟退火的多项式相位信号自适应时频变换方法,其特征在于,步骤6所述的对所述频率分布范围依次进行离散化具体为:
fis=[fimin fimin+Δf fimin+2Δf fimin+3Δf ... fimax]T
其中,fis为离散频率向量的瞬时频率值,Δf为设定的频率分辨率,具体的:
<mrow>
<mi>&Delta;</mi>
<mi>f</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mi>i</mi>
<mi>max</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>i</mi>
<mi>min</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mi>K</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mfrac>
</mrow>
其中,fimin为原始多项式相位信号s(t)的最小频率,fimax为原始多项式相位信号s(t)的最大频率,K为设定的离散频率维数,i为虚部单位。
10.如权利要求1所述的一种基于模拟退火的多项式相位信号自适应时频变换方法,其特征在于,步骤6所述的对所述频率分布范围进行保留主瓣响应具体计算方法为:
<mrow>
<msubsup>
<mi>f</mi>
<mi>c</mi>
<mo>&prime;</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>A</mi>
<mi>c</mi>
</msub>
<mo>&times;</mo>
<mi>sin</mi>
<mi>c</mi>
<mrow>
<mo>{</mo>
<mrow>
<mn>0.886</mn>
<mfrac>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mi>i</mi>
<mi>s</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>i</mi>
<mi>c</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>f</mi>
</mrow>
</mfrac>
</mrow>
<mo>}</mo>
</mrow>
<mo>&times;</mo>
<mrow>
<mo>{</mo>
<mrow>
<mo>-</mo>
<mn>1</mn>
<mo>&le;</mo>
<mn>0.886</mn>
<mfrac>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mi>i</mi>
<mi>s</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>i</mi>
<mi>c</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>f</mi>
</mrow>
</mfrac>
<mo>&le;</mo>
<mn>1</mn>
</mrow>
<mo>}</mo>
</mrow>
</mrow>
其中,Ac为第c个信号分量hc的频谱包络最大值处的强度复数取值,i为虚部单位,Δf为设定的频率分辨率,fis为离散频率向量fs的瞬时频率值、fic(t)为各信号分量hc对应的瞬时频率曲线fc(t)的瞬时频率值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710914700.5A CN107622035B (zh) | 2017-09-30 | 2017-09-30 | 一种基于模拟退火的多项式相位信号自适应时频变换方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710914700.5A CN107622035B (zh) | 2017-09-30 | 2017-09-30 | 一种基于模拟退火的多项式相位信号自适应时频变换方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107622035A true CN107622035A (zh) | 2018-01-23 |
CN107622035B CN107622035B (zh) | 2020-07-17 |
Family
ID=61091513
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710914700.5A Active CN107622035B (zh) | 2017-09-30 | 2017-09-30 | 一种基于模拟退火的多项式相位信号自适应时频变换方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107622035B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112130138A (zh) * | 2020-07-06 | 2020-12-25 | 电子科技大学 | 一种超视距雷达的电离层相位污染校正方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102023963A (zh) * | 2010-12-13 | 2011-04-20 | 西安电子科技大学 | 高速多模式时频域变换方法 |
US20130028297A1 (en) * | 2011-05-04 | 2013-01-31 | Casey Stephen D | Windowing methods and systems for use in time-frequency analysis |
CN104731762A (zh) * | 2015-03-25 | 2015-06-24 | 西安电子科技大学 | 基于循环移位的立方相位信号参数估计方法 |
CN105606894A (zh) * | 2016-01-28 | 2016-05-25 | 南京信息工程大学 | 基于模拟退火算法的瞬时频率估计方法 |
-
2017
- 2017-09-30 CN CN201710914700.5A patent/CN107622035B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102023963A (zh) * | 2010-12-13 | 2011-04-20 | 西安电子科技大学 | 高速多模式时频域变换方法 |
US20130028297A1 (en) * | 2011-05-04 | 2013-01-31 | Casey Stephen D | Windowing methods and systems for use in time-frequency analysis |
CN104731762A (zh) * | 2015-03-25 | 2015-06-24 | 西安电子科技大学 | 基于循环移位的立方相位信号参数估计方法 |
CN105606894A (zh) * | 2016-01-28 | 2016-05-25 | 南京信息工程大学 | 基于模拟退火算法的瞬时频率估计方法 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112130138A (zh) * | 2020-07-06 | 2020-12-25 | 电子科技大学 | 一种超视距雷达的电离层相位污染校正方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107622035B (zh) | 2020-07-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107729289A (zh) | 一种基于遗传优化的多项式相位信号自适应时频变换方法 | |
Kim et al. | Softflow: Probabilistic framework for normalizing flow on manifolds | |
CN109101910B (zh) | 一种基于噪声甄别的大地电磁信号去噪方法 | |
CN103778919B (zh) | 基于压缩感知和稀疏表示的语音编码方法 | |
CN107945809B (zh) | 一种复调音乐多音高估计方法 | |
CN107516065B (zh) | 经验模态分解结合字典学习的复杂信号去噪方法 | |
CN110334580A (zh) | 基于集成增量的动态权重组合的设备故障分类方法 | |
CN109524020A (zh) | 一种语音增强处理方法 | |
CN106537420A (zh) | 用于转换输入信号的方法 | |
CN104771163A (zh) | 基于csp和r-csp算法的脑电信号特征提取方法 | |
CN105677035A (zh) | 基于eemd和小波阈值的运动想象脑电信号消噪方法 | |
CN105760347A (zh) | 一种基于数据/极值联合对称延拓的hht端点效应抑制方法 | |
CN107729288A (zh) | 一种基于粒子群优化的多项式相位信号时频变换方法 | |
CN102609910A (zh) | 基于Ridgelet冗余字典的遗传进化图像重构方法 | |
CN104089774A (zh) | 一种基于并行多字典正交匹配的齿轮故障诊断方法 | |
CN108334872A (zh) | 基于改进hht变换的特征提取方法 | |
Bayram et al. | A penalty function promoting sparsity within and across groups | |
CN102012466A (zh) | 数字x射线成像***的噪声测量方法 | |
CN107622036A (zh) | 一种基于蚁群优化的多项式相位信号自适应时频变换方法 | |
CN107783939A (zh) | 一种模型驱动的多项式相位信号自适应时频变换方法 | |
CN102332268B (zh) | 基于自适应冗余字典的语音信号稀疏表示方法 | |
CN107622035A (zh) | 一种基于模拟退火的多项式相位信号自适应时频变换方法 | |
Tomassi et al. | Wavelet shrinkage using adaptive structured sparsity constraints | |
CN114067819A (zh) | 基于跨层相似性知识蒸馏的语音增强方法 | |
Siripatana et al. | Combining ensemble Kalman filter and multiresolution analysis for efficient assimilation into adaptive mesh models |
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 |