CN113288075B - 浮中沉脉位识别方法及系统 - Google Patents

浮中沉脉位识别方法及系统 Download PDF

Info

Publication number
CN113288075B
CN113288075B CN202110607165.5A CN202110607165A CN113288075B CN 113288075 B CN113288075 B CN 113288075B CN 202110607165 A CN202110607165 A CN 202110607165A CN 113288075 B CN113288075 B CN 113288075B
Authority
CN
China
Prior art keywords
pulse
voltage value
applied pressure
sample
value
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110607165.5A
Other languages
English (en)
Other versions
CN113288075A (zh
Inventor
不公告发明人
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanxi Hemei Internet Smart Traditional Chinese Medicine Industry Technology Research Institute Co ltd
Shanxi University of Chinese Mediciine
Original Assignee
Shanxi Hemei Internet Smart Traditional Chinese Medicine Industry Technology Research Institute Co ltd
Shanxi University of Chinese Mediciine
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shanxi Hemei Internet Smart Traditional Chinese Medicine Industry Technology Research Institute Co ltd, Shanxi University of Chinese Mediciine filed Critical Shanxi Hemei Internet Smart Traditional Chinese Medicine Industry Technology Research Institute Co ltd
Priority to CN202110607165.5A priority Critical patent/CN113288075B/zh
Publication of CN113288075A publication Critical patent/CN113288075A/zh
Application granted granted Critical
Publication of CN113288075B publication Critical patent/CN113288075B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/48Other medical applications
    • A61B5/4854Diagnosis based on concepts of traditional oriental medicine
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Surgery (AREA)
  • General Health & Medical Sciences (AREA)
  • Veterinary Medicine (AREA)
  • Physics & Mathematics (AREA)
  • Public Health (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Molecular Biology (AREA)
  • Physiology (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Psychiatry (AREA)
  • Signal Processing (AREA)
  • Cardiology (AREA)
  • Alternative & Traditional Medicine (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

本申请涉及浮中沉脉位识别方法及系统,涉及中医脉搏识别领域。浮中沉脉位识别方法,方法包括:根据预设第一算法,通过外加压力电压值、脉搏波电压值和测量时间计算待测样本下单周期脉搏信号的平方差指数;若平方差指数满足第一预设条件,则待测样本申请脉搏为平;若平方差指数不满足第一预设条件,根据预设第二算法,通过外加压力电压值和幅值电压值计算得到脉搏斜率指数,若脉搏斜率指数满足第二预设条件,则待测样本的脉搏为浮脉,若脉搏斜率指数不满足第二预设条件,则待测样本的脉搏为沉脉;本申请将原来需要中医获取的脉搏转化为可以通过设备获取,并且设备可以直接获取浮脉,中脉,沉脉。

Description

浮中沉脉位识别方法及系统
技术领域
本申请涉及中医脉搏识别领域,具体而言,涉及一种浮中沉脉位识别方法及系统。
背景技术
中医主要通过感知寸口脉的脉象来进行疾病诊断。传统中医学认为脉象可以反映人体系统功能的平衡状态,是中医判断人体疾病的重要依据。脉诊是传统中医最具特色的疾病诊察手段之一,其中,“位”是指脉位的“浮,中,沉”“寸,关,尺”,是对身体变化最敏感的脉象属性。
现有技术中中医通常是通过诊脉来判断患者的身体情况,而诊脉是一种技术手段,需要有经验中医使用手指对患者的脉搏进行号脉,并将脉象与医书上的病症或者医师的经验进行匹配,完成对患者疾病的诊断。
但是,由于中医号脉是一种技术,需要长期的学习锻炼,并且号脉的准确程度也很大程度的依靠医师的自身经验,使得号脉结果严重依靠人为经验,易于产生人为误差,并且有诊脉难度较大,不利于中医的推广。
发明内容
本发明的目的在于,针对上述现有技术中的不足,提供一种浮中沉脉位识别方法及系统,以解决现有技术中由于中医号脉是一种技术,需要长期的学习锻炼,并且号脉的准确程度也很大程度的依靠医师的自身经验,使得号脉结果严重依靠人为经验,易于产生人为误差,并且有诊脉难度较大,不利于中医的推广的问题。
为实现上述目的,本发明实施例采用的技术方案如下:
第一方面,本申请提供一种浮中沉脉位识别方法,方法包括:
获取待测样本的外加压力电压值、脉搏波幅值电压值、脉搏波电压值和测量时间;
根据预设第一算法,通过外加压力电压值、脉搏波电压值和测量时间计算待测样本下单周期脉搏信号的平方差指数;
若平方差指数满足第一预设条件,则待测样本的脉搏为平;若平方差指数不满足第一预设条件,根据预设第二算法,通过外加压力电压值和脉搏波幅值电压值计算得到脉搏斜率指数,若脉搏斜率指数满足第二预设条件,则待测样本的脉搏为浮脉,若脉搏斜率指数不满足第二预设条件,则待测样本的脉搏为沉脉。
可选地,该预设第一算法为:
Figure GDA0003688697420000031
其中,f1、f2以及Amax均代表不同外加压力电压值时脉搏波电压值,f1为外加压力从初始值至最大值范围内的任意值,f2为外加压力从最大值至最小值范围内的任意值,Amax为外加压力最大时的脉搏波电压值。
可选地,该根据预设第一算法,通过外加压力电压值、脉搏波电压值和测量时间计算待测样本下单周期脉搏信号的平方差指数的步骤之前还包括:
以待测样本的脉搏的波谷为起点,计算不同压力下待测样本的脉搏波的平方差之和;
根据预设条件选择主波振幅,并对主波振幅进行归一化处理。
可选地,该预设第二算法为:
Figure GDA0003688697420000032
其中,A1为一个周期内外加压力电压值上升段范围内,任意外加压力P1下的脉搏波幅值电压值,A2为一个周期内外加压力电压值下降段范围内,任意外加压力P2下的脉搏波幅值电压值,Amax为一个周期内在最大外加压力Pmax下获得的最大脉搏波幅值电压值。
可选地,该第一预设条件为小于10,第二预设条件为大于0。
第二方面,本申请提供一种浮中沉脉位识别系统,系统用于实现第一方面中任意一项的方法,系统包括:获取模块、计算模块和判断模块;
获取模块,用于获取待测样本的外加压力电压值、脉搏波幅值电压值、脉搏波电压值和测量时间;
计算模块,用于根据预设第一算法,通过外加压力电压值、脉搏波电压值和测量时间计算待测样本下单周期脉搏信号的平方差指数;
判断模块,用于若平方差指数满足第一预设条件,则待测样本的脉搏为平;若平方差指数不满足第一预设条件,根据预设第二算法,通过外加压力电压值和脉搏波幅值电压值计算得到脉搏斜率指数,若脉搏斜率指数满足第二预设条件,则待测样本的脉搏为浮脉,若脉搏斜率指数不满足第二预设条件,则待测样本的脉搏为沉脉。
可选地,该预设第一算法为:
Figure GDA0003688697420000051
其中,f1、f2以及Amax均代表不同外加压力电压值时脉搏波电压值,f1为外加压力从初始值至最大值范围内的任意值,f2为外加压力从最大值至最小值范围内的任意值,Amax为外加压力最大时的脉搏波电压值。
可选地,该系统还包括预处理模块,预处理模块用于以待测样本的脉搏的波谷为起点,计算不同压力下待测样本的脉搏波的平方差之和;根据预设条件选择主波振幅,并对主波振幅进行归一化处理。
可选地,该预设第二算法为:
Figure GDA0003688697420000052
其中,A1为一个周期内外加压力电压值上升段范围内,任意外加压力P1下的脉搏波幅值电压值,A2为一个周期内外加压力电压值下降段范围内,任意外加压力P2下的脉搏波幅值电压值,Amax为一个周期内在最大外加压力Pmax下获得的最大脉搏波幅值电压值。
可选地,该第一预设条件为小于10,第二预设条件为大于0。
第三方面,本申请提供一种电子设备,包括:存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现第一方面的浮中沉脉位识别方法。
第四方面,本申请提供一种计算机可读存储介质,所述计算机可读存储介质包括计算机程序,所述计算机程序运行时控制所述计算机可读存储介质所在电子设备执行第一方面的浮中沉脉位识别方法。
本发明的有益效果是:
浮中沉脉位识别方法,方法包括:获取待测样本的外加压力电压值、脉搏波幅值电压值、脉搏波电压值和测量时间;根据预设第一算法,通过外加压力电压值、脉搏波电压值和测量时间计算待测样本下单周期脉搏信号的平方差指数;若平方差指数满足第一预设条件,则待测样本的脉搏为平;若平方差指数不满足第一预设条件,根据预设第二算法,通过外加压力电压值和脉搏波幅值电压值计算得到脉搏斜率指数,若脉搏斜率指数满足第二预设条件,则待测样本的脉搏为浮脉,若脉搏斜率指数不满足第二预设条件,则待测样本的脉搏为沉脉;本申请将原来需要中医获取的脉搏转化为可以通过设备获取,并且设备可以直接获取浮脉,中脉,沉脉,减少人为检测的误差,并且不依靠医师自身经验,有利于中医和中医理论的推广。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,应当理解,以下附图仅示出了本发明的某些实施例,因此不应被看作是对范围的限定,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他相关的附图。
图1为本申请提供的一种应用场景示意图;
图2为本申请实施例提供的一种浮中沉脉位识别方法的流程示意图;
图3为本申请实施例提供的一种浮中沉脉位识别方法的脉搏波电压值与时间的关系图;
图4为本申请实施例提供的一种浮中沉脉位识别方法的复制电压值与外加压力电压值的关系图;
图5为本申请实施例提供的另一种浮中沉脉位识别方法的流程示意图;
图6为本申请实施例提供的一种浮中沉脉位识别方法的脉宽示意图;
图7为本申请实施例提供的另一种浮中沉脉位识别方法的脉宽示意图。
具体实施方式
下面将结合本申请中附图,对本申请中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。通常在此处附图中描述和示出的本申请的组件可以以各种不同的配置来布置和设计。因此,以下对在附图中提供的本申请的实施例的详细描述并非旨在限制要求保护的本申请的范围,而是仅仅表示本申请的选定实施例。基于本申请的实施例,本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本申请保护的范围。
应注意到:相似的标号和字母在下面的附图中表示类似项,因此,一旦某一项在一个附图中被定义,则在随后的附图中不需要对其进行进一步定义和解释。
图1为本申请提供的一种应用场景示意图,如图1所示,本申请中的方法可以应用于图1所示的电子设备10。如图1所示,电子设备10可以包括:存储器11、处理器12、网络模块13。
存储器11、处理器12、网络模块13相互之间直接或间接地电性连接,以实现数据的传输或交互。例如,这些元件相互之间可通过一条或多条通讯总线或信号线实现电性连接。存储器11中存储有至少一个可以软件或固件(firmware)的形式存储于所述存储器11中的功能模块,所述处理器12通过运行存储在存储器11内的软件或硬件形式的功能模块,从而执行各种功能应用以及数据处理,即实现本申请中由电子设备10执行的方法。
其中,所述存储器11可以是,但不限于,随机存取存储器(Random Access Memory,RAM),只读存储器(Read Only Memory,ROM),可编程只读存储器(Programmable Read-OnlyMemory,PROM),可擦除只读存储器(Erasable Programmable Read-Only Memory,EPROM),电可擦除只读存储器(Electric Erasable Programmable Read-Only Memory,EEPROM)、磁盘或固态硬盘等。其中,存储器11用于存储程序,所述处理器12在接收到执行指令后,执行所述程序。
所述处理器12可能是一种集成电路芯片,具有数据的处理能力。上述的处理器12可以是通用处理器,包括中央处理器(Central Processing Unit,CPU)、网络处理器(Network Processor,NP)等。可以实现或者执行本申请中的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等,一般的,该处理器12可以是PCB603C01的芯片,该芯片灵敏度为100mV/g,量程为±50g,频带为5-10kHz,适用温度为-54℃~+121℃。
网络模块13用于通过网络建立电子设备10与外部通信终端之间的通信连接,实现网络信号及数据的收发操作。上述网络信号可包括无线信号或者有线信号。
可以理解,图1所示的结构仅为示意,电子设备10还可包括比图1中所示更多或者更少的组件,或者具有与图1所示不同的配置。图1中所示的各组件可以采用硬件、软件或其组合实现,
本申请的电子设备10还包括获取装置14,该获取装置14与该处理器12电连接,方法的获取装置14为脉诊仪,本申请的脉诊仪执行本申请的方法,以实现对寸、关、尺的检测。
在上述基础上,本申请还提供一种计算机可读存储介质,所述计算机可读存储介质包括计算机程序,所述计算机程序运行时控制所述计算机可读存储介质所在电子设备10执行下述的方法。
为了使本发明的实施过程更加清楚,下面将会结合附图进行详细说明。
图2为本申请实施例提供的一种浮中沉脉位识别方法的流程示意图;如图2所示,本申请提供一种浮中沉脉位识别方法,方法包括:
S101、获取待测样本的外加压力电压值、脉搏波幅值电压值、脉搏波电压值和测量时间。
通过脉诊仪分别获取待测样本的外加压力电压值、脉搏波幅值电压值、脉搏波电压值和测量时间,根据该脉搏波电压值和测量时间的关系,绘制该脉搏波电压值和测量时间的曲线图,根据该外加压力电压值和脉搏波幅值电压值的关系,绘制有关该外加压力电压值和脉搏波幅值电压值的曲线图,其中,该脉搏波电压值和测量时间的曲线图在绘制的时候,先将该脉搏波电压值和测量时间的数据中,提取一个周期的数据,以脉搏波电压值信号为同一个起点,将本周期内的脉搏波电压值和测量时间的曲线对齐设置,在该脉搏波电压值和测量时间的曲线图中,该外加压力电压值为该纵坐标的最大值对应的外加压力电压值最高峰点纵坐标;在绘制外加压力电压值和脉搏波幅值电压值的关系时,使用贝塞尔函数对外加压力电压值和脉搏波幅值电压值的曲线进行拟合,得到该外加压力电压值和脉搏波幅值电压值的曲线图。
S102、根据预设第一算法,通过外加压力电压值、脉搏波电压值和测量时间计算待测样本下单周期脉搏信号的平方差指数。
可选地,该预设第一算法为:
Figure GDA0003688697420000121
其中,f1、f2以及Amax均代表不同外加压力电压值时脉搏波电压值,f1为外加压力从初始值至最大值范围内的任意值,f2为外加压力从最大值至最小值范围内的任意值,Amax为外加压力最大时的脉搏波电压值。
由于该Amax为一个周期内外加压力电压值最高峰点纵坐标,则Amax将该脉搏波电压值和测量时间的曲线图分为两部分,一部分为上升段,一部分为下降段,在实际应用中,该Amax的左半部分为上升段,右半部分为下降段,根据上述预设第一算法,将外加压力电压值和脉搏波幅值电压值中上升段脉搏信号的纵坐标,外加压力电压值和脉搏波幅值电压值中下降段脉搏信号的纵坐标,对应压力下脉搏信号的周期,一个周期内外加压力电压值最高峰点纵坐标,带入该预设第一算法中,得到计算结果为平方差指数,对该平方差指数的计算结果得到的平方差指数进行判断。
S103、若平方差指数满足第一预设条件,则待测样本的脉搏为平;若平方差指数不满足第一预设条件,根据预设第二算法,通过外加压力电压值和脉搏波幅值电压值计算得到脉搏斜率指数,若脉搏斜率指数满足第二预设条件,则待测样本的脉搏为浮脉,若脉搏斜率指数不满足第二预设条件,则待测样本的脉搏为沉脉。
可选地,该第一预设条件为小于10,第二预设条件为大于0。
可选地,该预设第二算法为:
Figure GDA0003688697420000131
其中,A1为一个周期内外加压力电压值上升段范围内,任意外加压力P1下的脉搏波幅值电压值,A2为一个周期内外加压力电压值下降段范围内,任意外加压力P2下的脉搏波幅值电压值,Amax为一个周期内在最大外加压力Pmax下获得的最大脉搏波幅值电压值。
将该预设第一算法的平方差指数与第一预设条件进行比较,在实际应用中,该第一预设条件为小于10,若通过该预设第一算法得到的平方差指数小于10,则待测样本的脉象为平脉,若通过该预设第一算法得到的平方差指数不小于10,则进行第二次判断,使得预设第二算法,将一个周期内外加压力电压值最高峰点纵坐标,一个周期内外加压力电压值上升段的纵坐标,一个周期内外加压力电压值下降段的纵坐标,一个周期内外加压力电压值最高峰点横坐标,一个周期内外加压力电压值上升段的横坐标,一个周期内外加压力电压值下降段的横坐标,带入到该预设第二算法中,得到计算结果为脉搏斜率指数,并该通过预设第二算法的计算的结果脉搏斜率指数与该第二预设条件进行比较,若该脉搏斜率指数满足第二预设条件,则待测样本的脉搏为浮脉,若脉搏斜率指数不满足第二预设条件,则待测样本的脉搏为沉脉,即若该脉搏斜率指数大于0,则待测样本的脉搏为浮脉,若脉搏斜率指数不不大于0,则待测样本的脉搏为沉脉。
图3为本申请实施例提供的一种浮中沉脉位识别方法的脉搏波电压值与时间的关系图;图4为本申请实施例提供的一种浮中沉脉位识别方法的复制电压值与外加压力电压值的关系图;如图3和图4所示,为了进一步说明,在此进行举例说明,图3中标有Amax的线为0.6347V,标有f1的线为0.5447V,标有f2的线为0.7277V,其中,图3和图4均为通过实际检测得到的数据,使用上述预设第一算法,将图3中的数据代入到该预设第一算法中,计算得到该平方差指数为19.84640243,由于通过计算得到的平方差指数大于10,则表示该待测样本的脉象为浮脉或者沉脉,之后将图4中的数据代入到该预设第二算法中,计算得到斜率指数为0.103,由于使用预设第二算法得到的斜率指数大于0,则表示该待测样本的脉象为浮脉。
图5为本申请实施例提供的另一种浮中沉脉位识别方法的流程示意图;如图5所示,可选地,该根据预设第一算法,通过外加压力电压值、脉搏波电压值和测量时间计算待测样本下单周期脉搏信号的平方差指数的步骤之前还包括:
S201、以待测样本的脉搏的波谷为起点,计算不同压力下待测样本的脉搏波的平方差之和。
S202、根据预设条件选择主波振幅,并对主波振幅进行归一化处理。
上述以待测样本的脉搏的波谷为起点的归一化处理的步骤,待测样本的脉搏波谷为起点进行归一化处理,将有量纲数值转化为无量纲的数值,提升不同压力梯度分辨浮中沉模型的精度,提高浮中沉脉位识别模型的准确率和可靠性。
浮中沉脉位识别方法,方法包括:获取待测样本的外加压力电压值、脉搏波幅值电压值、脉搏波电压值和测量时间;根据预设第一算法,通过外加压力电压值、脉搏波电压值和测量时间计算待测样本下单周期脉搏信号的平方差指数;若平方差指数满足第一预设条件,则待测样本的脉搏为平;若平方差指数不满足第一预设条件,根据预设第二算法,通过外加压力电压值和脉搏波幅值电压值计算得到脉搏斜率指数,若脉搏斜率指数满足第二预设条件,则待测样本的脉搏为浮脉,若脉搏斜率指数不满足第二预设条件,则待测样本的脉搏为沉脉;本申请将原来需要中医获取的脉搏转化为可以通过设备获取,并且设备可以直接获取浮脉,中脉,沉脉,减少人为检测的误差,并且不依靠医师自身经验,有利于中医和中医理论的推广,即本申请是通过梯度施加压力的方式来获取脉搏波,通过绘制脉搏幅值随取脉压力变化的曲线来实现浮中沉脉位的判别。传统判断脉位主要根据施加的压力大小来确定脉位,但是由于个体腕部肌肉的松紧程度不同,浮中沉脉位对应的取脉压力也不同,因此通过单纯通过压力的大小判断脉位不够准确,本申请相比于传统单纯根据取脉压力的大小来判断脉位更加准确。
图6为本申请实施例提供的一种浮中沉脉位识别方法的脉宽示意图;如图6所示,本申请提供一种浮中沉脉位识别系统,系统用于实现第一方面中任意一项的方法,系统包括:获取模块20、计算模块30和判断模块40;
获取模块20,用于获取待测样本的外加压力电压值、脉搏波幅值电压值、脉搏波电压值和测量时间;
计算模块30,用于根据预设第一算法,通过外加压力电压值、脉搏波电压值和测量时间计算待测样本下单周期脉搏信号的平方差指数;
判断模块40,用于若平方差指数满足第一预设条件,则待测样本的脉搏为平;若平方差指数不满足第一预设条件,根据预设第二算法,通过外加压力电压值和脉搏波幅值电压值计算得到脉搏斜率指数,若脉搏斜率指数满足第二预设条件,则待测样本的脉搏为浮脉,若脉搏斜率指数不满足第二预设条件,则待测样本的脉搏为沉脉。
可选地,该预设第一算法为:
Figure GDA0003688697420000171
其中,f1、f2以及Amax均代表不同外加压力电压值时脉搏波电压值,f1为外加压力从初始值至最大值范围内的任意值,f2为外加压力从最大值至最小值范围内的任意值,Amax为外加压力最大时的脉搏波电压值。
图7为本申请实施例提供的另一种浮中沉脉位识别方法的脉宽示意图;如图7所示,可选地,该系统还包括预处理模块50,预处理模块50用于以待测样本的脉搏的波谷为起点,计算不同压力下待测样本的脉搏波的平方差之和;根据预设条件选择主波振幅,并对主波振幅进行归一化处理。
可选地,该预设第二算法为:
Figure GDA0003688697420000181
其中,A1为一个周期内外加压力电压值上升段范围内,任意外加压力P1下的脉搏波幅值电压值,A2为一个周期内外加压力电压值下降段范围内,任意外加压力P2下的脉搏波幅值电压值,Amax为一个周期内在最大外加压力Pmax下获得的最大脉搏波幅值电压值。
可选地,该第一预设条件为小于10,第二预设条件为大于0。
本申请提供一种电子设备,包括:存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现第一方面的浮中沉脉位识别方法。
本申请提供一种计算机可读存储介质,所述计算机可读存储介质包括计算机程序,所述计算机程序运行时控制所述计算机可读存储介质所在电子设备执行第一方面的浮中沉脉位识别方法。
以上仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (4)

1.一种浮中沉脉位识别方法,其特征在于,所述方法包括:
获取待测样本的外加压力电压值、脉搏波幅值电压值、脉搏波电压值和测量时间;
根据预设第一算法,通过所述外加压力电压值、所述脉搏波电压值和所述测量时间计算所述待测样本下单周期脉搏信号的平方差指数;
若所述平方差指数满足第一预设条件,则所述待测样本的脉搏为平;若所述平方差指数不满足第一预设条件,根据预设第二算法,通过所述外加压力电压值和所述脉搏波幅值电压值计算得到所述脉搏斜率指数,若所述脉搏斜率指数满足第二预设条件,则所述待测样本的脉搏为浮脉,若所述脉搏斜率指数不满足第二预设条件,则所述待测样本的脉搏为沉脉;
所述预设第一算法为:
Figure FDA0003688697410000011
其中,f1、f2以及Amax均代表不同外加压力电压值时脉搏波电压值,f1为外加压力从初始值至最大值范围内的任意值,f2为外加压力从最大值至最小值范围内的任意值,Amax为外加压力最大时的脉搏波电压值;
所述预设第二算法为:
Figure FDA0003688697410000021
其中,A1为一个周期内外加压力电压值上升段范围内,任意外加压力P1下的脉搏波幅值电压值,A2为一个周期内外加压力电压值下降段范围内,任意外加压力P2下的脉搏波幅值电压值,Amax为一个周期内在最大外加压力Pmax下获得的最大脉搏波幅值电压值;
所述第一预设条件为小于10,所述第二预设条件为大于0。
2.根据权利要求1所述的浮中沉脉位识别方法,其特征在于,所述根据预设第一算法,通过所述外加压力电压值、所述脉搏波电压值和所述测量时间计算所述待测样本下单周期脉搏信号的平方差指数的步骤之前还包括:
以所述待测样本的脉搏的波谷为起点,计算不同压力下所述待测样本的脉搏波的平方差之和;
根据预设条件选择主波振幅,并对所述主波振幅进行归一化处理。
3.一种浮中沉脉位识别系统,其特征在于,所述系统用于实现权利要求1-2任意一项所述的方法,所述系统包括:获取模块、计算模块和判断模块;
所述获取模块,用于获取待测样本的外加压力电压值、脉搏波幅值电压值、脉搏波电压值和测量时间;
所述计算模块,用于根据预设第一算法,通过所述外加压力电压值、所述脉搏波电压值和所述测量时间计算所述待测样本下单周期脉搏信号的平方差指数;
所述判断模块,用于若所述平方差指数满足第一预设条件,则所述待测样本的脉搏为平;若所述平方差指数不满足第一预设条件,根据预设第二算法,通过所述外加压力电压值和所述脉搏波幅值电压值计算得到所述脉搏斜率指数,若所述脉搏斜率指数满足第二预设条件,则所述待测样本的脉搏为浮脉,若所述脉搏斜率指数不满足第二预设条件,则所述待测样本的脉搏为沉脉;
所述预设第一算法为:
Figure FDA0003688697410000031
其中,f1、f2以及Amax均代表不同外加压力电压值时脉搏波电压值,f1为外加压力从初始值至最大值范围内的任意值,f2为外加压力从最大值至最小值范围内的任意值,Amax为外加压力最大时的脉搏波电压值;
所述预设第二算法为:
Figure FDA0003688697410000041
其中,A1为一个周期内外加压力电压值上升段范围内,任意外加压力P1下的脉搏波幅值电压值,A2为一个周期内外加压力电压值下降段范围内,任意外加压力P2下的脉搏波幅值电压值,Amax为一个周期内在最大外加压力Pmax下获得的最大脉搏波幅值电压值;
所述第一预设条件为小于10,所述第二预设条件为大于0。
4.根据权利要求3所述的浮中沉脉位识别系统,其特征在于,所述系统还包括预处理模块,所述预处理模块用于以所述待测样本的脉搏的波谷为起点,计算不同压力下所述待测样本的脉搏波的平方差之和;根据预设条件选择主波振幅,并对所述主波振幅进行归一化处理。
CN202110607165.5A 2021-06-01 2021-06-01 浮中沉脉位识别方法及系统 Active CN113288075B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110607165.5A CN113288075B (zh) 2021-06-01 2021-06-01 浮中沉脉位识别方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110607165.5A CN113288075B (zh) 2021-06-01 2021-06-01 浮中沉脉位识别方法及系统

Publications (2)

Publication Number Publication Date
CN113288075A CN113288075A (zh) 2021-08-24
CN113288075B true CN113288075B (zh) 2022-07-19

Family

ID=77326649

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110607165.5A Active CN113288075B (zh) 2021-06-01 2021-06-01 浮中沉脉位识别方法及系统

Country Status (1)

Country Link
CN (1) CN113288075B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114431840A (zh) * 2022-04-08 2022-05-06 北京大学深圳研究生院 一种脉搏采集装置、脉搏采集方法和系统
CN114711733B (zh) * 2022-06-07 2023-07-28 北京大学深圳研究生院 一种脉搏信号提取方法、装置、电子设备及存储介质

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0299396A2 (en) * 1987-07-13 1989-01-18 ISCO, Inc. Chromatographic system
WO1999026529A1 (fr) * 1997-11-20 1999-06-03 Seiko Epson Corporation Appareil de diagnostic lie aux formes d'impulsion, dispositif de controle de pression arterielle, dispositif de controle de la configuration des formes d'impulsion et dispositif de controle d'effet pharmacologique
CN104997493A (zh) * 2015-06-23 2015-10-28 田胜利 一种基于脉搏传感器的脉象分析方法
CN109464134A (zh) * 2017-09-08 2019-03-15 纳智源科技(唐山)有限责任公司 脉搏监测装置及系统
CN110517229A (zh) * 2019-07-31 2019-11-29 平安科技(深圳)有限公司 一种脉搏检测方法、系统、电子装置及存储介质
CN110853755A (zh) * 2019-11-07 2020-02-28 山西中医药大学 一种更年期女性脉象特点的研究方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6510687B1 (en) * 1996-06-14 2003-01-28 Sharav Sluices Ltd. Renewable resource hydro/aero-power generation plant and method of generating hydro/aero-power
CN103027669B (zh) * 2011-09-30 2016-09-14 Ge医疗系统环球技术有限公司 用于判断脉象浮沉度的方法和设备
CN108670209A (zh) * 2018-03-29 2018-10-19 中国科学院微电子研究所 一种自动识别中医脉象的方法和系统

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0299396A2 (en) * 1987-07-13 1989-01-18 ISCO, Inc. Chromatographic system
WO1999026529A1 (fr) * 1997-11-20 1999-06-03 Seiko Epson Corporation Appareil de diagnostic lie aux formes d'impulsion, dispositif de controle de pression arterielle, dispositif de controle de la configuration des formes d'impulsion et dispositif de controle d'effet pharmacologique
CN104997493A (zh) * 2015-06-23 2015-10-28 田胜利 一种基于脉搏传感器的脉象分析方法
CN109464134A (zh) * 2017-09-08 2019-03-15 纳智源科技(唐山)有限责任公司 脉搏监测装置及系统
CN110517229A (zh) * 2019-07-31 2019-11-29 平安科技(深圳)有限公司 一种脉搏检测方法、系统、电子装置及存储介质
CN110853755A (zh) * 2019-11-07 2020-02-28 山西中医药大学 一种更年期女性脉象特点的研究方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Jang, DG ; Park, SH and Hahn, M.Enhancing the Pulse Contour Analysis-Based Arterial Stiffness Estimation Using a Novel Photoplethysmographic Parameter.《IEEE JOURNAL OF BIOMEDICAL AND HEALTH INFORMATICS》.2015, *
一种区分检测中医脉学脉波信息的新方法;陈先农, 罗致诚.;《中国生物医学工程学报》;19841231;全文 *
中医脉象仪研究进展;庞瑜丹,柳建敏,张理兵.;《北京生物医学工程》;20141231;101-107 *

Also Published As

Publication number Publication date
CN113288075A (zh) 2021-08-24

Similar Documents

Publication Publication Date Title
CN113288075B (zh) 浮中沉脉位识别方法及系统
CN108024740B (zh) 血压测量方法、血压测量装置以及终端
CN204909419U (zh) 一种人体生理指标检测装置
AU2016201728B2 (en) A system and method for determining psychological stress of a person
CN109730723B (zh) 确定脉搏传输时间的方法、动脉硬化检测设备及系统
CN113707308A (zh) 一种医疗数据的分析装置及计算机可读存储介质
JP2023517017A (ja) 体内挿入用センサの校正感度を計算する方法
Basjaruddin et al. Measurement device for stress level and vital sign based on sensor fusion
CN109009004A (zh) 一种基于中医脉象分析的体质检测方法
CN115831371A (zh) 一种儿童体质健康评估系统和设备
CN110353699A (zh) 传感器脱落检测方法、装置及存储介质
CN103519791A (zh) 一种人体体温的采集系统及方法
CN112786185B (zh) 血压健康状态获取方法、装置和系统
CN113509151B (zh) 肌张力等级的评估方法和装置
CN115770028A (zh) 血压检测方法、系统、装置及存储介质
CN113693604B (zh) 一种基于佩戴者疲劳状态对肌张力等级的评估方法和装置
CN115752800A (zh) 一种基于互联网的医用无线体温监测系统
CN109009009A (zh) 血管评估方法、装置及智能终端
CN114340483A (zh) 一种血压标定选择方法及其建模方法
CN109118522B (zh) 一种呼吸道纤毛运动特征分析方法和装置
CN114176532A (zh) 一种测定cfPWV参数的临床验证方法及其应用系统
CN108742538B (zh) 基于大数据与人工智能的体征测量方法和医疗机器人系统
RU2230485C2 (ru) Способ определения концентрации глюкозы в крови человека (варианты)
CN106214153A (zh) 一种红外线臂围测量仪
KR101094993B1 (ko) 맥파검사를 통한 건강 자가진단장치

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
TA01 Transfer of patent application right

Effective date of registration: 20211028

Address after: 030032 room 07, 5 / F, block F, Internet of things Industrial Park, Xuefu Industrial Park, transformation and comprehensive reform demonstration zone, Taiyuan, Shanxi Province

Applicant after: Shanxi Hemei Internet smart traditional Chinese Medicine Industry Technology Research Institute Co.,Ltd.

Applicant after: SHANXI University OF CHINESE MEDICINE

Address before: 030051 No.65, xihou street, Shanglan Street office, Jiancaoping District, Taiyuan City, Shanxi Province

Applicant before: Shanxi Hengchao Technology Co.,Ltd.

Applicant before: SHANXI University OF CHINESE MEDICINE

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant