CN101310676A - 中医脉诊的分析***与方法 - Google Patents
中医脉诊的分析***与方法 Download PDFInfo
- Publication number
- CN101310676A CN101310676A CNA2007101055016A CN200710105501A CN101310676A CN 101310676 A CN101310676 A CN 101310676A CN A2007101055016 A CNA2007101055016 A CN A2007101055016A CN 200710105501 A CN200710105501 A CN 200710105501A CN 101310676 A CN101310676 A CN 101310676A
- Authority
- CN
- China
- Prior art keywords
- pulse
- aforementioned
- pulse wave
- blood pressure
- parameter
- 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.)
- Pending
Links
Images
Landscapes
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
本发明是关于一种中医脉诊的分析***与方法，可在脉诊中对脉象变化进行定量分析，尤其是中医脉象构成要素的「位」、「数」、「形」与「势」进行定量分析，大幅提升脉诊的准确性。本发明的分析***包含一脉波信号收集装置，用以收集并产生血压信号与心电信号；以及一信号处理单元，接收并分析连续的前述血压信号与心电信号，前述信号处理单元的执行包含：建构所测量到前述血压信号的部位参数；建构在由前述血压信号或心电信号所运算的心率参数；建构在由前述血压信号所分析的血压波形参数；建构在动脉压变异性(APV)及心率变异性(HRV)的频谱分析参数；以及基于上述参数界定中医脉象构成要素的定量指标。
Description
技术领域
本发明是关于一种中医脉诊的分析***与方法，尤其是一种关于中医脉象构成要素的定量分析的中医脉诊的分析***与方法。
背景技术
脉诊是中医根据脉象的变化来判断人体脏腑的气血、阴阳、生理与病理的状况；而脉象***，如弦脉、滑脉、平脉、短脉等28脉象，并不是只有心跳而已，脉诊需考虑到病患脉搏的位、数、形、势的变化。中医师作脉诊时，以指按压饶骨动脉，利用手指感觉脉动，进行病患整体状态的诊断，就″位″而言，桡骨动脉处即可分为寸、关、尺三部份。古代周学海提出脉诊时要力求把位、数、形、势讲得真切：“夫脉有四科，位、数、形、势而已。位者，浮沉尺寸也。数者，迟数促结也。形者，长短、广狭、厚薄、粗细、刚柔，犹算学家之有线、面、体也。势者，敛舒伸缩进退起伏之有盛衰也”。上述的位、数、形、势是判断脉象的纲领，但脉诊因为没有一个统一的标准和简便而又科学的表述方法，必需依靠人手指感觉出微妙的触觉而诊断，临床上并不能很好地实行。
自20世纪50年代所发展的先前技术开始以西方传来的脉搏描绘技术进行脉波图形化的尝试。一直到20世纪80年代所发展的检测重复性好、长期测试稳定性佳、取脉压力的读取和标记方便、具有足够的灵敏度及动态范围的仪器。台湾地区则自20世纪70年代起所发展的先前技术，如脉波仪，是由脉搏感应器、压力转换器、多频道记录器，将脉波仪与计算机相结合，使脉波图与心电图同步显现。另一种先前技术则利用傅立叶转换分析由压力转换器取得的脉波，以不同频率共振波强度作为不同脏腑的健康状态的指标。上述先前技术价格昂贵，脉位的选取仍靠人为操作，且报告的分析仍不够详细客观，仍需靠中医师判读，且相关数据无法供西医界了解与讨论有关。
现代医学的研究，构成脉搏的因素，主要是心脏搏动(包括搏出量与搏出力)所产生的压力、动脉管壁弹性与末梢抵抗力、血液粘稠度三个条件。正常脉搏的成因是心跳频率、心脏活动节律、心脏射血功能、动脉壁弹性、小动脉紧张度、血管充盈度及神经、内分泌调节功能等多种因素综合反映。
另外，压力脉图取得的生理信息具有一定的局限性，动脉脉搏除发出压力搏动的信号之外，还有管腔容积、血流速度、脉管的三度空间运动等多种信息。脉象图的研究大多还局限于28脉的范围内，如弦脉、滑脉、疾脉等，某脉的脉象图实质上可能是某复合脉的脉搏波图，况且某些脉象的压力脉图之间没有明显的区别。因此，目前用压力脉图的指标，主要集中在″位″与″形″上的探讨，有其局限性，不能将脉象的脉位、脉宽、脉长、脉率、脉律、脉力、脉流、脉体等指标全部定量地反映出来，尚需要配合其它方法进行脉象研究。
发明内容
本发明的目的在于提供一种中医脉诊的分析***，通过信号收集定量分析中医脉象构成要素，尤其是对「形」、「数」、「势」进行定量分析。
本发明的又一目的在于提供一种中医脉诊的分析方法，参考连续的血压信号与心电信号，以界定中医脉象构成要素的定量指标。
本发明的另一目的在于提供一种中医脉诊的分析***与方法，从西医的角度来定义「形」、「数」、「势」分析。
为达到前述目的，本发明提供一种中医脉诊的分析***，包含一脉波信号收集装置，用以收集并产生血压信号与心电信号；以及一信号处理单元，接收并分析连续的前述血压信号与心电信号，前述信号处理单元的执行包含：建构所测量到前述血压信号的部位参数；建构在由前述血压信号或心电信号所运算的心率参数；建构在由前述血压信号所分析的血压波形参数；建构在心率变异性(HRV)及动脉压变异性(APV)的频谱分析参数；以及基于上述参数界定中医脉象构成要素的定量指标。
本发明又提供一种中医脉诊的分析方法，参考连续的血压信号与心电信号，该分析方法包含：建构所测量到前述血压信号的部位参数；建构在由前述血压信号或心电信号所运算的心率参数；建构在由前述血压信号所分析的血压波形参数；建构在动脉压变异性(APV)及心率变异性(HRV)的频谱分析参数；以及基于上述参数界定中医脉象构成要素的定量指标。
本发明的优点为把难以定量的脉象，通过自动化运算出定量指标，使得脉象判定标准上有一致性，避免人为诊断的误差，而提升脉诊的准确性，增进脉诊仪的效能，并与西医结合，以上皆为本发明优异的效能。
配合所附附图、实施例的详细说明及申请专利范围，对上述及本发明的其它目的与优点进行了详述。然而，当可了解所附附图纯为解说本发明的精神而设，不当视为本发明范畴的定义。有关本发明保护范围的定义，请参照所附的权利要求书。
附图说明
图1为本发明中医脉诊的分析***的***方块图；
图2为本发明中医脉诊的分析方法5秒钟原始血压波形图；
图3为本发明中医脉诊的分析方法5分钟原始血压波形图；
图4为图3血压信息运算得出平均血压(MSBP)及平均心跳间距(MPPI)的结果图。
图中，
10脉波信号收集装置
11血压收集器
111压力转换器 112放大器
113模拟数字转换器
12心电信号收集器
121电极 122放大器
123模拟数字转换器
20信号处理单元
具体实施方式
定义
本发明中的“脉诊”是指中医根据脉象的变化来判断人体脏腑的气血、阴阳、生理与病理的状况。
本发明中的“脉象”是指中医师以指按压病患脉搏，利用手指感觉脉动的变化。
本发明中的中医脉象构成要素的“位”是指量测脉搏的位置。
本发明中的中医脉象构成要素的“数”是指每单位时间内的脉搏次数。
本发明中的中医脉象构成要素的“形”是指脉的形状，脉相为波。
本发明中的中医脉象构成要素的“势”是指脉相的趋势(trend)，趋势即为改变的趋向(variation oftrend)。
本发明中的“QRS波群(QRS complex)”是指心电图中典型的QRS波群包括三个紧密相连的波，第一个向下的波称为Q波，继Q波后的一个高尖的直立波称为R波，R波后向下的波称为S波。因其紧密相连，且反映了心室电激动过程，故统称为QRS波群。这个波群反映了左、右两心室的除极过程。
首先，参阅图1，图1为本发明中医脉诊的分析***的***方块图。在本发明的较佳实施例中，本发明的中医脉诊的分析***，包含一脉波信号收集装置10，该脉波信号收集装置10进一步包含一血压收集器11和一心电信号收集器12，而该心电信号收集器12在本发明***中为选择性安装的装置，其中，血压信号收集器11从人体特定部位收集血压信号，例如桡骨动脉处的寸、关、尺位置，或是人体其它部位。收集到的血压信号，经由压力转换器111将压力信号转为电信号，再经放大器112将信号放大，放大后的血压信号再经由模拟数字转换器113(analog-digital converter)将信号转换为数字信号，取样频率为2048Hz，持续数分钟之后，取得并产生用数字信号表示的连续血压信号，以时间为横轴将血压信号逐点划出立即性的血压波形信号。
另一方面，心电信号收集器12可经由电极121收集心电信号，心电信号经放大器122将信号放大，放大后的心电信号再经由模拟数字转换器123将信号转换为数字信号，取样频率为256Hz，持续数分钟之后，取得并产生用数字信号表示的连续心电信号，以时间为横轴将血压信号逐点划出立即性的心电波形信号。
一脉波处理单元20接收并分析连续的前述血压信号和心电信号，该脉波处理单元20可为一个人计算机、PDA、行动电话内含的微电脑或手表内含的微电脑。该血压信号是从人体特定部位收集，以此建构所测量到前述血压信号的部位参数，界定该部位参数为中医脉象构成要素「位」。
脉波处理单元20接收到前述的血压波形信号，在方块201中，分析相对高点即为一收缩压，相对低点即为一舒张压，利用该收缩压数值及该舒张压数值经程序运算出一平均压(平均压＝1/3收缩压+2/3舒张压)，该收缩压数值减舒张压数值即为一脉搏压(脉搏压＝收缩压-舒张压)。在方块202中，脉波处理单元20利用上述血压数字信号做一次微分(dp/dt)后，每次脉搏的最高值减去最低值，便可得到心脏收缩力(contractility)数值。而心搏量(srtoke volume)是利用以下公式得出：
SV＝KPsa(1+Ts/Td)
K为一常数，Psa是单个脉博血压曲线的面积，Ts是收缩持续时间，Td是舒张持续时间。若将上述心搏量乘以心率即得一心输出量(心输出量＝心搏量×心率)。在方块203中，脉波处理单元20进行血压谐振波分析(harmonicanalysis)，是将前述得到的血压数字信号以快速傅立叶转换为频谱，并利用积分的方式来订量每一个特定频率区间的功率。每一频谱以两度空间直方图表达，或以数值方式输出，运算出动脉压变异性数值，其中随心跳倍率变化的超高频成分即为血压谐振波。而动脉顺应性(artery compliance)是利用以下公式得出：
Cp＝(SV/K)×(bexpbp/expbp s-expbp d)
另外，动脉阻力值(resistance ofartery)是利用以下公式得出：
Zc＝Pf/Qf
Zc即为动脉阻力，Pf为血压，Qf为血流。前述方块201中的收缩压、舒张压、平均压、及脉博压，是为脉波处理单元20接收到血压波形信号后，经前述运算得到的数值，此种经由血压波形分析而得到的参数，为血压波形参数，脉波处理单元20界定该血压波形参数为中医脉象构成要素的「形」，是以收缩压、舒张压、平均压、及脉博压为定量指标。而前述方块202、及方块203中心脏收缩力、心搏量、心输出、血压谐振波分析、动脉顺应性、及动脉阻力值同时为中医脉象构成要素「数」与「形」的定量指标。
另一方面，脉波处理单元20接收到前述的血压波形信号，经运算后同时得到心跳周期间距(PPI)的数值，即为方块210中的心率参数，或是利用心电波形信号，同样可得到方块210中的心率参数，脉波处理单元20界定该心率参数为中医脉象构成要素的「数」，是以心率为定量指标。
变异性分析是指任何一个参数经长时间收集后都可为其趋势或变动进行分析，如心率经长时间累积足够的数量后可进行心率变异性分析、及动脉压经长时间累积足够的数量后可进行动脉压变异性分析。
在信号处理单元20内，前述心率参数与血压波形参数经由方块220的频域分析法(Frequency domain analysis)分析，在动脉压变异性分析(APV)方面，其中血管交感神经活性(vascular sympathetic activity)是将前述得到的数字血压信号以快速傅立叶转换为频谱，并利用积分的方式来定量每一个特定频率区间的功率，得到方块230的动脉压的功率频谱。每一频谱以两度空间直方图表达，或以数值方式输出，运算出动脉压变异性数值，其中方块231的低频成分(BLF)可作为血管交感神经活性的指标。
在信号处理单元20内，前述心率参数与血压波形参数经由方块220的频域分析法(Frequency domain analysis)分析，在心率变异性分析(HRV)方面，是将前述运算得到的心电信号以快速傅立叶转换为频谱，并利用积分的方式来订量每一个特定频率区间的功率，得到方块240的心率功率频谱。每一频谱以两度空间直方图表达，或以数值方式输出，运算出心率变异性数值，并以方块241的心脏交感神经活性LF/HF为心脏交感神经活性(cardiac sympatheticactivity)的指标。方块242的心脏副交感神经活性(cardiac vagal activity)以HF为心脏迷走神经活性的指标，而LF则视为交感和副交感神经功能的统合指针。
信号处理单元20界定变异性的频谱分析参数为中医脉象构成要素的「势」，是以血管交感神经活性、心脏交感神经活性、心脏副交感神经活性、以及感压反射灵敏度分析为定量指标。
前述方块230和231中的动脉压的功率频谱和血管交感神经活性的数值，是为该脉波处理单元20接收到血压波形参数后，经过连续时间运算得到的数值，此种经由血压波形参数连续时间分析而得到的参数，是为动脉压变异性(APV)的频谱分析参数，脉波处理单元20界定动脉压变异性(APV)的频谱分析参数为中医脉象的「形之势」。
前述方块240、241和242中的心率的功率频谱、心脏交感神经活性、和心脏副交感神经活性的数值，是为该脉波处理单元20接收到心率参数后，经过连续时间运算得到的数值，此种经由心率参数连续时间分析而得到的参数，是为心率变异性(HRV)的频谱分析参数，脉波处理单元20界定心率变异性(HRV)的频谱分析参数为中医脉象的「数之势」。
由于心率变异与血压变异是息息相关；意即当血压有所变动时，心脏会调节其心率以维持血压的恒定性。于是「数之势」与「形之势」其实是环环相扣，互为消长的一套「势」的***；方块250的血压和心率的交叉频谱和方块251的感压反射灵敏度，就是在表达「数之势」与「形之势」分析之下的统合的「势」。以下进一步说明脉波处理单元20的频域分析工作内容。
动脉压变异性分析(Arterial Pressure Variability Analysis)
动脉压(APV)的收集，可由一脉波信号收集装置进行，是利用血压信号收集器如臂式血压计、腕式血压计、手指体积描记器(finger plethysmography)、笔式脉波传感器、压平眼压计(tonometry)、环式血压计等收集血压波形信号，血压波形信号经放大器放大10,000倍并滤波后，再经12位(bit)模拟数字转换器(PCL-818L，Advantech，Taiwan)连接至信号处理单元20，如IBM兼容计算机，并以1024Hz取样速率存录至信号处理单元20。
在一信号处理单元内将动脉压信号进行重新取样(resample)以维持其时间连续性。首先消除信号的基准线飘移(baseline drift)以防止低频带的干扰，原本连续的血压信号切分为每64秒(4096点)为一组(或称为一个分析窗)进行快速傅立叶转换(Fast Fourier transform，FFT)估算功率密度频谱(powerspectral density)，同时采用海明窗(Hamming window)运算以避免频谱中个别频率成份的互相渗漏(leakage)作用。
本发明以积分的方式定量功率密度频谱中各频带的功率，包括动脉压频谱中的低频功率(BLF，0.04-0.15Hz)、高频功率(BHF，0.15-0.4Hz)等量化参数，其中，动脉压的低频功率为血管的交感神经活性指标。
心率变异性分析(Heart Rate Variability Analysis)
心率变异性(HRV)的定义，即是心跳与心跳之间隔时间(interval)改变的变异度。心率的收集，可由本发明的脉波信号收集装置收集并产生的血压波形信号，经由信号处理单元接收后，运算得到心跳周期，即心率。
另一方面，心率的收集也可由各式心电信号收集器12记录至少5分钟信息。心电信号经放大一千倍，并通过滤波器(0.68-16Hz)之后，用一8位(bit)的模拟数字心电图信号转换器以256Hz采样记录。
心率变异频谱分析的程序首先以QRS尖峰检测程序将心电信号中的最高点(peaks)找出，视为每次心跳的QRS波群(QRS complex)。自每个QRS中，测量其振幅(amplitude)和持续时间(duration)等参数，并将各参数的平均值和标准差算出作为标准模版。
接下来，每个QRS都以此模版进行比对，如果其中QRS的比对结果落在标准模版三个标准差之外，将被认为是噪声或异位心跳(ectopic beat)而删除。将合格QRS波群的R点作为该心跳的时间点，而心跳与心跳间的时间差作为一次的心跳周期(RR interval)。接下来进行心跳周期的过滤程序，首先将所有心跳周期的平均值和标准差算出，再进行所有心跳周期的确认。如果某一心跳周期落在4个标准差之外，它会被认为是错误或不稳定(non-stationary)信号而删除。
心率估算方面，首先以尖峰检测程序将R波产生的时间点找出，由其前后相隔时间的倒数可算得每次心脏跳动的心率，并以取样-保值程序维持其时间连贯性，此取样-保值步骤的更新率为每秒16次。将前后延伸的心率及心电信号切分为每64秒(1024点)为一组(或称为一个分析窗)。在每一个分析窗中，前后并有50％的重迭，首先消除信号的直线飘移以防止低频带的干扰，亦采用海明(Hamming)运算以避免频谱中个别频率成份的互相干扰。而信号再以快速傅立叶转换(Fast Fourier transform)的方式转为功率频谱。
本发明以积分的方式定量心率频谱中2个频率成份，包括低频(low-frequency，LF，0.04-0.15Hz)和高频(high-frequency，HF，0.15-0.4Hz)成份。同时求出总功率(total power)和低频/高频比值(LF/HF)，其中，HF为心脏迷走神经活性的指标，并以LF/HF为心脏交感神经活性的指标，而LF则视为交感和副交感神经功能的统合指针。
动脉压/心率平均频谱(APV/HRV Average Periodogram Analysis)
在数据收集后进行动脉压/心率平均频谱及动脉压-心率分析转换函数分析。每次分析须取896秒的动脉压和心率数据，这些数据可分为55个分析窗(每个分析窗1024点)，前后并有50％的重迭。对于每一个分析窗，首先消除信号的基准线飘移以防止低频带的干扰，亦采用Hamming window运算以避免频谱中个别频率成份的互相渗漏(leakage)作用，再将这六个分析窗的频谱平均，则得到平均频谱(average periodogram)。平均之后可衰减随机的噪声而突显再生性高的频谱成分。
动脉压/心率交叉频谱(cross spectrogram)及转换函数分析
得到动脉压/心率的平均频谱后，接下来，进行交叉频谱分析，交叉频谱分析的来源同样是上述55组资料：
一、计算动脉压和心率的平均频谱线性相关程度(coherence)，运算定义如下：
κ2(f)＝|SHB(f)2/[SBB(f)×SHH(f)]
SBB(f)及SHH(f)分别代表动脉压及心率的功率频谱信号，而SHB(f)则为它们的交叉频谱。运算结果分布于0到1之间，表示两信号间各频率的线性相关程度，同时表示其相对转换函数真确性。若其值大于0.5表示相关性具有统计意义，方可信其转换函数。
二、转换函数(transfer function)：
H(f)＝SHB(f)/SBB(f)
振幅定义为：
{[HR(f)]2+[HI(f)]2}1/2
其中HR(f)和HI(f)为别为H(f)的实部及虚部的值，转换幅度单位为bpm/mmHg。
感压反射灵敏度评估
过去许多研究感压反射，需通过扰动血压后，使受试者的感压反射而造成心率的极大变动，而求得感压反射灵敏度。本发明非侵入性在线连续测量动脉压及心率，便能连续分析感压反射灵敏度，不会干扰受试者。
感压反射灵敏度以两种方式进行评估：
一、动脉压-心率转换函数法：根据动脉压-心率转换函数分析的理论，动脉压-心率转换函数的幅度是反应心率感压反射的灵敏度。依所量取的频率范围，此数值又分为低频(BrrLF)和高频(BrrHF)两者。
二、动脉压及心率的连续变化：体动脉压在自然情况下并非一成不变，而是随机产生一些微小的上升或下降，在此血压变化的瞬时，心率因感压反射也会产生相对的下降或上升。将动脉压持续上升且心率持续下降，超过三次的动脉压-心率数据组进行线性回归可得一斜率(BrrA)。反之，将动脉压持续降压且心率持续上升，超过三次的动脉压-心率数据组进行线性回归也可以得一斜率(BrrD)。(BrrA)与(BrrD)皆能反应当时的感压反射灵敏度。因此，前后连续三次以上的动脉压和心率变化，即可求得感压反射灵敏度。
根据前述本发明的分析***所实施的中医脉诊的分析方法，参考连续的血压信号与心电信号，该分析方法包含：建构所测量到前述血压信号的部位参数，并界定部位参数为中医脉象构成要素的「位」；建构在由前述血压信号或心电信号所运算的心率参数，并界定心率参数为中医脉象构成要素的「数」；建构在由前述血压信号所分析的血压波形参数，并界定血压波形参数为中医脉象构成要素的「形」；建构在动脉压变异性(APV)及心率变异性(HRV)的频谱分析参数，并界定变异性的频谱分析参数为中医脉象构成要素的「势」；以及基于上述参数界定中医脉象构成要素的定量指标。其中，以心率定量分析前述「数」，以平均压、收缩压、舒张压及脉搏压定量分析前述「形」，以心输出、心脏收缩力、心搏量、血压谐振波分析、动脉顺应性及动脉阻力定量分析前述「数」及「形」，以血管交感神经活性、心脏交感神经活性、心脏副交感神经活性、以及感压反射灵敏度分析定量分析前述「势」。
因为，中国医学利用阴阳虚实来辩证，诊断病人，将病人体质与疾病分为阴证、阳证、虚证及实证。再通过中医治疗方式，将病人阴阳虚实的偏颇调节，使之趋向平衡而达到健康。中医的阴阳虚实均与全身的各个***有相关，中医书籍更指出主要是与西方医学所谓的自主神经***最有关系，一直以来缺乏科学证据来证实。
本发明的分析***所实施的中医脉诊的分析方法，通过测量多项生理指针，如平均压、收缩压、舒张压、脉搏压、动脉顺应性、及动脉阻力定量分析等，提供了心率变异性分析和动脉压变异性分析的「势」的信息，突破中医脉象上定量的限制。并可建构一预设的脉象范围，如28脉象；以及判读上述参数的脉象构成要素以对应该脉象范围。
本发明实施例运算结果如图2、图3和图4。参阅图2，图2为5秒钟原始血压波形图，是显示本发明中医脉诊的分析***，经由8点平均1点的方式(bunching)缩减频率为32Hz。原本无限延伸的血压信号切分为32秒(1024点)为一个分析窗。以时间为横轴将血压信号逐点划出即得立即性的血压波型。图中的5秒钟原始血压波形(BP)，经运算后同时得到的收缩压(SBP)及心跳周期间距(PPI)的数值。其中原始血压波形及收缩压对应于脉诊的「形」，而心跳周期间距则对应于脉诊的「数」。
参阅图3，图3为5分钟原始血压波形(BP)图，经运算后同时得到的收缩压(SBP)及心跳周期间距(PPI)的数值。其中，原始血压波形及收缩压对应于脉诊的「形」，而心跳周期间距则对应于脉诊的「数」，而由收缩压所表现的波形变动，及心跳周期间距所表现的波形变动则对应于脉诊的「势」。
参阅图4，图4为图3中血压信息运算得出平均血压(MSBP)及平均心跳间距(MPPI)的结果图，并经频谱分析后得到心率功率密度(HPSD)及血压功率密度(BPSD)。并将特定频率分为极低频(BVLF、HVLF)、低频(BLF、HLF)及高频(BHF、HHF)成分。
根据上述图示结果，本发明安排有经验的中医师对病患进行评估，同时测量心率变异各种参数，已有中医阴阳虚实辩证及自主神经相关的科学证据。首先，将不同体质分组，接着，与各种心率变异分析参数与体质分数作线性回归分析，结果发现阴阳平衡与交感副交感平衡有很好的正相关，阳证体质者较阴证体质者有较低的副交感神经活性与较高的交感神经活性，阳实者比阳虚者有较高的交感神经活性。
经由中医的适当调理体质后，发现受试者除了改变辩证分数外，同时亦改变自主神经参数。由此可知，中医的阴阳虚实辩证确实与心脏自主神经有明确相关，其中的LF/HF可代表阳证强弱，交感副交感平衡可代表阳阴平衡，LF可代表实虚体质的强弱，且体质调理变化与自主神经改变也相关。BVLF、BLF、BHF、HVLF、HLF、HHF的数值不但可以了解交感神经、副交感神经的状态，更提供了中医师在阴阳虚实诊断上的参考。
Claims (25)
1.一种中医脉诊的分析***，其特征在于，包含：
一脉波信号收集装置，用以收集并产生血压信号；以及
一信号处理单元，接收并分析连续的前述血压信号，
其特征在于，前述信号处理单元的执行包含：
建构所测量到前述血压信号的部位参数；
建构在由前述血压信号所运算的心率参数；
建构在由前述血压信号所分析的血压波形参数；
建构在动脉压变异性及心率变异性的频谱分析参数；以及
基于上述参数界定中医脉象构成要素的定量指标。
2.根据权利要求1所述的中医脉诊的分析***，其特征在于，前述脉波信号收集装置能进一步收集并产生心电信号，且前述信号处理单元能接收并分析连续的心电信号。
3.根据权利要求2所述的中医脉诊的分析***，其特征在于，前述信号处理单元的执行，进一步包含建构在由前述心电信号所运算的心率参数。
4.根据权利要求1所述的中医脉诊的分析***，其特征在于，前述信号处理单元界定部位参数为中医脉象构成要素的位。
5.根据权利要求1或3所述的中医脉诊的分析***，其特征在于，前述信号处理单元界定心率参数为中医脉象构成要素的数。
6.根据权利要求1所述的中医脉诊的分析***，其特征在于，前述信号处理单元界定血压波形参数为中医脉象构成要素的形。
7.根据权利要求1所述的中医脉诊的分析***，其特征在于，前述信号处理单元界定变异性的频谱分析参数为中医脉象构成要素的势。
8.根据权利要求7所述的中医脉诊的分析***，其特征在于，前述信号处理单元界定动脉压变异性的频谱分析参数为中医脉象的形之势。
9.根据权利要求7所述的中医脉诊的分析***，其特征在于，前述信号处理单元界定心率变异性的频谱分析参数为中医脉象的数之势。
10.根据权利要求5所述的中医脉诊的分析***，其特征在于，前述数以心率为定量指标。
11.根据权利要求6所述的中医脉诊的分析***，其特征在于，前述形以平均压、收缩压、舒张压及脉搏压为定量指标。
12.根据权利要求10或11所述的中医脉诊的分析***，其特征在于，前述数与形皆以心输出、心脏收缩力、心搏量、血压谐振波分析、动脉顺应性及动脉阻力为定量指标。
13.根据权利要求7所述的中医脉诊的分析***，其特征在于，前述势以血管交感神经活性、心脏交感神经活性、心脏副交感神经活性、以及感压反射灵敏度分析为定量指标。
14.一种中医脉诊的分析方法，参考连续的血压信号，其特征在于，该分析方法包含：
建构所测量到前述血压信号的部位参数；
建构在由前述血压信号所运算的心率参数；
建构在由前述血压信号所分析的血压波形参数；
建构在动脉压变异性及心率变异性的频谱分析参数；以及
基于上述参数界定中医脉象构成要素的定量指标。
15.根据权利要求14所述的中医脉诊的分析方法，其特征在于，该方法进一步包含参考连续的心电信号。
16.根据权利要求15所述的中医脉诊的分析方法，其特征在于，该方法进一步包含建构在由前述心电信号所运算的心率参数。
17.根据权利要求14所述的中医脉诊的分析方法，其特征在于，该方法进一步包含：界定部位参数为中医脉象构成要素的位。
18.根据权利要求14所述的中医脉诊的分析方法，其特征在于，该方法进一步包含：界定心率参数为中医脉象构成要素的数。
19.根据权利要求14所述的中医脉诊的分析方法，其特征在于，该方法进一步包含：界定血压波形参数为中医脉象构成要素的形。
20.根据权利要求14所述的中医脉诊的分析方法，其特征在于，该方法进一步包含：界定变异性的频谱分析参数为中医脉象构成要素的势。
21.根据权利要求18所述的中医脉诊的分析方法，其特征在于，该方法进一步包含：以心率定量分析前述数。
22.根据权利要求19所述的中医脉诊的分析方法，其特征在于，该方法进一步包含：以平均压、收缩压、舒张压及脉搏压定量分析前述形。
23.根据权利要求21或22所述的中医脉诊的分析方法，其特征在于，该方法进一步包含：以心输出、心脏收缩力、心搏量、血压谐振波分析、动脉顺应性及动脉阻力定量分析前述数及形。
24.根据权利要求20所述的中医脉诊的分析方法，其特征在于，该方法进一步包含：以血管交感神经活性、心脏交感神经活性、心脏副交感神经活性、以及感压反射灵敏度分析定量分析前述势。
25.根据权利要求14所述的中医脉诊的分析方法，其特征在于，该方法进一步包含：建构一预设的脉象范围；以及判读上述参数的脉象构成要素以对应该脉象范围。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNA2007101055016A CN101310676A (zh) | 2007-05-24 | 2007-05-24 | 中医脉诊的分析***与方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNA2007101055016A CN101310676A (zh) | 2007-05-24 | 2007-05-24 | 中医脉诊的分析***与方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN101310676A true CN101310676A (zh) | 2008-11-26 |
Family
ID=40099652
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CNA2007101055016A Pending CN101310676A (zh) | 2007-05-24 | 2007-05-24 | 中医脉诊的分析***与方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN101310676A (zh) |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101972141A (zh) * | 2010-08-17 | 2011-02-16 | 中国科学院苏州纳米技术与纳米仿生研究所 | 面向基层医疗的心脑血管特征信号实时感知与监护*** |
| CN102319061A (zh) * | 2011-09-14 | 2012-01-18 | 魏伯卿 | 中医电子号脉机 |
| CN103027668A (zh) * | 2011-09-30 | 2013-04-10 | Ge医疗***环球技术有限公司 | 用于中医脉诊的脉搏信号的处理 |
| CN103049631A (zh) * | 2011-10-13 | 2013-04-17 | 客隆生医科技股份有限公司 | 计算机脉波诊断分析装置及其方法 |
| CN104688200A (zh) * | 2013-12-10 | 2015-06-10 | 张国源 | 心脏信息分析*** |
| CN105725983A (zh) * | 2016-01-07 | 2016-07-06 | 深圳市和来科技有限公司 | 一种周围动脉硬化的早期筛查方法及*** |
| CN107126203A (zh) * | 2016-02-28 | 2017-09-05 | 南京全时陪伴健康科技有限公司 | 人体健康信息实时测量显示*** |
| CN110575151A (zh) * | 2018-06-11 | 2019-12-17 | 当代汉云企业有限公司 | 脉诊测量装置及方法 |
| CN111212602A (zh) * | 2017-10-19 | 2020-05-29 | 东洋纺株式会社 | 生理信息显示***和训练方法 |
| CN113545760A (zh) * | 2021-07-29 | 2021-10-26 | 西安仙峒科技有限责任公司 | 基于模态分析技术的中医脉诊测试仪及其使用方法 |
-
2007
- 2007-05-24 CN CNA2007101055016A patent/CN101310676A/zh active Pending
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101972141A (zh) * | 2010-08-17 | 2011-02-16 | 中国科学院苏州纳米技术与纳米仿生研究所 | 面向基层医疗的心脑血管特征信号实时感知与监护*** |
| CN101972141B (zh) * | 2010-08-17 | 2012-07-04 | 中国科学院苏州纳米技术与纳米仿生研究所 | 面向基层医疗的心脑血管特征信号实时感知与监护*** |
| CN102319061A (zh) * | 2011-09-14 | 2012-01-18 | 魏伯卿 | 中医电子号脉机 |
| CN102319061B (zh) * | 2011-09-14 | 2013-12-04 | 魏伯卿 | 中医电子号脉机 |
| CN103027668A (zh) * | 2011-09-30 | 2013-04-10 | Ge医疗***环球技术有限公司 | 用于中医脉诊的脉搏信号的处理 |
| CN103049631A (zh) * | 2011-10-13 | 2013-04-17 | 客隆生医科技股份有限公司 | 计算机脉波诊断分析装置及其方法 |
| CN104688200A (zh) * | 2013-12-10 | 2015-06-10 | 张国源 | 心脏信息分析*** |
| CN105725983A (zh) * | 2016-01-07 | 2016-07-06 | 深圳市和来科技有限公司 | 一种周围动脉硬化的早期筛查方法及*** |
| CN105725983B (zh) * | 2016-01-07 | 2020-12-08 | 深圳市和来科技有限公司 | 一种周围动脉硬化的早期筛查方法及*** |
| CN107126203A (zh) * | 2016-02-28 | 2017-09-05 | 南京全时陪伴健康科技有限公司 | 人体健康信息实时测量显示*** |
| CN111212602A (zh) * | 2017-10-19 | 2020-05-29 | 东洋纺株式会社 | 生理信息显示***和训练方法 |
| CN110575151A (zh) * | 2018-06-11 | 2019-12-17 | 当代汉云企业有限公司 | 脉诊测量装置及方法 |
| CN113545760A (zh) * | 2021-07-29 | 2021-10-26 | 西安仙峒科技有限责任公司 | 基于模态分析技术的中医脉诊测试仪及其使用方法 |
| CN113545760B (zh) * | 2021-07-29 | 2023-08-18 | 西安仙峒科技有限责任公司 | 基于模态分析技术的中医脉诊测试仪及其使用方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101310676A (zh) | 中医脉诊的分析***与方法 | |
| Lu et al. | A comparison of photoplethysmography and ECG recording to analyse heart rate variability in healthy subjects | |
| Li et al. | Characters available in photoplethysmogram for blood pressure estimation: beyond the pulse transit time | |
| Bolanos et al. | Comparison of heart rate variability signal features derived from electrocardiography and photoplethysmography in healthy individuals | |
| CN104382571B (zh) | 一种基于桡动脉脉搏波传导时间的测量血压方法及装置 | |
| KR101210828B1 (ko) | 다중 생체 신호 측정을 이용하여 손목혈압의 정확도를 향상시키는 방법 및 장치 | |
| US20100125217A1 (en) | Method and Apparatus for Presenting Heart Rate Variability by Sound and/or Light | |
| Tadi et al. | Seismocardiography: Toward heart rate variability (HRV) estimation | |
| Lu et al. | Limitations of oximetry to measure heart rate variability measures | |
| US20090124914A1 (en) | Analysis system and a method for pulse diagnosis in chinese medicine | |
| JP2008295517A (ja) | 漢方医における脈診の分析システムと方法 | |
| Wang et al. | An advanced signal processing technique for impedance cardiography | |
| Singh et al. | Effects of RR segment duration on HRV spectrum estimation | |
| CN113160921A (zh) | 一种基于血流动力学的数字人体心血管***的构建方法和应用 | |
| CN110881967A (zh) | 一种无创多节段外周动脉血管弹性功能检测方法及其仪器 | |
| JP2008513073A (ja) | 心調律信号のシリーズ（ｒｒ）を処理するための方法、及び心調律の変動性を分析するための、特に生物の痛み又はストレスを評価するためのその使用 | |
| EP1982646A2 (en) | Analysis system and a method for pulse diagnostic in chinese medicine | |
| CN111839488A (zh) | 基于脉搏波的无创连续血压测量装置和方法 | |
| US20040181159A1 (en) | Method and apparatus for detecting yin-yang and asthenia-sthenia | |
| KR102171883B1 (ko) | 멀티스케일 진폭 인식 치환 엔트로피를 사용한 심박변이율 분석 방법 | |
| Ahmad et al. | A prototype of an integrated blood pressure and electrocardiogram device for multi-parameter physiologic monitoring | |
| Chou et al. | Comparison between heart rate variability and pulse rate variability for bradycardia and tachycardia subjects | |
| Firoozabadi et al. | Efficient noise-tolerant estimation of heart rate variability using single-channel photoplethysmography | |
| Amanipour et al. | The effects of blood glucose changes on frequency-domain measures of HRV signal in type 1 diabetes | |
| Sakajiri et al. | Non-contact capacitive ballistocardiogram measurements using in-bed fabric-sheet electrode for blood pressure estimation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Open date: 20081126 |