CN1030164C - 人体生物学年龄测极仪 - Google Patents
人体生物学年龄测极仪 Download PDFInfo
- Publication number
- CN1030164C CN1030164C CN 90102886 CN90102886A CN1030164C CN 1030164 C CN1030164 C CN 1030164C CN 90102886 CN90102886 CN 90102886 CN 90102886 A CN90102886 A CN 90102886A CN 1030164 C CN1030164 C CN 1030164C
- Authority
- CN
- China
- Prior art keywords
- pulse
- human
- parameter
- age
- human body
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
一种测报人的生物学年龄的仪器,由通用微机,专用硬件和软件组成,其特征是根据人体脉搏波形随年龄增长而变化的特点,统计大量不同年龄人体的脉搏波形图,定量地分析脉搏波上特征点参数之间的关系,建立一个能测定健康人生物学年龄的数学模型。本发明只要把脉搏传感器与受试者人的脉搏搏动点接触,很快就可以测报出受试者的生物学年龄,既精确又方便。
Description
一种用于测报人体生物学年龄的仪器,属诊断的其他方法或仪器类。
人体脉象计算机测报技术,国内已有许多研究,主要用来诊断和分析人体各种疾病,但是利用人体脉象测报人的生物学年龄却很少,而且测报不准确,测报复杂费时,不能在临床上应用。
本发明的目的是根据大量不同年龄健康人脉象的统计数据,建立一种数学模型并设计一种自动识别脉图特征点参数的方法,精确测报人体生物学年龄的仪器。
根据人体脉搏图象测报生物学年龄的原理是,由于人的动脉壁弹性随年龄增长而降低,动脉硬度增加,故中年以后脉象逐渐带弦,老年人多见弦脉,人的脉搏图重搏前波的相对位置随增龄而由低到高,主波的上1/3与下1/3交界处宽度越大,切脉感觉由平、缓、滑,向弦的方向发展,统计分析表明,脉图的有关指标值,随增龄而增高,因此,测量出不同年龄人体的脉图参数,可以计算出不同人的生物学年龄,在某种病理情况下(如高血压、冠心病等),生物学年龄可能显著大于时序年龄(也称日历年龄),从而可以推断人体可能有某种疾病及衰老的程度。由此可见,只要测量出人体脉搏图波形特征点的参数,就可推算出人的生物学年龄,关键是如何计算,本发明是以数百例20~70岁健康成年人的脉搏图的统计数据为依据建立了一种计算人的生物学年龄的通用数学模型,以此来测报人的生物学年龄,该数学模型的建立过程是,设变量Y与另外P个变量X1,X2,……Xp的内在联系是线性的,它的第α次试验数据是Yα;Xα1,Xα2…Xαp,其中α=1,2,…N,则这一组数据可以有如下的结构式:
Y1=β0+β1X11,+β2X12+……+βpX1p+ε1
Y2=β0+β1X21,+β2X22+……+βpX2p+ε2
┋
YN=β0+β1XN1+β2XN2+……+βpXNp+εN
其中β0,β1,β2,……βP是P+1个待估计参数,X1,X2,……Xp是P个可以精确测量或可控制的一般变量,ε1,ε2,……εN是N个相互独立且服从同一正态分布N(0,σ)的随机变量,这就是多元线性回归的数学模型,如果已知某个确定人的年龄Y1,及该人的脉搏图波形上的各参数X11,X12,……X1p(本发明设计一个专用程序如图4,计算此参数),则可介该联立方程,从而确定系数β0,β1,……βP。设b0,b1……bP分别是参数β0,β1,……βP的最小二乘估计,则该回归方程可写为Y=b0+b1X1+b2X2+……+bpXpb为
待定系数,本发明以N=103名20~70岁成年女子为例,在测量出她们的脉图参数Xp后(p=7),介出计算女性生物学年龄数学模型的系数b是,y=50.61647+0.3067497X1+23.54114X2-38.42367X3+48.23479X4-5.810995X5-7.702139X6-2.595945X7,其中,X1,X2,X3,X4,X5,X6,X7分别代表X1=M,X2=H1,X3=H2,X4=H3,X5=H5,X6=W/t,X7=t5/t4,y表示被测人的年龄。附图2是一个典型的人体脉搏波形,由附图所标出的符号,定义为M=t4/t0是与动脉硬度有关的参数,H1=(h2-h1)/h1,H2=h4/h1,H3=(h3-h1)/h1是与外周阻力有关的指标,H5=h5/h1也是与动脉硬度有关的指标,建立了数学模型之后,下一步便是如何利用计算机自动确定不同人脉搏波形上各特征点的参数。图3是几种常见的人的脉搏波形,呈三峰波形,也有双峰波形,本发明采用极值点和曲率混合结合经验阀值的方法,使用图4的脉图特征点自动识别的程序流程图,可以很方便地找出脉搏图上特征点A、B、C、D、E、F、G,流程图中的5是首先找出一个完整的脉图,并求各特征点一次二次导数及曲率,6是找脉图最高点MAX值(X,Y),7是判断X(横座标)是否大于170ms,如果是大于,则转到9最高点是D点,然后接下下去转到10,程序在A、D之间找B、C点。如果X不大于170ms,则转到8,确定最高点是B点,接下去程序转到11,在30~50%周期内找E、F点,12是判断D点是否已找到,未找到,则转到13,在B、E之间找C、D点,如已找到,则返回,例如要确定脉搏图中的C点,程序首先分析判断出C点在脉图上的某断区域,然后在该区域寻找曲率(绝对值)最大值的点C1,如果在C1点附近发现极小值点C2的话,那么C2点即为C点,否则,要在C1点附近逐渐向左、右搜索,直到找到的点满足已经设定的经验阀值时为止,此点即为要找的C点,本发明计算脉搏图曲率的公式是Ki=K1×{arctg〔K2×(Xi-Xi-3)〕-arctg〔K2×Ki+3-Xi)〕},K1.2为估计系数。
图1是本发明的方框图,图中1是脉搏传感器,它与人体脉搏占接触,并把脉搏跳动的压力信号变为电信号,该传感器有市售产品。2是滤波及放大,一般通用电子线路,3是A/D转换器,2,3做成的硬件,按放在通用计算机机箱内。A/D变换后的采样频率为200点/秒,4是通用微机。微机的显示器上可实时显示人的脉搏波形,该波形可用键盘控制,冻结在显示屏上,微机的打印机可以打印出受试者的生物学年龄。
本发明的特点是设备简单,测报人的生物学年龄准确。
Claims (2)
1、一种利用人体脉搏波形图象测报人体生物学年龄的仪器,由脉搏传感器(1)、电子线路(2)、A/D转换器(3)及微型计算机(4)组成,脉搏传感器采集人体的脉搏波形图象,经电子线路滤波放大和A/D转换后由微形计算机处理,其特征是根据大量不同年龄健康人脉搏波形图象参数的统计数据,建立用来测报人体生物学年龄的数学模型Y=b0+b1X1+b2X2+…+bpXp,其中X1=M1,X2=H1,X4=H3,X5=H5,X6=W/t,X7=t5/t4,表示人体脉搏波形图上特征点的参数,Y为被测人的年龄;
M和H5是与动脉硬度有关的参数,H1~H3是与动脉外周阻力有关的指标,微型计算机的专用软件计算出被测人体脉搏波形图上的参数X1~X7。
2、按权利要求1所述的仪器,其特征是微型计算机专用软件流程(5)~(13),找出由脉搏传感器测得的人体脉搏波形a~h上的特征点A、B、C、D、E、F、G并计算出参数X1~X7。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 90102886 CN1030164C (zh) | 1990-06-28 | 1990-06-28 | 人体生物学年龄测极仪 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 90102886 CN1030164C (zh) | 1990-06-28 | 1990-06-28 | 人体生物学年龄测极仪 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1057579A CN1057579A (zh) | 1992-01-08 |
CN1030164C true CN1030164C (zh) | 1995-11-01 |
Family
ID=4877766
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 90102886 Expired - Fee Related CN1030164C (zh) | 1990-06-28 | 1990-06-28 | 人体生物学年龄测极仪 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN1030164C (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2015039542A (ja) * | 2013-08-22 | 2015-03-02 | セイコーエプソン株式会社 | 脈波測定装置 |
-
1990
- 1990-06-28 CN CN 90102886 patent/CN1030164C/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN1057579A (zh) | 1992-01-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
DE69431575T2 (de) | Pulswellen- diagnosegerät | |
Wübbeler et al. | Verification of humans using the electrocardiogram | |
US7409242B2 (en) | Active muscle display device | |
EP3358485A1 (en) | General noninvasive blood glucose prediction method based on timing analysis | |
EP0699413B1 (en) | Apparatus and method for analyzing information relating to physical and mental condition | |
CN109730663B (zh) | 基于脉搏波传导速度非线性分析的血压评估方法 | |
CN101815466B (zh) | 用于对动脉nadi脉搏波形进行定量检测的非侵入式设备nadi tarangini | |
CN105193431A (zh) | 一种人体精神压力状态分析装置 | |
DE112009004271T5 (de) | Elektronisches Blutdruckmessgerät und Blutdruckmessverfahren | |
JP2009540915A (ja) | アレイ圧力センサーを利用した脈波分析装置 | |
KR20070046889A (ko) | 도플러 속도 스펙트럼 내에서 유속에 대한 조정가능한트레이싱 | |
CN111466882A (zh) | 一种智能中医手诊分析系统及方法 | |
CN113729648B (zh) | 一种基于多脉搏传感器的穿戴式脉诊手环系统 | |
CN113456033B (zh) | 生理指标特征值数据处理方法、系统及计算机设备 | |
CN111067500A (zh) | 一种基于ppg信号实现连续血压的监测系统 | |
JP2000166877A (ja) | 生体リズム検査装置及び生体リズム検査方法 | |
CN110200642A (zh) | 一种认知负荷与心理压力的测量方法及终端 | |
CN1030164C (zh) | 人体生物学年龄测极仪 | |
CN109567869A (zh) | 一种处理胎心率曲线上加速活动的方法及系统 | |
CN116942182A (zh) | 一种体质健康监控装置 | |
CN114176532B (zh) | 一种测定cfPWV参数的临床验证方法及其应用系统 | |
US20080097719A1 (en) | Service Providing System, Service Providing Method, and Program | |
CN115040126A (zh) | 一种孕产妇心理评估及辅助治疗系统 | |
Heimendinger et al. | Growth changes: Measuring the effect of an intervention | |
KR100768586B1 (ko) | 가속도맥파를 이용한 표준지표 산출 방법 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C19 | Lapse of patent right due to non-payment of the annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |