专利名称：信号时域波形极值点和周期的检测方法
技术领域：
本发明涉及数据和图形的处理方法，尤其是涉及基于有限状态机转换的信号时域波形极值点的检测和处理方法。
背景技术：
现有技术采集和分析数据，例如检测人体生理信息，包括血压，呼吸，心电，脑电和胎音信号，再又如收集和分析水文信号，大气压力和温度信号等，及其他气体检测中气体各成分浓度信号，以及工程上检测炼钢炉的炉温信号等，一般是将这类时域信号取固定的一段时间内的数据进行分析，找到该区段的极值点，并计算周期。这样分析数据的方法，缺点是容易漏检波形，而且不能做到实时分析。这种方法主要考虑的是在一段时间内可以检测出波形的最大值最小值和周期就可以了，并没有考虑如何使算法能实时准确地检测出波形的最大值和最小值。比如，目前在一些实时呼吸波形检测中使用的传统方法，就有漏检波形呼吸率计算不准确等问题。
发明内容
本发明要解决的技术问题在于避免上述现有技术的不足之处而提供一种信号时域波形极值点和周期的检测方法，从而能快速实时准确地检测出波形的最大值最小值和周期。
本发明解决上述技术问题所采用的技术方案是提供一种信号时域波形极值点和周期的检测方法，包括步骤1.一种信号时域波形极值点和周期的检测方法，包括步骤(1)在波形开始的一两个周期内，估计最大值和最小值的范围最大值Vmax的范围是开始一段时间内检测到的最大值H的1至2倍，即Vmax＝(1～2)·H最小值Vmin的范围是开始一段时间内检测到的最小值L的0.5至1倍，即Vmin＝(0.5～1)·L
(2)在开始搜索前，对最大值最小值和中值进行初始化，并给出门限滤波参数，对最大值初始值和最小值初始值分别赋于最小值和最大值，中值初始值为开始一段时间的算术平均值；(3)检测实时的信号值，比较它是否大于中值与滤波门限之和；(4)如果步骤(3)的结果为“是”，则进入“搜索最大值”的状态，并将当前检测到的信号值作为最大值，如果前一状态为“搜索最小值”状态，则用当前最大值和最小值的算术平均值来更新中值，如果前一状态也为“搜索最大值”状态，则维持中值不变，然后返回步骤(3)；(5)如果步骤(3)的结果为“否”，则比较该实时信号是否小于中值与滤波门限之差；如果步骤(5)的结果为“是”，则进入“搜索最小值”的状态，并将当前检测到的信号值作为最小值，如果前一状态为“搜索最大值”状态，则用当前最大值和最小值的算术平均值来更新中值，如果前一状态也为“搜索最小值”状态，则维持中值不变，然后返回步骤(3)。
同现有技术相比较，本发明的优点在于能够实时快速准确地检测各种类似正弦波的生理信号和非生理信号的周期和极值点；每次只来一个点就可以处理，只要有一个波形就检测能得到最大值最小值和周期，而不是要收集一段数据才可以开始处理，并且能够很好地跟上波形突变和波动的情况。
图1是基于有限状态机波形检测算法的示意图；图2是状态转移图其中X为实时信号的输入，M为上一个周期检测得到的中值，D为滤波门限；图3是医用二氧化碳浓度监测系统实现框图；图4是正常人的呼吸二氧化碳浓度波形曲线；图5是算法实现的简明流程图；图6是实际检测到的呼吸二氧化碳浓度曲线。
具体实施例方式
以下结合各附图所示之最佳实施例对本发明作进一步详述。
对于一个如图1所示形态类似正弦波的波形，在开始的一段时间内，通常是一至两个周期，只估计最大值最小值和中值的范围。最大值Vmax的范围是开始一段时间内检测到的最大值H的1至2倍，即Vmax＝(1～2)·H，最小值Vmin的范围是开始一段时间内检测到的最小值L的0.5至1倍，即Vmin＝(0.5～1)·L。在开始搜索之前，最大值MAX和最小值MIN都应该分别赋以估计的最小值和最大值，这样保证以后的搜索可以找到最大值和最小值。中值初始值使用的是开始一段时间内的算术平均值。得到中值之后即开始一个完整的波形检测过程。对于每个实时信号的输入X，都会和中值M较，大于中值M与滤波门限D之和，就进入搜索最大值状态，简称“0状态”，并搜索最大值MAX；小于中值M与滤波门限D之差就进入搜索最小值状态，简称“1状态”，并搜索最小值MIN。这两个状态就这样交替切换，如图2所示。在完成两次状态交替的时候就要用当前一个周期的数据，即最大值和最小值的算术平均值来更新中值M，最大值MAX和最小值MIN在一个周期末的时候也要复位。这样在一个周期之内就可以得到极值点和周期，而不需要用一段数据才能做波形识别。对正常情况下的波形，如图1所示，在两条竖线ab和cd之间就可以得到最大值和最小值从而算出周期，并且根据实际需要记录当前的最大值最小值和周期。当波形发生突变的时候，由于每个周期都更新了中值最大值和最小值，所以检测到的波形在一个周期之内就会跟上实时波形的变化，马上检测出波形的最大值和最小值。当波形发生突变的时候，实际波形的中值点在横线ef处，但是搜索最大值和最小值用的中值是原来上一个周期的中值，这样中值如果能跟上波形的变化，就仍然可以准确的找到波形的最大值最小值和周期。但是如果中值没有跟上波形的变化，那么就会漏检波形。如果每一个周期都更新了中值，中值有较快变化，就能很快跟上波形的变化，不会漏检波形。如果是非常大的突变，则用两个异常参数来判断以保证在一定的时间内恢复波形检测。首先是时间上的异常处理，如果在规定的时间范围内没有搜索到极值点，就会重新开始估计最大值最小值和中值的范围并开始搜索。其次是每次更新中值的时候，用当前的中值和原来的中值做比较，如果超过了正常变化范围，就认为这是一个突发的干扰，于是限制中值更新，让它仍然保持在正常的范围之内，所述中值的正常范围是指中值在当前中值减去波形峰值和当前中值加上波形峰值之间。
医用二氧化碳浓度监测装置是监测呼吸过程中二氧化碳浓度水平的仪器。它实现的框图如图3所示传感器得到的小信号经放大和采集电路进入计算机系统，在计算机系统中经过信号预处理之后，开始呼吸二氧化碳浓度波形特征点的检测。
正常人的呼吸二氧化碳浓度波形曲线如图4所示，正常人呼吸的二氧化碳浓度曲线由4部分组成第一部分为吸气基线，处于零位，是呼气的开始部分，即AB段；第二部分为呼气上升支，比较陡直，相当于BC段；第三部分为呼气平台，呈水平形，相当于图中的CD段；第四部分为吸气下降支，迅速而陡直下降到吸气开始，即相当于图中的DE段。
对于呼吸气体的二氧化碳浓度波形，一般观察基线、高度、形态和频率四个参数。基线代表吸气末二氧化碳浓度，一般比较接近于零；高度代表呼气末二氧化碳浓度；形态主要观察正常和不正常的波形；频率反应了每分钟呼吸次数。对于这样一个类似正弦波的二氧化碳浓度曲线，就可以利用上述的检测算法进行波形特征点的检测。
本发明算法实现的简明流程如图5所示。首先开始进行最大值和最小值范围估计，估计得到最大值和最小值范围之后给出一个最大值、最小值和中值的初始值，并给出门限滤波参数和时间异常参数。有了这些基本的参数，就可以开始一个波形检测的过程。每来一个新值就判断它是否大于中值与滤波门限之和，若大于就进入“搜索最大值”状态搜索最大值；否则将新值和中值与滤波门限之差比较，如果小于就进入“搜索最小值”状态搜索最小值；否则就判断是否时间异常，若非时间异常就继续下一个点的搜索；若为时间异常就更新一次中值，最大值和最小值复位，完成更新中值和最大值最小值复位之后开始新的搜索过程。在进入“搜索最大值”状态搜索到最大值之后将更新最大值，同样进入“搜索最小值”状态搜索到最小值之后将更新最小值。在开始下一个新的波形检测之前以搜索到的最大值和最小值的算术平均值来更新中值，同时将最大值和最小值复位到上一次波形检测前的状态。在更新中值时，要将当时的中值与原来的中值相比较，如果超过正常变化的范围，就限制中值更新，让它保持在正常的范围之内，所述中值的正常范围是指中值在当前中值减去波形峰值和当前中值加上波形峰值之间。
应用该算法检测的实际效果如图6所示图中双线为检测最大值的上升支，双线末端就是最大值所在位置。刷状线为检测最小值的下降支，刷状线末端就是最小值所在位置。可见该算法可以十分清楚检测到最大值和最小值，能准确定位最大值和最小值位置就可以计算得到准确的周期。
本发明还可用于其它形态类似正弦波的波形极值点和周期的检测。如呼吸波形、脉搏信号、心电图信号、肌电信号等生理信号；控制领域的炉温信号；水利领域的水文信号；电子通讯领域内类似正弦波的信号；大气监测的气压信号，温度信号，及其它气体检测中各种气体成分的浓度信号。
权利要求
1.一种信号时域波形极值点和周期的检测方法，包括步骤(1)在波形开始的一两个周期内， 估计最大值和最小值的范围最大值Vmax的范围是开始一段时间内检测到的最大值H的1至2倍，即Vmax＝(1～2)·H最小值Vmin的范围是开始一段时间内检测到的最小值L的0.5至1倍，即Vmin＝(0.5～1)·L(2)在开始搜索前，对最大值最小值和中值进行初始化，并给出门限滤波参数，对最大值初始值和最小值初始值分别赋于最小值和最大值，中值初始值为开始一段时间的算术平均值；(3)检测实时的信号值，比较它是否大于中值与滤波门限之和；(4)如果步骤(3)的结果为“是”，则进入“搜索最大值”的状态，并将当前检测到的信号值作为最大值，如果前一状态为“搜索最小值”状态，则用当前最大值和最小值的算术平均值来更新中值，如果前一状态也为“搜索最大值”状态，则维持中值不变，然后返回步骤(3)；(5)如果步骤(3)的结果为“否”，则比较该实时信号是否小于中值与滤波门限之差；(6)如果步骤(5)的结果为“是”，则进入“搜索最小值”的状态，并将当前检测到的信号值作为最小值，如果前一状态为“搜索最大值”状态，则用当前最大值和最小值的算术平均值来更新中值，如果前一状态也为“搜索最小值”状态，则维持中值不变，然后返回步骤(3)；
2.如权利要求1所述的信号时域波形极值点和周期的检测方法，其特征在于在所述步骤(2)中还给出时间参数，所述时间参数为周期的1至2倍，并再执行如下步骤(7)如果步骤(5)的结果为“否”，则用时间参数来判断是否为时间上的异常；(8)如果步骤(7)的结果为“是”，则返回到步骤(1)；(9)如果步骤(7)的结果为“否”，则直接返回到步骤(3)；
3.如权利要求1或2所述的信号时域波形极值点和周期的检测方法，其特征在于在所述步骤(4)和(6)中更新中值前，将当前的中值与原来的中值相比较，如果超过正常变化范围，就限制中值更新，让它仍然保持在正常范围之内；
4.如权利要求1或2中任一项所述的信号时域波形极值点和周期的检测方法，其特征在于所述门限滤波参数不大于波形幅值的百分之十；
5.如权利要求3所述的信号时域波形极值点和周期的检测方法，其特征在于所述中值的正常范围是指中值在当前中值减去波形峰值和当前中值加上波形峰值之间。
全文摘要
一种信号时域波形极值点和周期的检测方法，首先在进行最大值和最小值范围估计后给出一个最大值、最小值和中值的初始值，并给出门限滤波参数，然后开始一个波形检测的过程。每来一个新值就判断它是否大于中值与滤波门限之和，若大于就进入“搜索最大值”状态搜索最大值；否则将新值和中值与滤波门限之差比较，如果小于就进入“搜索最小值”状态搜索最小值。在开始下一个新的波形检测之前以搜索到的最大值和最小值的算术平均值来更新中值，同时将最大值和最小值复位到上一次波形检测前的状态。本发明能够实时快速准确地检测各种类似正弦波的生理信号和非生理信号的周期和极值点，能够很好地跟上波形突变和波动的情况。
文档编号A61B5/0205GK1525395SQ0311383
公开日2004年9月1日 申请日期2003年2月24日 优先权日2003年2月24日
发明者周慧玲, 叶继伦 申请人:深圳迈瑞生物医疗电子股份有限公司