CN106096303B - 一种自主神经系统的分析方法、服务器及系统 - Google Patents

一种自主神经系统的分析方法、服务器及系统 Download PDF

Info

Publication number
CN106096303B
CN106096303B CN201610460442.3A CN201610460442A CN106096303B CN 106096303 B CN106096303 B CN 106096303B CN 201610460442 A CN201610460442 A CN 201610460442A CN 106096303 B CN106096303 B CN 106096303B
Authority
CN
China
Prior art keywords
measured
user
electrocardiosignal
peak value
heart rate
Prior art date
Application number
CN201610460442.3A
Other languages
English (en)
Other versions
CN106096303A (zh
Inventor
杨荣骞
张磊
吕瑞雪
付小婷
陈秀文
王之辉
郭加成
宋传旭
Original Assignee
深圳市是源医学科技有限公司
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 深圳市是源医学科技有限公司 filed Critical 深圳市是源医学科技有限公司
Priority to CN201610460442.3A priority Critical patent/CN106096303B/zh
Publication of CN106096303A publication Critical patent/CN106096303A/zh
Application granted granted Critical
Publication of CN106096303B publication Critical patent/CN106096303B/zh

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F19/00Digital computing or data processing equipment or methods, specially adapted for specific applications
    • G06F19/30Medical informatics, i.e. computer-based analysis or dissemination of patient or disease data
    • G06F19/34Computer-assisted medical diagnosis or treatment, e.g. computerised prescription or delivery of medication or diets, computerised local control of medical devices, medical expert systems or telemedicine
    • G06F19/3418Telemedicine, e.g. remote diagnosis, remote control of instruments or remote monitoring of patient carried devices
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/20ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L67/00Network-specific arrangements or communication protocols supporting networked applications
    • H04L67/12Network-specific arrangements or communication protocols supporting networked applications adapted for proprietary or special purpose networking environments, e.g. medical networks, sensor networks, networks in a car or remote metering networks

Abstract

本发明属于信号处理技术领域,提供了一种自主神经系统的分析方法、服务器及系统,所述方法包括:接收移动终端发送的预处理后的待测用户的心电信号;根据所述预处理后的心电信号,获得所述待测用户的心率信号和呼吸信号;对所述心率信号进行频谱分析,以获得所述心率信号的功率谱曲线图;以及对所述呼吸信号进行频谱分析,以获得所述呼吸信号的功率谱曲线图;根据所述心率信号的功率谱曲线图和所述呼吸信号的功率谱曲线图,分析所述待测用户的自主神经系统,并将分析结果发送给所述移动终端,以使得所述移动终端显示所述分析结果。通过本发明可提高分析自主神经系统的准确性。

Description

一种自主神经系统的分析方法、服务器及系统

技术领域

[0001] 本发明属于信号处理技术领域,尤其涉及一种自主神经系统的分析方法、服务器 及系统。

背景技术

[0002] 自主神经系统又称为植物神经系统,包括交感神经和副交感神经,可以支配和调 节机体内脏、血管、平滑肌、心肌和腺体等的活动。由于慢性疾病伴随着自主神经紊乱的发 生,可以通过维持自主神经的平衡状态有效降低患病率和死亡率。因此,实时检测自主神经 的状态在预防疾病领域具有重要的意义,而心电信号的频域分析在分析自主神经系统时尤 为重要。

[0003] 目前,通常是对心电信号的RR间期进行傅里叶变换功率谱分析,其中,具有固定带 宽0.04〜0.15Hz的低频功率LF代表交感神经,具有固定带宽0.15〜0.4Hz的高频功率HF代 表副交感神经,通过LF、HF和LF/HF这三个值来分析自主神经系统。然而,现有技术并没有考 虑呼吸信号对自主神经系统的影响,进而导致分析自主神经系统的准确性不高。

[0004] 故,有必要提出一种新的技术方案,以解决上述技术问题。

发明内容

[0005] 鉴于此,本发明实施例提供一种自主神经系统的分析方法、服务器及系统,以提高 分析自主神经系统的准确性。

[0006] 本发明实施例的第一方面,提供一种自主神经系统的分析方法,所述方法包括:

[0007] 接收移动终端发送的预处理后的待测用户的心电信号;

[0008] 根据所述预处理后的心电信号,获得所述待测用户的心率信号和呼吸信号;

[0009] 对所述心率信号进行频谱分析,以获得所述心率信号的功率谱曲线图;以及对所 述呼吸信号进行频谱分析,以获得所述呼吸信号的功率谱曲线图;

[0010] 根据所述心率信号的功率谱曲线图和所述呼吸信号的功率谱曲线图,分析所述待 测用户的自主神经系统,并将分析结果发送给所述移动终端,以使得所述移动终端显示所 述分析结果。

[0011] 本发明实施例的第二方面,提供一种服务器,所述服务器包括:

[0012] 接收模块,用于接收移动终端发送的预处理后的待测用户的心电信号;

[0013] 获得模块,用于根据所述预处理后的心电信号,获得所述待测用户的心率信号和 呼吸信号;

[0014] 第一分析模块,用于对所述心率信号进行频谱分析,以获得所述心率信号的功率 谱曲线图;以及对所述呼吸信号进行频谱分析,以获得所述呼吸信号的功率谱曲线图;

[0015] 第二分析模块,用于根据所述心率信号的功率谱曲线图和所述呼吸信号的功率谱 曲线图,分析所述待测用户的自主神经系统,并将分析结果发送给所述移动终端,以使得所 述移动终端显示所述分析结果。

[0016] 本发明实施例的第三方面,提供一种自主神经系统的分析系统,所述系统包括:

[0017] 心电信号采集盒、移动终端以及服务器;

[0018] 所述心电信号采集盒,用于采集待测用户的心电信号,并将所述心电信号发送给 所述移动终έ而;

[0019] 所述移动终端,用于对接收到的所述心电信号进行预处理,并将预处理后的心电 信号发送给所述服务器;

[0020] 所述服务器,用于接收所述移动终端发送的所述预处理后的心电信号;根据所述 预处理后的心电信号,获得所述待测用户的心率信号和呼吸信号;对所述心率信号进行频 谱分析,以获得所述心率信号的功率谱曲线图;以及对所述呼吸信号进行频谱分析,以获得 所述呼吸信号的功率谱曲线图;根据所述心率信号的功率谱曲线图和所述呼吸信号的功率 谱曲线图,分析所述待测用户的自主神经系统,并将分析结果发送给所述移动终端,以使得 所述移动终端显示所述分析结果。

[0021] 本发明实施例与现有技术相比存在的有益效果是:本发明实施例通过将待测用户 的心率信号的功率谱曲线图和呼吸信号的功率谱曲线图相结合,共同作用分析待测用户的 自主神经系统,从而使得在分析待测用户的自主神经系统时,考虑呼吸信号对自主神经系 统的影响,提高了分析待测用户的自主神经系统的准确性。

附图说明

[0022] 为了更清楚地说明本发明实施例中的技术方案,下面将对实施例或现有技术描述 中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些 实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些 附图获得其他的附图。

[0023] 图1是本发明实施例一提供的自主神经系统的分析方法的实现流程图;

[0024] 图2是本发明实施例二提供的服务器的组成示意图;

[0025] 图3是本发明实施例三提供的自主神经系统的分析系统的组成示意图。

具体实施方式

[0026] 为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对 本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并 不用于限定本发明。

[0027] 实施例一:

[0028] 图1示出了本发明实施例一提供的自主神经系统的分析方法的实现流程,所述实 现流程详述如下:

[0029] 在步骤SlOl中,接收移动终端发送的预处理后的待测用户的心电信号;

[0030] 在本发明实施例中,服务器接收移动终端发送的预处理后的待测用户的心电信 号,所述预处理包括但不限于去基线和滤波处理,所述去基线处理用于去除心电信号中的 基线漂移,所述滤波处理用于滤除心电信号中的噪声。

[0031] 在本发明实施例中,服务器可以将移动终端发送的预处理后的待测用户的心电信 号以日期或者姓名的方式保存在数据库中,以便于学者或专家等进行统计学分析,深入挖 掘心电信号隐藏的其他特征。

[0032] 在步骤S102中,根据所述预处理后的心电信号,获得所述待测用户的心率信号和 呼吸信号;

[0033] 在本发明实施例中,服务器可以通过预处理后的心电信号的RR间期,获得心率数 值,再对心率数值进行插值运算可获得所述待测用户的心率信号,所述RR间期为预处理后 的心电信号的两个连续QRS波中R波之间的时间,由于心率为心脏每分钟跳动的次数,如果 RR间期的单位为分,则RR间期的倒数即为心率数值,如果RR间期的单位为秒,则60除以RR间 期即为心率数值,所述插值运算包括但不限于三次样条插值,插值频率可以为4Hz。

[0034] 示例性的,RR间期为0.8秒,则对应的心率数值为75次/分。

[0035] 进一步的,所述根据所述预处理后的心电信号,获得所述待测用户的呼吸信号可 以包括:

[0036] 步骤一,提取所述预处理后的心电信号的峰值,并对所述预处理后的心电信号的 峰值进行标号;

[0037] 在本发明实施例中,提取所述预处理后的心电信号的峰值及每个峰值所对应的时 间,并对每个峰值按顺序进行标号,即所述预处理后的心电信号的第一个峰值标记为1,第 二个峰值标记为2,以此类推,直到标记完所述预处理后的心电信号的所有峰值。

[0038] 步骤二,根据所述预处理后的心电信号的峰值,计算第i个峰值的前向平均值

Figure CN106096303BD00061

其中,R (j)为第j个峰值,i和j为大于零的整数; ,

[0039] 步骤三,根据所述第i个峰值和所述第i个峰值的前向平均值R⑴avg,计算与所述 第i个峰值对应的呼吸数值

Figure CN106096303BD00062

[0040] 在本发明实施例中,服务器根据所述第i个峰值和所述第i个峰值的前向平均值R (i)avg,计算与所述第i个峰值所对应的时间对应的呼吸数值

Figure CN106096303BD00063

,即所述呼吸数 值

Figure CN106096303BD00064

所对应的时间与所述第i个峰值所对应的时间相同。

[0041] 步骤四,重复执行步骤二和步骤三,直到计算出与所述预处理后的心电信号的所 有峰值对应的所有呼吸数值;

[0042] 在本发明实施例中,服务器可以先根据步骤二和步骤三,计算与所述预处理后的 心电信号的第一个峰值对应的呼吸数值,再根据步骤二和步骤三,计算与所述预处理后的 心电信号的第二个峰值对应的呼吸数值,根据步骤二和步骤三,计算与所述预处理后的心 电信号的第三个峰值对应的呼吸数值,以此类推,直到计算与所述预处理后的心电信号的 最后一个峰值对应的呼吸数值,即可计算出与所述预处理后的心电信号的所有峰值对应的 所有呼吸数值。

[0043] 步骤五,对所述所有呼吸数值进行插值运算,以获得所述待测用户的呼吸信号。

[0044] 在本发明实施例中,所述插值运算包括但不限于三次样条插值,插值频率可以为 4Hz〇

[0045] 在步骤S103中,对所述心率信号进行频谱分析,以获得所述心率信号的功率谱曲 线图;以及对所述呼吸信号进行频谱分析,以获得所述呼吸信号的功率谱曲线图;

[0046] 在本发明实施例中,服务器可以通过小波变换对所述心率信号进行频谱分析,获 得所述心率信号的功率谱曲线图,对所述呼吸信号进行频谱分析,获得所述呼吸信号的功 率谱曲线图,所述小波变换的小波基包括但不限于mor let小波,其表达式为:

Figure CN106096303BD00071

其中,f c为小波基的中心频率, 可取值为〇. 8Hz,t为小波基宽度,a为小波基的尺度,U用于调节小波基带通滤波器的宽度,U 值越小,带通滤波器的宽度就越窄,i为虚数单位,Pi为圆周率I

[0047] 在步骤S104中,根据所述心率信号的功率谱曲线图和所述呼吸信号的功率谱曲线 图,分析所述待测用户的自主神经系统,并将分析结果发送给所述移动终端,以使得所述移 动终端显示所述分析结果。

[0048] 在本发明实施例中,服务器将分析结果发送给所述移动终端,所述移动终端向用 户显示所述分析结果,可以使用户实时了解自己的自主神经状态。

[0049] 进一步的,所述根据所述心率信号的功率谱曲线图和所述呼吸信号的功率谱曲线 图,分析所述待测用户的自主神经系统可以包括:

[0050] 从所述呼吸信号的功率谱曲线图中选取出功率谱最大值所对应的频率,并将该频 率作为所述待测用户的呼吸频率frf,其中,frf大于零;

[0051] 根据所述呼吸频率frf,确定所述待测用户的交感神经的带宽和副交感神经的带 宽;

[0052] 根据所述心率信号的功率谱曲线图和所述交感神经的带宽,计算所述待测用户的 交感神经张力;以及根据所述心率信号的功率谱曲线图和所述副交感神经的带宽,计算所 述待测用户的副交感神经张力;

[0053] 在本发明实施例中,将所述交感神经的带宽对应到所述心率信号的功率谱曲线图 中,获得所述交感神经的带宽在所述心率信号的功率谱曲线图中对应的功率谱,所述功率 谱与所述交感神经的带宽在所述心率信号的功率谱曲线图中所围面积即为所述待测用户 的交感神经张力LFa;将所述副交感神经的带宽对应到所述心率信号的功率谱曲线图中,获 得所述副交感神经的带宽在所述心率信号的功率谱曲线图中对应的功率谱,所述功率谱与 所述副交感神经的带宽在所述心率信号的功率谱曲线图中所围面积即为所述待测用户的 副交感神经张力RFa。

[0054] 根据所述交感神经张力和所述副交感神经张力,分析所述待测用户的自主神经系 统。

[0055] 进一步的,所述根据所述呼吸频率frf,确定所述待测用户的交感神经的带宽和副 交感神经的带宽可以包括:

[0056] 根据所述呼吸频率f rf,确定所述待测用户的交感神经的带宽为Ή〜f rf *X2,其中, X1、X2大于零;

[0057] 较佳的,X1SO. 04,X2为0.65,交感神经的带宽单位为Hz。

[0058] 根据所述呼吸频率frf,确定所述待测用户的副交感神经的带宽为frf*X2〜frf* X3,其中,X3大于零。

[0059] 较佳的,X3为1.35, g[J交感神经的带宽单位为Hz。

[0060] 进一步的,所述根据所述交感神经张力和所述副交感神经张力,分析所述待测用 户的自主神经系统可以包括:

[0061] 当所述交感神经张力和所述副交感神经张力符合预设条件时,确定所述待测用户 处于自主神经平衡状态;

[0062] 在本发明实施例中,所述预设条件为LFa和RFa的取值范围为0.5〜lObpm2,且LFa/ RFa比值范围为0.4〜3,所述bpm2为所述交感神经张力和所述副交感神经张力的单位,可以 表示为单位时间内的心率。其中,所述单位时间为1分钟,所述心率为心脏每分钟跳动的次 数。

[0063] 当所述交感神经张力和所述副交感神经张力不符合预设条件时,确定所述待测用 户处于自主神经紊乱状态。

[0064] 本发明实施例通过将待测用户的心率信号的功率谱曲线图和呼吸信号的功率谱 曲线图相结合,共同作用分析待测用户的自主神经系统,从而使得在分析待测用户的自主 神经系统时,考虑呼吸信号对自主神经系统的影响,提高了分析待测用户的自主神经系统 的准确性。

[0065] 实施例二:

[0066] 图2示出了本发明实施例二提供的服务器的组成示意图,为了便于说明,仅示出了 与本发明实施例相关的部分,详述如下:

[0067] 接收模块21,用于接收移动终端发送的预处理后的待测用户的心电信号;

[0068] 获得模块22,用于根据所述预处理后的心电信号,获得所述待测用户的心率信号 和呼吸信号;

[0069] 进一步的,所述获得模块22具体包括:

[0070] 提取单元221,用于提取所述预处理后的心电信号的峰值,并对所述预处理后的心 电信号的峰值进行标号;

[0071] 第一计算单元222,用于根据所述预处理后的心电信号的峰值,计算第i个峰值的 前向平均值

Figure CN106096303BD00081

其中,R (j)为第j个峰值,i和j为大于零的整数;

[0072] 第二计算单元223,用于根据所述第i个峰值和所述第i个峰值的前向平均值R ⑴avg,计算与所述第i个峰值对应的呼吸数值

Figure CN106096303BD00082

[0073] 第三计算单元224,用于重复执行所述第一计算单元和所述第二计算单元,直到计 算出与所述预处理后的心电信号的所有峰值对应的所有呼吸数值;

[0074] 插值单元225,用于对所述所有呼吸数值进行插值运算,以获得所述待测用户的呼 吸信号。

[0075] 第一分析模块23,用于对所述心率信号进行频谱分析,以获得所述心率信号的功 率谱曲线图;以及对所述呼吸信号进行频谱分析,以获得所述呼吸信号的功率谱曲线图;

[0076] 第二分析模块24,用于根据所述心率信号的功率谱曲线图和所述呼吸信号的功率 谱曲线图,分析所述待测用户的自主神经系统,并将分析结果发送给所述移动终端,以使得 所述移动终端显示所述分析结果。

[0077] 进一步的,所述第二分析模块24具体包括:

[0078] 选取单元241,用于从所述呼吸信号的功率谱曲线图中选取出功率谱最大值所对 应的频率,并将该频率作为所述待测用户的呼吸频率frf,其中,frf大于零;

[0079] 确定单元242,用于根据所述呼吸频率frf,确定所述待测用户的交感神经的带宽 和副交感神经的带宽;

[0080] 计算单元243,用于根据所述心率信号的功率谱曲线图和所述交感神经的带宽,计 算所述待测用户的交感神经张力;以及根据所述心率信号的功率谱曲线图和所述副交感神 经的带宽,计算所述待测用户的副交感神经张力;

[0081] 分析单元244,用于根据所述交感神经张力和所述副交感神经张力,分析所述待测 用户的自主神经系统。

[0082] 进一步的,所述确定单元242具体包括:

[0083] 第一带宽确定子单元,用于根据所述呼吸频率frf,确定所述待测用户的交感神经 的带宽为〜frf*x2,其中,χι、χ2大于零;

[0084] 第二带宽确定子单元,用于根据所述呼吸频率frf,确定所述待测用户的副交感神 经的带宽为frf*X2〜frf*X3,其中,X3大于零。

[0085] 进一步的,所述分析单元244具体包括:

[0086] 第一张力确定子单元,用于当所述交感神经张力和所述副交感神经张力符合预设 条件时,确定所述待测用户处于自主神经平衡状态;

[0087] 第二张力确定子单元,用于当所述交感神经张力和所述副交感神经张力不符合预 设条件时,确定所述待测用户处于自主神经紊乱状态。

[0088] 本发明实施例提供的服务器可以使用在前述对应的方法实施例一中,详情参见上 述实施例一的描述,在此不再赘述。

[0089] 本发明实施例通过将待测用户的心率信号的功率谱曲线图和呼吸信号的功率谱 曲线图相结合,共同作用分析待测用户的自主神经系统,从而使得在分析待测用户的自主 神经系统时,考虑呼吸信号对自主神经系统的影响,提高了分析待测用户的自主神经系统 的准确性。

[0090] 实施例三:

[0091] 图3示出了本发明实施例三提供的自主神经系统的分析系统的组成示意图,为了 便于说明,仅示出了与本发明实施例相关的部分,详述如下:

[0092] 所述系统包括心电信号采集盒31、移动终端32以及服务器33;

[0093] 所述心电信号采集盒31,用于采集待测用户的心电信号,并将所述心电信号发送 给所述移动终纟而32;

[0094] 在本发明实施例中,所述心电信号采集盒31为便携式佩戴设备,可以通过心电传 感器实时采集待测用户的心电信号,并将所述心电信号通过无线方式发送给所述移动终端 32,所述无线方式包括但不限于蓝牙、WiFi等,所述待测用户为所述心电信号采集盒31佩戴 者。

[0095] 所述移动终端32,用于对接收到的所述心电信号进行预处理,并将预处理后的心 电信号发送给所述服务器33;

[0096] 在本发明实施例中,所述移动终端32为具有显示功能的终端,如手机、平板电脑 等,可以实时显示所述心电信号采集盒31发送的心电信号,并保存所述心电信号,以便于用 户的多次使用。

[0097] 在本发明实施例中,所述预处理包括但不限于去基线和滤波处理,可以采用0.8Hz 二阶巴特沃斯高通滤波器进行去基线处理,所述去基线处理用于去除心电信号中的基线漂 移,可以采用30Hz四阶巴特沃斯低通斯滤波器进行滤波处理,所述滤波处理用于滤除心电 信号中的噪声。所述移动终端32可以通过无线方式将预处理后的心电信号发送给所述服务 器33,所述无线方式包括但不限于蓝牙、WiFi等。

[0098] 所述服务器33,用于接收所述移动终端32发送的所述预处理后的心电信号;根据 所述预处理后的心电信号,获得所述待测用户的心率信号和呼吸信号;对所述心率信号进 行频谱分析,以获得所述心率信号的功率谱曲线图;以及对所述呼吸信号进行频谱分析,以 获得所述呼吸信号的功率谱曲线图;根据所述心率信号的功率谱曲线图和所述呼吸信号的 功率谱曲线图,分析所述待测用户的自主神经系统,并将分析结果发送给所述移动终端32, 以使得所述移动终端32显示所述分析结果。

[0099] 本发明实施例中的服务器33与实施例二中的服务器相同,详情参见上述实施例二 的描述,在此不再赘述。

[0100] 所述领域的技术人员可以清楚地了解到,为描述的方便和简洁,仅以上述各功能 模块的划分进行举例说明,实际应用中,可以根据需要而将上述功能分配由不同的功能模 块完成,即所述服务器的内部结构划分成不同的功能模块,上述功能模块既可以采用硬件 的形式实现,也可以采用软件的形式实现。另外,各功能模块的具体名称也只是为了便于相 互区别,并不用于限制本申请的保护范围。

[0101] 综上所述,本发明实施例将待测用户的心率信号的功率谱曲线图和呼吸信号的功 率谱曲线图相结合,共同作用分析待测用户的自主神经系统,从而使得在分析待测用户的 自主神经系统时,考虑呼吸信号对自主神经系统的影响,提高了分析待测用户的自主神经 系统的准确性。

[0102] 本领域普通技术人员还可以理解,实现上述实施例方法中的全部或部分步骤是可 以通过程序来指令相关的硬件来完成,所述的程序可以在存储于一计算机可读取存储介质 中,所述的存储介质,包括R0M/RAM、磁盘、光盘等。

[0103] 以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精 神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (4)

1. 一种服务器,其特征在于,所述服务器包括: 接收模块,用于接收移动终端发送的预处理后的待测用户的心电信号; 获得模块,用于根据所述预处理后的心电信号,获得所述待测用户的心率信号和呼吸 信号; 第一分析模块,用于对所述心率信号进行频谱分析,以获得所述心率信号的功率谱曲 线图;以及对所述呼吸信号进行频谱分析,以获得所述呼吸信号的功率谱曲线图; 第二分析模块,用于根据所述心率信号的功率谱曲线图和所述呼吸信号的功率谱曲线 图,分析所述待测用户的自主神经系统,并将分析结果发送给所述移动终端,以使得所述移 动终端显示所述分析结果; 所述获得模块具体包括: 提取单元,用于提取所述预处理后的心电信号的峰值,并对所述预处理后的心电信号 的峰值进行标号; 第一计算单元,用于根据所述预处理后的心电信号的峰值,计算第i个峰值的前向平均 值
Figure CN106096303BC00021
其中,R (j)为第j个峰值,i和j为大于零的整数; 第二计算单元,用于根据所述第i个峰值和所述第i个峰值的前向平均值R (i)avg,计算 与所述第i个峰值对应的呼吸数值
Figure CN106096303BC00022
第三计算单元,用于重复执行所述第一计算单元和所述第二计算单元,直到计算出与 所述预处理后的心电信号的所有峰值对应的所有呼吸数值; 插值单元,用于对所述所有呼吸数值进行插值运算,以获得所述待测用户的呼吸信号。
2. 根据权利要求1所述的服务器,其特征在于,所述第二分析模块具体包括: 选取单元,用于从所述呼吸信号的功率谱曲线图中选取出功率谱最大值所对应的频 率,并将该频率作为所述待测用户的呼吸频率frf,其中,frf大于零; 确定单元,用于根据所述呼吸频率frf,确定所述待测用户的交感神经的带宽和副交感 神经的带宽; 计算单元,用于根据所述心率信号的功率谱曲线图和所述交感神经的带宽,计算所述 待测用户的交感神经张力;以及根据所述心率信号的功率谱曲线图和所述副交感神经的带 宽,计算所述待测用户的副交感神经张力; 分析单元,用于根据所述交感神经张力和所述副交感神经张力,分析所述待测用户的 自主神经系统。
3. 根据权利要求2所述的服务器,其特征在于,所述确定单元具体包括: 第一带宽确定子单元,用于根据所述呼吸频率frf,确定所述待测用户的交感神经的带 宽为Xl〜frf*X2,其中,X1、X2大于零; 第二带宽确定子单元,用于根据所述呼吸频率frf,确定所述待测用户的副交感神经的 带宽为frf*X2〜frf*X3,其中,X3大于零; 所述分析单元具体包括: 第一张力确定子单元,用于当所述交感神经张力和所述副交感神经张力符合预设条件 时,确定所述待测用户处于自主神经平衡状态; 第二张力确定子单元,用于当所述交感神经张力和所述副交感神经张力不符合预设条 件时,确定所述待测用户处于自主神经紊乱状态。
4. 一种自主神经系统的分析系统,其特征在于,所述系统包括: 心电信号采集盒、移动终端以及服务器; 所述心电信号采集盒,用于采集待测用户的心电信号,并将所述心电信号发送给所述 移动终端; 所述移动终端,用于对接收到的所述心电信号进行预处理,并将预处理后的心电信号 发送给所述服务器; 所述服务器,用于接收所述移动终端发送的所述预处理后的心电信号;根据所述预处 理后的心电信号,获得所述待测用户的心率信号和呼吸信号;对所述心率信号进行频谱分 析,以获得所述心率信号的功率谱曲线图;以及对所述呼吸信号进行频谱分析,以获得所述 呼吸信号的功率谱曲线图;根据所述心率信号的功率谱曲线图和所述呼吸信号的功率谱曲 线图,分析所述待测用户的自主神经系统,并将分析结果发送给所述移动终端,以使得所述 移动终端显示所述分析结果; 所述服务器,具体用于执行步骤一,提取所述预处理后的心电信号的峰值,并对所述预 处理后的心电信号的峰值进行标号;步骤二,根据所述预处理后的心电信号的峰值,计算第
Figure CN106096303BC00031
i个峰值的前向平均值 其中,R (j)为第j个峰值,i和j为大于零的整数;步 9 骤三,根据所述第i个峰值和所述第i个峰值的前向平均值R⑴avg,计算与所述第i个峰值对 应的呼吸数彳1
Figure CN106096303BC00032
;步骤四,重复执行步骤二和步骤三,直到计算出与所述预处理 后的心电信号的所有峰值对应的所有呼吸数值;步骤五,对所述所有呼吸数值进行插值运 算,以获得所述待测用户的呼吸信号。
CN201610460442.3A 2016-06-22 2016-06-22 一种自主神经系统的分析方法、服务器及系统 CN106096303B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610460442.3A CN106096303B (zh) 2016-06-22 2016-06-22 一种自主神经系统的分析方法、服务器及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610460442.3A CN106096303B (zh) 2016-06-22 2016-06-22 一种自主神经系统的分析方法、服务器及系统

Publications (2)

Publication Number Publication Date
CN106096303A CN106096303A (zh) 2016-11-09
CN106096303B true CN106096303B (zh) 2018-12-18

Family

ID=57253139

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610460442.3A CN106096303B (zh) 2016-06-22 2016-06-22 一种自主神经系统的分析方法、服务器及系统

Country Status (1)

Country Link
CN (1) CN106096303B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107115657A (zh) * 2017-05-16 2017-09-01 广东乐源数字技术有限公司 一种运动指导系统
CN107212480A (zh) * 2017-05-16 2017-09-29 广东乐源数字技术有限公司 一种监测运动呼吸率的智能服装及运动指导系统

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1527235A (zh) * 2003-09-22 2004-09-08 南京大学 获得自主神经功能信号的方法
CN101637383A (zh) * 2008-08-01 2010-02-03 深圳市新元素医疗技术开发有限公司 一种心电信号的监测方法、系统及移动终端
CN101642369A (zh) * 2008-08-04 2010-02-10 南京大学 自主神经功能生物反馈方法和系统
CN104127193A (zh) * 2014-07-14 2014-11-05 华南理工大学 一种抑郁症程度量化的评估系统及其评估方法
CN204293140U (zh) * 2014-10-12 2015-04-29 吴健康 一种心率变异性生物反馈康复装置
CN204745287U (zh) * 2015-06-05 2015-11-11 华南理工大学 一种基于心率变异性的催眠治疗抑郁症效果评价系统
CN204931634U (zh) * 2015-07-30 2016-01-06 华南理工大学 基于生理信息的抑郁症评估系统
CN105380618A (zh) * 2014-08-27 2016-03-09 精工爱普生株式会社 生物体信息解析系统及生物体信息处理系统

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1527235A (zh) * 2003-09-22 2004-09-08 南京大学 获得自主神经功能信号的方法
CN101637383A (zh) * 2008-08-01 2010-02-03 深圳市新元素医疗技术开发有限公司 一种心电信号的监测方法、系统及移动终端
CN101642369A (zh) * 2008-08-04 2010-02-10 南京大学 自主神经功能生物反馈方法和系统
CN104127193A (zh) * 2014-07-14 2014-11-05 华南理工大学 一种抑郁症程度量化的评估系统及其评估方法
CN105380618A (zh) * 2014-08-27 2016-03-09 精工爱普生株式会社 生物体信息解析系统及生物体信息处理系统
CN204293140U (zh) * 2014-10-12 2015-04-29 吴健康 一种心率变异性生物反馈康复装置
CN204745287U (zh) * 2015-06-05 2015-11-11 华南理工大学 一种基于心率变异性的催眠治疗抑郁症效果评价系统
CN204931634U (zh) * 2015-07-30 2016-01-06 华南理工大学 基于生理信息的抑郁症评估系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于心电信号中提取呼吸信号的探讨;赵素文等;《桂林航天工业高等专科学校学报》;20120315(第1期);20-24 *
基于心电信号提取呼吸信号的算法;郭永丛等;《吉林大学学报( 信息科学版)》;20160515;第34卷(第3期);327-333 *

Also Published As

Publication number Publication date
CN106096303A (zh) 2016-11-09

Similar Documents

Publication Publication Date Title
Charlton et al. An assessment of algorithms to estimate respiratory rate from the electrocardiogram and photoplethysmogram
Wöstmann et al. Spatiotemporal dynamics of auditory attention synchronize with speech
Paalasmaa et al. Adaptive heartbeat modeling for beat-to-beat heart rate measurement in ballistocardiograms
CN105147248B (zh) 基于生理信息的抑郁症评估系统及其评估方法
Lauritzen et al. Top–down flow of visual spatial attention signals from parietal to occipital cortex
Chu et al. Impulsive noise suppression and background normalization of electrocardiogram signals using morphological operators
Taelman et al. Instantaneous changes in heart rate regulation due to mental load in simulated office work
Barry et al. Muscle sounds are emitted at the resonant frequencies of skeletal muscle
Gysels et al. Phase synchronization for the recognition of mental tasks in a brain-computer interface
Ivanov et al. Levels of complexity in scale-invariant neural signals
Blechert et al. ANSLAB: Integrated multichannel peripheral biosignal processing in psychophysiological science
US7460901B2 (en) Procedure for deriving reliable information on respiratory activity from heart period measurement
KR101386649B1 (ko) 사용자의 상태 적용 게임 장치 및 그 게임 제공 방법
Sato et al. A novel extraction method of fetal electrocardiogram from the composite abdominal signal
Rosenblum et al. Phase synchronization: from theory to data analysis
Lohman et al. A digital signal processor for Doppler radar sensing of vital signs
Towle et al. Electrocorticographic coherence patterns
Salehizadeh et al. A novel time-varying spectral filtering algorithm for reconstruction of motion artifact corrupted heart rate signals during intense physical activities using a wearable photoplethysmogram sensor
Xu et al. Pattern recognition of motor imagery EEG using wavelet transform
Lutz et al. Guiding the study of brain dynamics by using first-person data: synchrony patterns correlate with ongoing conscious states during a simple visual task
US20160242665A1 (en) Devices and methods for real-time denoising of electrocardiograms
Wang et al. Fractal dimension based neurofeedback in serious games
CN103584840B (zh) 基于脑电、心率变异及两者相干性的自动睡眠分期方法
Schäfer et al. Estimation of breathing rate from respiratory sinus arrhythmia: comparison of various methods
US9675282B2 (en) Method and apparatus for determining a respiration signal

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant