使用小波熵来证明正念练习如何增加不规则脑与心脏活动之间的协调

摘要

该手稿描述了如何使用小波熵指数来分析高密度脑电图（EEG）和心电图（ECG）数据。我们表明，在基于正念的减压实践中，脑和心脏活动的不规则性变得更加协调。

摘要

在东西方，传统教义说，心灵与心灵有着密切的联系，特别是在精神实践中。客观地证明这一点的一个难点在于脑和心脏活动的本质是完全不同的。在本文中，我们提出了一种使用小波熵来测量脑电图（EEG）和心电图（ECG）数据的混沌水平的方法，并展示了如何在不同的实验条件下探索心灵与心脏之间的潜在协调。此外，使用统计参数映射（SPM）来鉴定EEG​​小波熵受实验条件影响最大的脑区域。作为一个例证，在8周标准的基于正念的减轻压力（MBSR）训练课程（预测试）开始时，脑电图和心电图记录在两种不同的条件（正常休息和心灵呼吸）下，之后c我们（后测）。使用所提出的方法，结果一致地表明，与封闭眼睛静止状态期间相比，在脑脊髓脑缺血呼吸状态期间，脑EEG的小波熵降低。类似地，在MBSR心灵呼吸期间发现较低的心率的小波熵。然而，在预测试和后测试之间发现了MBSR注意呼吸期间小波熵的差异。在所有参与者的正常休息期间，脑波熵与心率的熵之间没有观察到相关性，而在MBSR心灵呼吸期间观察到显着的相关性。此外，最相关的脑区域位于大脑的中心区域。这项研究提供了建立正义实践（ 即注意呼吸）可能增加心灵和心脏活动之间协调的证据的方法。

引言

随着神经科学和医学的进步¹ ，我们在了解正念和冥想对精神和身体疾病的非药理学益处方面处于更好的地位。许多研究表明，基于心灵的压力减轻（MBSR）训练是一种精神锻炼，可以有效减少心灵和身体的不必要的压力。在冥想的实践中，身体和心灵可能会变得更加夹带。这些问题吸引了学术界和大众的兴趣，并促使我们设计出一种在不同条件下探索心灵与身体（或更具体地，大脑与心脏之间）潜在协调的方法来评估效果的MBSR培训和实践。

脑电图不均匀性（EEG）在常规EEG谱分析中通常被认为是噪声，b它可能包含有用信息，可以通过信号熵来测量。 EEG熵被用作癫痫，阿尔茨海默病和麻醉等临床病症的指标，一些证据表明，脑电图的熵可能比传统的EEG频谱分析更为敏感。熵长期以来被用于热力学，作为描述系统混沌状态的指标。由脑电测量的脑电波由于大脑网络在秩序和混沌之间的运作而闻名于其动态活动。在各种熵中，小波熵提供了一种在不同观测尺度下计算能量熵的方法，非常适合于分析不同尺度的具有高度变化特征的脑电信号。对于完全随机的信号，不同分辨率级别的相对能量将是相同，熵将达到最大值。因此，小波熵作为脑电混沌活动的量度，可以作为区分不同心态的指标。

以前的心电图（ECG）研究提供了与冥想状态相比心悸变异性的冥想相关变化的一致证据⁷ 。心跳由自主神经（包括交感神经和副交感神经）神经系统调节。自主神经调节由神经节之间的相互作用调节;中耳节迷走神经和高阶脑结构，如延髓，岛内和边缘系统。这些结构形成一个分层系统，其非线性活动可以通过熵⁸进行测量。例如，高阶脑结构中的神经突发可以影响心跳节律⁹ 。在临床情况下，脑紊乱如癫痫发作会扰乱心律¹⁰ 。神经心脏学研究还表明心脏和大脑之间的紧密相互作用¹¹ 。

检查脑和心脏信号之间的关联的一个困难在于它们在信号内容和时间尺度上是完全不同的。因此，一个主要的挑战是确定一个不仅适用于脑和心脏信号的常见措施，而且对于身心协调的背景下的解释也是有意义的。在本文中，我们提出使用小波熵来检测脑和心脏活动不规则之间的可能协调关系，并将脑电信号与心跳信号相关联，然后通过MBSR实验进行评估。

因为MBSR涉及身体和心灵的正念实践，我们假设MBSR正念呼吸动作可能不仅影响脑活动，还可能影响心跳。以前的研究主要根据不同方面分别研究正念训练和/或实践对脑或心脏的影响;这两个紧密相关的系统的同时变化仍然在很大程度上未被探索。只有少数研究报告，正念训练后身体和心灵的协调性更好¹² 。使用提出的方法，代表性的结果和以前的研究¹³显示，MBSR正念训练可能会增加身体和心灵在其混沌活动方面的协调，这可能会提供新的洞察心理训练对中央和外围环境的影响神经系统。

研究方案

本研究由香港大学机构审查委员会批准。所有参与者在参与本研究之前提供书面知情同意书。

实验设计

表1. MBSR课程概述。参加者8周的MBSR训练课程。

1. 招募参加者
   注意:来自当地MBSR课程的11名健康成年人（5名女性）参加了这项脑电图研究，年龄为28-52岁。基于Beck抑郁症患者抑郁症的参与者被排除在研究之外。课程遵循标准的MBSR培训（ 表1 ），参与者必须承诺培训结构。
2. 脑电录音系统
   1. 使用128通道EEG在安静的房间内进行脑电图记录ystem由EEG帽，放大器，流浆箱和台式电脑组成。

EEG和ECG数据采集

注意:EEG数据收集遵循类似于之前描述的方法¹⁴的标准程序。要求每个参与者在进入实验地点之前清洗头发和头皮。通知参与者实验程序，主要是他/她将正常休息10分钟，并在记录EEG和ECG设备的同时进行10分钟的呼吸。

1. 用酒精拭子清洁参与者的脸部区域和乳突。
2. 使用卷尺，测量参与者的头围，然后选择适当大小的帽子。从nasion进行一次测量，并在耳朵顶部和头皮上进行另一个测量。标记顶点（nasion之间的中间点两只耳朵之间的距离和中距离）用柔软的标记笔。
3. 根据10-5电极系统¹⁵设置电极位置。定位帽子，使得Cz电极在顶点之上，Nz电极处于nasion，Lz电极处于该位置，RM电极处于右乳突，LM电极处于左乳突。
   1. 将ECG电极置于左右锁骨下窝。
   2. 使用钝点注射器将电极固定器用凝胶填充。
   3. 将每个电极的阻抗保持在20kΩ以下。通过调节电极位置来减少阻抗，增加与头皮的接触;必要时加入更多的凝胶。
   4. 拧紧下腹部，请求参与者保持静止。
   5. 将EEG设备的采样速率设置为1,000 Hz。使用左乳突作为原始参考点。
4. 记录EEG数据（步骤1.2。1）在MBSR课程开始时（2周内）。
   1. 让参与者进行简短的身体扫描以放松身体。参与者呼吸呼吸时要注意呼吸。在EEG数据采集期间，每个参与者进行10分钟的MBSR心灵呼吸（MBSR正念状态）和10分钟的正常休息（控制条件）。这产生了具有两个条件的前MBSR训练数据集。
      注意:平衡呼吸暂停的顺序和参与者之间的正常休息。
5. 参加者完成MBSR课程约1个月后再次记录脑电数据。
   1. 请求每位参与者在脑电数据采集期间进行10分钟的正常呼吸和10分钟的正常休息。这产生了具有两个条件的后MBSR训练数据集。

EEG和ECG数据分析

注意:前四个步骤是用于数据预处理，其余步骤用于小波熵计算和相关分析。在这里，一个名为EEGLAB ¹⁶的开源软件用于EEG分析，但是不同软件的操作应该是相似的。以下将以EEGLAB为基础的操作进行说明。有关详细信息，请参阅EEGLAB手册（https://sccn.ucsd.edu/wiki/EEGLAB_Wiki#EEGLAB_Tutorial）。

1. 使用EEG软件通过选择"工具">"更改采样率"，以250 Hz重新采样数据。使用有限脉冲响应（FIR）滤波器，通过选择工具>滤波数据>基本FIR滤波器，通过0.5 - 100 Hz通带滤波。
   1. 为了减少由于国家特定的市电交流电而产生的噪音，请使用短路非线性无限脉冲响应（IIR）滤波器，通过选择"工具>滤波数据>短非线性IIR过滤器。这覆盖了频谱分析中的所有相关频率，同时消除了不必要的噪音。
2. 使用EEG软件通过选择Plot>通道数据（滚动）来可视地滚动并检查EEG信号。
   1. 识别和删除包含明显肌肉噪音和任何其他奇怪事件的脑电图段。左键单击并将鼠标拖放到不良段以突出显示，然后使用EEG软件"删除"或"拒绝"段。
      注意:肌肉活动如吞咽或面部表情常常产生不连续和不规则的信号。奇怪的事件通常由头部或身体动作引起，这通常导致EEG信号的显着偏移（低频）。
   2. 确定是否有任何不良通道。通过选择"工具">"插值"通道，使用球面插值方法重建每个不良通道。
3. 使用EEG软件执行In依赖组件分析¹⁷ （ICA）对数据的选择工具>运行ICA。有经验丰富的EEG操作员通过选择工具>拒绝数据使用ICA>拒绝组件通过地图，目视识别和丢弃眼睛运动和眨眼，肌肉运动和其他可能的噪音的组成部分。工具>删除组件。 EEG软件将使用保留的组件自动重建数据。
4. 在进一步分析之前，使用EEG软件将数据重新引导至所有通道的平均值，方法是选择"工具">"重新引用"。
5. 频谱分析。
   1. 使用EEG软件命令行功能"spectopo"使用Welch的功率谱密度估计¹⁸计算EEG的光谱，并获得δ（1-4Hz），θ（4-8Hz），α（ 8-12Hz），β（12-30Hz）和γ（30-80Hz）波。
   2. 计算相对po每种类型的波¹⁸的频率（ 即相对于总频谱功率的各个频带的功率比例）。
6. 使用统计参数映射（SPM）方法进行脑电图源分析。
   1. 使用最小范数估计¹⁹将头皮电极的脑电信号变换成3D脑源信号。
      注意:这里，原来的频道号是122，并被转换为分布在三维脑模型中的8,196个频道。显着性水平设定为p = 0.01。请参阅SPM手册中的源分析（www.fil.ion.ucl.ac.uk/spm/doc/spm8_manual.pdf）。
7. 使用以下步骤计算每个EEG通道的小波熵（以及估计的三维脑源信号）:
   1. 给了一个信号 ，计算小波系数为 ，哪里是Haar小波族的正交基础表示分解水平（本研究中N = 10）。
   2. 将相对能量定义为 ，与 。
   3. 计算小波熵为 。平均所有信道的小波熵获得每个参与者的EEG小波熵。
8. 使用软件包对ECG数据进行峰值分析，以获得心跳信号，其由ECG波峰的数量和峰之间的间隔组成。
   注意:假设在正常休息和正常呼吸状态下，正常心跳60-100次/分钟，10分钟持续时间的心跳情况N的数量应相对较小，满足600 < N <1,000。
   1. 为了从心跳信号获得可靠的小波熵，使用500点的滑动窗口和10点的步进增量来计算一系列小波熵，采用与EEG相同的程序，详见3.7节。使用平均值作为最终的小波熵。
9. 使用具有Pearson相关函数的任何统计软件，分析心电图的小波熵与头皮每个通道的脑电图脑的脑和心脏相关性。
   1. 使用EEG软件命令行功能"topoplot"来构建相关性的统计头皮图。绘制线性相关图与心率熵和平均EEG熵只有那些nnels与心率熵显着相关（见代表性结果部分）。

结果

频谱分析

