本發(fā)明涉及分娩監(jiān)護(hù)領(lǐng)域，特別涉及一種基于分段濾波處理的子宮肌電強(qiáng)度信息提取方法。

背景技術(shù)：

目前，分娩力描記法(tocodynamometer，TOCO)作為一種檢測(cè)子宮收縮活動(dòng)的方法，被廣泛用于產(chǎn)前和產(chǎn)中的子宮收縮壓力檢測(cè)，但該方法在測(cè)量準(zhǔn)確性和靈敏度上存在一定局限性。子宮肌電信號(hào)作為一種最有可能替代傳統(tǒng)子宮收縮活動(dòng)監(jiān)護(hù)方法進(jìn)行子宮收縮狀態(tài)監(jiān)護(hù)。研究表明，該方法能夠同時(shí)滿足宮縮檢測(cè)的高靈敏度和非侵入性要求，用來(lái)監(jiān)測(cè)分娩過(guò)程中的子宮收縮壓力。

子宮肌電信號(hào)，即子宮平滑肌產(chǎn)生的電信號(hào)。目前，子宮肌電信號(hào)檢測(cè)主要通過(guò)放置在孕產(chǎn)婦腹部的電極進(jìn)行腹部子宮肌電信號(hào)的采集。宮縮，即子宮有規(guī)律的收縮，是臨產(chǎn)的一個(gè)重要特征。每一次宮縮可分為3個(gè)階段：上升期，峰值期和下降期，宮縮結(jié)束后進(jìn)入間歇期。臨床研究表明，宮縮壓力與子宮肌電信號(hào)強(qiáng)度具有較高的對(duì)應(yīng)性，宮縮期子宮肌電信號(hào)會(huì)出現(xiàn)明顯的爆發(fā)波，而進(jìn)入間歇期后子宮電信號(hào)趨于平緩。此外，通過(guò)對(duì)比研究宮縮壓力測(cè)量“金標(biāo)準(zhǔn)”—子宮內(nèi)壓力導(dǎo)管法(internal uterine pressure catheter，IUPC)檢測(cè)的子宮壓力與子宮肌電得到的子宮肌電強(qiáng)度，結(jié)果表明兩者具有較高的一致性。由于通過(guò)子宮肌電信號(hào)得到的子宮肌電信號(hào)包絡(luò)反映了子宮肌電強(qiáng)度信號(hào)，因此，利用子宮肌電信號(hào)包絡(luò)來(lái)預(yù)測(cè)宮縮壓力具有臨床應(yīng)用意義。

目前，在眾多子宮肌電強(qiáng)度信號(hào)提取方法中，對(duì)于子宮肌電信號(hào)中沒(méi)有明顯脈沖性噪聲的情況下，這些方法都能夠得到較好的子宮肌電強(qiáng)度信號(hào)，但當(dāng)子宮肌電信號(hào)中混疊有明顯的脈沖性噪聲時(shí)，利用這些方法得到的信號(hào)會(huì)出現(xiàn)由脈沖噪聲導(dǎo)致的虛假波峰，從而影響最后的檢測(cè)結(jié)果。為了克服上述問(wèn)題，本方法提出在識(shí)別出子宮肌電信號(hào)爆發(fā)波段與非爆發(fā)波段的基礎(chǔ)上對(duì)信號(hào)的不同波段采用不同的濾波處理方法，從而保證在濾除子宮肌電信號(hào)中混疊的脈沖性噪聲的同時(shí)，又能夠保留有效的子宮肌電信號(hào)成分，進(jìn)而得到準(zhǔn)確的子宮肌電強(qiáng)度信息，為臨床判斷提供可靠的參考，便于實(shí)際臨床應(yīng)用。

技術(shù)實(shí)現(xiàn)要素：

本發(fā)明的主要目的在于克服現(xiàn)有技術(shù)的缺點(diǎn)與不足，提供一種基于分段濾波處理的子宮肌電強(qiáng)度信息提取方法，濾除子宮肌電信號(hào)中混疊的脈沖性噪聲，從而得到更準(zhǔn)確的子宮肌電強(qiáng)度信息。

