参数化时频分析理论方法及应用_2016_西北工大.pptxVIP

西北工业大学第三届杰青论坛参数化时频分析理论、方法及应用报告人: 彭志科上海交通大学机械系统与振动国家重点实验室2016.3.26-27非平稳信号是实际应用中的常见信号，对其分析的准确性是解决实际工程问题系统的前提与基础，也是信号处理领域中的热点。体波、表面波超声波效应生命医学信号旋转机械振动信号地震波多普勒效应风力发电机脑电信号兰姆波心电图微多普勒效应水轮机研究背景研究背景Fourier级数Fourier变换 - King of transforms Jean Baptiste Joseph Fourier (1768 - 1830) 本质思想：线性空间的正交基分解和重构适用对象：平稳信号(频率不随时间变化)研究背景非平稳信号 特点：非平稳信号的频率常随时间变化例1：(0≤ t ≤15)采样频率120HzFourier变换能反映信号的频率范围不能反映频率随时间变化的变化规律非平稳信号 例2：x2x1Fourier变换准确反映信号所含频率分量不能反映频率分量存在的时间段时频分析方法 短时傅立叶变换 定义：a=1窗函数本质：加窗傅立叶变换适用对象：分段平稳信号Freq/Hz例2- x1例2- x2t/sect/sec连续小波变换 定义：ω0 = 5母波函数本质：变分辨率带通滤波适用对象：局部奇异性信号示例转子碰摩故障典型时频特征a = 0.5a = 1Wigner-Ville分布 定义：本质：瞬时相关函数的傅立叶变换适用对象：单分量信号交叉项示例单分量信号x1单分量信号x2多分量信号x线调频小波变换 定义：调频窗函数本质：加调频窗的傅立叶变换STFT适用对象：线性调频信号示例脉冲反射波信号简化模型反射波Chirplet时频分析方法的不足 例1信号真实时频曲线STFTCWT不足不能正确反映分量的瞬时幅值集中性较差Chirplet(α = 3)WVD非线性调频分量 特点：频率是时间的非线性函数(A)(B)(C)Opt Exp, 19(2011) 26174J. Ac. Soc. Am. 107(2000), Pt.1 App. Ac. 71 (2010) 1070–1080(E)(D)激光脉冲信号Lamb波信号鲸鱼声波水轮机停机振动信号战斗机机翼测试信号J S. Vib. 330 (2011)1225–1243J. G. Con. DY 21(1998)375-382参数化时频分析-原理 线调频信号正弦信号STFTSTFT参数化时频分析-原理 Chirplet 定义新表达Chirplet 工作原理?平移旋转算子旋转平移算子t0参数化时频分析-原理 例3：Chirp Rate = 5π / s α = 0 α = 5π α = 2.5π α = -2.5π 参数化时频分析-原理 信号模型:瞬时频率:参数化时频分析原理图解旋转算子平移算子参数化时频分析-方法 参数化频率时延分析 工作原理图解定义参数化时频分析-方法 性质定义线性可加时移不变频移不变参数化时频分析-方法 多项式调频小波变换 (PCT)旋转算子平移算子思想：闭区间上的连续函数可 用多项式函数一致逼近STFTPCT示例参数化时频分析-方法 示例- PCTCWTSTFTPCTChirplet(α = 3)WVD参数化时频分析-方法 样条调频小波变换 (SCT)旋转算子平移算子思想：分段多项式逼近，避免Runge现象PCTSCT示例参数化时频分析-方法 示例-SCTSCTCWTSTFTPCTWVD参数化时频分析-方法 广义Warblet变换 (GWT)旋转算子平移算子思想：任何可积函数都可 用三角函数展开示例STFTGWT参数化时频分析-方法 示例- GWTSTFTCWTWVDGWT参数化时频分析-方法 如何确定变换核参数多项式调频小波变换样条调频小波变换广义Warblet变换参数化时频分析-方法 变换核参数估计-PCT最小二乘法拟合++++++++++++++++++++++++++++++++STFTPCT最小二乘法拟合参数化时频分析-方法 变换核参数估计-GWTFourier变换STFTGWTFourier变换参数化时频分析-方法 如何分析多频率分量信号？示例- 多频率分量信号IF1IF2参数化时频分析-方法 多频率分量信号分析-时频融合法ψ1IF1IF2ψ2参数化时频分析-方法 多频率分量信号分析-时频融合法高通滤波高通滤波参数化时频分析-方法 时频融合法示例蝙蝠回波定位信号STFTPWVDPCTChirplet参数化时频分析-方法 多频率分量信号分析-信号分解法信号模型:问题:Fourier级数参数化时频分析-方法 多频率分量信号分析-信号分解法示例信号A-B-C-Ind:(A)-0.1007 ; (B)-0.4650; (C)-0.6060 参数化时频分析-方