一种基于同步压缩变换重构的幅值校正方法
【技术领域】
[0001] 本发明属于信号处理技术领域，具体设及一种基于同步压缩变换重构的幅值校正 方法。
【背景技术】
[0002] 时频分析方法是处理非平稳信号的有力工具，目前常用的有经典的时频方法，如 短时傅里叶变换和小波变换，近几年有学者通过运些经典的方法，得出了一种新的时频分 析方法：同步压缩变换。该方法类似经验模式分解(EMD)，可W进行信号分解，得到单成分信 号的时域图，即对信号进行重构。
[0003] 根据文献【latsenko D,Mcclintock P V E,Stefanovska A.Linear and synchrosqueezed time-frequency representations revisited:Overview,standards of use.resolution.reconstruction , concentration,and algorithms[J].Digital Signal Processing，2015，42.】中可W总结出一般的同步压缩变换，但是一般的同步压缩 变换中使用了瞬时频率的近似值，该近似值大大降低了时域重构信号的精度，所得到的信 号的幅值相比原始信号具有很大的差别。
【发明内容】
[0004] 为了克服上述现有技术的缺点，本发明的目的在于提供一种基于同步压缩变换重 构的幅值校正方法，校正所得到信号的幅值，提高时域重构信号的精度。
[000引为了达到上述目的，本发明采取的技术方案为：
[0006] -种基于同步压缩变换重构的幅值校正方法，包括W下步骤：
[0007] 步骤一，对线性调频信号s(t)进行短时傅里叶变换（STFT)，得到短时傅里叶变换 系数Vs，其中信号的形式为s(t)=Asin(2JT(i) (t))，A为信号的幅值，t为时间参数，Φ (t)为 信号的相位函数；
[000引步骤二，用得到的短时傅里叶变换系数Vs估算瞬时频率，即
得 到估计瞬时频率氣，在该估计瞬时频率的基础上进行改进，用ω =? +0 0旅一了)来计 算准确瞬时频率S，其中Φ"(τ)为相位函数的二阶导数，τ为过渡参数，不影响结果；
[0009] 步骤Ξ，将信号的所需要分析的频率范围设定为[Wmin，Wmax]，将该范围划分为多 个Wwk为中屯、频率，ΚΔ ω为间隔的频率段，通过c〇k=化-k〇)A ω来计算中屯、频率，其中
,Δ ω在不加大计算量的前提下尽可能的 取小值；
[0010] 步骤四，整合前Ξ步计算出来的参数V S、窃、ω k和Δ ω，通过公式
完成同步压缩变换，即将所有对应同一 频率段的准确瞬时频率的短时傅里叶变换系数转化为同步压缩系数，如果某几个准确瞬时 频率对应第k个频率段，那么就将相应的所有短时傅里叶变换系数转化为和该频率段相对 应的第k个同步压缩系数，从而得到同步压缩变换系数，公式中j是计数参数，在该公式中只 有计数作用，不影响结果，％为通过步骤二计算出来的第j个准确瞬时频率，Vs(t，W^为在 时间为t和频率为ω川寸的短时傅里叶变换系数，而ω J是属于所需分析的频率范围，为已知 量；
[0011] 步骤五，提取同步压缩变换的系数脊线，并将脊线附近的系数进行整合，得到相应 的信号成分。
[0012] 所述的步骤一中线性调频信号s(t)为单成分信号或多成分信号。
[001引所述的步骤ミ中Δ ω为0.巧化.2之间。
[0014] 本发明的有益效果为：由于本发明在步骤二中对瞬时频率的计算进行了改进，能 够得到准确的瞬时频率，因此提高了时域重构信号的精度。
【附图说明】
[0015] 图1为使用本发明方法对单成分信号进行信号重构，其中（a)、（b)和（C)图分别为 仿真出的四阶相位信号、Ξ阶相位信号和二阶相位信号。
[0016] 图2为使用一般的同步压缩变换对单成分信号进行信号重构，其中（a)、（b)和(C) 图分别为仿真出的四阶相位信号、Ξ阶相位信号和二阶相位信号。
[0017] 图3为使用本发明方法对多成分信号进行分解得到的各个成分信号与原始信号的 对比图。
[0018] 图4为使用一般的同步压缩变换对多成分信号进行分解得到的各个成分信号与原 始信号的对比图。
【具体实施方式】
[0019] 下面结合附图和实施例对本发明作进一步详细阐述。
[0020] 实施例一
[0021] -种基于同步压缩变换重构的幅值校正方法，包括W下步骤：
[0022] 步骤一，分别仿真不同阶相位的单成分信号，其中单成分信号分别为四阶相位信 号 s(t) = sin(23i(250t+100t4))，s阶相位信号 s(t) = sin(23i(250t+100t3))和二阶相位信 号s(t) = sin(化(250t+l(K)t2))，对线性调频信号s(t)进行短时傅里叶变换(STFT)，得到短 时傅里叶变换系数Vs;
[0023] 步骤二，用得到的短时傅里叶变换系数Vs估算瞬时频率，即
得 到估计瞬时频率耸，在该估计瞬时频率的基础上进行改进，用
农计 算准确瞬时频率谅，其中Φ"(τ)为相位函数的二阶导数，τ为过渡参数，不影响结果；
[0024] 步骤Ξ，将信号的所需要分析的频率范围设定为[Wmin，Wmax]，将该范围划分为多 个Wwk为中屯、频率，ΚΔ ω为间隔的频率段，通过c〇k=化-ko)A ω来计算中屯、频率，其中
，.Δ ω -般应该在不加大计算量的前提下尽 可能的取小值，本实施例中为0.135;
[002引步骤四，整合前Ξ步计算出来的参数Vs、?、Wk和Δ ω，通过公式
完成同步压缩变换，即将所有对应同一 频率段的准确瞬时频率的短时傅里叶变换系数转化为同步压缩系数，即如果某几个准确瞬 时频率对应第k个频率段，那么就将相应的所有短时傅里叶变换系数转化为和该频率段相 对应的第k个同步压缩系数，从而得到同步压缩变换系数，如图1中（a)、（b)、（c)中的SWFT 图；其中j是计数参数，在该公式中只有计数作用，不影响结果，鸣为通过步骤二计算出来的 第j个准确瞬时频率，Vs(t，Wj)为在时间为t和频率为Wj时的短时傅里叶变换系数，而COj 是属于所需分析的频率范围，为已知量；
[0026] 步骤五，提取同步压缩变换的系数脊线，并将脊线附近的系数进行整合，得到相应 的信号成分，如图1中(a)、(b)、（c)中的重构与原始对比图。
[0027] 由一般的同步压缩变换所得到的结果如图2所示，对比图1和图2,可W看出本发明 方法相对一般的方法得到的重构信号与原始信号更加一致，幅值精度得到了极大提高。
[002引实施例二
[0029] -种基于同步压缩变换重构的幅值校正方法，包括W下步骤：
[0030] 步骤一，仿真一个双成分信号 s = sl+s2,其中 sl = sin(化(250t+100t3))，s2 = sin (2JT(l〇〇t+l(K)t2))，其中对信号s(t)进行短时傅里叶变换（STFT)，得到短时傅里叶变换系 数Vs;
[0031] 步骤二，用得到的短时傅里叶变换系数Vs估算瞬时频率，即
得到 估计瞬时频率4，在该估计瞬时频率的基础上进行改进，用
来计算 准确瞬时频率ffl，:其中φ"(τ)为相位函数的二阶导数，τ为过渡参数，不影响结果；
[0032] 步骤Ξ，将信号的所需要分析的频率范围设定为[Wmin，Wmax]，将该范围划分为多 个Wwk为中屯、频率，ΚΔ ω为间隔的频率段，通过c〇k=化-k〇)A ω来计算中屯、频率，其中
，Δ ω-般应该在不加大计算量的前提下尽 可能的取小值，本实施例中取为0.135;
[0033] 步骤四，整合前Ξ步计算出来的参数Vs、孩、Wk和Δ ω，通过公式
完成同步压缩变换，即将所有对应同一 频率段的准确瞬时频率的短时傅里叶变换系数转化为同步压缩系数，即如果某几个准确瞬 时频率对应第k个频率段，那么就将相应的所有短时傅里叶变换系数转化为和该频率段相 对应的第k个同步压缩系数，从而得到同步压缩变换系数;其中j是计数参数，在该公式中只 有计数作用，不影响结果，兩;为通过步骤二计算出来的第j个准确瞬时频率，Vs(t，W^为在 时间为t和频率为ω川寸的短时傅里叶变换系数，而ω J是属于所需分析的频率范围，为已知 量；
[0034] 步骤五，提取同步压缩变换的系数脊线，并将脊线附近的系数进行整合，得到相应 的信号成分，如图3所示。
[0035] 由一般的同步压缩变换所得到的结果如图4所示，对比图3和图4,可W看出本发明 方法相对一般的方法得到的重构信号与原始信号更加一致，幅值精度得到了极大提高。
【主权项】
1. 一种基于同步压缩变换重构的幅值校正方法，其特征在于，包括以下步骤： 步骤一，对线性调频信号S(t)进行短时傅里叶变换（STFT)，得到短时傅里叶变换系数 Vs，其中信号的形式为s(t)=A Sin(23i(Ht))，A为信号的幅值，t为时间参数，Φ(〇为信号 的相位函数； 步骤二，用得到的短时傅里叶变换系数1估算瞬时频率，即，得到估计 瞬时频率?,,在该估计瞬时频率的基础上进行改进，用丨来计算准确瞬 时频率其中φ〃（τ)为相位函数的二阶导数，τ为过渡参数，不影响结果； 步骤三，将信号的所需要分析的频率范围设定为[?min，comax]，将该范围划分为多个以 wk为中心频率，以Δ ω为间隔的频率段，通过cok= (k-k〇) Δ ω来计算中心频率，其中，△ω在不加大计算量的前提下尽可能的 取小值； 步骤四，整合前三步计算出来的参数Vs、孩:、wk和Δ ω，通过公式完成同步压缩变换，即将所有对应同 一频率段的准确瞬时频率的短时傅里叶变换系数转化为同步压缩系数，如果某几个准确瞬 时频率对应第k个频率段，那么就将相应的所有短时傅里叶变换系数转化为和该频率段相 对应的第k个同步压缩系数，从而得到同步压缩变换系数，公式中j是计数参数，在该公式中 只有计数作用，不影响结果，％为通过步骤二计算出来的第j个准确瞬时频率，V s(t，〇^)为 在时间为t和频率为时的短时傅里叶变换系数，而是属于所需分析的频率范围，为已 知量； 步骤五，提取同步压缩变换的系数脊线，并将脊线附近的系数进行整合，得到相应的信 号成分。2. 根据权利要求1所述的一种基于同步压缩变换重构的幅值校正方法，其特征在于:所 述的步骤一中线性调频信号s(t)为单成分信号或多成分信号。3. 根据权利要求1所述的一种基于同步压缩变换重构的幅值校正方法，其特征在于:所 述的步骤三中Δ ω为0.1到0.2之间。
【专利摘要】一种基于同步压缩变换重构的幅值校正方法，对线性调频信号进行短时傅里叶变换，得到短时傅里叶变换系数，然后估算瞬时频率，在该估计瞬时频率的基础上进行改进，计算准确瞬时频率，再计算中心频率，然后完成同步压缩变换，即将所有对应同一频率段的准确瞬时频率的短时傅里叶变换系数转化为同步压缩系数，最后提取同步压缩变换的系数脊线，并将脊线附近的系数进行整合，得到相应的信号成分，由于本发明对瞬时频率的计算进行了改进，能够得到准确的瞬时频率，因此提高了时域重构信号的精度。
【IPC分类】G06F17/14
【公开号】CN105446940
【申请号】CN201510821373
【发明人】廖与禾, 李帅, 陈鹏, 孙鹏, 林京, 雷亚国
【申请人】西安交通大学
【公开日】2016年3月30日
【申请日】2015年11月23日