為了達(dá)到上述目的，本發(fā)明采用以下技術(shù)方案：

本發(fā)明的一種基于分段濾波處理的子宮肌電強(qiáng)度信息提取方法，包括以下步驟：

步驟A1、采集孕產(chǎn)婦的子宮肌電信號(hào)，記為原始子宮肌電信號(hào)；

步驟A2、識(shí)別上述子宮肌電信號(hào)的爆發(fā)波段與非爆發(fā)波段；

步驟A3、對(duì)所述子宮肌電信號(hào)不同信號(hào)波段進(jìn)行不同方式的濾波處理，即將子宮肌電信號(hào)劃分為爆發(fā)波段與非爆發(fā)波段，并分段進(jìn)行不同的濾波處理，得到濾波處理后的子宮肌電信號(hào)消除子宮肌電信號(hào)中的脈沖性噪聲干擾；

步驟A4、利用包絡(luò)提取算法提取所述子宮肌電強(qiáng)度信息。

作為優(yōu)選的技術(shù)方案，所述步驟A1中，所述子宮肌電信號(hào)是通過(guò)貼附在孕產(chǎn)婦腹部的Ag-AgCl表面電極采集得到。

作為優(yōu)選的技術(shù)方案，所述步驟A2中，識(shí)別爆發(fā)波段與非爆發(fā)波段的方法，包括以下步驟：

步驟B1、對(duì)所述子宮肌電信號(hào)進(jìn)行預(yù)處理；

步驟B2、由預(yù)處理后子宮肌電信號(hào)特點(diǎn)將其分為爆發(fā)波段與非爆發(fā)波段。

作為優(yōu)選的技術(shù)方案，所述步驟B1中，對(duì)于采集到的原始子宮肌電信號(hào)預(yù)處理包括如下步驟：

步驟C1、對(duì)采集到的子宮肌電信號(hào)進(jìn)行降采樣；

步驟C2、計(jì)算上述降采樣后信號(hào)中待求樣本點(diǎn)所在窗口的過(guò)零率。

作為優(yōu)選的技術(shù)方案，所述步驟B2中，子宮肌電信號(hào)爆發(fā)波與非爆發(fā)波段的判斷方法包括：

步驟D1、設(shè)定子宮肌電信號(hào)爆發(fā)波段過(guò)零率閾值；

步驟D2、比較權(quán)利要求3中計(jì)算待求樣本點(diǎn)所在窗口的過(guò)零率，若其過(guò)零率大于等于所述爆發(fā)波過(guò)零率設(shè)定閾值，則判斷該點(diǎn)處于爆發(fā)波段；反之處于非爆發(fā)波段。

作為優(yōu)選的技術(shù)方案，所述子宮肌電信號(hào)爆發(fā)波段與非爆發(fā)波段的判斷方法還包括：

步驟E1、計(jì)算子宮肌電信號(hào)基線；

步驟E2、通過(guò)判斷子宮肌電信號(hào)波形與子宮肌電信號(hào)基線之間的位置關(guān)系得出爆發(fā)波段與非爆發(fā)波段信號(hào)，當(dāng)子宮肌電信號(hào)波形處于該信號(hào)基線上方時(shí)，則認(rèn)為該信號(hào)處于爆發(fā)波段；反之處于非爆發(fā)波段。

作為優(yōu)選的技術(shù)方案，所述步驟A3中，對(duì)子宮肌電信號(hào)不同信號(hào)波段進(jìn)行不同方式的濾波處理，方法如下：

在相同濾波函數(shù)的條件下，對(duì)子宮肌電信號(hào)的爆發(fā)波段與非爆發(fā)波段進(jìn)行不同長(zhǎng)度的濾波窗口處理，即實(shí)現(xiàn)對(duì)子宮肌電信號(hào)進(jìn)行分段濾波處理，得到濾波處理后的子宮肌電信號(hào)。

