CN104274211B - 一种基于延时估计的多普勒胎儿心音信号的提取方法 - Google Patents

Abstract

本发明公开一种基于延时估计的多普勒胎儿心音信号的提取方法，包括以下步骤：采集母亲腹部超声多普勒回波信号；分别取多普勒胎儿心音信号、亲子宫杂音信号和母亲腹部主动脉音信号的延时估计的先验信息；对先验信息做k阶信号差分运算，k为正整数；进行线性无关化处理；求解最优化问题，求得权值向量；分别建立范德蒙方程组并对其求解，确定各信号的延时；取母亲干扰信号的延时；采用基于斜投影的单通道盲分离方法将多普勒胎儿心音从采集的信号中提取出来，获得纯净的多普勒胎儿心音信号。本发明适用于胎儿健康状况检查中的多普勒胎儿心音信号的提取，能够有效分离出母亲干扰信号。

Description

一种基于延时估计的多普勒胎儿心音信号的提取方法

技术领域

本发明涉及产科胎儿监护技术领域，更具体地，涉及一种基于延时估计的多普勒胎儿心音信号的提取方法。

背景技术

胎儿心率监护是产科胎儿监护的重要方法，是避免胎儿窘迫，保障母亲和胎儿安全的重要手段。常见的胎儿心率监护方法有胎儿心电监护方法和胎儿心音监护方法。

胎儿心电监护法提取和分析胎儿心电现象的方法。具体的检测方法可分为两种：一是腹壁置电极的间接法。腹壁胎儿心电法是由母体腹壁置三个电极诱导胎儿心电的方法，此法会同时检测出母体心电信号和微弱的胎儿心电信号，二者混合在一起，不利于处理；二是电极直接置胎体的方法。直接胎儿心电法必须等到胎膜破裂或人工破膜后才能使用，不利于胎儿的前期监护和长期监护。

胎儿心音监护方法又可以分为采用心音换能器检测胎儿心音信号的方法和采用超声多普勒技术检测胎儿心音信号的方法。采用心音换能器检测胎儿心音信号的方法是基于胎儿心脏收缩、舒张的机械现象。胎儿心脏的一次搏动由心脏的收缩和舒张形成第一音和第二音。两个第一音或两个第二音的时间间隔是一个完整的心跳周期。但是，胎儿的心音信号很弱，杂音多，诸如胎动音、血管跳动音、母亲说话及开门声、脚步声等都会很轻易混入。往往由于心音探头位置稍有不当，心音间隔就不是从两个第一音或两个第二音检测出来，而是以前一次搏动的第一音和后一次搏动的第二音或前一次搏动的第二音和后一次搏动的第一音为心音间隔，从而导致测量的瞬时心率具有很大误差。基于超声多普勒技术的胎儿心音检测方法是利用超声探头产生的超声波束直接作用于胎儿心脏，由于胎儿心脏的运动会导致反射的超声波束产生多普勒效应，检测反射超声波的多普勒频移就可获得与胎儿心脏活动有关的心音信号。由于超声波具有很好的方向性,且信号的强度可控制,故采用超声探头检测来的胎心音信号强，杂音小，受宫缩影响小，并适用妊娠各阶段。超声信号的回波主要由以下几种成分组成：

(1)频率在0.01％至0.02％的载波频率处的胎儿心脏壁室的肌肉运动；

(2)频率在0.02％至0.05％的载波频率处的胎儿心脏瓣膜信号；

(3)频率在0.05至0.1％的载波频率处的胎儿动脉血流信号；

(4)子宫杂音；

(5)孕妇腹部主动脉音；

其中前三项构成了超声回波的主要部分和有用部分，称为多普勒胎儿心音信号。第四项和第五项与孕妇的脉搏同步，是超声回波信号的噪声成分。因此需要对超声回波信号进行一定的技术处理，尽量去除噪声信号成分，从而获得清晰的胎儿心音信号用于医学分析。然而，由于子宫杂音和孕妇腹部主动脉音基本与胎儿的多普勒心音频带重叠，无法用常规的频域滤波技术消除干扰，而仅通过调整传感器的位置避免噪声成分过强。

