心音信号短时傅立叶变换分析.pdfVIP

维普资讯 2004年 8月 重 庆 大 学 学 报 Aug．2004 第27卷第8期 JournalofChongqingUniversity Vo1．27 No．8 文章编号：1000—582X(2004)08—0083—03 ． 心音信号的短时傅立叶变换分析 ‘ 朱 冰 莲，杨 磊 (重庆大学 通信工程学院，重庆 400030) 摘 要：心音信号是一种典型的非平稳信号，短时傅立叶变换( ，又称窗口傅立叶变换)是用于 对非平稳信号进行时频分析的有效工具。在用STFT分析心音信号时，窗口宽度的选择是非常重要的。 一 方面，要使用尽可能窄的窗口来保证信号的局部平稳性，另一方面，要选用较宽的窗口来提高频率分 辨率。如何处理这一矛盾就成了问题的关键。通过调整滑动时间窗的宽度，达到了比较满意的效果。 首先用窄时间窗进行分析，频谱图具有高时间分辨率，得到了心动周期等时域特征参数。进而逐渐加宽 时间窗，最后得到高频率分辨率的频谱图，可以看出各心音成分的频谱特征。实验结果表明，不同时频 尺度的STFr分析可以很好地描述正常心音信号的时域和频域特征。 关键词：心音信号；短时傅立叶变换；时频分析 中图分类号：TN911．72 文献标识码：A 人体的心脏在一个心动周期会产生4个心音成 为了在时频领域表征心音成分sl和S2，采用短 分，在心音图上分别标识为Sl，S2，S3，．s4。由于S3和 时傅立叶变换(S，IF-I．)对一个正常的心音样本信号进 ．s4的强度较低，通常情况下仅能察觉到Sl和 ，如图 行了分析。实验结果表明，通过适当地调整 STFr的 l所示。Sl的频率范围一般在90Hz至 180Hz之间， 参数(窗口的类型和宽度)，可以比较准确地得到心音 它是在心室收缩早期，房室瓣关闭到半月瓣开放期间， 信号的时频表示。当选用窄窗进行分析时，时间分辨 由于瓣膜和心肌的振动以及血流的加速而形成的。 率高，可得到心音信号的时域能量包络图，由此可以准 s2是在心室的舒张期，半月瓣关闭到房室瓣开放期间 确地量度心音成分的时域特征参数。而选用宽窗进行 形成的，频率可能达到200Hzn1引。 分析时，便得到每个心音成分(sl和S2)的频谱，由此 } S 可以以高的频率分辨率识别出它们的组成。 1 短时傅立叶变换 (STFT) ■■_I … 。． … ． ．．．． 一 一 … ～ ． … JJ …．一． …’…一 l一… …… 。。lr’…… 短时傅立叶变换 (STFT)是建立在传统傅立叶变 换基础之上的，其基本思想是引入一个具有时频局域 性的窗函数 (t)，让它沿着信号s(t)滑动，对每一段 被窗口截取的部分作傅立叶变换 4【】。由于滑动窗的 图 1 正常心音图样本 位置引入了时间的信息，这样就可以得到一个时变的 因为心音各成分的持续时间很短且频率较低，人 频域分析结果。可以用下面的表达式来说明SrF-I： 耳是不能分辨出它们的时间和频率特征的。然而，尽 S(t =∑s(r)w(一t)e书 (1) 管心音信号具有时变的特性，仍然可以用现代数字信