作為優(yōu)選的技術(shù)方案，所述步驟A3中，對(duì)子宮肌電信號(hào)不同信號(hào)波段進(jìn)行不同方式的濾波處理，方法還包括：

對(duì)子宮肌電信號(hào)的爆發(fā)波段與非爆發(fā)波段分別使用不同長(zhǎng)度的濾波窗口與不同濾波函數(shù)的處理，得到濾波處理后的子宮肌電信號(hào)。

作為優(yōu)選的技術(shù)方案，步驟A4中，使用包絡(luò)提取算法從分段濾波處理后的子宮肌電信號(hào)中提取子宮肌電強(qiáng)度信號(hào)，具體為：

步驟F1、使用包絡(luò)提取算法提取子宮肌電強(qiáng)度信號(hào)；

步驟F2、為了抵消包絡(luò)提取過(guò)程中可能出現(xiàn)的頻譜泄露，對(duì)所述子宮肌電信號(hào)包絡(luò)進(jìn)行漢寧窗處理；

所述漢寧窗計(jì)算公式如下：

其中x_i為輸入信號(hào)，即第i點(diǎn)的子宮肌電強(qiáng)度信號(hào)幅值，y_i為當(dāng)前點(diǎn)對(duì)應(yīng)的輸出信號(hào)值，N為窗口長(zhǎng)度。

作為優(yōu)選的技術(shù)方案，所述步驟F1中，采用RMS法提取子宮肌電強(qiáng)度信號(hào)，在計(jì)算子宮肌電強(qiáng)度信號(hào)時(shí)需要選擇合適的窗口長(zhǎng)度L_RMS，且在信號(hào)的不同位置，其窗口長(zhǎng)度選擇也會(huì)存在差異。

本發(fā)明與現(xiàn)有技術(shù)相比，具有如下優(yōu)點(diǎn)和有益效果：

1、本發(fā)明通過(guò)分段濾波處理的子宮肌電強(qiáng)度信息提取方法，能夠有效去除脈沖性噪聲對(duì)于子宮肌電強(qiáng)度信號(hào)提取造成的影響，抑制子宮肌電強(qiáng)度信號(hào)中由脈沖性噪聲造成的虛假波峰，且得到的信號(hào)波形更加平滑，從而使得臨床檢測(cè)的準(zhǔn)確性和波形的可觀察性得到大大改善，降低臨床誤判率，進(jìn)而使得利用子宮肌電強(qiáng)度信息進(jìn)行宮縮活動(dòng)監(jiān)護(hù)成為可能，便于實(shí)際臨床應(yīng)用。

2、本發(fā)明有助于推動(dòng)子宮肌電強(qiáng)度信號(hào)檢測(cè)方法在臨床上的推廣使用，從而取代TOCO用以宮縮壓力檢測(cè)，為臨床上宮縮檢測(cè)提供一種更加準(zhǔn)確，舒適的宮縮檢測(cè)方法，保障分娩監(jiān)護(hù)正常，可靠的進(jìn)行。

附圖說(shuō)明

圖1是本發(fā)明子宮肌電強(qiáng)度信息特征提取算法流程圖；

圖2是本發(fā)明子宮肌電信號(hào)分為爆發(fā)波段與非爆發(fā)波段的示意圖；

圖3是本發(fā)明采集到的子宮肌電信號(hào)中混有脈沖性噪聲的示意圖；

圖4是含有虛假波峰的子宮肌電強(qiáng)度信號(hào)示意圖；

圖5是在濾波函數(shù)為中值濾波函數(shù)，不同長(zhǎng)度的濾波窗口下得到的子宮肌電信號(hào)與子宮肌電強(qiáng)度信號(hào)示意圖；

圖6是在濾波函數(shù)為均值濾波函數(shù)，不同長(zhǎng)度的濾波窗口下得到的子宮肌電信號(hào)與子宮肌電強(qiáng)度信號(hào)示意圖；