在EEG数据的频谱分析中，与正常休息相比，在MBSR心灵呼吸期间，增强了α（8-12Hz）和β（12-30Hz）和减小的δ（1-4Hz）波。 α波的增加在全球有显着性，特别是在额叶和枕叶，而β波的增量主要在额叶。在中部顶部地区注意到减少的三角洲波动（ 图1 ）。然而，我们在MBSR培训之前和之后没有发现重大变化。

图1. MBSR意识呼吸和正常休息条件的频谱分析。频谱分析表明，MBSR实践可以改变大脑活动，如不同之处所示脑电图谱。原始来源:参考文献¹³ 。 请点击此处查看此图的较大版本。

图2示出了所定义的感兴趣区域:枕叶（通道N19-21，N41-46，N67-72，N96-100和N119），中额叶（通道N53 -61和N79-83），以及中间顶叶（通道N48-50，N64-66和N74-776）。 表2显示了每个频谱的功率和MBSR心灵呼吸和正常休息条件下的熵，以及每个感兴趣区域的两个条件之间的差异的t检验的p值。我们结合了前后MBSR训练EEG数据，因为两个阶段之间没有差异。

图2.定义的兴趣区域 EEG通道的定义区域代表中额叶，中叶和枕叶。 请点击此处查看此图的较大版本。

		三角洲	西塔	Α	Beta版	伽玛	熵
中前额叶	MBSR	18.96±3.10	22.29±5.86	30.90±7.76	0.12±43.77	-467.25±79.48	0.753±0.060
	休息	22.40±6.56	22.76±5.98	26.91±7.14	-7.11±42.27	-449.76±102.92	0.785±0.066
	p值	0.0243	0.6555	0.0085	0.0114	0.4419	0.0084
枕叶	MBSR	15.39±3.30	17.91±6.53	33.38±6.50	-8.03±40.72	-466.23±61.46	0.719±0.048
	休息	19.82±6.32	20.38±8.11	29.73±5.30	-11.62±40.22	-439.03±102.81	0.763±0.055
	p值	0.0134	0.1213	0.0125	0.0796	0.8936	0.0098
顶叶	MBSR	17.95±3.60	19.46±6.54	32.84±6.68	2.68±38.23	-487.72±104.13	0.738±0.072
	休息	21.24±6.37	21.16±8.00	29.79±6.60	-1.27±34.42	-490.01±123.83	0.764±0.075
	p值	0.0157	0.0963	0.0177	0.1507	0.2878	0.0368

表2.脑区的频谱和熵。 MBSR注意呼吸和正常休息条件在三个预定义的感兴趣区域（ROI）中进行了比较。

小波熵分析

脑电图的小波熵分析显示，与正常休息相比，MBSR心灵呼吸期间的脑电熵减少与MBSR训练前和后期均有关。因为两个阶段之间没有显着差异，所以它们被合并以产生一个平均的图（ 如图3的第三行所示 ）。颅骨下降的主要区域为额叶和顶叶。

图3.小波熵分析。休息1和MBSR 1表示MBSR前培训，而休息2和MBSR 2表示MBSR后培训。小波熵分析表明，MBSR实践可以减少不规则性脑电子活动。原始来源:参考文献¹³ 。 请点击此处查看此图的较大版本。

小波熵的源分析

EEG信号的源分析可以通过将头皮脑电图解卷积到皮质表面^20上的电活动来提高空间分辨率。分析表明，受MBSR正念训练影响的主要脑区位于左中枕叶，前叶，上颞叶和左梭梭（ 图4 ）。在表3中 ，对于从源分析获得的四个解剖标签，我们提供了两个不同状态的熵差的t检验结果，即MBSR正念状态与正常休息。

图4.源分析。来源分析显示，在MBSR心灵呼吸状态期间，各种脑区域（红色突出显示）的熵减少。 请点击此处查看此图的较大版本。

解剖标签	素	x，y，z	t值
中间颞回	1728	-44 -60 16	3.77
原生动物	1324	2 -56 28	3.63
枕骨L	749	-2 -102-12	3.61
梭形L	142	-102	3.51

