TWI766783B - 訊號的波峰與波谷之偵測方法 - Google Patents

Abstract

本發明係一種訊號的波峰與波谷之偵測方法，包括下列步驟，波形中當前的採樣時間點之訊號值做為第一訊號值，前一個採樣時間點的訊號值作為第二訊號值，後一個採樣時間的訊號值作為第三訊號值，根據第一訊號值、第二訊號值及第三訊號值，將第一訊號值分類為波峰點、波谷點、第一假性波峰點、第二假性波峰點、第一假性波谷點、第二假性波谷點，其中第一假性波峰點或第一假性波谷點，再與後續的第二假性波峰點或第二假性波谷點確認出波峰點或波谷點，而且最終只有波峰點或波谷點被儲存在儲存單元中。

Description

訊號的波峰與波谷之偵測方法

本發明係有關於訊號處理中，尤指一種從訊號的波形中偵測出波峰與波谷的方法。

人體生理參數監測包括對心電、血壓、心率、血氧飽和度、呼吸頻率、體溫、呼吸中的二氧化碳濃度等，並應用於臨床醫學、可穿戴設備、移動終端等領域。例如，目前成熟的心率測量方法有脈搏測量法、心電信號法Electrocardiogram, ECG)、光電容積脈搏波描記法(Photoplethysmography，PPG)等。

在傳統的心電信號法的檢測過程中，通常採用標準的十二級導程形式、單級導程形式以及其他改進的導形式，用以量測竇房結有節律地控制心臟收縮舒張，從而向軀幹泵血所產生的電信號。這個電信號會逐漸擴散到體表，可以在皮膚通過電極測量。目前醫院使用的心電儀就是採用這個原理。這個電信號也被稱為心電信號，心電信號提供很多資訊可以提供給醫生診斷，而此心電信號在進行分析的過程中，有些心電信號的特徵需要偵測局部最低與局部最高之信號轉折點，也就是要找出心電信號之波形中的波峰及波谷的位置。

除了心電信號以外，許多的波形分析當中，也是需要找出波形中的波峰及波谷的位置，而傳統尋找波形中的波峰或波谷的方法而言，例如：美國第US5987392號專利案(發明名稱：Wave form peak detector，以下稱美國專利前案)，峰值檢測是通過將連續波形轉換為離散樣本來進行的。評估離散樣本以識別連續波形峰值附近的三個點。在識別出峰值附近的三個點後，將它們存儲在離散域數據結構中，該數據結構與連續變換函數生成的表格結合使用，以識別準確的峰值幅度和時間。

又，例如中國第CN105637360B號發明專利案(發明名稱：波形中的峰值端點檢測方法及檢測裝置，以下稱中國專利前案)，係一種波形中的峰值端點檢測方法及檢測裝置。其中的檢測方法包括如下步驟：a)獲取峰值檢測對象波形的拐點提取波形，其中所述拐點提取波形為峰值檢測對象波形的二階微分波形；b)獲取所述峰值檢測物件波形的峰頂位置；c)在所述峰值檢測物件波形的所述峰頂位置的左右側中的任一側，檢測所述拐點提取波形的極大值及獲得所述極大值的橫軸的位置；d)所述極大值乘以具有規定的0至1之間的值的相對閾值率而計算相對閾值，其中所述相對閾值率是基於根據所述峰值檢測物件波形的峰值強度、峰值寬度計算出的高斯函數計算而得到；e)以獲得所述極大值的橫軸的位置為基點，在遠離與所述峰頂位置相對應的橫軸的位置的方向上，檢測出所述拐點提取波形的值下降至所述相對閾值為止的最早的點作為峰值端點的橫軸位置，並檢測出與所述峰值端點的橫軸位置相對應的所述峰值檢測對象波形上的點作為峰值端點。

前述的美國專利前案及中國專利前案都需要複雜的數學計算過程，使得波形的峰值(波峰或波谷)的偵測方法顯得相當的複雜且困難，而且這些數學計算過程需要使用複雜的電路才可以完成。

然而，波形中的峰值也是可以從簡單的峰值偵測方法取得局部最低與局部最高之信號轉折點，也就是要找出波形中的波峰及波谷的位置，例如：當波形中的被採樣的訊號值，同時大於前一個採樣時間、下一個採樣時間上前後被採樣的訊號值，即為波形中的波峰位置，當波形中的被採樣的訊號值，同時小於前個採樣時間、下一個採樣時間上前後被採樣的訊號值，即為波形中的波谷位置。

前述的波形中存在以下的問題，例如：當前採樣的訊號值，若僅大於(或小於)前一個採樣時間點之訊號值，但是等於下一個採樣時間點之訊號值時，即需持續監測直至與下一個採樣時間點之訊號值的大小不同為止，方能決定是否為坡峰或波谷。但是當在即時偵測時，因無法預測會有多少個採樣時間點之訊號值是相同大小，因此有可能需要無窮大之儲存空間，來存放無法預期長度之連續採樣時間點之訊號值。故，我們需要解決前述的問題，以便於使用有限的儲存空間完成波形的峰值偵測。

有鑑於先前技術的問題，本發明之一目的，係於為了更有效率的使用儲存空間，完成波形的峰值與波谷偵測，尤其是解決在即時偵測時，無論多少個採樣時間點之訊號值是相同大小，都只要使用少許固定的儲存空間即可完成波形的峰值偵測。

根據本發明之一目的，係提供一種訊號的波峰與波谷之偵測方法，包括下列步驟，波形中當前的採樣時間點之訊號值做為第一訊號值，前一個採樣時間點的訊號值作為第二訊號值，後一個採樣時間的訊號值作為第三訊號值，當第一訊號值同時大於第二訊號值與第三訊號值，則第一訊號值為波峰點，並儲存在儲存單元中，當第一訊號值同時小於第二訊號值與第三訊號值，則第一訊號值作為波谷點，並儲存在儲存單元中，或當第一訊號值大於第二訊號值且等於第三訊號值，則將第一訊號值做為第一假性波峰點，或當第一訊號值等於第二訊號值且大於第三訊號值，則將第一訊號值做為第二假性波峰點，或當第一訊號值小於第二訊號值且等於第三訊號值，則將第一訊號值做為第一假性波谷點，或當第一訊號值等於第二訊號值且小於第三訊號值，則將第一訊號值做為第二假性波谷點。

其中，當第一訊號值為第一假性波峰點時，且在儲存單元最近一次儲存為波谷點資訊，則將第一假性波峰點以第一假性波峰點資訊儲存在儲存單元，並且後續的第一訊號值等於第二訊號值及第三訊號值時，所有的第一訊號值不予儲存在儲存單元，直至取得第二假性波峰點時，則從儲存單元取出最後儲存之第一假性波峰點資訊的第一假性波峰點，並選擇第一假性波峰點或第二假性波峰點的其中之一個作為波峰點進行合併，再作為波峰點資訊儲存於儲存單元。

其中，當第一訊號值為第一假性波峰點時，且在儲存單元最近一次儲存為波谷點資訊，則將第一假性波峰點以第一假性波峰點資訊儲存在儲存單元，並且後續的第一訊號值等於第二訊號值及第三訊號值時，所有的第一訊號值不予儲存在儲存單元，直至取得第二假性波谷點時，則從儲存單元刪除最後儲存之第一假性波峰點資訊。

其中，當第一訊號值為第一假性波谷點時，且在儲存單元最近一次儲存為波峰點資訊，則將第一假性波谷點以第一假性波谷資訊儲存在儲存單元，並且後續的第一訊號值等於第二訊號值及第三訊號值時，所有的第一訊號值不予儲存在儲存單元，直至取得第二假性波谷點時，則從儲存單元取出最後儲存之第一假性波谷點資訊的第一假性波谷點，並選擇第一假性波谷點或第二假性波谷點的其中之一個作為波谷點進行合併，再做為波谷點資訊儲存於儲存單元。

其中，當第一訊號值為第一假性波谷點時，且在儲存單元最近一次儲存的第一訊號值係為波峰點，則將第一假性波谷點作為一第一假性波谷點資訊儲存在儲存單元，並且後續的第一訊號值等於第二訊號值及第三訊號值時，所有的第一訊號值不予儲存在儲存單元，直至取得第二假性波峰點時，則從儲存單元刪除最後儲存之第一假性波谷點。

綜上所述，儲存單元中最終只儲存波峰點資訊或者波谷點資訊，第一訊號值等於第二訊號值及第三訊號值時，並不會有任何訊號值被儲存在儲存單元中，所有的第一假性波峰點資訊或第一假性波谷點資訊僅暫時儲存在儲存單元中，並不會佔用到儲存單元的最終儲存空間，且第二假性波峰點或第二假性波谷點，完全不會被暫時儲存在儲存單元中，使得儲存單元以最節省儲存空間的方式，或者是較有效率的方式記錄波形中所有的波峰點或者波谷點。

為了使本發明的目的、技術方案及優點更加清楚明白，下面結合附圖及實施例，對本發明進行進一步詳細說明。應當理解，此處所描述的具體實施例僅用以解釋本發明，但並不用於限定本發明。

請參閱圖1~圖6所示，本發明係一種訊號的波峰與波谷之偵測方法，包括下列步驟： 波形1中當前的採樣時間點之訊號值做為第一訊號值10； 前一個採樣時間點的訊號值作為第二訊號值12； 後一個採樣時間的訊號值作為第三訊號值14； 再利用第一訊號值10、第二訊號值12及第三訊號值14進行下列的分類： 當第一訊號值10同時大於第二訊號值12與第三訊號值14，則確認第一訊號值10為波峰點2(如圖1所示)； 當第一訊號值10同時小於第二訊號值12與第三訊號值14，則確認第一訊號值10為波谷點3(如圖2所示)； 當第一訊號值10大於第二訊號值12且等於第三訊號值14，則將第一訊號值10做為第一假性波峰點4(如圖3所示)； 當第一訊號值10等於第二訊號值12且大於第三訊號值14，則將第一訊號值10做為第二假性波峰點5(如圖4所示)； 當第一訊號值10小於第二訊號值12且等於第三訊號值14，則將第一訊號值10做為第一假性波谷點6(如圖5所示)； 當第一訊號值10等於第二訊號值12且小於第三訊號值14，則將第一訊號值10做為第二假性波谷點7(如圖6所示)。

前述的各判斷流程並無先後順序，只要能由第一訊號值10比對大於、小於或等於第二訊號值12與第三訊號值14的關係，而確定出第一訊號值10為波峰點2、波谷點3、第一假性波峰點4、第二假性波峰點5、第一假性波谷點6或第二假性波谷點7的其中之一者即可，並確定出第一訊號值10為波峰點2或波谷點3後，依照其採樣時間點先後順序以波峰點資訊80或波谷點資訊81依序儲存在儲存單元8中而第一假性波峰點4及第一假性波谷點6則以作為第一假性波峰點資訊82或第一假性波谷點資訊83依照其採樣時間點先後順序暫存在儲存單元8。第二假性波峰點5及第二假性波谷點7不會被儲存在儲存單元8中，但是會被存到暫存器(圖中未示)，以便後續跟第一假性波峰點4或第一假性波谷點6進行辨識判斷，以便與該第一假性波峰點4或第一假性波谷點6合併為波峰點2或波谷點3，並將之儲存於儲存單元8中，或是將該第一假性波峰點4或第一假性波谷點6由儲存單元8中刪除。儲存單元8係可為堆疊式儲存結構或陣列式儲存結構，但不以此為限。

在本發明中，請參閱圖7所示，當第一訊號值10為第一假性波峰點4時，且在儲存單元8最近一次儲存的第一訊號值10係為波谷點資訊81，則將第一假性波峰點4以第一假性波峰點資訊82儲存在儲存單元8，並且後續的第一訊號值10等於第二訊號值12及第三訊號值14時，所有的第一訊號值10不予儲存在儲存單元8，直至取得第二假性波峰點5時，則從儲存單元8取出第一假性波峰點資訊82的第一假性波峰點4，並將該第一假性波峰點4與第二假性波峰點5合併為波峰點2，選擇該第一假性波峰點4或第二假性波峰點5之資訊做為波峰點資訊80儲存於儲存單元8。

在本發明中，請參閱圖8所示，當第一訊號值10為第一假性波峰點4時，且在儲存單元8最近一次儲存的第一訊號值10係為波谷點資訊81，則將第一假性波峰點4以第一假性波峰點資訊82儲存在儲存單元8，並且後續的第一訊號值10等於第二訊號值12及第三訊號值14時，所有的第一訊號值10不予儲存在儲存單元8，直至取得第二假性波谷點7時，則從儲存單元8刪除第一假性波峰點資訊82。

在本發明中，請參閱圖9所示，當第一訊號值10為第一假性波谷點6時，且在儲存單元8最近一次儲存的第一訊號值10係為波峰點資訊80，則將第一假性波谷點6以第一假性波谷點資訊83儲存在儲存單元8，並且後續的第一訊號值10等於第二訊號值12及第三訊號值14時，所有的第一訊號值10不予儲存在儲存單元8，直至取得第二假性波谷點7時，則從儲存單元8取出第一假性波谷點資訊83的第一假性波谷點6，並將該第一假性波谷點6與第二假性波谷點7合併為波谷點3，選擇該第一假性波谷點6或第二假性波谷點7之資訊作為波谷點資訊81儲存於儲存單元8。

