信号时域波形极值点和周期的检测方法技术

一种信号时域波形极值点和周期的检测方法，首先在进行最大值和最小值范围估计后给出一个最大值、最小值和中值的初始值，并给出门限滤波参数，然后开始一个波形检测的过程。每来一个新值就判断它是否大于中值与滤波门限之和，若大于就进入“搜索最大值”状态搜索最大值；否则将新值和中值与滤波门限之差比较，如果小于就进入“搜索最小值”状态搜索最小值。在开始下一个新的波形检测之前以搜索到的最大值和最小值的算术平均值来更新中值，同时将最大值和最小值复位到上一次波形检测前的状态。本发明专利技术能够实时快速准确地检测各种类似正弦波的生理信号和非生理信号的周期和极值点，能够很好地跟上波形突变和波动的情况。（*该技术在2023年保护过期，可自由使用*）

【技术实现步骤摘要】

本专利技术涉及数据和图形的处理方法，尤其是涉及基于有限状态机转换的信号时域波形极值点的检测和处理方法。
技术介绍
现有技术采集和分析数据，例如检测人体生理信息，包括血压，呼吸，心电，脑电和胎音信号，再又如收集和分析水文信号，大气压力和温度信号等，及其他气体检测中气体各成分浓度信号，以及工程上检测炼钢炉的炉温信号等，一般是将这类时域信号取固定的一段时间内的数据进行分析，找到该区段的极值点，并计算周期。这样分析数据的方法，缺点是容易漏检波形，而且不能做到实时分析。这种方法主要考虑的是在一段时间内可以检测出波形的最大值最小值和周期就可以了，并没有考虑如何使算法能实时准确地检测出波形的最大值和最小值。比如，目前在一些实时呼吸波形检测中使用的传统方法，就有漏检波形呼吸率计算不准确等问题。
技术实现思路
本专利技术要解决的技术问题在于避免上述现有技术的不足之处而提供一种，从而能快速实时准确地检测出波形的最大值最小值和周期。本专利技术解决上述技术问题所采用的技术方案是提供一种，包括步骤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所示。首先开始进行最大值和最小值范围估计，估计得到最大值和最小值范围之后给出一个最大值、最小值和中值的初始值，并给出门限滤波参数和时间异常参数。有了这些基本的参数，就可以开始一个波形检测的过程。每来一个新值就判断它是否大于中值与滤波门限之和，若大于就进入“搜索最大值”状态搜索最大值；否则将新值和中值与滤波门限之差比较，如果小于就进入“搜索最小值”状态搜索最小值；否则就判断是否时间异常，若非时间异常就继续下一个点的搜索；若为时间异常就更新一次中值，最大值和最小值复位，完成更新中值和最大值最小值复位之后开始新的搜索过程。在进入“搜索最大值”状态搜索到最大值之后将更新最大值，同样进入“搜索最小值”状态搜索到最小值之本文档来自技高网...

【技术保护点】
一种信号时域波形极值点和周期的检测方法，包括步骤：　　　　（１）在波形开始的一两个周期内，估计最大值和最小值的范围：　　　　最大值Ｖｍａｘ的范围是开始一段时间内检测到的最大值Ｈ的１至２倍，即Ｖｍａｘ＝（１～２）.Ｈ　　　　最小值Ｖｍｉｎ的范围是开始一段时间内检测到的最小值Ｌ的０．５至１倍，即Ｖｍｉｎ＝（０．５～１）.Ｌ　　　　（２）在开始搜索前，对最大值、最小值和中值进行初始化，并给出门限滤波参数，对最大值初始值和最小值初始值分别赋于最小值和最大值，中值初始值为开始一段时间的算术平均值；　　　　（３）检测实时的信号值，比较它是否大于中值与滤波门限之和；　　　　（４）如果步骤（３）的结果为“是”，则进入“搜索最大值”的状态，并将当前检测到的信号值作为最大值，如果前一状态为“搜索最小值”状态，则用当前最大值和最小值的算术平均值来更新中值，如果前一状态也为“搜索最大值”状态，则维持中值不变，然后返回步骤（３）；　　　　（５）如果步骤（３）的结果为“否”，则比较该实时信号是否小于中值与滤波门限之差；　　　　（６）如果步骤（５）的结果为“是”，则进入“搜索最小值”的状态，并将当前检测到的信号值作为最小值，如果前一状态为“搜索最大值”状态，则用当前最大值和最小值的算术平均值来更新中值，如果前一状态也为“搜索最小值”状态，则维持中值不变，然后返回步骤（３）。...

【技术特征摘要】

【专利技术属性】
技术研发人员：叶继伦，周慧玲，
申请(专利权)人：深圳迈瑞生物医疗电子股份有限公司，
类型：发明
国别省市：94[中国|深圳]