表3.源分析报告。解剖脑区域与MBSR正气呼吸与正常休息条件之间的熵差异显着。 L，大脑左侧。 R，大脑右侧（ p <0.001， t检验，未校正）。

心率小波熵

分析显示，平均心率没有显着性差异，但MBSR注意呼吸（MBSR训练前后）的心率熵较低（ 表4 ）。

	Rest1	MBSR1	Rest2	MBSR2
心率（每分钟节拍）	68.2±9.5	67.7±9.3	71.8±8.1	70.7±8.4
心率小波熵	0.89±0.05	0.79±0.11 *	0.89±0.07	0.80±0.12 ＃

表4. MBSR心灵呼吸和正常休息条件的心率和心率小波熵。在MBSR训练前和后期比较两个条件。 *休息1和MBSR 1之间的显着差异（ p <0.05）。＃2与MBSR 2之间的显着性差异（ p <0.05）。

关联心与脑活动之间

鉴于大脑和心脏之间的潜在联系，我们分析了全脑脑电图熵（ 即单个受试者的所有脑电信道上EEG熵的平均值）与MBSR心灵呼吸和正常休息状态之间的心跳熵的相关性。大脑和心脏的熵在MBSR心脏呼吸期间显着相关，但在正常休息期间并不显着。 如图5的上半部分所示 ，每个通道的脑电图熵与心率熵之间的相关性的意义在脑中央部位最突出。下图显示了中心区域（与心率熵显着相关的虚线通道的平均值）和心率熵之间的EEG熵的相关性。

图5.心与脑电子活动的小波熵之间的相关性。大脑活动的小波熵是中心区域重要通道的平均值，如右上图所示。心脏的小波熵是基于心跳间隔。在所有受试者之间计算MBSR训练前和后期训练结果的相关性（下图）。原始来源:参考文献¹³ 。 请点击此处查看此图的较大版本。

讨论

这项EEG研究提供了具体的证据，与正常休息相比，MBSR心脏呼吸增强了正面和枕叶的α波。这与大多数关于冥想中增加α波的报道一致，特别是初学者²¹ 。阿尔法波浪通常表现出唤醒和警惕的状态，不受任何具体任务的束缚。这样的精神状态是冥想的目的之一;也就是说，保持清醒和清晰，无任何幻想的航班或入睡。在MBSR实践中增强的α波同步也可以增强心灵与身体²²之间的协调。

大脑中大量神经元的电子活动形成非线性分层系统，其不规则性可以通过熵⁸很好地测量。我们发现在MBSR实践中，脑电的不规则大部分头皮脑电图，尤其是额叶和双侧枕叶脑区域的活动减少。小波熵的进一步来源分析表明，MBSR实践中枕骨和脑脊骨活动的不规则性降低。涉及的其他领域包括右中扣带和上颞叶。这符合另一项研究，表明在MBSR培训²³之后，灰质增加。高级冥想从业者通常在非快速眼动睡眠期间具有较高的枕头脑电图伽马波¹ 。 MBSR训练的一个重要特征是不考虑脑/脑信息流。这可能潜在地减少大脑活动的不规则。因为大脑中的主要信息来源是视觉皮质，视觉相关的大脑区域比其他区域更受MBSR实践的影响是合理的纳秒。这使得心灵对不相关信息的反应较少，并且更注意注意呼吸。

小心熵分析显示，与正常休息状态相比，注意呼吸期间心脏心脏的不规则性较少，这意味着由于副交感神经和交感神经紧张之间的相对不稳定的相互作用，自主神经系统的精确平衡。注意呼吸也可能影响心脏收缩，因为肺感觉传入对中脑中的迷走神经运动神经元发出抑制性投射，其通过迷走神经²⁴抑制心动过速。心脑血管相互作用由脑干中的脑干和边缘系统介导²⁵ 。鉴于中枢神经系统与心脏活动之间的联系，似乎合理的是，通过训练心灵，MBSR或其他冥想练习也可以影响其他身体功能，如重置压力反射灵敏度，提高肺部气体交换效率，平衡自主神经系统²⁶ 。