在本發明中，請參閱圖10所示，當第一訊號值10為第一假性波谷點6時，且在儲存單元8最近一次儲存的第一訊號值10係為波峰點資訊80，則將第一假性波谷點6以第一假性波谷點資訊83儲存在儲存單元8，並且後續的第一訊號值10等於第二訊號值12及第三訊號值14時，所有的第一訊號值10不予儲存在儲存單元8，直至取得第二假性波峰點5時，則從儲存單元8刪除第一假性波谷點資訊83。

綜上所述，儲存單元8中只儲存波峰點2或者波谷點3，當第一訊號值10等於第二訊號值12及第三訊號值14時，並不會有任何訊號值被儲存在儲存單元8中，而且當儲存單元8最近一次儲存的第一訊號值10係為波峰點資訊80或者波谷點資訊81，而儲存單元8後續儲存的內容為第一假性波峰點資訊82或第一假性波谷點資訊83會被取出進行比較判斷，當是第一假性波峰點4後續遇到第二假性波峰點5，則合併成為波峰點2，並選擇以第一假性波峰點4或第二假性波峰點5之資訊作為波峰點資訊80儲存在儲存單元8，當是第一假性波峰點4後續遇到第二假性波谷點7，則刪除被儲存的第一假性波峰點資訊82，另外，若是第一假性波谷點6後續遇到第二假性波谷點7，則合併成為波谷點3，並選擇以第一假性波谷點6或第二假性波谷點7之資訊作為波谷點資訊81儲存在儲存單元8，若是第一假性波谷點6後續遇到第二假性波峰點5，則刪除被儲存的第一假性波谷點資訊83。

另外需要在此說明的是，本發明中不會出現第二假性波峰點5與第二假性波谷點7之資訊被暫存在儲存單元8的情況，其原因在於第二假性波峰點5會與第一假性波峰點4合併成波峰點，或者第二假性波峰點5會與第一假性波谷點6相互抵消並將該第一假性波谷點6資訊83刪除，所以第二假性波峰點5不會有機會存入到儲存單元8，同時暫存之第一假性波峰點4資訊82或是第一假性波谷點6資訊83會被合併或取消，也不會出現於最終之儲存單元8。而第二假性波谷點7與第一假性波谷點6合併成波谷點，或者第二假性波谷點7會與第一假性波峰點4相互抵消刪除，所以第二假性波谷點7不會有機會存入到儲存單元8。意即，最終之儲存單元8中只會出現波峰點資訊80與波谷點資訊81，而第一假性波峰點資訊82或第一假性波谷點資訊83僅會於偵測過程中一次性的暫存在儲存單元8。

據上所述，本發明不會將第一訊號值10、第二訊號值12及第三訊號值14不斷地儲存在儲存單元8，而只在出現波峰點2、波谷點3才會作為波峰點資訊80、波谷點資訊81儲存在儲存單元8，另外第一假性波峰點資訊82或第一假性波谷點資訊83是一次性的暫存在儲存單元8，而第二假性波峰點5與第二假性波谷點7完全不會有機會存入儲存單元8，所以在偵測訊號的波峰或波谷的過程中儲存單元8將是以較有效率的方式進行資料儲存，解決了傳統偵測方法可能需要無窮大之儲存空間之問題，或者是至少解決傳統偵測方法較無效率的使用儲存空間的問題。

上列詳細說明係針對本發明的可行實施例之具體說明，惟前述的實施例並非用以限制本發明之專利範圍，凡未脫離本發明技藝精神所為之等效實施或變更，均應包含於本案之專利範圍中。

1:波形

10:第一訊號值

12:第二訊號值

14:第三訊號值

2:波峰點

3:波谷點

4:第一假性波峰點

5:第二假性波峰點

6:第一假性波谷點

7:第二假性波谷點

8:儲存單元

80:波峰點資訊

81:波谷點資訊

82:第一假性波峰點資訊

83:第一假性波谷點資訊