圖7是經(jīng)過(guò)不同長(zhǎng)度的濾波窗口和不同的濾波函數(shù)處理得到的子宮肌電信號(hào)和子宮肌電強(qiáng)度信號(hào)示意圖。

具體實(shí)施方式

下面結(jié)合實(shí)施例及附圖對(duì)本發(fā)明作進(jìn)一步詳細(xì)的描述，但本發(fā)明的實(shí)施方式不限于此。

實(shí)施例1

基于分段濾波處理的子宮肌電強(qiáng)度信息提取方法，如圖1所示，主要包括以下步驟：

步驟A1、采集孕產(chǎn)婦的子宮肌電信號(hào)，記為原始子宮肌電信號(hào)；

需要指出的是，這里的子宮肌電信號(hào)是通過(guò)貼附在孕產(chǎn)婦腹部的Ag-AgCl表面電極采集得到。

步驟A2、識(shí)別上述子宮肌電信號(hào)的爆發(fā)波段與非爆發(fā)波段；

步驟A3、對(duì)所述子宮肌電信號(hào)不同信號(hào)波段進(jìn)行不同方式的濾波處理，即將子宮肌電信號(hào)劃分為爆發(fā)波段與非爆發(fā)波段，并分段進(jìn)行不同的濾波處理，得到濾波處理后的子宮肌電信號(hào)消除子宮肌電信號(hào)中的脈沖性噪聲干擾；

步驟A4、利用包絡(luò)提取算法提取所述子宮肌電信息，比如其信號(hào)包絡(luò)，即子宮肌電強(qiáng)度信號(hào)。

作為本發(fā)明進(jìn)一步的技術(shù)方案，所述步驟A2包括以下步驟：

步驟B1、對(duì)所述子宮肌電信號(hào)進(jìn)行預(yù)處理；

作為本發(fā)明進(jìn)一步的技術(shù)方案，所述步驟B1包括以下步驟：

步驟C1、對(duì)采集到的子宮肌電信號(hào)進(jìn)行降采樣；

具體實(shí)施方案為：為了減小計(jì)算量，提高信號(hào)處理的速度，這里對(duì)采集到的原始子宮肌電信號(hào)進(jìn)行降采樣。降采樣的比值控制在一定范圍，從而不影響信號(hào)本身特征，這里降采樣比值設(shè)為1/100，這樣可以保證降采樣后的信號(hào)能夠反映信號(hào)的原有性質(zhì)。將所述子宮肌電信號(hào)在原有的幅值基礎(chǔ)上加上或減去一個(gè)微小正值P_val，即將子宮肌電信號(hào)整體向上抬高或降低，從而得到位置調(diào)整后的子宮肌電信號(hào)波形。

優(yōu)先的，調(diào)整幅值P_val的取值范圍為0.004～0.006mV，在本實(shí)施例中，幅值P_val為0.005mV。

步驟C2、計(jì)算上述降采樣后信號(hào)中待求樣本點(diǎn)所在窗口的過(guò)零率。

具體實(shí)施方案為：

首先，計(jì)算當(dāng)前樣本點(diǎn)所在窗口的過(guò)零點(diǎn)個(gè)數(shù)，計(jì)算公式為：

其中，x(i)為當(dāng)前窗口中第i個(gè)樣本點(diǎn)的幅值，M為窗口長(zhǎng)度，ZC為所述當(dāng)前窗口的過(guò)零點(diǎn)個(gè)數(shù)；

然后，計(jì)算當(dāng)前樣本點(diǎn)所在窗口的過(guò)零率ZCR，計(jì)算公式為：

ZCR＝ZC/M

步驟B2、根據(jù)預(yù)處理后子宮肌電信號(hào)特點(diǎn)將其分為爆發(fā)波段與非爆發(fā)波段。作為本發(fā)明進(jìn)一步的技術(shù)方案，所述步驟B2包括以下步驟：

步驟D1、設(shè)定子宮肌電信號(hào)爆發(fā)波段過(guò)零率閾值；