发明内容

本发明为克服上述现有技术所述的至少一种缺陷(不足)，提供一种能够有效分离母亲干扰信号的基于延时估计的多普勒胎儿心音信号的提取方法。

为解决上述技术问题，本发明的技术方案如下：

一种基于延时估计的多普勒胎儿心音信号的提取方法，包括以下步骤：

S1：采集母亲腹部超声多普勒回波信号，记为x；

S2：分别取一个心动周期的常规多普勒胎儿心音离散信号、一个心动周期的常规母亲子宫杂音离散信号和母亲腹部主动脉音离散信号，作为延时估计的先验信息，记为F(n)、U(n)和V(n)；

S3：分别对F(n)、U(n)和V(n)做k阶信号差分运算，k为正整数，并分别与F(n)、U(n)和V(n)一起构成一个矢量的集合，记为F₀、U₀和V₀；

S4：分别对集合F₀、U₀和V₀中的元素进行线性无关化处理，处理结果记为F、U和V；

S5：求解最优化问题min||x-α·F-β·U-γ·V||₂，求得权值向量α，β，γ，其中α＝[α₁，α₂，…,α_k+1]，β＝[β₁,β₂,…,β_k+1]，γ＝[γ₁，γ₂,…，γ_k+1]；

S6：分别建立范德蒙方程组并对其求解，确定各个胎儿心动周期的多普勒胎儿心音信号及各个母亲心动周期的母亲子宫杂音信号、母亲腹部主动脉音信号相对于信号起始处的延时；

S7：对于各个母亲心动周期的母亲子宫杂音信号和母亲腹部主动脉音信号相对于信号起始处的延时，取二者的算术平均作为母亲干扰信号的延时；

S8：根据多普勒胎儿心音信号和母亲干扰信号的延时，采用基于斜投影的单通道盲分离方法将多普勒胎儿心音从采集的信号中提取出来，获得纯净的多普勒胎儿心音信号；

基于斜投影的单通道盲分离方法参考专利CN102525452A。

在一种优选的方案中，步骤S3中，分别对F(n)、U(n)和V(n)做k阶信号差分运算的表达式为：

F₀＝{F(n)，F′(n),F″(n),F⁽³⁾(n),…,F^(k)(n)}

U₀＝{U(n),U′(n)，U″(n)，U⁽³⁾(n)，…,U^(k)(n)}

V₀＝{V(n),V′(n)，V″(n),V⁽³⁾(n),…，V^(k)(n)}

其中n＞＞k。

在一种优选的方案中，步骤S4中，分别对集合F₀、U₀和V₀中的元素进行线性无关化处理，其具体方法为：

F＝{F₁(n)，F₂(n)，…，F_k+1(n)}

＝{F(n)+aE₁(n),F′(n)+aE₂(n)，F″(n)+aE₃(n),F⁽³⁾(n)

+aE₄(n),…,F^(k)(n)+aE_k+1(n)}

U＝{U₁(n),U₂(n),…,U_k+1(n)}

＝{U(n)+bE₁(n),U′(n)+bE₂(n),U″(n)+bE₃(n),U⁽³⁾(n)

+bE₄(n)，…，U^(k)(n)+bE_k+1(n)}

V＝{V₁(n)，V₂(n)，…，V_k+1(n)}

＝{V(n)+cE₁(n),V′(n)+cE₂(n)，V″(n)+cE₃(n)，V⁽³⁾(n)

+cE₄(n)，…，V^(k)(n)+cE_k+1(n)}

其中，a，b，c为互不相等常实数，即a≠b，b≠c，c≠a，|a|,|b|,|c|分别小于集合F₀、U₀和V₀中元素绝对值的最大值的千分之一，即