圖1係本發明之確認第一訊號值為波峰點之示意圖。 圖2係本發明之確認第一訊號值為波谷點之示意圖。 圖3係本發明之確認第一訊號值為第一假性波峰點之示意圖。 圖4係本發明之確認第一訊號值為第二假性波峰點之示意圖。 圖5係本發明之確認第一訊號值為第一假性波谷點之示意圖。 圖6係本發明之確認第一訊號值為第二假性波谷點之示意圖。 圖7係本發明之第一假性波峰點與第二假性波峰點組合合併形成波峰點之示意圖。 圖8係本發明之第一假性波峰點與第二假性波谷點組合刪除儲存單元中的第一假性波峰點之示意圖。 圖9係本發明之第一假性波谷點與第二假性波谷點組合合併形成波谷點之示意圖。 圖10係本發明之第一假性波谷點與第二假性波峰點組合刪除儲存單元中的第一假性波谷點之示意圖。

1:波形

10:第一訊號值

12:第二訊號值

14:第三訊號值

4:第一假性波峰點

5:第二假性波峰點

8:儲存單元

80:波峰點資訊

81:波谷點資訊

82:第一假性波峰點資訊

Claims (4)

1. 一種訊號的波峰與波谷之偵測方法，包括：將波形中當前的採樣時間點之訊號值做為一第一訊號值；前一個採樣時間點的訊號值作為一第二訊號值；後一個採樣時間的訊號值作為一第三訊號值：再依據該第一訊號值、該第二訊號值及該第三訊號值，對該第一訊號值進行下列的分類；當該第一訊號值同時大於該第二訊號值與該第三訊號值，則確認該第一訊號值為一波峰點；當該第一訊號值同時小於該第二訊號值與該第三訊號值，則確認該第一訊號值為一波谷點；當該第一訊號值大於該第二訊號值且等於該第三訊號值，則將該第一訊號值做為一第一假性波峰點；當該第一訊號值等於該第二訊號值且大於該第三訊號值，則將該第一訊號值做為一第二假性波峰點；當該第一訊號值小於該第二訊號值且等於該第三訊號值，則將該第一訊號值做為一第一假性波谷點；當該第一訊號值等於該第二訊號值且小於該第三訊號值，則將該第一訊號值做為一第二假性波谷點；其中，該波峰點與該波谷點係依照各該採樣時間點的先後順序作為一波峰點資訊與一波谷點資訊被儲存在一儲存單元中； 其中當該第一訊號值為該第一假性波峰點時，且在該儲存單元最近一次儲存為該波谷點資訊，則將該第一假性波峰點作為一第一假性波峰點資訊儲存在該儲存單元，並且後續的該第一訊號值等於該第二訊號值及該第三訊號值時，所有的該第一訊號值不予儲存在該儲存單元，直至取得該第二假性波峰點時，則從該儲存單元取出該第一假性波峰點資訊的該第一假性波峰點，將該第一假性波峰點與該第二假性波峰點進行合併，即選擇該第一假性波峰點或該第二假性波峰點的其中之一個的資訊作為該波峰點資訊，再將該波峰點資訊儲存於該儲存單元。
2. 如請求項1所述的訊號的波峰與波谷之偵測方法，其中當該第一訊號值為該第一假性波峰點時，且在該儲存單元最近一次儲存為該波谷點資訊，則將該第一假性波峰點作為一第一假性波峰點資訊儲存在該儲存單元，並且後續的該第一訊號值等於該第二訊號值及該第三訊號值時，所有的該第一訊號值不予儲存在該儲存單元，直至取得該第二假性波谷點時，則從該儲存單元刪除該第一假性波峰點資訊。
3. 如請求項1所述的訊號的波峰與波谷之偵測方法，其中當該第一訊號值為該第一假性波谷點時，且在該儲存單元最近一次儲存為該波峰點資訊，則將該第一假性波谷點以一第一假性波谷資訊儲存在該儲存單元，並且後續的該第一訊號值等於該第二訊號值及該第三訊號值時，所有的該第一訊號值不予儲存在該儲存單元，直至取得該第二假性波谷點時，則從該儲存單元取出該第一假性波谷點資訊的該第一假性波谷點，將該第一假性波谷點與該第二假性波谷點進行合併，即選擇該第一假性波谷點或該第二假性波谷點的其中之一個的資訊作為該波谷點資訊，再將該波谷點資訊儲存於該儲存單元。
4. 如請求項1所述的訊號的波峰與波谷之偵測方法，其中當該第一訊號值為該第一假性波谷點時，且在該儲存單元最近一次儲存為該波峰點資訊，則將該第一假性波谷點作為一第一假性波谷點資訊儲存在該儲存單元，並且後續的該第一訊號值等於該第二訊號值及第三訊號值時，所有的該第一訊號值不予儲存在該儲存單元，直至取得該第二假性波峰點時，則從該儲存單元刪除該第一假性波谷點。

Images