步驟D2、計(jì)算待求樣本點(diǎn)所在窗口的過(guò)零率，若其過(guò)零率大于等于所設(shè)閾值，則判斷其處于爆發(fā)波段，反之則處于非爆發(fā)波段；

所述閾值為：經(jīng)過(guò)大量實(shí)驗(yàn)統(tǒng)計(jì)得到的過(guò)零點(diǎn)閾值ZC_thd。

優(yōu)先的，閾值ZC_thd的取值范圍為0.05～0.06。在本實(shí)施例中，閾值ZC_thd為0.055。

如圖2所示，將子宮肌電信號(hào)分為爆發(fā)波段與非爆發(fā)波段。從圖中可以看出爆發(fā)波段的過(guò)零點(diǎn)個(gè)數(shù)明顯高于非爆發(fā)波段，且爆發(fā)波段的信號(hào)幅值要遠(yuǎn)高于非爆發(fā)波段的值。

作為本發(fā)明進(jìn)一步的技術(shù)方案，所述步驟A3包括以下步驟：

對(duì)所述子宮肌電信號(hào)爆發(fā)波段與非爆發(fā)波段信號(hào)使用相同的濾波函數(shù)，但兩種波段的濾波窗口長(zhǎng)度不同，方法如下：

步驟E1、對(duì)子宮肌電信號(hào)爆發(fā)波段使用指定長(zhǎng)度L₁濾波窗口進(jìn)行濾波處理；

步驟E2、對(duì)子宮肌電信號(hào)非爆發(fā)波段使用指定長(zhǎng)度L₂濾波窗口進(jìn)行濾波處理。

需要指出的是，本實(shí)施例采用中值濾波函數(shù)進(jìn)行不同信號(hào)波段的濾波處理，從而得到經(jīng)過(guò)分段濾波處理的子宮肌電信號(hào)。

如圖3所示，在采集到的子宮肌電信號(hào)中混有脈沖性噪聲，其分布于整個(gè)信號(hào)中；如圖5中頂圖所示，經(jīng)過(guò)分段濾波后，非爆發(fā)波段的脈沖基本上被消除；爆發(fā)波段在濾除脈沖性噪聲的同時(shí)，保留較為完整的爆發(fā)波有效成分。

具體實(shí)施方案為：對(duì)于子宮肌電信號(hào)中的爆發(fā)波段，采用窗口長(zhǎng)度為L(zhǎng)₁的中值濾波函數(shù)；對(duì)于子宮肌電信號(hào)中的非爆發(fā)波段，采用窗口長(zhǎng)度為L(zhǎng)₂的中值濾波函數(shù)。其中，當(dāng)前樣本點(diǎn)位于窗口中部，設(shè)窗口寬度為W(奇數(shù))時(shí)，則取該窗口當(dāng)前樣本點(diǎn)之前的(W-1)/2個(gè)點(diǎn)和之后的(W-1)/2個(gè)點(diǎn)。計(jì)算該窗口的中值作為當(dāng)前點(diǎn)的幅值。

優(yōu)先的，窗口L₁的范圍為9～17，窗口L₂的范圍為61～85。在本實(shí)施例中，窗口L₁為11，窗口L₂為71。

作為本發(fā)明進(jìn)一步的技術(shù)方案，所述步驟A3中采用的濾波函數(shù)還包括均值濾波函數(shù)等多種濾波函數(shù)，這里以均值濾波函數(shù)為例加以說(shuō)明。

該濾波函數(shù)處理與上述采用中值濾波函數(shù)操作基本相同，不同之處在于，這里采用的濾波函數(shù)為均值濾波函數(shù)。對(duì)子宮肌電信號(hào)的爆發(fā)波段與非爆發(fā)波段進(jìn)行不同窗口長(zhǎng)度的均值濾波，得到分段濾波處理的子宮肌電信號(hào)，如圖6中頂圖所示。對(duì)比圖5與圖6，我們可以看到使用兩種濾波函數(shù)均能夠有效濾波子宮肌電信號(hào)中混疊的脈沖性噪聲。