$\left|a\right|<\frac{1}{1000}max\left\{\right|F\left(n\right)|,|F^{\&prime;}\left(n\right)|,|F^{\&prime;\&prime;}\left(n\right)|,|F^{\left(3\right)}\left(n\right)|,...,|F^{\left(k\right)}\left(n\right)\left|\right\}$

$\left|b\right|<\frac{1}{1000}max\left\{\right|U\left(n\right)|,|U^{\&prime;}\left(n\right)|,|U^{\&prime;\&prime;}\left(n\right)|,|U^{\left(3\right)}\left(n\right)|,...,|U^{\left(k\right)}\left(n\right)\left|\right\}$

$\left|c\right|<\frac{1}{1000}max\left\{\right|V\left(n\right)|,|V^{\&prime;}\left(n\right)|,|V^{\&prime;\&prime;}\left(n\right)|,|V^{\left(3\right)}\left(n\right)|,...,|V^{\left(k\right)}\left(n\right)\left|\right\}$

E_i(n),i＝1,2,…,k+1表示：n个数据组成的E_i(n)中第i个元素为1，其余均为0。

在一种优选的方案中，步骤S6中，求解多普勒胎儿心音延时的范德蒙方程组表示为：

$\left\{\begin{array}{c}{\mathrm{\&tau;}}_1+{\mathrm{\&tau;}}_2+...+{\mathrm{\&tau;}}_m={\mathrm{\&alpha;}}_1\\ {\mathrm{\&tau;}}_1^2+{\mathrm{\&tau;}}_2^2+...+{\mathrm{\&tau;}}_m^2=2!{\mathrm{\&alpha;}}_2\\ \begin{array}{c}.\\ .\\ .\end{array}\\ {\mathrm{\&tau;}}_1^{k+1}+{\mathrm{\&tau;}}_2^{k+1}+...+{\mathrm{\&tau;}}_m^{k+1}=\left(k+1\right)!{\mathrm{\&alpha;}}_{k+1}\end{array}\right.$

其中τ₁，τ₂，…，τ_m表示采集的母亲腹部超声多普勒回波信号中每一个胎儿心音波形的延时,m为正整数且m≤k+1；

求解母亲子宫杂音延时的范德蒙方程组表示为：

$\left\{\begin{array}{c}{\mathrm{\&delta;}}_1+{\mathrm{\&delta;}}_2+...+{\mathrm{\&delta;}}_p={\mathrm{\&beta;}}_1\\ {\mathrm{\&delta;}}_1^2+{\mathrm{\&delta;}}_2^2+...+{\mathrm{\&delta;}}_p^2=2!{\mathrm{\&beta;}}_2\\ \begin{array}{c}.\\ .\\ .\end{array}\\ {\mathrm{\&delta;}}_1^{k+1}+{\mathrm{\&delta;}}_2^{k+1}+...+{\mathrm{\&delta;}}_p^{k+1}=(k+1)!{\mathrm{\&beta;}}_{k+1}\end{array}\right.$

其中δ₁,δ₂,…,δ_p表示采集的超声多普勒回波信号中每一个母亲子宫杂音波形的延时,p为正整数且p≤k+1；

求解母亲腹部主动脉音延时的范德蒙方程组表示为：

$\left\{\begin{array}{c}{\mathrm{\&sigma;}}_1+{\mathrm{\&sigma;}}_2+...+{\mathrm{\&sigma;}}_p={\mathrm{\&gamma;}}_1\\ {\mathrm{\&sigma;}}_1^2+{\mathrm{\&sigma;}}_2^2+...+{\mathrm{\&sigma;}}_p^2=2!{\mathrm{\&gamma;}}_2\\ \begin{array}{c}.\\ .\\ .\end{array}\\ {\mathrm{\&sigma;}}_1^{k+1}+{\mathrm{\&sigma;}}_2^{k+1}+...+{\mathrm{\&sigma;}}_p^{k+1}=(k+1)!{\mathrm{\&gamma;}}_{k+1}\end{array}\right.$

其中σ₁,σ₂，…，σ_p表示采集的母亲腹部超声多普勒回波信号中每一个母亲腹主动脉音波形的延时，由于胎儿心律比母亲心律快，因此，超声多普勒回波信号中胎心音信号的数量比母亲干扰信号多，上述公式中，p＜m≤k+1。

在一种优选的方案中，步骤S7中，母亲干扰信号的延时表示为：

${\&Element;}_i=\frac{1}{2}({\mathrm{\&delta;}}_i+{\mathrm{\&sigma;}}_i),i=1,2,...,p.$

与现有技术相比，本发明技术方案的有益效果是：本发明利用多普勒胎儿心音先验信息和母亲干扰信号的先验信息，采用延时估计技术分别在采集信号中对多普勒胎儿心音信号和母亲干扰信号进行定位，即延时估计，并利用单通道盲分离技术将多普勒胎儿心音信号和母亲干扰信号进行分离，该方法简单且有效地将母亲干扰信号和多普勒胎儿心音信号进行分离，从而得到纯净的多普勒胎儿心音信号，提高了胎儿医学检测效率和准确率。

附图说明

图1为本发明基于延时估计的多普勒胎儿心音信号的提取方法的处理流程图。

图2为本发明利用延时估计技术确定超声多普勒胎儿心音信号、母亲子宫杂音信号和母亲腹部主动脉音信号延时的处理流程图。

图3为本发明通过超声多普勒探头采集的带有母亲干扰信号的胎心音信号。

图4为本发明所述方法提取的纯净超声多普勒胎心音信号。

图5为本发明所述方法分离的母亲子宫杂音和腹主动脉音混合噪声信号。

具体实施方式

附图仅用于示例性说明，不能理解为对本专利的限制；对于本领域技术人员来说，附图中某些公知结构及其说明可能省略是可以理解的。

下面结合附图和实施例对本发明的技术方案做进一步的说明。

实施例1

本发明为克服上述现有技术所述的至少一种缺陷(不足)，提供一种能够有效分离母亲干扰信号的基于延时估计的多普勒胎儿心音信号的提取方法。

为解决上述技术问题，本发明的技术方案如下：

一种基于延时估计的多普勒胎儿心音信号的提取方法，如图1所示，包括以下步骤：

S1：通过超声多普勒探头采集母亲腹部超声多普勒回波信号，记为x；

S2：分别取一个心动周期的常规多普勒胎儿心音离散信号、一个心动周期的常规母亲子宫杂音离散信号和母亲腹部主动脉音离散信号，作为延时估计的先验信息，记为F(n)、U(n)和V(n)；

S3：分别对F(n)、U(n)和V(n)做k阶信号差分运算，k为正整数，并分别与F(n)、U(n)和V(n)一起构成一个矢量的集合，记为F₀、U₀和V₀；

S4：分别对集合F₀、U₀和V₀中的元素进行线性无关化处理，处理结果记为F、U和V；

S5：求解最优化问题min||x-α·F-β·U-γ·V||₂，求得权值向量α，β，γ，其中α＝[α₁，α₂，…,α_k+1]，β＝[β₁,β₂，…，β_k+1]，γ＝[γ₁，γ₂,…,γ_k+1]；

S6：分别建立范德蒙方程组并对其求解，确定各个胎儿心动周期的多普勒胎儿心音信号及各个母亲心动周期的母亲子宫杂音信号、母亲腹部主动脉音信号相对于信号起始处的延时；

S7：对于各个母亲心动周期的母亲子宫杂音信号和母亲腹部主动脉音信号相对于信号起始处的延时，取二者的算术平均作为母亲干扰信号的延时；

S8：根据多普勒胎儿心音信号和母亲干扰信号的延时，采用基于斜投影的单通道盲分离方法将多普勒胎儿心音从采集的信号中提取出来，获得纯净的多普勒胎儿心音信号；

基于斜投影的单通道盲分离方法参考专利CN102525452A。

其中，如图2所示，步骤S1～S7为利用多普勒胎儿心音先验信息和母亲干扰信号的先验信息，采用延时估计技术分别在采集信号中对多普勒胎儿心音信号和母亲干扰信号进行延时估计，即定位；步骤S8为利用单通道盲分离技术将多普勒胎儿心音信号和母亲干扰信号进行分离。

在具体实施过程中，步骤S3中，分别对F(n)、U(n)和V(n)做k阶信号差分运算的表达式为：

F₀＝{F(n)，F′(n)，F″(n)，F⁽³⁾(n),…,F^(k)(n)}

U₀＝{U(n)，U^′(n)，U″(n),U⁽³⁾(n),…，U^(k)(n)}

V₀＝{V(n),V′(n),V″(n),V⁽³⁾(n),…,V^(k)(n)}

其中n＞＞k。

在具体实施过程中，步骤S4中，分别对集合F₀、U₀和V₀中的元素进行线性无关化处理，其具体方法为：

F＝{F₁(n)，F₂(n),…,F_k+1(n)}

＝{F(n)+aE₁(n),F′(n)+aE₂(n),F″(n)+aE₃(n),F⁽³⁾(n)

+aE₄(n),…,F^(k)(n)+aE_k+1(n)}

U＝{U₁(n)，U₂(n),…,U_k+1(n)}

＝{U(n)+bE₁(n)，U′(n)+bE₂(n)，U″(n)+bE₃(n)，U⁽³⁾(n)

+bE₄(n),…,U^(k)(n)+bE_k+1(n)}

V＝{V₁(n),V₂(n),…,V_k+1(n)}

＝{V(n)+cE₁(n),V′(n)+cE₂(n),V″(n)+cE₃(n),V⁽³⁾(n)

+cE₄(n),…,V^(k)(n)+cE_k+1(n)}

其中，a，b，c为互不相等常实数，即a≠b，b≠c，c≠a，|a|,|b|,|c|分别小于集合F₀、U₀和V₀中元素绝对值的最大值的千分之一，即

$\left|a\right|<\frac{1}{1000}max\left\{\right|F\left(n\right)|,|F^{\&prime;}\left(n\right)|,|F^{\&prime;\&prime;}\left(n\right)|,|F^{\left(3\right)}\left(n\right)|,...,|F^{\left(k\right)}\left(n\right)\left|\right\}$

$\left|b\right|<\frac{1}{1000}max\left\{\right|U\left(n\right)|,|U^{\&prime;}\left(n\right)|,|U^{\&prime;\&prime;}\left(n\right)|,|U^{\left(3\right)}\left(n\right)|,...,|U^{\left(k\right)}\left(n\right)\left|\right\}$

$\left|c\right|<\frac{1}{1000}max\left\{\right|V\left(n\right)|,|V^{\&prime;}\left(n\right)|,|V^{\&prime;\&prime;}\left(n\right)|,|V^{\left(3\right)}\left(n\right)|,...,|V^{\left(k\right)}\left(n\right)\left|\right\}$

E_i(n),i＝1,2,…,k+1表示：n个数据组成的E_i(n)中第i个元素为1，其余均为0。

在具体实施过程中，步骤S6中，求解多普勒胎儿心音延时的范德蒙方程组表示为：

$\left\{\begin{array}{c}{\mathrm{\&tau;}}_1+{\mathrm{\&tau;}}_2+...+{\mathrm{\&tau;}}_m={\mathrm{\&alpha;}}_1\\ {\mathrm{\&tau;}}_1^2+{\mathrm{\&tau;}}_2^2+...+{\mathrm{\&tau;}}_m^2=2!{\mathrm{\&alpha;}}_2\\ \begin{array}{c}.\\ .\\ .\end{array}\\ {\mathrm{\&tau;}}_1^{k+1}+{\mathrm{\&tau;}}_2^{k+1}+...+{\mathrm{\&tau;}}_m^{k+1}=\left(k+1\right)!{\mathrm{\&alpha;}}_{k+1}\end{array}\right.$

其中τ₁，τ₂，…，τ_m表示采集的母亲腹部超声多普勒回波信号中每一个胎儿心音波形的延时,m为正整数且m≤k+1；

求解母亲子宫杂音延时的范德蒙方程组表示为：

$\left\{\begin{array}{c}{\mathrm{\&delta;}}_1+{\mathrm{\&delta;}}_2+...+{\mathrm{\&delta;}}_p={\mathrm{\&beta;}}_1\\ {\mathrm{\&delta;}}_1^2+{\mathrm{\&delta;}}_2^2+...+{\mathrm{\&delta;}}_p^2=2!{\mathrm{\&beta;}}_2\\ \begin{array}{c}.\\ .\\ .\end{array}\\ {\mathrm{\&delta;}}_1^{k+1}+{\mathrm{\&delta;}}_2^{k+1}+...+{\mathrm{\&delta;}}_p^{k+1}=(k+1)!{\mathrm{\&beta;}}_{k+1}\end{array}\right.$

其中δ₁,δ₂，…，δ_p表示采集的超声多普勒回波信号中每一个母亲子宫杂音波形的延时,p为正整数且p≤k+1；

求解母亲腹部主动脉音延时的范德蒙方程组表示为：

$\left\{\begin{array}{c}{\mathrm{\&sigma;}}_1+{\mathrm{\&sigma;}}_2+...+{\mathrm{\&sigma;}}_p={\mathrm{\&gamma;}}_1\\ {\mathrm{\&sigma;}}_1^2+{\mathrm{\&sigma;}}_2^2+...+{\mathrm{\&sigma;}}_p^2=2!{\mathrm{\&gamma;}}_2\\ \begin{array}{c}.\\ .\\ .\end{array}\\ {\mathrm{\&sigma;}}_1^{k+1}+{\mathrm{\&sigma;}}_2^{k+1}+...+{\mathrm{\&sigma;}}_p^{k+1}=(k+1)!{\mathrm{\&gamma;}}_{k+1}\end{array}\right.$

其中σ₁，σ₂,…，σ_p表示采集的母亲腹部超声多普勒回波信号中每一个母亲腹主动脉音波形的延时，由于胎儿心律比母亲心律快，因此，超声多普勒回波信号中胎心音信号的数量比母亲干扰信号多，上述公式中，p＜m≤k+1。

在具体实施过程中，步骤S7中，母亲干扰信号的延时表示为：

${\&Element;}_i=\frac{1}{2}({\mathrm{\&delta;}}_i+{\mathrm{\&sigma;}}_i),i=1,2,...,p.$

本实施例通过超声多普勒探头采集的带有母亲干扰信号的胎心音信号如图3所示；本实施例提取的纯净超声多普勒胎儿心音信号如图4所示；本实施例分离的母亲子宫杂音和腹部主动脉音混合噪声信号如图5所示。

本实施例利用多普勒胎儿心音先验信息和母亲干扰信号的先验信息，采用延时估计技术分别在采集信号中对多普勒胎儿心音信号和母亲干扰信号进行定位，即延时估计，并利用单通道盲分离技术将多普勒胎儿心音信号和母亲干扰信号进行分离，该方法简单且有效地将母亲干扰信号和多普勒胎儿心音信号进行分离，从而得到纯净的多普勒胎儿心音信号，提高了胎儿医学检测效率和准确率。

显然，本发明的上述实施例仅仅是为清楚地说明本发明所作的举例，而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明权利要求的保护范围之内。

Claims (5)

1.一种基于延时估计的多普勒胎儿心音信号的提取方法，其特征在于，包括以下步骤：

S1：采集母亲腹部超声多普勒回波信号，记为x；

S2：分别取一个心动周期的常规多普勒胎儿心音离散信号、一个心动周期的常规母亲子宫杂音离散信号和母亲腹部主动脉音离散信号，作为延时估计的先验信息，记为F(n)、U(n)和V(n)；

S3：分别对F(n)、U(n)和V(n)做k阶信号差分运算，k为正整数，并分别与F(n)、U(n)和V(n)一起构成一个矢量的集合，记为F₀、U₀和V₀；

S4：分别对集合F₀、U₀和V₀中的元素进行线性无关化处理，处理结果记为F、U和V；

S5：求解最优化问题min‖x-α·F-β·U-γ·V‖₂，求得权值向量α,β,γ，其中α＝[α₁,α₂,…,α_k+1]，β＝[β₁,β₂,…,β_k+1]，γ＝[γ₁,γ₂,…,γ_k+1]；

S6：分别建立范德蒙方程组并对其求解，确定各个胎儿心动周期的多普勒胎儿心音信号及各个母亲心动周期的母亲子宫杂音信号、母亲腹部主动脉音信号相对于信号起始处的延时；

S7：对于各个母亲心动周期的母亲子宫杂音信号和母亲腹部主动脉音信号相对于信号起始处的延时，取二者的算术平均作为母亲干扰信号的延时；

S8：根据多普勒胎儿心音信号和母亲干扰信号的延时，采用基于斜投影的单通道盲分离方法将多普勒胎儿心音从采集的信号中提取出来，获得纯净的多普勒胎儿心音信号。

2.根据权利要求1所述的基于延时估计的多普勒胎儿心音信号的提取方法，其特征在于，步骤S3中，分别对F(n)、U(n)和V(n)做k阶信号差分运算的表达式为：

F₀＝{F(n),F′(n),F″(n),F⁽³⁾(n),…,F^(k)(n)}

U₀＝{U(n),U′(n),U″(n),U⁽³⁾(n),…,U^(k)(n)}

V₀＝{V(n),V′(n),V″(n),V⁽³⁾(n),…,V^(k)(n)}

其中n＞＞k。

3.根据权利要求2所述的基于延时估计的多普勒胎儿心音信号的提取方法，其特征在于，步骤S4中，分别对集合F₀、U₀和V₀中的元素进行线性无关化处理，其具体方法为：

F＝{F₁(n),F₂(n),…,F_k+1(n)}

＝{F(n)+aE₁(n),F′(n)+aE₂(n),F″(n)+aE₃(n),F⁽³⁾(n)

+aE₄(n),…,F^(k)(n)+aE_k+1(n)}

U＝{U₁(n),U₂(n),…,U_k+1(n)}

＝{U(n)+bE₁(n),U′(n)+bE₂(n),U″(n)+bE₃(n),U⁽³⁾(n)

+bE₄(n),…,U^(k)(n)+bE_k+1(n)}

V＝{V₁(n),V₂(n),…,V_k+1(n)}

＝{V(n)+cE₁(n),V′(n)+cE₂(n),V″(n)+cE₃(n),V⁽³⁾(n)

+cE₄(n),…,V^(k)(n)+cE_k+1(n)}

其中，a，b，c为互不相等常实数，|a|,|b|,|c|分别小于集合F₀、U₀和V₀中元素绝对值的最大值的千分之一，即

$\left|a\right|<\frac{1}{1000}max\left\{\right|F\left(n\right)|,|F^{\&prime;}\left(n\right)|,|F^{\&prime;\&prime;}\left(n\right)|,|F^{\left(3\right)}\left(n\right)|,...,|F^{\left(k\right)}\left(n\right)\left|\right\}$

$\left|b\right|<\frac{1}{1000}max\left\{\right|U\left(n\right)|,|U^{\&prime;}\left(n\right)|,|U^{\&prime;\&prime;}\left(n\right)|,|U^{\left(3\right)}\left(n\right)|,...,|U^{\left(k\right)}\left(n\right)\left|\right\}$

$\left|c\right|<\frac{1}{1000}max\left\{\right|V\left(n\right)|,|V^{\&prime;}\left(n\right)|,|V^{\&prime;\&prime;}\left(n\right)|,|V^{\left(3\right)}\left(n\right)|,...,|V^{\left(k\right)}\left(n\right)\left|\right\}$

E_i(n),i＝1,2,…,k+1表示：n个数据组成的E_i(n)中第i个元素为1，其余均为0。

4.根据权利要求1所述的基于延时估计的多普勒胎儿心音信号的提取方法，其特征在于，步骤S6中，求解多普勒胎儿心音延时的范德蒙方程组表示为：

$\left\{\begin{array}{c}{\mathrm{\&tau;}}_1+{\mathrm{\&tau;}}_2+...+{\mathrm{\&tau;}}_m={\mathrm{\&alpha;}}_1\\ {\mathrm{\&tau;}}_1^2+{\mathrm{\&tau;}}_2^2+...+{\mathrm{\&tau;}}_m^2=2!{\mathrm{\&alpha;}}_2\\ .\\ .\\ .\\ {\mathrm{\&tau;}}_1^{k+1}+{\mathrm{\&tau;}}_2^{k+1}+...+{\mathrm{\&tau;}}_m^{k+1}=(k+1)!{\mathrm{\&alpha;}}_{k+1}\end{array}\right.$

其中τ₁,τ₂,…,τ_m表示采集的母亲腹部超声多普勒回波信号中每一个胎儿心音波形的延时,m为正整数且m≤k+1；

求解母亲子宫杂音延时的范德蒙方程组表示为：

$\left\{\begin{array}{c}{\mathrm{\&delta;}}_1+{\mathrm{\&delta;}}_2+...+{\mathrm{\&delta;}}_p={\mathrm{\&beta;}}_1\\ {\mathrm{\&delta;}}_1^2+{\mathrm{\&delta;}}_2^2+...+{\mathrm{\&delta;}}_p^2=2!{\mathrm{\&beta;}}_2\\ .\\ .\\ .\\ {\mathrm{\&delta;}}_1^{k+1}+{\mathrm{\&delta;}}_2^{k+1}+...+{\mathrm{\&delta;}}_p^{k+1}=(k+1)!{\mathrm{\&beta;}}_{k+1}\end{array}\right.$

其中δ₁,δ₂,…,δ_p表示采集的超声多普勒回波信号中每一个母亲子宫杂音波形的延时,p为正整数且p<m≤k+1；

求解母亲腹部主动脉音延时的范德蒙方程组表示为：

$\left\{\begin{array}{c}{\mathrm{\&sigma;}}_1+{\mathrm{\&sigma;}}_2+...+{\mathrm{\&sigma;}}_p={\mathrm{\&gamma;}}_1\\ {\mathrm{\&sigma;}}_1^2+{\mathrm{\&sigma;}}_2^2+...+{\mathrm{\&sigma;}}_p^2=2!{\mathrm{\&gamma;}}_2\\ .\\ .\\ .\\ {\mathrm{\&sigma;}}_1^{k+1}+{\mathrm{\&sigma;}}_2^{k+1}+...+{\mathrm{\&sigma;}}_p^{k+1}=(k+1)!{\mathrm{\&gamma;}}_{k+1}\end{array}\right.$

其中σ₁,σ₂,…,σ_p表示采集的母亲腹部超声多普勒回波信号中每一个母亲腹主动脉音波形的延时。

5.根据权利要求4所述的基于延时估计的多普勒胎儿心音信号的提取方法，其特征在于，步骤S7中，母亲干扰信号的延时表示为：

${\&Element;}_i=\frac{1}{2}({\mathrm{\&delta;}}_i+{\mathrm{\&sigma;}}_i),i=1,2,...,p.$