[发明专利]一种电子听诊器的摩擦音干扰自动消除方法

权利要求书

1.一种电子听诊器的摩擦音干扰自动消除方法，其特征在于，包括：

步骤一：读取缓存中时长为N的听诊信号采样序列x(n),n＝1,2,...,N，并将其表示为向量形式x；

步骤二：利用梅尔倒谱系数以及支持向量机，判定x上受摩擦音干扰的K个不交迭的干扰时间区间，按时间先后顺序依次为：

[n_K,begin,n_K,end],[n_K-1,begin,n_K-1,end],...,[n_1,begin,n_1,end]；

其中，K≥0；

步骤三：判断K的大小：若K＝0，执行步骤九；若K＞0，执行步骤四；

步骤四：将K个所述干扰时间区间，按后出现区间先进入后退出的方式形成堆栈；

步骤五：对x中第K个干扰时间区间上的数据x_K进行经验模态分解得到M_K个本征模函数分量：

步骤六：对IMF₂与IMF₃作低通滤波，然后分帧计算IMF₂与IMF₃的相关系数，确定相关系数大于预设阈值Th_c的P个不交迭的数据段；

按先后顺序将各所述数据段对应的数据点区间分别记为：对IMF₂与IMF₃上对应这些数据点区间以外的数据全置零；

步骤七：计算并依次对P个数据点区间的边缘数据段通过三次样条插值赋以更新值；

步骤八：更新x第K个干扰时间区间上的数据为：x(n_K，begin：n_k，end)＝x_K，移除堆栈顶部的时间区间[n_K,begin,n_K,end]，令K＝K-1，并返回步骤三；

步骤九：输出消除摩擦音干扰的听诊信号x。

2.根据权利要求1所述的一种电子听诊器的摩擦音干扰自动消除方法，其特征在于，步骤二中判定x上受摩擦音干扰的K个不交迭的干扰时间区间具体包括如下：

首先，把数据x以0.2s长度作为一个数据单元进行分割，且数据单元之间交迭0.1s，计算各数据单元的梅尔倒谱系数；

然后，将各数据单元上计算所得的MFCC向量c除以max(|c|)，并代入线性支持向量机f(c)＝sign(w^Tc+b)，其中w和b分别为该线性支持向量机的法向量和截距，当f(c)＞0时判断为该数据单元有摩擦音干扰；当f(c)＜0时判断为该数据单元无摩擦音干扰；

最后，若相近的两个检测出摩擦音干扰的数据单元的距离不大于0.1s，则将它们归并到同一个不间断的受摩擦音干扰的干扰时间区间，最终得到x上受摩擦音干扰的K个不交迭的干扰时间区间，按时间先后顺序依次为：[n_K,begin,n_K,end],[n_K-1,begin,n_K-1,end],...,[n_1,begin,n_1,end]。