作為本發(fā)明進(jìn)一步的技術(shù)方案，所述步驟A4中，使用包絡(luò)提取算法從分段濾波處理后的子宮肌電信號(hào)中提取子宮肌電強(qiáng)度信號(hào)，步驟如下：

步驟F1、使用包絡(luò)提取算法提取子宮肌電強(qiáng)度信號(hào)；

需要指出的是，本實(shí)施例選擇子宮肌電強(qiáng)度信號(hào)提取方法為較常用的RMS法。詳細(xì)方法可參見(jiàn)：Horoba K,Jezewski J,Wrobel J,et al.Algorithm for detection of uterine contractions from electrohysterogram[C]//Engineering in Medicine and Biology Society,2001.Proceedings of the,International Conference of the IEEE.IEEE,2001:2161-2164vol.3.。對(duì)于利用RMS法提取子宮肌電強(qiáng)度信號(hào)，窗口長(zhǎng)度的選取尤為重要，長(zhǎng)度太大，則計(jì)算量迅速增加；長(zhǎng)度太小，則不能準(zhǔn)確反應(yīng)當(dāng)前點(diǎn)的幅值強(qiáng)度。因此，在計(jì)算子宮肌電強(qiáng)度信號(hào)時(shí)需要選擇合適的窗口長(zhǎng)度L_RMS，且在信號(hào)的不同位置，其窗口長(zhǎng)度選擇也會(huì)存在差異。

具體實(shí)施方案為：對(duì)于信號(hào)的前50個(gè)點(diǎn)而言，其窗口長(zhǎng)度可選為50；前50～100個(gè)點(diǎn)時(shí)，窗口長(zhǎng)度選為100；101至信號(hào)最后300個(gè)點(diǎn)之前的樣本點(diǎn)窗口長(zhǎng)度選為W_F，后300個(gè)點(diǎn)窗口為100。

優(yōu)先的，窗口長(zhǎng)度W_F范圍為360～450。在本實(shí)施例中，W_F為400。

步驟F2、為了抵消包絡(luò)提取過(guò)程中可能出現(xiàn)的頻譜泄露，對(duì)所述子宮肌電信號(hào)包絡(luò)進(jìn)行漢寧窗處理。

所述漢寧窗計(jì)算公式如下：

其中x_i為輸入信號(hào)，即第i點(diǎn)的子宮肌電強(qiáng)度信號(hào)幅值，y_i為當(dāng)前點(diǎn)對(duì)應(yīng)的輸出信號(hào)值，N為窗口長(zhǎng)度。

由分段濾波處理的子宮肌電信號(hào)提取的子宮肌電強(qiáng)度信號(hào)，如圖5中底圖所示。對(duì)比圖4與圖5可以發(fā)現(xiàn)，經(jīng)過(guò)分段濾波的子宮肌電強(qiáng)度信號(hào)能夠避免由脈沖性噪聲造成的虛假波峰，且波形更為平滑。

實(shí)施例2

本實(shí)施例與實(shí)施例1基本相同，不同之處僅在于待求樣本點(diǎn)所在窗口過(guò)零率的計(jì)算方法還包括：

通過(guò)判斷子宮肌電信號(hào)相鄰兩點(diǎn)之間幅值乘積的符號(hào)，同樣可以得到所要計(jì)算樣本點(diǎn)所在窗口的過(guò)零點(diǎn)數(shù)目及過(guò)零率，公式如下：

F＝S(n)S(n+1)

其中S(n)與S(n+1)為所要計(jì)算樣本點(diǎn)所在窗口相鄰兩個(gè)信號(hào)點(diǎn)幅值，為了計(jì)算當(dāng)前樣本點(diǎn)窗口的過(guò)零點(diǎn)數(shù)目，這里只需判斷F的正負(fù)號(hào)即可，如F為負(fù)，則將過(guò)零點(diǎn)數(shù)目N加1，最后通過(guò)將當(dāng)前窗口的過(guò)零點(diǎn)數(shù)目除以窗口總長(zhǎng)得到當(dāng)前樣本點(diǎn)所在窗口的過(guò)零率，這里同樣設(shè)窗口長(zhǎng)度為M，則過(guò)零率為N/M。