许多东方传统都强调了心灵的联系。在这方面，我们的研究提供了熵的第一个客观证据。在MBSR实践中，脑和心脏活动的不规则变得更加协调，双边体感区域的协调最大。这个地区直接负责人体的感官输入和运动，这种增加的协调意味着在MBSR实践中更多的身体意识。发现正念训练在功能MRI研究中增加了内脏区域的活动，其具有比EEG ²⁷更好的空间分辨率。我们的研究结果表明，在MBSR实践中，大脑和心脏可能变得更加夹带，从而在系统中卸载不必要的熵 up class ="xref"> 28。因为心跳是由自主神经系统调节的，其中自主神经调节由神经节之间的相互作用调节;中耳节迷走神经和高阶脑结构，如延髓，胰岛和边缘系统，自主神经和中枢神经系统在MBSR实践中可能变得更加协调。

在我们这样的研究中的一个重要步骤是选择正确的熵指数来测量脑电图和心电图的不规则性。与使用近似熵^29的 Renu Madhavi 等人不同，我们选择使用小波熵来集中在不同尺度的能量之间的复杂性。由于脑电子活动源自于不同尺度和频谱具有不同特征的大量神经元，所以小波熵应更适合于测量脑电图的不规则性> 30。类似地，瞬时心率在不同的时间尺度下也具有不同的特征，因此对心电图³¹应用小波熵是合理的。此外，香农的熵和近似熵对采样频率敏感，近似熵需要大量的计算时间，与数据长度呈指数增长，除了源空间中的计算外。

目前研究的几个限制值得注意。首先，头皮脑电图的空间分辨率相对较差，尽管其时间分辨率高。为了解决这个问题，我们的研究中采用了源分析。 EEG源重建是一个不成问题，另一种方法是从以前的功能磁共振成像（fMRI）研究中获得空间先验，这有助于提高重构精度^32,33 。另一个限制是t帽子之间的脑和心脏活动之间的协调被评估受试者之间;课程内协调将留给未来的研究。另外，心电图的熵不是直接从ECG信号中得到的，而是来自心跳间隔，因为心电图的振幅可能受到电极位置和皮肤阻抗的很大影响，因此在受试者之间是不可比的。

以前的研究通常对EEG信号进行频谱分析，将EEG中的不规则视为噪声。相比之下，我们建议使用小波熵作为更敏感的指标来测量大脑的混沌电子活动，作为思维状态的反映。所提出的方法提供了客观证据，表明思维和心灵在冥想中更加协调的传统假设。

因为小波熵可以很敏感地监视冥想状态，所以可以用于心理训练正如在临床实践中使用熵来测量手术期间的麻醉水平一样。使用类似于小波熵的索引，可以使用无线EEG设备来方便地监视个人的日常冥想练习或用于教育目的的其他心理训练。然而，由于无线EEG通常只能稳定地测量前额中的通道，因此我们需要探讨这些通道是否足够敏感以监测心理状态。

该过程的关键步骤包括将每个电极的皮肤阻抗保持为低，以获得更好的信号质量和更低的噪声;适当进行MBSR培训;在实验之前确保一个好的，不疲劳的精神状态，因为冥想需要很大的精力;选择适当的熵指数（小波熵）;并熟悉EEG数据的SPM源分析。

材料

Name	Company	Catalog Number	Comments
Synamps 2/RT 128-channel EEG/EP/ERP	Compumedics NeuroScan	Model 8050	Amplifiers and EEG acquisition
Quik-CapEEG 128 electrodes	Compumedics NeuroScan		EEG caps
Quik-Gel	Compumedics NeuroScan		conductive electrode gel
EEGLAB	Swartz Center for Computational Neuroscience		EEG analysis software
Matlab R2011a	MathWorks		EEGLAB is based on Matlab
SPM	Wellcome Trust Centre for Neuroimaging		EEG source analysis software