1.一种时域自适应加窗的胎心音降噪方法,其特征在于,包括:
步骤A,对第i秒原始胎心音数据di(n)进行滤波,滤除胎心音频段外的干扰信号,得到胎心音滤波数据fi(n),其中,n=1,2,…,Ⅳ,Ⅳ为1秒总采样点数;
步骤B,对滤波数据fi(n)取包络,得到包络数据ei(n);
步骤C,根据前M秒胎心率值判定是否开启时域自适应加窗,若开启加窗,则进入步骤D,若不开启加窗,则跳至步骤E;
步骤D,在包络数据ei(n)上进行自适应加窗,输出新包络数据:
步骤E,根据包络峰值间隔计算当前1秒心率值。
2.根据权利要求1所述的方法,其特征在于,所述步骤C开启自适应加窗的条件是:前M秒连续计算出在正常范围内的胎心率。
3.根据权利要求1所述的方法,其特征在于,所述步骤D具体包括:
子步骤D1:根据上1秒胎心率值,生成当前1秒两个矩形窗之间的参考间距L及间距的变化长度dH和dL;
子步骤D2:令矩形窗宽度为w,在当前1秒的包络上移动,移动步长为v,计算每个起始位置为lk处的矩形窗内包络的总能量pk,其中,k=1,2,…,K,K为总位置点数;
子步骤D3:将当前1秒包络数据划分成Q段,从pk中找出每段中能量最大值pq及其相应位置lq,q=1,2,…,Q;
子步骤D4:从Q段最大值pq中找出满足条件的两个最大值pmax1和pmax2,依次记录其矩形窗的起始位置lmax1和lmax2;
子步骤D5:以lmax1和lmax2为加矩形窗的起始点,分别在包络上加宽度为w的矩形窗,然后将在[lmax1+w,lmax2]范围内的包络置零,以及lmax1位置前和在lmax2+w位置后在lmax2-lmax1-w距离内的包络置零。
4.根据权利要求3所述的方法,其特征在于,所述子步骤D1中参考间距L、间距的变化长度dH和dL分别为:
其中,Ⅳ为1秒总采样点数,f0为上一秒胎心率,胎心率为f0表示每分钟f0次心跳,fv为胎心率变化值。
5.根据权利要求3所述的方法,其特征在于,所述子步骤D4选择两个最大值的必需满足以下三个条件:
条件a:lmax2-lmax1的值在[L-dL,L+dH]范围内;
条件b:lmax1在设定的起始范围内;
条件c:pmax1与pmax2的相互比值不超过2。