[发明专利]一种心音信号医学指标数据的获取方法

摘要

本发明公开了一种心音信号医学指标数据的获取方法，使用双阈值限定来进行心音信号分段，解决了单一阈值分段存在的抗干扰能力差的问题，提高了心音信号分段的准确性，为后续可以准确的提取出心音信号的医学指标数据提供了保障。同时，在分段前对采集的心音信号进行尖峰信号消除，采用取对数，用线性的低通滤波器进行低通滤波得获得一个非常平滑的心音信号包络，这在对心音信号进行分段，确定心音中的各项成分是非常有帮助的。最后，根据心音信号的分段，得到心音信号的第一心音、第二心音、收缩期以及舒张期的持续时间，然后根据这些持续时间，就可以提取出心音信号的各项医学指标数据用作后续诊断使用。

主权项

1.一种心音信号医学指标数据的获取方法，其特征在于，包括以下步骤：/n(1)、心音信号采集与预处理/n首先使用电子听诊器以采样频率F_s采集人体心音信号，得到离散的心音信号x₀(n)，然后对所采集到的离散心音信号x₀(n)进行预处理：/n1.1)、分别使用截止频率为25Hz四阶巴特沃斯高通滤波器和截止频率为400Hz四阶巴特沃斯低通滤波器，对离散心音信号x₀(n)进行滤波处理，得到滤波后的心音x(n)；/n1.2)、对滤波后的心音信号x(n)进行尖峰消除处理，消除心音信号x(n)中的尖峰；/n1.3)、对消除尖峰的心音信号x(n)进行归一化处理，得到归一化的心音信号x(n)；/n(2)、心音信号包络特征提取/n对步骤(1)获得的归一化心音信号x(n)的绝对值取对数，得到心音信号ln|x(n)|，并通过一个线性的低通滤波器对心音信号ln|x(n)|进行低通滤波，得到除去心音信号中高频成分的信号L[ln|x(n)|]，将信号L[ln|x(n)|]作为e的指数，得到心音信号的同态包络(简称心音信号包络)Env_homomorphic(n)；/n(3)、心音信号分段/n3.1)、设定两个阈值thr₁与thr₂，其中：/nthr₁＝para₁*max(Env_homomorphic(n)),para₁∈[0.1,0.4]/nthr₂＝para₂*mean(Env_homomorphic(n)),para₂∈[0,0.1]/nmax()函数为求最大值函数，mean()函数为求平均值函数；/n3.2)、对步骤(2)提取到的心音信号包络Env_homomorphic(n)，按照时间顺序找出所有幅值大于阈值thr₁包络段的第一个点，并依次存入数组LargerThanThr₁中；/n3.3)、从数组LargerThanThr₁中依次逐个取点，从该点前后两个方向，心音信号包络Env_homomorphic(n)找出小于thr₂的第一个点，前向找到的点依次存入数组thr2end，后向找到的点依次存入数组thr2start；/n3.4)、计算时间间隔T₁、T₂：/nT1＝thr2start[2]-thr2end[1]/nT2＝thr2start[3]-thr2end[2]/n其中，thr2start[2]为数组thr2start的第二个点对应的时刻点，thr2start[3]为数组thr2start的第二个点对应的时刻点，thr2end[1]为数组thr2end的第一个点对应的时刻点，thr2end[2]为数组thr2end的第二个点对应的时刻点；/n如果T1<T2，则数组thr2start的奇数项为第一心音的起始点，偶数项为第二心音的起始点，数组thr2end数组奇数项为第一心音的终止点，偶数项为第二心音的终止点；/n如果T1>T2，则数组thr2start的奇数项为第二心音的起始点，偶数项为第一心音的起始点，数组thr2end奇数项为第二心音的终止点，偶数项为第一心音的终止点；/n3.5)、记录下第一、第二心音的起始与终止点，并根据第一、第二心音的起始与终止点得到心音信号中第一心音、第二心音以及收缩期、舒张期的确切位置，从而完成心音信号的分段工作；/n(4)、心音信号医学指标数据提取/n根据心音信号的分段，得到心音信号的第一心音、第二心音、收缩期以及舒张期的持续时间。/n