[发明专利]一种脉搏波的分析方法及装置有效
申请号: | 201811119099.1 | 申请日: | 2018-09-25 |
公开(公告)号: | CN109171678B | 公开(公告)日: | 2021-06-22 |
发明(设计)人: | 韦海成;李云琴;吴官胜;肖明霞;潘俊涛;许亚杰 | 申请(专利权)人: | 北方民族大学 |
主分类号: | A61B5/02 | 分类号: | A61B5/02;A61B5/00 |
代理公司: | 北京睿智保诚专利代理事务所(普通合伙) 11732 | 代理人: | 周新楣 |
地址: | 750021 宁夏回族*** | 国省代码: | 宁夏;64 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 一种 脉搏 分析 方法 装置 | ||
1.一种脉搏波的分析方法,其特征在于,包括以下步骤:
步骤1,接收待分析的原始脉搏波;采用红外传感器对脉搏信号进行采集;
步骤2,对所述原始脉搏波进行分解,从原始脉搏波中提取出若干个本征模态函数;
步骤3,选取与原始脉搏波最接近的一个本征模态函数,基于该本征模态函数进行希尔伯特变换,得到原始脉搏波的边际谱密度函数,以该边际谱密度函数最大能量的瞬时频率作为所述原始脉搏波的最大频率输出;
所述对所述原始脉搏波进行分解,从原始脉搏波中提取出若干个本征模态函数,包括步骤:
步骤2-1,提取脉搏波x(t)中所有的局部极大值和局部极小值;若某个时刻的值既大于前一时刻的值也大于后一时刻的值,则该时刻的值为局部极大值;若某个时刻的值既小于前一时刻的值也小于后一时刻的值,则该时刻的值为局部极小值;
步骤2-2,使用三次样条函数进行拟合,获得上包络线xmax(t)和下包络线xmin(t),并计算上包络线、下包络线的均值,m(t)=[xmax(t)+xmin(t)]/2;
步骤2-3,用脉搏波信号x(t)减去均值m(t),得到一个组件h(t)=x(t)-m(t);
步骤2-4,判断h(t)是否为本征模态函数,若是,则进入步骤2-5,否则进入步骤2-6;
步骤2-5,令I1(t)=h(t),r(t)=x(t)-I1(t),并判断r(t)是否为单调序列或常值序列,若是则结束,否则以r(t)代替x(t),并返回步骤2-1;
步骤2-6,以h(t)代替x(t),并返回步骤2-1。
2.根据权利要求1所述的方法,其特征在于,所述步骤2-6用以下步骤代替:
判断标准偏差SD是否大于设定值,若是,则以h(t)代替x(t),并返回步骤2-1,否则进入步骤2-5;其中,h1(k-1)(t)和h1k(t)分别表示两个连续的筛选结果。
3.根据权利要求1所述的方法,其特征在于,所述基于该本征模态函数进行希尔伯特变换,得到原始脉搏波的边际谱密度函数,包括步骤:
步骤3-1,对该本征模态函数做希尔伯特变换,得到希尔伯特变换函数步骤3-2,构造解析信号zi(t),Zi(t)=Ii(t)+H[Ii(t)]=αi(t)exp[jφi(t)],式中,αi(t)为幅值,φi(t)为相位,则有瞬时频率
步骤3-3,将原始脉搏波信号改写为,则有把H(ω,t)对时间积分,得到Hilbert边际谱函数为T表示总的数据长度。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于北方民族大学,未经北方民族大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201811119099.1/1.html,转载请声明来源钻瓜专利网。
- 上一篇:一种血流动力学参数测量装置及方法
- 下一篇:一种心内科医生用胸腔监测急救装置