脉搏波数据处理方法及相关设备与流程

1.本发明涉及信号处理领域，尤其涉及一种脉搏波数据处理方法及相关设备。


背景技术：

2.脉搏波频域特征分析又称为脉搏波谐频特征分析，是血液动力学研究分析重要的一环，所得到的脉搏波谐频的振幅和相位是人体血液循环系统中重要的生理参数，反映了动脉系统状态以及传递出心血管疾病病理变化的重要信息。因此，脉搏波谐频的振幅和相位特征也被用于中药对循环系统影响的一个量化测试方式，应用于药物动力学中，被称之为脉象药性动力学研究。所以，准确的获取稳定脉搏波频域特征有非常重要的意义。中医的经典《黄帝内经》中，素问三部九候论篇是中国医家先贤从临床脉象分析的总结，说明了不同动脉间的脉搏波波型特征的变化是人体状态的重要脉象表征。然而，一个中草药的完整特性包含了许多复合挥发性成份、亲水性成分、不可溶成份等等，因此，由于短煮久煮、大火小火，炮制、精油提纯等得到部分复合成份药性皆不相同，中草药的特性无法进行量化分析。


技术实现要素：

3.本发明实施例旨在提供一种脉搏波数据处理方法，能够在获取指定目标在目标药态下的脉搏波信号后，量化目标药态下脉搏波频域特征在动脉系统上的变量分布，旨在解决现有技术提供的脉搏波频域特征由于一个中草药的完整特性包含了许多复合挥发性成份、亲水性成分、不可溶成份等等，导致短煮久煮、大火小火，炮制、精油提纯等得到部分复合成份药性皆不相同，中草药的完整特性无法进行量化分析的问题，从而将中草药在不同药态下的脉搏波转换成可量测、可再现、直观的量化分析数据。
4.第一方面，本发明实施例提供一种脉搏波数据处理方法，所述脉搏波数据处理方法包括：
5.获取指定目标在目标药态下的脉搏波信号；
6.对所述脉搏波信号进行谐频特征提取，得到所述脉搏波信号对应的谐频特征数据；
7.计算所述谐频特征数据的变量，得到所述脉搏波信号的变量分布数据。
8.可选的，所述目标药态包括纯露、汤液、固态物对应的三项药态，获取指定目标在目标药态下的脉搏波信号，包括：
9.获取所述指定目标分别处于纯露药态、汤液药态、固态物药态下对应的脉搏波信号。
10.可选的，对所述脉搏波信号进行谐频特征提取，得到所述脉搏波信号对应的谐频特征数据，包括：
11.将所述脉搏波信号按信号周期进行切分，得到所述脉搏波信号的切分序列，每个切分序列包括与所述脉搏波信号周期数相同的切分信号；
12.对所述切分序列中的各个切分信号进行复利叶级数转换处理，得到所述脉搏波信号对应的谐频特征数据。
13.可选的，所述对所述切分序列中的各个所述切分信号进行复利叶级数转换处理，得到所述脉搏波信号对应的谐频特征数据，包括：
14.对所述切分序列中的各个所述切分信号进行复利叶级数转换处理，得到所述切分信号的谐频数据；
15.对所述切分信号的谐频数据进行正规化，得到所述脉搏波信号对应的谐频特征数据。
16.可选的，所述对所述切分序列中的各个所述切分信号进行复利叶级数转换处理，得到所述切分信号的谐频数据，包括：
17.计算所述脉搏波信号的平均振幅和平均相位；
18.根据所述脉搏波信号的平均振幅和平均相位，对所述切分信号的谐频数据进行正规化，得到所述脉搏波信号对应的谐频特征数据。
19.可选的，所述计算所述谐频特征数据的变量，得到所述脉搏波信号的变量分布数据，包括：
20.计算所述谐频特征数据的百分比变量，得到所述脉搏波信号的百分比变量分布数据。
21.可选的，所述计算所述谐频特征数据的变量，得到所述脉搏波信号的变量分布数据，包括：
22.计算所述谐频特征数据的相位差变量，得到所述脉搏波信号的相位差变量分布数据。
23.第二方面，本发明实施例提供一种脉搏波数据处理，所述装置包括：
24.获取模块，用于获取指定目标在目标药态下的脉搏波信号；
25.提取模块，用于对所述脉搏波信号进行谐频特征提取，得到所述脉搏波信号对应的谐频特征数据；
26.计算模块，用于计算所述谐频特征数据的变量，得到所述脉搏波信号的变量分布数据。
27.第三方面，本发明实施例提供一种电子设备，包括：存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现本发明实施例提供的脉搏波数据处理方法中的步骤。
28.第四方面，本发明实施例提供一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时实现发明实施例提供的脉搏波数据处理方法中的步骤。
29.本发明实施例中，获取指定目标在目标药态下的脉搏波信号；对所述脉搏波信号进行谐频特征提取，得到所述脉搏波信号对应的谐频特征数据；计算所述谐频特征数据的变量，得到所述脉搏波信号的变量分布数据。能够在获取在指定目标在目标药态下的脉搏波信号后，量化目标药态下脉搏波频域特征在动脉系统上的变量分布，旨在解决现有技术提供的脉搏波频域特征由于一个中草药的完整特性包含了许多复合挥发性成份、亲水性成分、不可溶成份等等，导致短煮久煮、大火小火，炮制、精油提纯等得到部分复合成份药性皆
不相同，中草药的完整特性无法进行量化分析的问题，从而将中草药在不同药态下的脉搏波转换成可量测、可再现、直观的量化分析数据。
附图说明
30.为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
31.图1是本发明实施例提供的一种脉搏波数据处理系统的架构图；
32.图2是本发明实施例提供的一种药材蒸馏装置的结构示意图；
33.图3是本发明实施例提供的一种脉搏波数据处理方法的流程示意图；
34.图4是本发明实施例提供的一种服食中草药前后的脉搏波示意图；
35.图5是本发明实施例提供一种谐频振幅平均值变化示意图；
36.图6是本发明实施例提供一种谐频振幅变异系数变化示意图；
37.图7是本发明实施例提供一种谐频相位变化示意图；
38.图8是本发明实施例提供一种标准相位偏差变化示意图；
39.图9是本发明实施例提供的一种脉搏波数据处理装置的结构示意图；
40.图10是本发明实施例提供的一种电子设备的结构示意图。
具体实施方式
41.下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
42.请参见图1，图1是本发明实施例提供的一种脉搏波数据处理系统的架构图，如图1所示，该脉搏波数据处理系统包括：药材蒸馏装置、压力信号采集模块、脉搏谐频分析单元以及谐频特征变异量分析单元。
43.其中，药材蒸馏装置包括恒温定时加热锅炉、风扇、冷凝管、纯露收集瓶；压力信号采集模块包括压力传感器和直流、交流脉搏信号获取电路；压力信号采集模块与脉搏谐频分析单元信号连接，压力信号采集模块用于将采集到的脉搏信号传输到脉搏波谐频分析单元中；脉搏波谐频分析单元与谐频特征变异量分析单元信号连接。
44.具体的，请参考图2，图2是本发明实施例提供的一种药材蒸馏装置的结构示意图，如图2所示，冷凝管设置在恒温定时加热锅炉的上方，风扇设置在冷凝管的上方，冷凝管的出口设置有集液嘴，集液嘴的下方设置有纯露收集瓶。更具体的，可以通过下述步骤进行中草药蒸馏：
45.d1.准备原材料：将目标中草药进行清洗、适当切片、切碎、压扁，便于烹煮提纯。其中，目标中草药为新鲜中草药。
46.d2.将目标中草药投入恒温定时加热锅炉中，以水或液态溶剂充分浸泡盖过目标中草药，将恒温定时加热锅炉加热至特定恒温进行蒸馏，通过恒温定时加热锅炉上方的风
扇、冷凝管和纯露收集瓶收集冷凝后的液态纯露，上述纯露可以包含精油、水、低沸点挥发物。蒸馏至固定时间后开盖冷却，再将锅炉内蒸馏后的汤液和固态物通过过滤的方式分开。
47.在一种可能的实施例中，上述步骤d2中的恒温加热的温度可以维持在80-120摄氏度，上述步骤d2中的蒸馏时间可以为30-420分钟，上述步骤d2中的液态溶剂可以为亲水性溶剂，上述步骤d2中的液态溶剂体积可以为药材1倍至10倍，上述步骤d2中可以使用低压低温蒸馏法，其中减压蒸馏温度在35-80摄氏度。
48.在本发明实施例中，获取在指定目标在目标药态下的脉搏波信号，指定目标可以为哺乳动物，比如白鼠、白兔等哺乳动物，其中，上述的目标药态包括服食纯露前后、服食汤液前后、服食固态物前后等三种药态。通过对中草药进行蒸馏，得到中草药的纯露、汤液和固态物，记录哺乳动物服食前、服食纯露、汤液、固态物后的脉搏波信号，将采集到脉搏波信号作为数据样本进行分析。
49.在本发明实施例中，可以选用三组哺乳动物，第一组哺乳动物记录哺乳动物服食前和服食纯露后的脉搏波信号，第二组哺乳动物记录哺乳动物服食前和服食汤液后的脉搏波信号，第三组哺乳动物记录哺乳动物服食前和服食固态物后的脉搏波信号。
50.具体的，可以通过压力信号采集模块采集哺乳动物服食前、服食纯露、汤液、固态物后的脉搏波信号，通过脉搏谐频分析单元对采集到的脉搏波信号进行谐频处理，通过谐频特征变异量分析单元对经过谐频处理后的脉搏波信号进行变量分析。
51.上述信号连接可以是有线连接方式或者无线连接方式，上述无线连接方式可以包括但不限于3g/4g/5g连接、wifi连接、蓝牙连接、wimax连接、zigbee连接、uwb(ultra wideband)连接、以及其他现在已知或将来开发的无线连接方式。
52.请参见图3，图3是本发明实施例提供的一种脉搏波数据处理方法的流程图，如图3所示，该脉搏波数据处理方法具体包括以下步骤：
53.301、获取在指定目标在目标药态下的脉搏波信号。
54.在本发明实施例中，获取在指定目标在目标药态下的脉搏波信号，指定目标可以为哺乳动物，比如白鼠、白兔等哺乳动物，其中，上述的目标药态包括服食纯露前后、服食汤液前后、服食固态物前后等三种药态。通过对中草药进行蒸馏，得到中草药的纯露、汤液和固态物，记录哺乳动物服食前、服食纯露、汤液、固态物后的脉搏波信号，将采集到脉搏波信号作为数据样本进行分析。
55.在本发明实施例中，可以选用三组哺乳动物，第一组哺乳动物记录哺乳动物服食前和服食纯露后的脉搏波信号，第二组哺乳动物记录哺乳动物服食前和服食汤液后的脉搏波信号，第三组哺乳动物记录哺乳动物服食前和服食固态物后的脉搏波信号。
56.上述脉搏波信号可以是压力信号采集模块采集到的交流压力信号后，经过adc采集、数字滤波去除高频噪声后得到的脉搏波信号。
57.可选的，在获取在指定目标在目标药态下的脉搏波信号的步骤中，可以保持任意两个脉搏波信号的平均静压压力读值差在小于预设值的条件下，对多个指定人体部位进行脉搏波信号的采集，以得到多个指定人体部位的脉搏波信号。
58.具体来说，可以在脉搏波信号的采集过程中，通过定压控制模块将压力信号采集模块中的袖带气囊维持10-60秒恒定的压力读值，压力值范围在60-100毫米汞柱区间。任意两个袖戴气囊的平均静压压力读值差《3毫米汞柱。在将采集到的一系列单周期信号从频域
方面进行稳定特征提取时，当任意单个脉搏波脉率》平均脉率110％或单个脉搏波脉率《平均脉率90％时该单个脉搏波将不列入后续计算。
59.可选的，目标药态包括纯露、汤液、固态物对应的三项药态，具体的，在获取在指定目标在目标药态下的脉搏波信号的步骤中，上述的目标药态可以包括服食纯露前后、服食汤液前后、服食固态物前后等三种药态。上述脉搏波信号包括服食纯露前后、服食汤液前后、服食固态物前后等六个类型的脉搏波信号。
60.302、对脉搏波信号进行谐频特征提取，得到脉搏波信号对应的谐频特征数据。
61.在本发明实施例中，需要说明的是，脉搏波频域特征又称为脉搏波谐频特征，初始的脉搏波信号是属于时域信号，可以对脉搏波信号进行谐频特征提取，提取得到属于频域信号的谐频特征数据。
62.可选的，可以将脉搏波信号按信号周期进行切分，得到脉搏波信号的切分序列，每个切分序列包括与脉搏波信号周期数相同的切分信号；对切分序列中的各个切分信号进行复利叶级数转换处理，得到脉搏波信号对应的谐频特征数据。
63.具体的，可以对连续的脉搏波按最低点位切分成逐个脉搏波，循序对每一个部位的脉搏波以复利叶级数转换得到谐频特征数据cn。
64.可选的，在对切分序列中的各个切分信号进行复利叶级数转换处理，得到脉搏波信号对应的谐频特征数据的步骤中，还可以对切分序列中的各个切分信号进行复利叶级数转换处理，得到切分信号的谐频数据；对切分信号的谐频数据进行正规化，得到脉搏波信号对应的谐频特征数据。
65.具体的，可以对连续的脉搏波按最低点位切分成逐个脉搏波，循序对每一个部位的脉搏波以复利叶级数转换与正规化手段计算谐频特征数据cn。
66.可选的，在对切分序列中的各个切分信号进行复利叶级数转换处理，得到切分信号的谐频数据的步骤中，可以计算脉搏波信号的平均振幅和平均相位；根据脉搏波信号的平均振幅和平均相位，对切分信号的谐频数据进行正规化，得到脉搏波信号对应的谐频特征数据。
67.在本发明实施例中，谐频数据包括可以振幅数据和相位数据，可以对切分信号的振幅数据进行正规化，得到脉搏波信号对应于振幅数据的谐频特征数据，具体的，可以计算脉搏波信号的平均振幅；根据脉搏波信号的平均振幅，对切分信号的振幅数据进行正规化，得到脉搏波信号对应于振幅数据的谐频特征数据。
68.进一步的，上述对应于振幅数据的谐频特征数据还包括经过正规化的振幅标准偏差cnsd，正规化谐频振幅变异系数cncv。
69.可以对切分信号的相位数据进行正规化，得到脉搏波信号对应于相位数据的谐频特征数据，具体的，可以计算脉搏波信号的平均相位；根据脉搏波信号的平均相位，对切分信号的相位数据进行正规化，得到脉搏波信号对应于相位数据的谐频特征数据，即谐频相位平均值。
70.进一步的，上述对应于相位数据的谐频特征数据还包括经过正规化相位的标准偏差pnsd。
71.为进一步进行说明，本发明实施例以12个周期以上的连续脉搏波信号为例进行说明，对连续脉搏波信号以每个脉搏波波的最低点为划分点进行周期划分；每个脉搏波波形
含有各个谐频的振幅an和相位θn特征，经过复利叶级数的特征计算可得出每个脉搏波的振幅an和相位θn；其中，每个脉搏波波形的通式表达示如下述式1所示：
[0072][0073]
其中y(k)为单个脉搏波信号离散数字化的值，an为第n个复利叶级数的振幅，θn为第n个复利叶级数的相位，n代表第n个谐频。n代表数据点总数。f0为脉搏波信号的主频率(基频)，上述复利叶级数的计算表达式可以如下述式2所示：
[0074][0075]
在振幅an的基础上除以脉搏波波形的平均值a0，进而得到正规化的谐频特征数据cn(也可以称为谐频振幅)，在得到每个脉搏波的振幅an后，可以通过下述式3和式4计算一段连续脉搏波经过正规化的谐频振幅平均值cnavg，式3和式4如下所示：
[0076][0077][0078]
其中，m为一段连续时间内脉搏波的总共个数，an，q为第q个脉搏波计算后第n个复利叶级数的振幅，a0，q为第q个脉搏波的平均值。cnavg为一段连续时间内第n个谐频的经过正规化的谐频振幅平均值。
[0079]
在得到一段连续脉搏波的经过正规化的谐频振幅平均值cnavg后，可以通过下述式5和式6计算一段连续时间内第n个谐频的经过正规化的标准振幅偏差cnsd和一段连续时间内第n个谐频的经过正规化的谐频振幅变异系数cncv，式5和式6如下所示：
[0080][0081][0082]
其中，cncv为一段连续时间内第n个谐频经过正规化振幅的变异系数。
[0083]
在得到每个脉搏波的相位θn后，可以通过下述式7和式8计算一段连续脉搏波经过正规化的谐频相位平均值pn和经过正规化的标准相位偏差pnsd，式7和式8如下所示：
[0084]
[0085][0086]
其中，m为一段连续时间内脉搏波的总共个数，θn，q为第q个脉搏波计算后第n个复利叶级数的相位，n代表第n个谐频。pn为一段连续时间内第n个谐频的相位平均值。pnsd为一段连续时间内第n个谐频经过正规化的标准相位偏差。
[0087]
303、计算谐频特征数据的变量，得到脉搏波信号的变量分布数据。
[0088]
在本发明实施例中，上述谐频特征数据的变量可以包括谐频特征数据的百分比变量和谐频特征数据的相位差变量。上述谐频特征数据的变量可以是服食中草药前后的变量。
[0089]
通过谐频特征数据的变量，得到在服食中草药前后的脉搏波信号的变量分布数据，进而可以得到中草药在不同药态下对于脉搏波信号的影响。具体的，通过谐频特征数据的百分比变量和谐频特征数据的相位差变量，可以得到中草药在不同药态下对于脉搏波信号的影响。
[0090]
可选的，在计算谐频特征数据的变量，得到脉搏波信号的变量分布数据的步骤中，可以计算谐频特征数据的百分比变量，得到脉搏波信号的百分比变量分布数据。
[0091]
具体的，谐频特征数据的百分比变量包括目标药态下的谐频振幅平均值cnavg百分比变量、标准相位偏差pnsd百分比变量和谐频振幅变异系数cncv百分比变量。
[0092]
更具体的，谐频特征数据的百分比变量包括服食纯露前后的谐频振幅平均值cnavg百分比变量、标准相位偏差pnsd百分比变量和谐频振幅变异系数cncv百分比变量；服食汤液前后的谐频振幅平均值cnavg百分比变量、标准相位偏差pnsd百分比变量和谐频振幅变异系数cncv百分比变量；服食固态物前后的谐频振幅平均值cnavg百分比变量、标准相位偏差pnsd百分比变量和谐频振幅变异系数cncv百分比变量。
[0093]
具体的，上述谐频特征数据的百分比变量可以通过下述式9至式11进行计算：
[0094][0095][0096][0097]
上述cnavg0代表服食中草药前脉搏波信号的谐频振幅平均值、cncv0表服食中草药前脉搏波信号的谐频振幅变异系数、pnsd0代表服食中草药前脉搏波信号的标准相位偏差；cnavg1代表服食中草药后脉搏波信号的谐频振幅平均值、cncv1代表服食中草药后脉搏波信号的谐频振幅变异系数、pnsd1代表服食中草药后脉搏波信号的标准相位偏差。
△
cnavg(％)表示在目标药态下的谐频振幅平均值cnavg百分比变量，
△
cncv(％)表示在目标药态下的谐频振幅变异系数cncv百分比变量，
△
pnsd(％)表示在目标药态下的标准相位偏差
pnsd百分比变量。
[0098]
可选的，在计算谐频特征数据的变量，得到脉搏波信号的变量分布数据的步骤中，可以计算谐频特征数据的相位差变量，得到脉搏波信号的相位差变量分布数据。
[0099]
具体的，谐频特征数据的百分比变量包括目标药态下的谐频相位平均值pn变量。上述相位差变量分布数据为谐频相位平均值pn变量。
[0100]
更具体的，谐频特征数据的百分比变量包括服食纯露前后的谐频相位平均值pn变量；服食汤液前后的谐频相位平均值pn变量；服食固态物前后的谐频相位平均值pn变量。
[0101]
上述谐频相位平均值pn变量可以通过下述式12进行计算：
[0102]
δpn(％)＝pn
1-pn0ꢀꢀꢀ
(12)
[0103]
上述pn0代表服食中草药前脉搏波信号的谐频相位平均值，pn1代表服食中草药后脉搏波信号的谐频相位平均值。
△
pn(％)表示在目标药态下的谐频相位平均值的相位偏差值。可以通过相位偏差值来呈现服食中草药前后脉象谐频特征的变化。
[0104]
本发明实施例中，获取在指定目标在目标药态下的脉搏波信号；对所述脉搏波信号进行谐频特征提取，得到所述脉搏波信号对应的谐频特征数据；计算所述谐频特征数据的变量，得到所述脉搏波信号的变量分布数据。能够在获取在指定目标在目标药态下的脉搏波信号后，量化目标药态下脉搏波频域特征在动脉系统上的变量分布，旨在解决现有技术提供的脉搏波频域特征由于一个中草药的完整特性包含了许多复合挥发性成份、亲水性成分、不可溶成份等等，导致短煮久煮、大火小火，炮制、精油提纯等得到部分复合成份药性皆不相同，中草药的完整特性无法进行量化分析的问题，从而将中草药在不同药态下的脉搏波转换成可量测、可再现、直观的量化分析数据。
[0105]
本发明实施例能够透过脉搏波信号的变量分布数据反应服食中医药造成的循环系统状态的改变，应用于血流动力学和中草药药性动力学的鉴定与分析相关的研究，或是作为中医三部九候量化数据的基石。能够看脉象或是血流动力学状态上更为细致、多元的变化特征，对于中医脉诊的科学量化研究有很大的正向促进作用。能够直观的查看中草药的药性影响数据，在血流动力学状态和中草药药性动力学鉴定与分析上更具有标志性的特征，能够很完整的将一味中草药对脉象的影响，以及其来源成分的主要物化特性掌握。同时，亦能观察这些分离的复合成份，随时间代谢的特性，对中草药药性动力学的研究，是非常重要的数据量化的基石。
[0106]
请参考图4，图4是本发明实施例提供的一种服食中草药前后的脉搏波示意图，其中，(a)部分是服食中草药前的脉搏波，(b)部分是服食中草药后的脉搏波，可以看出，服食中草药前的脉搏波与服食中草药后的脉搏波存在波形上的差异，实施本发明实施例，可以将中草药在不同药态下的脉搏波转换成可量测、可再现、直观的量化分析数据。
[0107]
具体的，以15年生长白山人参脉象药性定量全图谱分析为例，通过本发明实施例提供的中草药蒸馏方法分离得到为人参纯露、人参汤液、人参固态物，采集健康人或健康哺乳动物在服食中草药前的脉搏波信号，在健康人或健康哺乳动物分别服食人参纯露、人参汤液、人参固态物30分钟后，采集服食中草药后的脉搏波信号。
[0108]
首先，可以根据本发明实施例提供的脉搏波数据处理方法得到谐频振幅及变化，如图5所示，图5是本发明实施例提供的一种谐频振幅平均值变化示意图。如图5所示，可以看出，受测者分别服食15年生长白山人参“纯露、汤液、固态物”30分钟后，脉搏波谐频振幅
的改变。n代表复利叶级数的第n谐频。具体的，可以看出，饮用15年生长白山人参分离出的汤液对谐频振幅影响最大，其中又以第三谐频的谐频振幅增加最多，从而可以得出“人参汤液”对血液动力影响较大的结论。相对而言，“纯露、固态物”对谐频振幅特征影响较小。
[0109]
其次，可以根据本发明实施例提供的脉搏波数据处理方法得到谐频振幅平均值变化，如图6所示，图6是本发明实施例提供的一种谐频振幅变异系数变化示意图。如图6所示，可以看出，受测者分别服食15年生长白山人参“纯露、汤液、固态物”30分钟后，脉搏波谐频振幅变异系数的改变。n代表复利叶级数的第n谐频。具体的，可以看出，对脉搏波谐频振幅变异系数cncv的改变，仍以长白山人参“汤液”影响最为直接，而“纯露、固态物”的影响较小。
[0110]
然后，可以根据可以根据本发明实施例提供的脉搏波数据处理方法得到谐频相位变化如图7所示，图7是本发明实施例提供的一种谐频相位变化示意图。如图7所示，可以看出，受测者分别服食15年生长白山人参“纯露、汤液、固态物”30分钟后，脉搏波谐频振幅变异系数cncv的改变。n代表复利叶级数的第n谐频。具体的，可以看出，对脉搏波谐频相位平均值pn的改变，仍以长白山人参“汤液、纯露”影响特征的方向相同，“汤液”影响特征的变化强度较高，而“固态物”影响pn的方向和前两者相反，呈现拮抗作用，且对9-11高频谐频影响较大。
[0111]
最后，可以根据可以根据本发明实施例提供的脉搏波数据处理方法得到标准相位偏差如图8所示，图8是本发明实施例提供的一种标准相位偏差变化示意图。如图8所示，可以看出，受测者分别服食15年生长白山人参“纯露、汤液、固态物”30分钟后，脉搏波标准相位偏差的改变。n代表复利叶级数的第n谐频。具体的，可以看出，对脉搏波谐频相位标准偏差pnsd的改变，仍以长白山人参"汤液"影响最大，增幅量在第3谐频、第4谐频皆超过300％。表示"人参汤液"对脾经、肺经相位标准偏差(血分)影响较大，而人参纯露、固态物的功效远不如汤液。
[0112]
综上，实施例本发明实施例，可以得到谐频振幅平均值cnavg(谐频振幅平均值也可以称为谐频振幅)、谐频振幅变异系数cncv、谐频相位平均值pn(谐频相位平均值也可以称为谐频相位)、相位标准偏差pnsd完整的量化过程，可以通过量化后的数据更全面的理解15年长白山人参的挥发成份、水溶成份、固态成分分别对循环系统影响的每个谐频特征面相，说明了人参对脉象改变中的主要影响成份为水溶性汤液、其次是挥发的纯露，而固态物的循环系统的影响较小。
[0113]
需要说明的是，本发明实施例提供的脉搏波数据处理方法可以应用于可以进行脉搏波数据处理的智能手机、电脑、服务器等设备。
[0114]
可选的，请参见图9，图9是本发明实施例提供的一种脉搏波数据处理装置的结构示意图，如图9所示，所述装置包括：
[0115]
获取模块901，用于获取在指定目标在目标药态下的脉搏波信号；
[0116]
提取模块902，用于对所述脉搏波信号进行谐频特征提取，得到所述脉搏波信号对应的谐频特征数据；
[0117]
计算模块903，用于计算谐频特征数据的变量，得到脉搏波信号的变量分布数据。
[0118]
可选的，所述目标药态包括纯露、汤液、固态物对应的三项药态，获取模块901，包括：
[0119]
获取子模块，用于获取所述指定目标分别处于纯露药态、汤液药态、固态物药态下对应的脉搏波信号。
[0120]
可选的，提取模块902，包括：
[0121]
切分子模块，用于将所述脉搏波信号按信号周期进行切分，得到所述脉搏波信号的切分序列，每个切分序列包括与所述脉搏波信号周期数相同的切分信号；
[0122]
处理子模块，用于对所述切分序列中的各个切分信号进行复利叶级数转换处理，得到所述脉搏波信号对应的谐频特征数据。
[0123]
可选的，所述处理子模块，包括：
[0124]
第一处理单元，用于对所述切分序列中的各个所述切分信号进行复利叶级数转换处理，得到所述切分信号的谐频数据；
[0125]
第二处理单元，用于对所述切分信号的谐频数据进行正规化，得到所述脉搏波信号对应的谐频特征数据。
[0126]
可选的，所述第一处理单元，包括：
[0127]
计算子单元，用于计算所述脉搏波信号的平均振幅和平均相位；
[0128]
处理子单元，用于根据所述脉搏波信号的平均振幅和平均相位，对所述切分信号的谐频数据进行正规化，得到所述脉搏波信号对应的谐频特征数据。
[0129]
可选的，所述计算模块903，包括：
[0130]
第一计算子模块，用于计算所述谐频特征数据的百分比变量，得到所述脉搏波信号的百分比变量分布数据。
[0131]
可选的，所述计算模块903，包括：
[0132]
第二计算子模块，用于计算所述谐频特征数据的相位差变量，得到所述脉搏波信号的相位差变量分布数据。
[0133]
需要说明的是，本发明实施例提供的数据中心设备可以应用于可以进行脉搏波数据处理的智能手机、电脑、服务器等设备。
[0134]
本发明实施例提供的数据中心设备能够实现上述方法实施例中脉搏波数据处理方法实现的各个过程，且可以达到相同的有益效果。为避免重复，这里不再赘述。
[0135]
参见图10，图10是本发明实施例提供的一种电子设备的结构示意图，如图10所示，包括：存储器1002、处理器1001及存储在所述存储器1002上并可在所述处理器1001上运行的脉搏波数据处理方法的计算机程序，其中：
[0136]
电子设备应用于数据中心，处理器1001用于调用存储器1002存储的计算机程序，执行如下步骤：
[0137]
获取指定目标在目标药态下的脉搏波信号；
[0138]
对所述脉搏波信号进行谐频特征提取，得到所述脉搏波信号对应的谐频特征数据；
[0139]
计算所述谐频特征数据的变量，得到所述脉搏波信号的变量分布数据。
[0140]
可选的，处理器1001执行的所述目标药态包括纯露、汤液、固态物对应的三项药态，获取指定目标在目标药态下的脉搏波信号，包括：
[0141]
获取所述指定目标分别处于纯露药态、汤液药态、固态物药态下对应的脉搏波信号。
[0142]
可选的，处理器1001执行的对所述脉搏波信号进行谐频特征提取，得到所述脉搏波信号对应的谐频特征数据，包括：
[0143]
将所述脉搏波信号按信号周期进行切分，得到所述脉搏波信号的切分序列，每个切分序列包括与所述脉搏波信号周期数相同的切分信号；
[0144]
对所述切分序列中的各个切分信号进行复利叶级数转换处理，得到所述脉搏波信号对应的谐频特征数据。
[0145]
可选的，处理器1001执行的所述对所述切分序列中的各个所述切分信号进行复利叶级数转换处理，得到所述脉搏波信号对应的谐频特征数据，包括：
[0146]
对所述切分序列中的各个所述切分信号进行复利叶级数转换处理，得到所述切分信号的谐频数据；
[0147]
对所述切分信号的谐频数据进行正规化，得到所述脉搏波信号对应的谐频特征数据。
[0148]
可选的，处理器1001执行的所述对所述切分序列中的各个所述切分信号进行复利叶级数转换处理，得到所述切分信号的谐频数据，包括：
[0149]
计算所述脉搏波信号的平均振幅和平均相位；
[0150]
根据所述脉搏波信号的平均振幅和平均相位，对所述切分信号的谐频数据进行正规化，得到所述脉搏波信号对应的谐频特征数据。
[0151]
可选的，处理器1001执行的所述计算所述谐频特征数据的变量，得到所述脉搏波信号的变量分布数据，包括：
[0152]
计算所述谐频特征数据的百分比变量，得到所述脉搏波信号的百分比变量分布数据。
[0153]
可选的，处理器1001执行的所述计算所述谐频特征数据的变量，得到所述脉搏波信号的变量分布数据，包括：
[0154]
计算所述谐频特征数据的相位差变量，得到所述脉搏波信号的相位差变量分布数据。
[0155]
需要说明的是，本发明实施例提供的电子设备可以应用于可以进行脉搏波数据处理的智能手机、电脑、服务器等设备。
[0156]
本发明实施例提供的电子设备能够实现上述方法实施例中脉搏波数据处理方法实现的各个过程，且可以达到相同的有益效果。为避免重复，这里不再赘述。
[0157]
本发明实施例还提供一种计算机可读存储介质，计算机可读存储介质上存储有计算机程序，该计算机程序被处理器执行时实现本发明实施例提供的脉搏波数据处理方法或应用端脉搏波数据处理方法的各个过程，且能达到相同的技术效果，为避免重复，这里不再赘述。
[0158]
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于一计算机可读取存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，所述的存储介质可为磁碟、光盘、只读存储记忆体(read-only memory，rom)或随机存取存储器(random access memory，简称ram)等。
[0159]
以上所揭露的仅为本发明较佳实施例而已，当然不能以此来限定本发明之权利范
围，因此依本发明权利要求所作的等同变化，仍属本发明所涵盖的范围。