[发明专利]基于时频分析的多普勒胎心音的瞬时心率检测方法及装置

权利要求书

1.一种基于时频分析的多普勒胎心音的瞬时心率检测方法，其特征在于，包括如下步骤：

S1.信号预处理：对采集到的多普勒胎心音进行带通滤波，带通滤波器的通带频率为f_L到f_H；

S2.时频分析：对于经过步骤S1预处理的多普勒胎心音进行时频分析，得到时频图；

S3.选取特征频段和模板：从时频图中选取信号中的特征频段，特征频段为f_cL到f_cH；从时频图选取包含S1音和S2音特征的时频块，所述时频块的时间间隔为0.2(秒)＜t₀＜0.5(秒)；

S4.求互相关函数：对特征频段和模板求互相关函数，根据互相关函数的结果绘出相关曲线；

S5.用阈值检测法求互相关曲线的峰值；

S6.对峰值做差分，求瞬时心率值；根据瞬时心率值绘制出瞬时心率图。

2.根据权利要求1所述的基于时频分析的多普勒胎心音的瞬时心率检测方法，其特征在于，所述步骤S1中f_L为50Hz，f_H为250Hz，即采用带通频带为50-250Hz的带通滤波器。

3.根据权利要求1所述的基于时频分析的多普勒胎心音的瞬时心率检测方法，其特征在于，所述步骤S2采用短时傅里叶变换STFT进行时频分析，短时傅里叶变换STFT定义如下：

s(w,t)=12π∫-∞+∞e-iwx(τ)h(τ-t)dτ---(1)]]>

其中h(t)是窗函数，x(τ)是信号，τ为信号自变量，t是时间变量，w为频率自变量，即沿时间轴移动分析窗，得到二维的时频图由s(w,t)表示。

4.根据权利要求1所述的基于时频分析的多普勒胎心音的瞬时心率检测方法，其特征在于，所述步骤S3中特征频段取200-400Hz，即f_cL为200Hz，f_cH为400Hz。