實(shí)施例3

本實(shí)施例與實(shí)施例1大致相同，不同之處在于步驟A2中識(shí)別爆發(fā)波段與非爆發(fā)波段的方法還包括：

通過(guò)計(jì)算所述子宮肌電信號(hào)的基線，再根據(jù)子宮肌電信號(hào)與其基線的位置關(guān)系，識(shí)別出爆發(fā)波段與非爆發(fā)波段。

子宮肌電信號(hào)基線的計(jì)算步驟如下：

a、選擇希爾伯特變換對(duì)所述子宮肌電信號(hào)進(jìn)行信號(hào)幅度轉(zhuǎn)變；

b、根據(jù)滑動(dòng)間隔的局部最小值選擇時(shí)間窗口；

c、計(jì)算對(duì)應(yīng)時(shí)間窗口中數(shù)據(jù)的中值；

d、將這些中值進(jìn)行擬合，得到所述子宮肌電信號(hào)的基線。

子宮肌電信號(hào)爆發(fā)波段與非爆發(fā)波段的識(shí)別如下：

若子宮肌電信號(hào)位于其基線的上部，則將其判定為爆發(fā)波段；反之為非爆發(fā)波段。

實(shí)施例4

本實(shí)施例與實(shí)施例1～2大致相同，不同之處僅在于步驟A3中對(duì)子宮肌電信號(hào)中不同波段的濾波處理方法還包括：

對(duì)所述子宮肌電信號(hào)爆發(fā)波段與非爆發(fā)波段使用不同長(zhǎng)度的濾波窗口和不同的濾波函數(shù)，步驟如下：

步驟G1、對(duì)子宮肌電信號(hào)爆發(fā)波段使用長(zhǎng)度為L(zhǎng)₃的濾波窗口，F(xiàn)un₁的濾波函數(shù)進(jìn)行濾波處理；

步驟G2、對(duì)子宮肌電信號(hào)非爆發(fā)波段使用長(zhǎng)度為L(zhǎng)₄的濾波窗口，F(xiàn)un₂的濾波函數(shù)進(jìn)行濾波處理。

需要指出的是，濾波窗口長(zhǎng)度和濾波函數(shù)的選擇要根據(jù)爆發(fā)波段與非爆發(fā)波段的特點(diǎn)，從而保證在濾除脈沖性噪聲的同時(shí)，又能夠最大化地保留子宮肌電信號(hào)有效成分。

優(yōu)先的，爆發(fā)波段濾波窗口長(zhǎng)度L₃范圍為9～17，非爆發(fā)波段濾波窗口長(zhǎng)度L₄范圍為61～85，濾波函數(shù)能夠有效濾除脈沖性噪聲，并最大程度保留子宮肌電信號(hào)成分。在本實(shí)施例中，L₃為11，L₄為71，F(xiàn)un₁為中值濾波函數(shù)，F(xiàn)un₂為均值濾波函數(shù)。

如圖7所示，經(jīng)過(guò)不同長(zhǎng)度濾波窗口和不同濾波函數(shù)處理的子宮肌電信號(hào)，其混疊的脈沖性噪聲基本被去除，得到準(zhǔn)確的子宮肌電強(qiáng)度信號(hào)，且較為平滑。

上述實(shí)施例為本發(fā)明較佳的實(shí)施方式，但本發(fā)明的實(shí)施方式并不受上述實(shí)施例的限制，其他的任何未背離本發(fā)明的精神實(shí)質(zhì)與原理下所作的改變、修飾、替代、組合、簡(jiǎn)化，均應(yīng)為等效的置換方式，都包含在本發(fā)明的保護(hù)范圍之內(nèi)。