CN105446940B - 一种基于同步压缩变换重构的幅值校正方法 - Google Patents
一种基于同步压缩变换重构的幅值校正方法 Download PDFInfo
- Publication number
- CN105446940B CN105446940B CN201510821373.XA CN201510821373A CN105446940B CN 105446940 B CN105446940 B CN 105446940B CN 201510821373 A CN201510821373 A CN 201510821373A CN 105446940 B CN105446940 B CN 105446940B
- Authority
- CN
- China
- Prior art keywords
- frequency
- signal
- synchronous compression
- coefficient
- fourier transform
- 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.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
Abstract
一种基于同步压缩变换重构的幅值校正方法,对线性调频信号进行短时傅里叶变换,得到短时傅里叶变换系数,然后估算瞬时频率,在该估计瞬时频率的基础上进行改进,计算准确瞬时频率,再计算中心频率,然后完成同步压缩变换,即将所有对应同一频率段的准确瞬时频率的短时傅里叶变换系数转化为同步压缩系数,最后提取同步压缩变换的系数脊线,并将脊线附近的系数进行整合,得到相应的信号成分,由于本发明对瞬时频率的计算进行了改进,能够得到准确的瞬时频率,因此提高了时域重构信号的精度。
Description
技术领域
本发明属于信号处理技术领域,具体涉及一种基于同步压缩变换重构的幅值校正方法。
背景技术
时频分析方法是处理非平稳信号的有力工具,目前常用的有经典的时频方法,如短时傅里叶变换和小波变换,近几年有学者通过这些经典的方法,得出了一种新的时频分析方法:同步压缩变换。该方法类似经验模式分解(EMD),可以进行信号分解,得到单成分信号的时域图,即对信号进行重构。
根据文献【Iatsenko D,Mcclintock P V E,Stefanovska A.Linear andsynchrosqueezed time-frequency representations revisited:Overview,standardsof use,resolution,reconstruction,concentration,and algorithms[J].DigitalSignal Processing,2015,42.】中可以总结出一般的同步压缩变换,但是一般的同步压缩变换中使用了瞬时频率的近似值,该近似值大大降低了时域重构信号的精度,所得到的信号的幅值相比原始信号具有很大的差别。
发明内容
为了克服上述现有技术的缺点,本发明的目的在于提供一种基于同步压缩变换重构的幅值校正方法,校正所得到信号的幅值,提高时域重构信号的精度。
为了达到上述目的,本发明采取的技术方案为:
一种基于同步压缩变换重构的幅值校正方法,包括以下步骤:
步骤一,对线性调频信号s(t)进行短时傅里叶变换(STFT),得到短时傅里叶变换系数Vs,其中信号的形式为s(t)=Asin(2πφ(t)),A为信号的幅值,t为时间参数,φ(t)为信号的相位函数;
步骤二,用得到的短时傅里叶变换系数Vs估算瞬时频率,即得到估计瞬时频率在该估计瞬时频率的基础上进行改进,用来计算准确瞬时频率其中φ”(τ)为相位函数的二阶导数,τ为过渡参数,不影响结果;
步骤三,将信号的所需要分析的频率范围设定为[ωmin,ωmax],将该范围划分为多个以ωk为中心频率,以Δω为间隔的频率段,通过ωk=(k-k0)Δω来计算中心频率,其中 Δω在不加大计算量的前提下尽可能的取小值;
步骤四,整合前三步计算出来的参数Vs、ωk和Δω,通过公式
完成同步压缩变换,即将所有对应同一频率段的准确瞬时频率的短时傅里叶变换系数转化为同步压缩系数,如果某几个准确瞬时频率对应第k个频率段,那么就将相应的所有短时傅里叶变换系数转化为和该频率段相对应的第k个同步压缩系数,从而得到同步压缩变换系数,公式中j是计数参数,在该公式中只有计数作用,不影响结果,为通过步骤二计算出来的第j个准确瞬时频率,Vs(t,ωj)为在时间为t和频率为ωj时的短时傅里叶变换系数,而ωj是属于所需分析的频率范围,为已知量;
步骤五,提取同步压缩变换的系数脊线,并将脊线附近的系数进行整合,得到相应的信号成分。
所述的步骤一中线性调频信号s(t)为单成分信号或多成分信号。
所述的步骤三中Δω的取值范围为0.1≤Δω<0.2。
本发明的有益效果为:由于本发明在步骤二中对瞬时频率的计算进行了改进,能够得到准确的瞬时频率,因此提高了时域重构信号的精度。
附图说明
图1为使用本发明方法对单成分信号进行信号重构,其中(a)、(b)和(c)图分别为仿真出的四阶相位信号、三阶相位信号和二阶相位信号。
图2为使用一般的同步压缩变换对单成分信号进行信号重构,其中(a)、(b)和(c)图分别为仿真出的四阶相位信号、三阶相位信号和二阶相位信号。
图3为使用本发明方法对多成分信号进行分解得到的各个成分信号与原始信号的对比图。
图4为使用一般的同步压缩变换对多成分信号进行分解得到的各个成分信号与原始信号的对比图。
具体实施方式
下面结合附图和实施例对本发明作进一步详细阐述。
实施例一
一种基于同步压缩变换重构的幅值校正方法,包括以下步骤:
步骤一,分别仿真不同阶相位的单成分信号,其中单成分信号分别为四阶相位信号s(t)=sin(2π(250t+100t4)),三阶相位信号s(t)=sin(2π(250t+100t3))和二阶相位信号s(t)=sin(2π(250t+100t2)),对线性调频信号s(t)进行短时傅里叶变换(STFT),得到短时傅里叶变换系数Vs;
步骤二,用得到的短时傅里叶变换系数Vs估算瞬时频率,即得到估计瞬时频率在该估计瞬时频率的基础上进行改进,用来计算准确瞬时频率其中φ”(τ)为相位函数的二阶导数,τ为过渡参数,不影响结果;
步骤三,将信号的所需要分析的频率范围设定为[ωmin,ωmax],将该范围划分为多个以ωk为中心频率,以Δω为间隔的频率段,通过ωk=(k-k0)Δω来计算中心频率,其中 Δω一般应该在不加大计算量的前提下尽可能的取小值,本实施例中为0.135;
步骤四,整合前三步计算出来的参数Vs、ωk和Δω,通过公式
完成同步压缩变换,即将所有对应同一频率段的准确瞬时频率的短时傅里叶变换系数转化为同步压缩系数,即如果某几个准确瞬时频率对应第k个频率段,那么就将相应的所有短时傅里叶变换系数转化为和该频率段相对应的第k个同步压缩系数,从而得到同步压缩变换系数,如图1中(a)、(b)、(c)中的SWFT图;其中j是计数参数,在该公式中只有计数作用,不影响结果,为通过步骤二计算出来的第j个准确瞬时频率,Vs(t,ωj)为在时间为t和频率为ωj时的短时傅里叶变换系数,而ωj是属于所需分析的频率范围,为已知量;
步骤五,提取同步压缩变换的系数脊线,并将脊线附近的系数进行整合,得到相应的信号成分,如图1中(a)、(b)、(c)中的重构与原始对比图。
由一般的同步压缩变换所得到的结果如图2所示,对比图1和图2,可以看出本发明方法相对一般的方法得到的重构信号与原始信号更加一致,幅值精度得到了极大提高。
实施例二
一种基于同步压缩变换重构的幅值校正方法,包括以下步骤:
步骤一,仿真一个双成分信号s=s 1+s 2,其中s1=sin(2π(250t+100t3)),s2=sin(2π(100t+100t2)),其中对信号s(t)进行短时傅里叶变换(STFT),得到短时傅里叶变换系数Vs;
步骤二,用得到的短时傅里叶变换系数Vs估算瞬时频率,即得到估计瞬时频率在该估计瞬时频率的基础上进行改进,用来计算准确瞬时频率其中φ”(τ)为相位函数的二阶导数,τ为过渡参数,不影响结果;
步骤三,将信号的所需要分析的频率范围设定为[ωmin,ωmax],将该范围划分为多个以ωk为中心频率,以Δω为间隔的频率段,通过ωk=(k-k0)Δω来计算中心频率,其中 Δω一般应该在不加大计算量的前提下尽可能的取小值,本实施例中取为0.135;
步骤四,整合前三步计算出来的参数Vs、ωk和Δω,通过公式
完成同步压缩变换,即将所有对应同一频率段的准确瞬时频率的短时傅里叶变换系数转化为同步压缩系数,即如果某几个准确瞬时频率对应第k个频率段,那么就将相应的所有短时傅里叶变换系数转化为和该频率段相对应的第k个同步压缩系数,从而得到同步压缩变换系数;其中j是计数参数,在该公式中只有计数作用,不影响结果,为通过步骤二计算出来的第j个准确瞬时频率,Vs(t,ωj)为在时间为t和频率为ωj时的短时傅里叶变换系数,而ωj是属于所需分析的频率范围,为已知量;
步骤五,提取同步压缩变换的系数脊线,并将脊线附近的系数进行整合,得到相应的信号成分,如图3所示。
由一般的同步压缩变换所得到的结果如图4所示,对比图3和图4,可以看出本发明方法相对一般的方法得到的重构信号与原始信号更加一致,幅值精度得到了极大提高。
Claims (3)
1.一种基于同步压缩变换重构的幅值校正方法,其特征在于,包括以下步骤:
步骤一,对线性调频信号s(t)进行短时傅里叶变换(STFT),得到短时傅里叶变换系数Vs,其中信号的形式为s(t)=Asin(2πφ(t)),A为信号的幅值,t为时间参数,φ(t)为信号的相位函数;
步骤二,用得到的短时傅里叶变换系数Vs估算瞬时频率,即得到估计瞬时频率在该估计瞬时频率的基础上进行改进,用来计算准确瞬时频率其中φ”(τ)为相位函数的二阶导数,τ为过渡参数,不影响结果;
步骤三,将信号的所需要分析的频率范围设定为[ωmin,ωmax],将该范围划分为多个以ωk为中心频率,以Δω为间隔的频率段,通过ωk=(k-k0)Δω来计算中心频率,其中 Δω在不加大计算量的前提下尽可能的取小值;
步骤四,整合前三步计算出来的参数Vs、ωk和Δω,通过公式
完成同步压缩变换,即将所有对应同一频率段的准确瞬时频率的短时傅里叶变换系数转化为同步压缩系数,如果某几个准确瞬时频率对应第k个频率段,那么就将相应的所有短时傅里叶变换系数转化为和该频率段相对应的第k个同步压缩系数,从而得到同步压缩变换系数,公式中j是计数参数,在该公式中只有计数作用,不影响结果,为通过步骤二计算出来的第j个准确瞬时频率,Vs(t,ωj)为在时间为t和频率为ωj时的短时傅里叶变换系数,而ωj是属于所需分析的频率范围,为已知量;
步骤五,提取同步压缩变换的系数脊线,并将脊线附近的系数进行整合,得到相应的信号成分。
2.根据权利要求1所述的一种基于同步压缩变换重构的幅值校正方法,其特征在于:所述的步骤一中线性调频信号s(t)为单成分信号或多成分信号。
3.根据权利要求1所述的一种基于同步压缩变换重构的幅值校正方法,其特征在于:所述的步骤三中Δω的取值范围为0.1≤Δω<0.2。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510821373.XA CN105446940B (zh) | 2015-11-23 | 2015-11-23 | 一种基于同步压缩变换重构的幅值校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510821373.XA CN105446940B (zh) | 2015-11-23 | 2015-11-23 | 一种基于同步压缩变换重构的幅值校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105446940A CN105446940A (zh) | 2016-03-30 |
CN105446940B true CN105446940B (zh) | 2018-07-03 |
Family
ID=55557157
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510821373.XA Expired - Fee Related CN105446940B (zh) | 2015-11-23 | 2015-11-23 | 一种基于同步压缩变换重构的幅值校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105446940B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109034042B (zh) * | 2018-07-20 | 2021-10-01 | 东北大学 | 基于广义线性调频双同步提取变换的非平稳信号处理方法 |
CN110617964A (zh) * | 2019-07-29 | 2019-12-27 | 中国铁道科学研究院集团有限公司城市轨道交通中心 | 用于滚动轴承故障诊断的同步压缩变换阶比分析法 |
CN111273345B (zh) * | 2020-03-05 | 2021-06-29 | 西南石油大学 | 一种基于高精度时频瞬时相位的地震数据时频谱处理方法 |
CN114674410A (zh) * | 2022-03-24 | 2022-06-28 | 安徽大学 | 一种分量数时变的水声信号瞬时频率估计方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101709997A (zh) * | 2009-12-11 | 2010-05-19 | 南京航空航天大学 | 一种振动信号处理的谐波窗函数 |
CN103455470A (zh) * | 2013-09-03 | 2013-12-18 | 上海交通大学 | 一种瞬时频率含交叉点的信号时频分解方法 |
CN104239702A (zh) * | 2014-09-03 | 2014-12-24 | 西安交通大学 | 基于克隆选择算法与改进快速s变换获取谐波参数的方法 |
CN104468427A (zh) * | 2014-10-25 | 2015-03-25 | 哈尔滨工业大学 | 基于fri时频域综合分析的信号高效采样及信号重构方法 |
CN104897962A (zh) * | 2015-06-19 | 2015-09-09 | 天津大学 | 基于互素感知的单频信号短样本高精度测频方法及其装置 |
-
2015
- 2015-11-23 CN CN201510821373.XA patent/CN105446940B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101709997A (zh) * | 2009-12-11 | 2010-05-19 | 南京航空航天大学 | 一种振动信号处理的谐波窗函数 |
CN103455470A (zh) * | 2013-09-03 | 2013-12-18 | 上海交通大学 | 一种瞬时频率含交叉点的信号时频分解方法 |
CN104239702A (zh) * | 2014-09-03 | 2014-12-24 | 西安交通大学 | 基于克隆选择算法与改进快速s变换获取谐波参数的方法 |
CN104468427A (zh) * | 2014-10-25 | 2015-03-25 | 哈尔滨工业大学 | 基于fri时频域综合分析的信号高效采样及信号重构方法 |
CN104897962A (zh) * | 2015-06-19 | 2015-09-09 | 天津大学 | 基于互素感知的单频信号短样本高精度测频方法及其装置 |
Non-Patent Citations (1)
Title |
---|
"基于自适应短时Chirp-Fourier变换的瞬时转速估计及应用";赵明等;《机械工程学报》;20150706;第51卷(第14期);第9-14页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105446940A (zh) | 2016-03-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105446940B (zh) | 一种基于同步压缩变换重构的幅值校正方法 | |
CN105783974B (zh) | 一种线性调频信号的检测、参数估计方法及系统 | |
CN107576943B (zh) | 基于瑞利熵的自适应时频同步压缩方法 | |
CN109946512B (zh) | 一种改进频域插值的动态功率分析方法 | |
CN105488341B (zh) | 一种基于混合经验模态分解的去噪方法 | |
CN102054480B (zh) | 一种基于分数阶傅立叶变换的单声道混叠语音分离方法 | |
CN104464725B (zh) | 一种唱歌模仿的方法与装置 | |
CN103575984A (zh) | 基于凯塞窗双谱线插值fft的谐波分析方法 | |
CN108694392A (zh) | 一种高精度同步提取广义s变换时频分析方法 | |
CN107608935A (zh) | 基于时间重排压缩变换的冲击类信号时频分析与重构方法 | |
CN101655520B (zh) | 对电力系统中雷击信号与暂态谐波信号的提取方法 | |
CN105137180B (zh) | 基于六项余弦窗四谱线插值的高精度谐波分析方法 | |
CN107478896A (zh) | 一种基于级联广义积分器的频率自适应谐波电流检测方法 | |
CN103018555B (zh) | 一种高精度的电力参数软件同步采样方法 | |
CN105845149B (zh) | 声音信号中主音高的获取方法及系统 | |
CN103529436A (zh) | 基于hht的无接触生命探测中呼吸和心跳信号的分离及时频分析方法 | |
CN103207319A (zh) | 数字化变电站电力信号非同步采样条件下的谐波测量方法 | |
CN105720983A (zh) | 用于时间交织模数转换系统的误差估计方法和装置 | |
CN101701985B (zh) | 定频变点电网谐波检测方法及其测量仪 | |
CN107204840A (zh) | 基于dft和迭代校正的正弦信号频率估计方法 | |
CN103399204A (zh) | 一种基于Rife-Vincent(II)窗插值FFT的谐波与间谐波检测方法 | |
CN105137181A (zh) | 基于Nuttall-Kaiser组合窗双谱线插值的谐波分析方法 | |
CN109655665A (zh) | 基于布莱克曼窗的全相位傅里叶谐波分析方法 | |
CN102930863B (zh) | 一种基于简化自适应内插加权谱模型的语音转换及重构方法 | |
CN110367986A (zh) | 基于ewt的呼吸信号近似熵特征提取方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180703 Termination date: 20211123 |
|
CF01 | Termination of patent right due to non-payment of annual fee |