２９卷５期２０１０年ｌＯ月
中国生物医学工程学报ｏｆＢｉｏｍｅｄｉｃａｌ
Ｅｎｇｉｎｅｅｒｉｎｇ
Ｖ０１．２９Ｏｅｔｏｂｅｒ
Ｎｏ．５２０１０
ＣｈｉｎｅｓｅＪｏｕｒｎａｌ
小波多重分形在脑电信号分析中的应用
赵大庆１
王
俊２‘
１（南京邮电大学，图像处理与图像通信江苏省重点实验室，通信与信息工程学院，南京２１０００３）２（南京邮电大学，图像处理与图像通信江苏省重点实验室，地理与生物信息学院，南京２１０００３）
摘要：脑电图（ＥＥＧ）作为客观反映大脑机能状态的一个重要方面，蕴涵着丰富的生理、心理和病理信息。分形和小波相结合的方法，发挥了小波多尺度分析和分形标度不变性的两个特点。利用小波模极大值多重分形的方法，对青老年人群的脑电图进行分析。研究发现，脑电信号具有奇异性质，脑电信号奇异谱的跨度会随着年龄的增大而减小，可以用来检测目标个体是否出现脑衰退迹象，从而辅助进一步的临床干预措施。关键词：多重分形分析；脑电图；小波模极大值
中图分类号Ｒ３１８．５
文献标识码Ｄ
文章编号０２５８—８０２１（２０１０）０５－０７９３－０４
Ａｐｐｌｉｃａｔｉｏｎｏｆ
ＷａｖｅｌｅｔＭｕｌｔｉｆｒａｃｔａｌ
Ｓｐｅｃｔｒｕｍｔｏｔｈｅ
ＡｎａｌｙｓｉｓｏｆＥｌｅｃｔｒｏｅｎｃｅｐｈａｌｏｇｒａｍＳｉｇｎａｌｓ
ＺＨＡＯ
１
Ｄａ．Ｑｉｎ９１
ＷＡＮＧＪｕｎ２’
ｏｆＴｅｌｅｃｏｍｍ．＆Ｉｎｆｏｒｍａｔｉｏｎ
ｌＩｍａｇｅＰｒｏｃｅｓｓｉｎｇａｎｄＩｍａｇｅＣｏｍｍｕｎｉｃａｔｉｏｎＫｅｙＬａｂ。Ｃｏｌｌｅｇｅ
Ｅｎｇｉｎｅｅｒｉｎｇ。Ｎ８咏ｎｇＵｎｉｔｅｒｓ幼ｏｆＰｏｓｔｓ＆Ｔｅｌｅｃｏｍｍｕｎｉｃａｔｉｏｎ。ＮⅡｎｊｉｎｇ２１０００３。Ｃｈｉｎａ、
２（Ｉｍａｇｅ
Ｐｒｏｃｅｓｓｉｎｇ
ａｎｄＩｍａｇｅ
ＣｏｍｍｕｎｉｃａｔｉｏｎＫｅｙＬａｂ。Ｃｏｌｌｅｇｅｏｆ
Ｇｅｏｇｒａｐｈｉｃａｌ＆Ｂｉｏｌｏｇｉｃａｌ
Ｉｎｆｏｒｍａｔｉｏｎ。
ＮａｎｊｉｎｇＵｎｉｔｅｒｓｉｔｙｏｆＰｏｓｔｓ＆Ｔｅｌｅｃｏｍｍｕｎｉｃａｔｉｏｎ。Ｎａｎｊｉｎｇ２１０００３－Ｃｈｉｎａ）
Ｋｅｙｗｏｒｄｓ：ｍｕｌｔｉｆｒａｃｔａｌａｎａｌｙｓｉｓ；ｅｌｅｃｔｒｏｅｎｃｅｐｈａｌｏｇｒａｍ；ｗａｖｅｌｅｔｔｒａｎｓｆｏｒｍｍａｘｉｍｕｍｍｏｄｕｌｕｓ
许多领域的奇异性分析中有很好的应用。因此，将
引言
脑电图（ｅｌｅｃｔｒｏｅｎｃｅｐｈａｌｏｇｒａｍ，ＥＥＧ）作为客观反映大脑机能状态的一个重要方面，蕴涵着丰富的生理、心理和病理信息，通过对ＥＥＧ信号的分析和处理，可以在临床上对一些脑疾病的诊断和治疗起到辅助作用。近年来，脑电信号的非线性动力学参数，如关联维、李雅普诺夫指数和近似熵等，已成功应用于脑功能状态及病理变化等方面的研究，证明了ＥＥＧ信号的混沌性¨。２１。
小波用于对脑电信号的分析已有很多研究¨“ｏ，分形用于对脑电信号的分析也有很多文献报道”。６１，但是小波和分形结合用来对脑电的分析还比较少。分形和小波相结合的方法发挥了小波多尺度分析和分形标度不变性两个特点，该方法在
ｄｏｉ：１０．３９６９／ｊ．ｉｓｓｎ．０２５８．．８０２１．２０１０．０５．０２４
分形和小波相结合的方法应用于脑电图的研究，可更加准确而深入地了解脑电信号的内在特征。
小波变换极大模值多重分形谱的方法一经提出…，已经在很多领域得到应用ｍ１０１。本研究采用小波变换极大模值多重分形谱的方法，对健康青年人和老年人的脑电信号进行分析。研究发现，脑电信号具有奇异性质，脑电信号奇异谱的跨度会随着年龄的增大而减小。这可以用来检测目标个体是否出现脑衰退迹象，从而辅以进一步的临床干预措施。１
理论和方法
１．１多重分形谱
分形是对没有特征长度但具有一定意义的自相似性图形和结构的总称。多重分形是许多个单
收稿日期：２０１０－０６－１８，修回日期：２０１０－０８．１３・通讯作者。
Ｅ－ｍａｉｌ：ｗａｎｇｊ＠ｎｊｕｐｔ．ｅｄｕ．ｃｎ
万方数据
中国生物医学工程学报２９卷
一分形在空间上的相互缠结、镶嵌，是单一分形的推广，单一分形可以看做是多重分形的一种特例。多重分形与单一分形一样，也是自相似的，与尺度无关。由每个测度的奇异性强度指数ｄ和奇异谱ｆＣｑ）构成了描述多重分形局部特性的奇异谱，八ｑ）也表示奇异性强度指数ａ所对应测度的分形维数。多重分形测度包括不同的奇异性强度，八ｑ）通常是分布于［ａ。；。，ａ…］之间的１个单峰图形，ａ…相对于较强的奇异性，ａ。；。相对于较弱的奇异性。计算多重分形谱可使用直接计算，（ｑ）的方法，即标准盒计数法¨川，这基于八ｑ）是一个特别测度的理论测度支架的维数。用尺寸为Ｌ的盒子来覆盖实验的测度，并且计算测度在每一个盒子中的概率，具体步骤如下：
１）定义归一化概率为
Ｐ；＝｝
ｒ．
（１）
∑ｔ
式中Ｚ是某测度被分成Ⅳ等份的时间序列中的第ｉ份。
那么，构建在盒子中的概率重整化测度的单参量簇可表示为
肛。（ｇ，￡）；半业（２）
耳
ｙ［Ｐ。（￡）］９
在一般维数的定义中，参量ｑ提供了对奇异测度不同区域的微观探测。对ｑ＞１，／．ｔ（ｑ）放大测度的强奇异区域；对ｑ＜１，弘（ｑ）强调测度的弱奇异区域；而对ｑ＝１，ｐ（ｑ）则复制原来的测度。
２）
计算ｐ（ｑ）理论测度支架的Ｈａｕｓｄｏｒｆｆ维数
以ｇ），有
以ｑ）＝一＿】｝ｉｍ赤．／ｚｐｉ（口，Ｌ）Ｉｎ［ｐ如，￡）］＝
ｌｉｍ墅堕笔型
（３）
３）计算／．ｔ（ｑ）奇异强度ａ；＝Ｉｎ（Ｐ。）／ｌｎＬ的平
均值ａ（ｇ），有
Ⅳ
ｌ
ａ（ｇ）＝一恕赤萎肛ｒ（ｇ，Ｌ）ｌ“Ｐｒ（Ｌ）］＝
躲巫拦掣ｎ型
￡＿＋０
ｌＬ
（４）
对应不同ｑ，求出Ｒ－ｆ（ａ）的单峰图像。１．２小波模极大值方法的多重分形谱
小波变换模极大值定义：在某一尺度疋。下，如
万方数据
果存在一点（Ｘｏ，Ｙ０），使得型生鸳尘型
上有一个模极大值（过零）点。如果对Ｙ。的某一邻（Ｘｏ，Ｙｏ）是局部极值点，且掣
＝０，则称点
ｄｙ
＝０在Ｙ＝Ｙｏ
域内的任意点Ｙ，有ｗｆ（茗。，Ｙ。）≤盯（茗。，Ｙ。），则称（茗。，Ｙ。）为小波变换模极大值点。小波变换的极大值的位置和大小包含了信号的主要信息川。
考虑用小波模极大值的方法求多重分形谱¨。。因为用小波模极大方法计算的多重分形谱能全面地反映信号的奇异性分布，而且，选择适当的小波函数还可以消除信号中平滑部分的影响，因此可以相信这种方法应该优于直接对信号求多重分形谱。小波模极大值法的多重分形谱计算过程如下：
１）对每一个尺度Ｓ计算出形，（Ⅱ，Ｓ）及其模极
大值，随．ｓ的变化连接小波模极大值，得出模极大值
序列。
２）然后用标准盒计数法对其模极大值序列进行分析，得出其多重分形谱。２
实验与结果分析
使用Ｎｅｕｒｏｅｌｅｃｔｒｉｃ
ａｎｄＭｙｏｅｌｅｅｔｒｉｃＤａｔａｂａｓｅｓ¨“，
由于该数据库提供的总样本数限制，用了９位健康年轻人（年轻人脑电图见图１）以及９位健康老年人的脑电信号（老年人脑电信号图见图２）。从图１和图２两种人群的脑电信号图上可以看到，它们的幅度和形态都差不多，在时域上不能得到明显的区分。因此，需要寻找其在表面上还显现不出来的内在特征，寻找甄别不同脑功能状态的敏感特征参量。
截取长度相同的数据对其进行分析，为了视图
图１一位的年轻人的脑电信号图
Ｆｉｇ．１
ＡｎＥＥＧｔａｋｅｎｆｒｏｍ
ｏｎｅ
ｙｏｕｎｇｓｕｂｊｅｃｔ
５期
赵大庆等：小波多重分形在脑电信号分析中的应用
图２一位老年人的脑电信号图
Ｆｉｇ．２
ＡｎＥＥＧｔａｋｅｎｆｒｏｍ
ｏｎｅ
ｏｌｄ
ｓｕｂｊｅｃｔ
清晰，只把其中３位年轻人和３位老年人的多重分形谱画出来（见图３）。从图３中可以看出，正常青年、老年ＥＥＧ信号的奇异谱曲线都呈单峰拱形，这正是多分形信号的一个重要特征，说明正常青年、老年ＥＥＧ信号都是多分形信号，都具有奇异性。
图３
３位青年人和３位老年人脑电信号多重分形谱对比Ｆｉｇ．３ＭｕｌｔｉｆｒａｃｔａｌｓｐｅｃｔｒｕｍｏｆＥＥＧｔａｋｅｎｆｒｏｍ３ｙｏｕｎｇ
ｓｕｂｊｅｃｔｓａｎｄ３ｏｌｄｓｕｂｊｅｃｔｓ
表１和表２分别列出９个年轻人和９个老年人脑电信号的小波模极大多重分形谱奇异强度最大值、最小值和奇异强度区间大小详细数值，并且计算了两组研究人群的奇异强度的均值和标准差。
根据多重分形理论，每种尺度上都有最强波动和最弱波动，它们的标度指数ａ岫和ａ一分别是ｏｔ（ｇ）中ｇ＿＋∞和ｇ－＋一∞的结果。根据分形谱的宽度公式△ａ＝ａ…一ａ嘶，分别计算出两种信号的多重分形谱分布范围值，并对其进行样本均数比较的ｔ检验，得到ｔ＝２．３４４３。查表可得，ｔｏ‘０５（９）＝２．２２８。因为ｔ＝２．３４４３＞ｔｏ．０５（９），所以Ｐ＜０．０５差异有统计学意义。可以认为，年轻人ＥＥＧ信号的
万方数据
△ａ平均值明显高于老年人ＥＥＧ信号的△口平均值，从而可以区分开来这两种信号，因此ｄｏｔ可以作为区分和识别这两种信号的一个标准。
表１年轻人脑电信号多重分形谱
Ｔａｂ．１
Ｔｈｅｍｕｌｔｉｆｒａｃｔａｌｓｐｅｃｔｒｕｍｏｆｙｏｕｎｇｓｉｇｎａｌｓ
表２老年人脑电信号多重分形谱
Ｔａｂ．２
Ｔｈｅｍｕｌｔｉｆｒａｃｔａｌｓｐｅｃｔｒｕｍｏｆ
ｏｌｄｅｒｓｉｇｎａｌｓ
３讨论和结论
由于大脑生理的复杂性，大脑功能体现出复杂的非线性动力学特性。本研究基于脑电信号的奇异特性，以及分形和小波相结合的方法，针对小波多尺度分析和分形标度不变性两个特点，用小波变换模极大值的多重分形方法，对健康青年人和健康老年人的ＥＥＧ信号进行了分析和比较。研究发现，青年、老年人群的ＥＥＧ信号的奇异度多重分形趋向于单重分形，预示着随着年龄的增长，大脑功能内在体现出的系统复杂度在趋于衰退。
随着社会老龄化的加剧，老年人群中患老年痴呆的比例在不断增大，给社会和家庭带来新的挑战
７９６
中国生物医学工程学报
２９卷
和难题，及早地预测和定量大脑的衰退程度显得尤为重要。一般来说，脑电信号奇异谱的跨度会随着年龄的增大而减小，因此医生可以用目标个体的脑电信号的奇异度区间，与同年龄段大部分人的脑电信号的奇异度区间统计平均值正常波动范围进行对比，预先检测出目标个体是否出现脑衰退迹象，从而辅以进一步的临床干预措施，这对提高高危老年痴呆人群的生存质鼍具有重要的意义。
下一步的工作将致力于增大每个目标人群样本量的搜集，在大样本量的基础上，给出目标人群的小波模极大多重分形谱区间的变化范围。目前，对于脑衰退程度的定量研究还没有被充分重视，因此，深入进行这方面的研究，制定出科学的标准，对于诊断及临床干预和治疗脑衰退具有重要的意义。
参考文献
［１］
ＡｃｈａｒｙａＵＲ，ＦａｕｓｔＯ。ＫａｎｎａｔｈａｌＮ，ｅｔａ１．Ｎｏｎ－ｌｉｎｅａｒａｎａｌｙｓｉｓｏｆ
ＥＥＧ
ｓｉｇｎａｌｓ
ａｔ
ｖａｒｉｏｕｓ
ｓｌｅｅｐ
ｓｔａｇｅｓ［Ｊ］．Ｃｏｍｐｕｔ
Ｍｅｔｈｏｄｓ
Ｐｒｏｇｒａｍｓ
Ｂｉｏｍｅｄ．２００５，舯（１）：３７—４５．
［２］江潮晖，冯焕清，刘大路，等．睡眠脑电的关联维数和近似熵分析［Ｊ］．生物医学工程学杂志．２００５，２２（４）：６４９—６５３．
［３］徐宝国，宋爱国．基于小波包变换和聚类分析的脑电信号识别方法［Ｊ］．仪器仪表学报，２００９，３０（１）：２５—２８．
［４】吴捷，张宁，杨卓，等．小波相干分析及其在听觉与震动刺激事件相关诱发脑电处理中的应用［Ｊ］．生物物理学报，２００７，
万方数据
２３（６）：４８２—４８７．
［５］
Ｐｏｐｉｖａｎｏｖ
Ｄ，ＪｉｖｋｏｖａＳ，ＳｔｏｍｏｎｙａｋｏｖＶ，眦ａ１．Ｅｆｆｅｃｔ
ｏｆ
ｉｎｄｅｐｅｎｄｅｎｔｃｏｍｐｏｎｅｎｔ
ａｎａｌｙｓｉｓ
ｏｎ
ｍｕｈｉｆｒａｃｔａｌｉｔｙｏｆＥＥＧｄｕｒｉｎｇ
ｖｉｓｕａｌ．ｍｏｔｏｒｔａｓｋ［Ｊ］．ＳｉｇｎａｌＰｒｏｃｅｓｓｉｎｇ，２００５，８５（１１）：２１１２
—２１２３．
［６］ＷａｎｇＷｅｉ，ＮｉｎｇＢａｏ，Ｗａｎｇ
Ｊｕｎ，ｅｔａ１．Ｉｎｔｅｒｌｅａｖｉｎｇｄｉｓｔｒｉｂｕｔｉｏｎｏｆ
ｍｕｌｔｉｆｒａｃｔａｌ
ｓｔｒｅｎｇｔｈ
ｏｆ
１６－ｃｈａｎｎｅｌ
ＥＥＧ
ｓｉｇｎａｌｓ［Ｊ］．Ｃｈｉｎｅｓｅ
ＳｃｉｅｎｃｅＢｕｌｌｅｔｉｎ４８，２００３，４８（１６），１７００—１７０３．［７］Ｍｕｚｙ
ＪＦ，Ｂａｃｒｙ
Ｅ，Ａｍｅｏｄｏ
Ａ．Ｍｕｈｉｆｒａｃｔａｌｆｏｒｍａｌｉｓｍ
ｆｏｒｆｒａｃｔａｌ
ｓｉｇｎａｌｓ：Ｔｈｅｓｔｒｕｃｔｕｒｅ—ｆｕｎｃｔｉｏｎａｐｐｒｏａｃｈ
ｖｅｒｓｕｓ
ｔｈｅｗａｖｅｌｅｔ—
ｔｒａｎｓｆｏｒｍ
ｍｏｄｕｌｕｓ－ｍａｘｉｍａ
ｍｅｔｈｏｄ［Ｊ］．Ｐｈｙｓｉｃａｌ
Ｒｅｖｉｅｗ
Ｅ，
１９９３。４７（２）：８７５—８８４．
［８］ＫｅｓｔｅｎｅｒＰ，Ａｒｕｅｏｄｏ
Ａ．Ｇｅｎｅｒａｌｉｚｉｎｇｔｈｅｗａｖｅｌｅｔ－ｂａｓｅｄ
ｍｕｌｔｉｆｒａｃｔａｌ
ｆｏｒｍａｌｉｓｍ
ｔｏ
ｒａｎｄｏｍ
ｖｅｃｔｏｒ
ｆｉｅｌｄｓ：ａｐｐｌｉｃａｔｉｏｎｔｏ
ｔｈｒｅｅ—ｄｉｍｅｎｓｉｏｎａｌｔｕｒｂｕｌｅｎｃｅｖｅｌｏｃｉｔｙ
ａｎｄ
ｖｏｒｔｉｃｉｔｙ
ｄａｔａ［Ｊ］．
ＰｈｙｓｉｃａｌＲｅｖｉｅｗ
Ｌｅｔｔｅｒｓ，２００４，９３（４）：０４４５０１—１—０４４５０１—４．
［９］苟学强。张义军，董万胜，等．基于小波的地闪首次回击辐射场的多重分形分析［Ｊ］．地球物理学报，２００７，５０（１）：１０１
—１０５．
［１０］芩为，杨世峰，薛蓉，等．基于小波变换模极大法的聚乙烯催化剂表面分形分析［Ｊ］．中国科学Ｂ辑，２００７，３７（４）：３７４
—３７９．
［１１］
ＣｈｈａｂｒａＡＢ，ＭｃｎｅｖｅａｕＣ，Ｊｅｎｓｅｎ
ＲＶ．Ｄｉｒｅｃｔ
ｏｆｔｈｅ
ｆ（ａ）
ｓｉｎｇｕｌａｒｉｔｙ
ｓｐｅｃｔｒｕｍ
ａｎｄｉｔｓ
ａｐｐｌｉｃａｔｉｏｎｔＯ
ｆｕｌｌｙ
ｄｅｖｅｌｏｐｅｄ
ｔｕｒｂｕｌｅｎｃｅ［Ｊ］．Ｐｈｙｓ
Ｒｅｖ
Ａ。１９８９，４０（９）：５２８４—５２９４．
［１２］
Ｎａｔｉｏｎａｌ
Ｉｎｓｔｉｔｕｔｅｓｏｆ
Ｈｅａｌｔｈ．ＰｈｙｓｉｏＮｅｔ［ＤＢ／ＯＬ］．ｈｔｔｐ：／／ｗｗｗ．
ｐｈｙｓｉｏｎｅｔ．ｏｒｇ／ｐｈｙｓｉｏｂａｎｋ／ｄａｔａｂａｓｅ／ｎｅｓｆｄｂ／。２００７－０９旬１／２０１０—０６．１０．
小波多重分形在脑电信号分析中的应用
作者：作者单位：
赵大庆， 王俊， ZHAO Da-Qing， WANG Jun
赵大庆,ZHAO Da-Qing(南京邮电大学,图像处理与图像通信江苏省重点实验室,通信与信息工程学院,南京,210003)， 王俊,WANG Jun(南京邮电大学,图像处理与图像通信江苏省重点实验室,地理与生物信息学院,南京,210003)中国生物医学工程学报
CHINESE JOURNAL OF BIOMEDICAL ENGINEERING2010,29(5)
刊名：英文刊名：年，卷(期)：
参考文献(12条)
1.Acharya UR;Faust O;Kannathal N Non-linear analysis of EEG signals at various sleep stages[外文期刊] 2005(01)
2.江潮晖;冯焕清;刘大路 睡眠脑电的关联维数和近似熵分析[期刊论文]-生物医学工程学杂志 2005(04)3.徐宝国;宋爱国 基于小波包变换和聚类分析的脑电信号识别方法[期刊论文]-仪器仪表学报 2009(01)
4.吴捷;张宁;杨卓 小波相干分析及其在听觉与震动刺激事件相关诱发脑电处理中的应用[期刊论文]-生物物理学报 2007(06)
5.Popivanov D;Jivkova S;Stomonyakov V Effect of independent component analysis on multifractality ofEEG during visual-motor task[外文期刊] 2005(11)
6.Wang Wei;Ning Bao;Wang Jun Interleaving distribution of multifractal strength of 16-channel EEGsignals[期刊论文]-Chinese Science Bulletin 48 2003(16)
7.Muzy JF;Bacry E;Arneodo A Multifraetal formalism for fractal signals:The structure-functionapproach versus the wavelettransform modulus-maxima method 1993(02)
8.Kestener P;Arueodo A Generalizing the wavelet-based multifractal formalism to random vectorfields:application to three-dimensional turbulence velocity and vorticity data[外文期刊] 2004(04)9.苟学强;张义军;董万胜 基于小波的地闪首次回击辐射场的多重分形分析[期刊论文]-地球物理学报 2007(01)10.芩为;杨世峰;薛蓉 基于小波变换模极大法的聚乙烯催化剂表面分形分析[期刊论文]-中国科学B辑 2007(04)11.Chhabra AB;Meneveau C;Jensen RV Direct of the f(a) singularity spectrum and its application tofully developed turbulence 1989(09)
12.National Institutes of Health PhysioNet 2010
本文链接：http://d.wanfangdata.com.cn/Periodical_zgswyxgcxb201005024.aspx
２９卷５期２０１０年ｌＯ月
中国生物医学工程学报ｏｆＢｉｏｍｅｄｉｃａｌ
Ｅｎｇｉｎｅｅｒｉｎｇ
Ｖ０１．２９Ｏｅｔｏｂｅｒ
Ｎｏ．５２０１０
ＣｈｉｎｅｓｅＪｏｕｒｎａｌ
小波多重分形在脑电信号分析中的应用
赵大庆１
王
俊２‘
１（南京邮电大学，图像处理与图像通信江苏省重点实验室，通信与信息工程学院，南京２１０００３）２（南京邮电大学，图像处理与图像通信江苏省重点实验室，地理与生物信息学院，南京２１０００３）
摘要：脑电图（ＥＥＧ）作为客观反映大脑机能状态的一个重要方面，蕴涵着丰富的生理、心理和病理信息。分形和小波相结合的方法，发挥了小波多尺度分析和分形标度不变性的两个特点。利用小波模极大值多重分形的方法，对青老年人群的脑电图进行分析。研究发现，脑电信号具有奇异性质，脑电信号奇异谱的跨度会随着年龄的增大而减小，可以用来检测目标个体是否出现脑衰退迹象，从而辅助进一步的临床干预措施。关键词：多重分形分析；脑电图；小波模极大值
中图分类号Ｒ３１８．５
文献标识码Ｄ
文章编号０２５８—８０２１（２０１０）０５－０７９３－０４
Ａｐｐｌｉｃａｔｉｏｎｏｆ
ＷａｖｅｌｅｔＭｕｌｔｉｆｒａｃｔａｌ
Ｓｐｅｃｔｒｕｍｔｏｔｈｅ
ＡｎａｌｙｓｉｓｏｆＥｌｅｃｔｒｏｅｎｃｅｐｈａｌｏｇｒａｍＳｉｇｎａｌｓ
ＺＨＡＯ
１
Ｄａ．Ｑｉｎ９１
ＷＡＮＧＪｕｎ２’
ｏｆＴｅｌｅｃｏｍｍ．＆Ｉｎｆｏｒｍａｔｉｏｎ
ｌＩｍａｇｅＰｒｏｃｅｓｓｉｎｇａｎｄＩｍａｇｅＣｏｍｍｕｎｉｃａｔｉｏｎＫｅｙＬａｂ。Ｃｏｌｌｅｇｅ
Ｅｎｇｉｎｅｅｒｉｎｇ。Ｎ８咏ｎｇＵｎｉｔｅｒｓ幼ｏｆＰｏｓｔｓ＆Ｔｅｌｅｃｏｍｍｕｎｉｃａｔｉｏｎ。ＮⅡｎｊｉｎｇ２１０００３。Ｃｈｉｎａ、
２（Ｉｍａｇｅ
Ｐｒｏｃｅｓｓｉｎｇ
ａｎｄＩｍａｇｅ
ＣｏｍｍｕｎｉｃａｔｉｏｎＫｅｙＬａｂ。Ｃｏｌｌｅｇｅｏｆ
Ｇｅｏｇｒａｐｈｉｃａｌ＆Ｂｉｏｌｏｇｉｃａｌ
Ｉｎｆｏｒｍａｔｉｏｎ。
ＮａｎｊｉｎｇＵｎｉｔｅｒｓｉｔｙｏｆＰｏｓｔｓ＆Ｔｅｌｅｃｏｍｍｕｎｉｃａｔｉｏｎ。Ｎａｎｊｉｎｇ２１０００３－Ｃｈｉｎａ）
Ｋｅｙｗｏｒｄｓ：ｍｕｌｔｉｆｒａｃｔａｌａｎａｌｙｓｉｓ；ｅｌｅｃｔｒｏｅｎｃｅｐｈａｌｏｇｒａｍ；ｗａｖｅｌｅｔｔｒａｎｓｆｏｒｍｍａｘｉｍｕｍｍｏｄｕｌｕｓ
许多领域的奇异性分析中有很好的应用。因此，将
引言
脑电图（ｅｌｅｃｔｒｏｅｎｃｅｐｈａｌｏｇｒａｍ，ＥＥＧ）作为客观反映大脑机能状态的一个重要方面，蕴涵着丰富的生理、心理和病理信息，通过对ＥＥＧ信号的分析和处理，可以在临床上对一些脑疾病的诊断和治疗起到辅助作用。近年来，脑电信号的非线性动力学参数，如关联维、李雅普诺夫指数和近似熵等，已成功应用于脑功能状态及病理变化等方面的研究，证明了ＥＥＧ信号的混沌性¨。２１。
小波用于对脑电信号的分析已有很多研究¨“ｏ，分形用于对脑电信号的分析也有很多文献报道”。６１，但是小波和分形结合用来对脑电的分析还比较少。分形和小波相结合的方法发挥了小波多尺度分析和分形标度不变性两个特点，该方法在
ｄｏｉ：１０．３９６９／ｊ．ｉｓｓｎ．０２５８．．８０２１．２０１０．０５．０２４
分形和小波相结合的方法应用于脑电图的研究，可更加准确而深入地了解脑电信号的内在特征。
小波变换极大模值多重分形谱的方法一经提出…，已经在很多领域得到应用ｍ１０１。本研究采用小波变换极大模值多重分形谱的方法，对健康青年人和老年人的脑电信号进行分析。研究发现，脑电信号具有奇异性质，脑电信号奇异谱的跨度会随着年龄的增大而减小。这可以用来检测目标个体是否出现脑衰退迹象，从而辅以进一步的临床干预措施。１
理论和方法
１．１多重分形谱
分形是对没有特征长度但具有一定意义的自相似性图形和结构的总称。多重分形是许多个单
收稿日期：２０１０－０６－１８，修回日期：２０１０－０８．１３・通讯作者。
Ｅ－ｍａｉｌ：ｗａｎｇｊ＠ｎｊｕｐｔ．ｅｄｕ．ｃｎ
万方数据
中国生物医学工程学报２９卷
一分形在空间上的相互缠结、镶嵌，是单一分形的推广，单一分形可以看做是多重分形的一种特例。多重分形与单一分形一样，也是自相似的，与尺度无关。由每个测度的奇异性强度指数ｄ和奇异谱ｆＣｑ）构成了描述多重分形局部特性的奇异谱，八ｑ）也表示奇异性强度指数ａ所对应测度的分形维数。多重分形测度包括不同的奇异性强度，八ｑ）通常是分布于［ａ。；。，ａ…］之间的１个单峰图形，ａ…相对于较强的奇异性，ａ。；。相对于较弱的奇异性。计算多重分形谱可使用直接计算，（ｑ）的方法，即标准盒计数法¨川，这基于八ｑ）是一个特别测度的理论测度支架的维数。用尺寸为Ｌ的盒子来覆盖实验的测度，并且计算测度在每一个盒子中的概率，具体步骤如下：
１）定义归一化概率为
Ｐ；＝｝
ｒ．
（１）
∑ｔ
式中Ｚ是某测度被分成Ⅳ等份的时间序列中的第ｉ份。
那么，构建在盒子中的概率重整化测度的单参量簇可表示为
肛。（ｇ，￡）；半业（２）
耳
ｙ［Ｐ。（￡）］９
在一般维数的定义中，参量ｑ提供了对奇异测度不同区域的微观探测。对ｑ＞１，／．ｔ（ｑ）放大测度的强奇异区域；对ｑ＜１，弘（ｑ）强调测度的弱奇异区域；而对ｑ＝１，ｐ（ｑ）则复制原来的测度。
２）
计算ｐ（ｑ）理论测度支架的Ｈａｕｓｄｏｒｆｆ维数
以ｇ），有
以ｑ）＝一＿】｝ｉｍ赤．／ｚｐｉ（口，Ｌ）Ｉｎ［ｐ如，￡）］＝
ｌｉｍ墅堕笔型
（３）
３）计算／．ｔ（ｑ）奇异强度ａ；＝Ｉｎ（Ｐ。）／ｌｎＬ的平
均值ａ（ｇ），有
Ⅳ
ｌ
ａ（ｇ）＝一恕赤萎肛ｒ（ｇ，Ｌ）ｌ“Ｐｒ（Ｌ）］＝
躲巫拦掣ｎ型
￡＿＋０
ｌＬ
（４）
对应不同ｑ，求出Ｒ－ｆ（ａ）的单峰图像。１．２小波模极大值方法的多重分形谱
小波变换模极大值定义：在某一尺度疋。下，如
万方数据
果存在一点（Ｘｏ，Ｙ０），使得型生鸳尘型
上有一个模极大值（过零）点。如果对Ｙ。的某一邻（Ｘｏ，Ｙｏ）是局部极值点，且掣
＝０，则称点
ｄｙ
＝０在Ｙ＝Ｙｏ
域内的任意点Ｙ，有ｗｆ（茗。，Ｙ。）≤盯（茗。，Ｙ。），则称（茗。，Ｙ。）为小波变换模极大值点。小波变换的极大值的位置和大小包含了信号的主要信息川。
考虑用小波模极大值的方法求多重分形谱¨。。因为用小波模极大方法计算的多重分形谱能全面地反映信号的奇异性分布，而且，选择适当的小波函数还可以消除信号中平滑部分的影响，因此可以相信这种方法应该优于直接对信号求多重分形谱。小波模极大值法的多重分形谱计算过程如下：
１）对每一个尺度Ｓ计算出形，（Ⅱ，Ｓ）及其模极
大值，随．ｓ的变化连接小波模极大值，得出模极大值
序列。
２）然后用标准盒计数法对其模极大值序列进行分析，得出其多重分形谱。２
实验与结果分析
使用Ｎｅｕｒｏｅｌｅｃｔｒｉｃ
ａｎｄＭｙｏｅｌｅｅｔｒｉｃＤａｔａｂａｓｅｓ¨“，
由于该数据库提供的总样本数限制，用了９位健康年轻人（年轻人脑电图见图１）以及９位健康老年人的脑电信号（老年人脑电信号图见图２）。从图１和图２两种人群的脑电信号图上可以看到，它们的幅度和形态都差不多，在时域上不能得到明显的区分。因此，需要寻找其在表面上还显现不出来的内在特征，寻找甄别不同脑功能状态的敏感特征参量。
截取长度相同的数据对其进行分析，为了视图
图１一位的年轻人的脑电信号图
Ｆｉｇ．１
ＡｎＥＥＧｔａｋｅｎｆｒｏｍ
ｏｎｅ
ｙｏｕｎｇｓｕｂｊｅｃｔ
５期
赵大庆等：小波多重分形在脑电信号分析中的应用
图２一位老年人的脑电信号图
Ｆｉｇ．２
ＡｎＥＥＧｔａｋｅｎｆｒｏｍ
ｏｎｅ
ｏｌｄ
ｓｕｂｊｅｃｔ
清晰，只把其中３位年轻人和３位老年人的多重分形谱画出来（见图３）。从图３中可以看出，正常青年、老年ＥＥＧ信号的奇异谱曲线都呈单峰拱形，这正是多分形信号的一个重要特征，说明正常青年、老年ＥＥＧ信号都是多分形信号，都具有奇异性。
图３
３位青年人和３位老年人脑电信号多重分形谱对比Ｆｉｇ．３ＭｕｌｔｉｆｒａｃｔａｌｓｐｅｃｔｒｕｍｏｆＥＥＧｔａｋｅｎｆｒｏｍ３ｙｏｕｎｇ
ｓｕｂｊｅｃｔｓａｎｄ３ｏｌｄｓｕｂｊｅｃｔｓ
表１和表２分别列出９个年轻人和９个老年人脑电信号的小波模极大多重分形谱奇异强度最大值、最小值和奇异强度区间大小详细数值，并且计算了两组研究人群的奇异强度的均值和标准差。
根据多重分形理论，每种尺度上都有最强波动和最弱波动，它们的标度指数ａ岫和ａ一分别是ｏｔ（ｇ）中ｇ＿＋∞和ｇ－＋一∞的结果。根据分形谱的宽度公式△ａ＝ａ…一ａ嘶，分别计算出两种信号的多重分形谱分布范围值，并对其进行样本均数比较的ｔ检验，得到ｔ＝２．３４４３。查表可得，ｔｏ‘０５（９）＝２．２２８。因为ｔ＝２．３４４３＞ｔｏ．０５（９），所以Ｐ＜０．０５差异有统计学意义。可以认为，年轻人ＥＥＧ信号的
万方数据
△ａ平均值明显高于老年人ＥＥＧ信号的△口平均值，从而可以区分开来这两种信号，因此ｄｏｔ可以作为区分和识别这两种信号的一个标准。
表１年轻人脑电信号多重分形谱
Ｔａｂ．１
Ｔｈｅｍｕｌｔｉｆｒａｃｔａｌｓｐｅｃｔｒｕｍｏｆｙｏｕｎｇｓｉｇｎａｌｓ
表２老年人脑电信号多重分形谱
Ｔａｂ．２
Ｔｈｅｍｕｌｔｉｆｒａｃｔａｌｓｐｅｃｔｒｕｍｏｆ
ｏｌｄｅｒｓｉｇｎａｌｓ
３讨论和结论
由于大脑生理的复杂性，大脑功能体现出复杂的非线性动力学特性。本研究基于脑电信号的奇异特性，以及分形和小波相结合的方法，针对小波多尺度分析和分形标度不变性两个特点，用小波变换模极大值的多重分形方法，对健康青年人和健康老年人的ＥＥＧ信号进行了分析和比较。研究发现，青年、老年人群的ＥＥＧ信号的奇异度多重分形趋向于单重分形，预示着随着年龄的增长，大脑功能内在体现出的系统复杂度在趋于衰退。
随着社会老龄化的加剧，老年人群中患老年痴呆的比例在不断增大，给社会和家庭带来新的挑战
７９６
中国生物医学工程学报
２９卷
和难题，及早地预测和定量大脑的衰退程度显得尤为重要。一般来说，脑电信号奇异谱的跨度会随着年龄的增大而减小，因此医生可以用目标个体的脑电信号的奇异度区间，与同年龄段大部分人的脑电信号的奇异度区间统计平均值正常波动范围进行对比，预先检测出目标个体是否出现脑衰退迹象，从而辅以进一步的临床干预措施，这对提高高危老年痴呆人群的生存质鼍具有重要的意义。
下一步的工作将致力于增大每个目标人群样本量的搜集，在大样本量的基础上，给出目标人群的小波模极大多重分形谱区间的变化范围。目前，对于脑衰退程度的定量研究还没有被充分重视，因此，深入进行这方面的研究，制定出科学的标准，对于诊断及临床干预和治疗脑衰退具有重要的意义。
参考文献
［１］
ＡｃｈａｒｙａＵＲ，ＦａｕｓｔＯ。ＫａｎｎａｔｈａｌＮ，ｅｔａ１．Ｎｏｎ－ｌｉｎｅａｒａｎａｌｙｓｉｓｏｆ
ＥＥＧ
ｓｉｇｎａｌｓ
ａｔ
ｖａｒｉｏｕｓ
ｓｌｅｅｐ
ｓｔａｇｅｓ［Ｊ］．Ｃｏｍｐｕｔ
Ｍｅｔｈｏｄｓ
Ｐｒｏｇｒａｍｓ
Ｂｉｏｍｅｄ．２００５，舯（１）：３７—４５．
［２］江潮晖，冯焕清，刘大路，等．睡眠脑电的关联维数和近似熵分析［Ｊ］．生物医学工程学杂志．２００５，２２（４）：６４９—６５３．
［３］徐宝国，宋爱国．基于小波包变换和聚类分析的脑电信号识别方法［Ｊ］．仪器仪表学报，２００９，３０（１）：２５—２８．
［４】吴捷，张宁，杨卓，等．小波相干分析及其在听觉与震动刺激事件相关诱发脑电处理中的应用［Ｊ］．生物物理学报，２００７，
万方数据
２３（６）：４８２—４８７．
［５］
Ｐｏｐｉｖａｎｏｖ
Ｄ，ＪｉｖｋｏｖａＳ，ＳｔｏｍｏｎｙａｋｏｖＶ，眦ａ１．Ｅｆｆｅｃｔ
ｏｆ
ｉｎｄｅｐｅｎｄｅｎｔｃｏｍｐｏｎｅｎｔ
ａｎａｌｙｓｉｓ
ｏｎ
ｍｕｈｉｆｒａｃｔａｌｉｔｙｏｆＥＥＧｄｕｒｉｎｇ
ｖｉｓｕａｌ．ｍｏｔｏｒｔａｓｋ［Ｊ］．ＳｉｇｎａｌＰｒｏｃｅｓｓｉｎｇ，２００５，８５（１１）：２１１２
—２１２３．
［６］ＷａｎｇＷｅｉ，ＮｉｎｇＢａｏ，Ｗａｎｇ
Ｊｕｎ，ｅｔａ１．Ｉｎｔｅｒｌｅａｖｉｎｇｄｉｓｔｒｉｂｕｔｉｏｎｏｆ
ｍｕｌｔｉｆｒａｃｔａｌ
ｓｔｒｅｎｇｔｈ
ｏｆ
１６－ｃｈａｎｎｅｌ
ＥＥＧ
ｓｉｇｎａｌｓ［Ｊ］．Ｃｈｉｎｅｓｅ
ＳｃｉｅｎｃｅＢｕｌｌｅｔｉｎ４８，２００３，４８（１６），１７００—１７０３．［７］Ｍｕｚｙ
ＪＦ，Ｂａｃｒｙ
Ｅ，Ａｍｅｏｄｏ
Ａ．Ｍｕｈｉｆｒａｃｔａｌｆｏｒｍａｌｉｓｍ
ｆｏｒｆｒａｃｔａｌ
ｓｉｇｎａｌｓ：Ｔｈｅｓｔｒｕｃｔｕｒｅ—ｆｕｎｃｔｉｏｎａｐｐｒｏａｃｈ
ｖｅｒｓｕｓ
ｔｈｅｗａｖｅｌｅｔ—
ｔｒａｎｓｆｏｒｍ
ｍｏｄｕｌｕｓ－ｍａｘｉｍａ
ｍｅｔｈｏｄ［Ｊ］．Ｐｈｙｓｉｃａｌ
Ｒｅｖｉｅｗ
Ｅ，
１９９３。４７（２）：８７５—８８４．
［８］ＫｅｓｔｅｎｅｒＰ，Ａｒｕｅｏｄｏ
Ａ．Ｇｅｎｅｒａｌｉｚｉｎｇｔｈｅｗａｖｅｌｅｔ－ｂａｓｅｄ
ｍｕｌｔｉｆｒａｃｔａｌ
ｆｏｒｍａｌｉｓｍ
ｔｏ
ｒａｎｄｏｍ
ｖｅｃｔｏｒ
ｆｉｅｌｄｓ：ａｐｐｌｉｃａｔｉｏｎｔｏ
ｔｈｒｅｅ—ｄｉｍｅｎｓｉｏｎａｌｔｕｒｂｕｌｅｎｃｅｖｅｌｏｃｉｔｙ
ａｎｄ
ｖｏｒｔｉｃｉｔｙ
ｄａｔａ［Ｊ］．
ＰｈｙｓｉｃａｌＲｅｖｉｅｗ
Ｌｅｔｔｅｒｓ，２００４，９３（４）：０４４５０１—１—０４４５０１—４．
［９］苟学强。张义军，董万胜，等．基于小波的地闪首次回击辐射场的多重分形分析［Ｊ］．地球物理学报，２００７，５０（１）：１０１
—１０５．
［１０］芩为，杨世峰，薛蓉，等．基于小波变换模极大法的聚乙烯催化剂表面分形分析［Ｊ］．中国科学Ｂ辑，２００７，３７（４）：３７４
—３７９．
［１１］
ＣｈｈａｂｒａＡＢ，ＭｃｎｅｖｅａｕＣ，Ｊｅｎｓｅｎ
ＲＶ．Ｄｉｒｅｃｔ
ｏｆｔｈｅ
ｆ（ａ）
ｓｉｎｇｕｌａｒｉｔｙ
ｓｐｅｃｔｒｕｍ
ａｎｄｉｔｓ
ａｐｐｌｉｃａｔｉｏｎｔＯ
ｆｕｌｌｙ
ｄｅｖｅｌｏｐｅｄ
ｔｕｒｂｕｌｅｎｃｅ［Ｊ］．Ｐｈｙｓ
Ｒｅｖ
Ａ。１９８９，４０（９）：５２８４—５２９４．
［１２］
Ｎａｔｉｏｎａｌ
Ｉｎｓｔｉｔｕｔｅｓｏｆ
Ｈｅａｌｔｈ．ＰｈｙｓｉｏＮｅｔ［ＤＢ／ＯＬ］．ｈｔｔｐ：／／ｗｗｗ．
ｐｈｙｓｉｏｎｅｔ．ｏｒｇ／ｐｈｙｓｉｏｂａｎｋ／ｄａｔａｂａｓｅ／ｎｅｓｆｄｂ／。２００７－０９旬１／２０１０—０６．１０．
小波多重分形在脑电信号分析中的应用
作者：作者单位：
赵大庆， 王俊， ZHAO Da-Qing， WANG Jun
赵大庆,ZHAO Da-Qing(南京邮电大学,图像处理与图像通信江苏省重点实验室,通信与信息工程学院,南京,210003)， 王俊,WANG Jun(南京邮电大学,图像处理与图像通信江苏省重点实验室,地理与生物信息学院,南京,210003)中国生物医学工程学报
CHINESE JOURNAL OF BIOMEDICAL ENGINEERING2010,29(5)
刊名：英文刊名：年，卷(期)：
参考文献(12条)
1.Acharya UR;Faust O;Kannathal N Non-linear analysis of EEG signals at various sleep stages[外文期刊] 2005(01)
2.江潮晖;冯焕清;刘大路 睡眠脑电的关联维数和近似熵分析[期刊论文]-生物医学工程学杂志 2005(04)3.徐宝国;宋爱国 基于小波包变换和聚类分析的脑电信号识别方法[期刊论文]-仪器仪表学报 2009(01)
4.吴捷;张宁;杨卓 小波相干分析及其在听觉与震动刺激事件相关诱发脑电处理中的应用[期刊论文]-生物物理学报 2007(06)
5.Popivanov D;Jivkova S;Stomonyakov V Effect of independent component analysis on multifractality ofEEG during visual-motor task[外文期刊] 2005(11)
6.Wang Wei;Ning Bao;Wang Jun Interleaving distribution of multifractal strength of 16-channel EEGsignals[期刊论文]-Chinese Science Bulletin 48 2003(16)
7.Muzy JF;Bacry E;Arneodo A Multifraetal formalism for fractal signals:The structure-functionapproach versus the wavelettransform modulus-maxima method 1993(02)
8.Kestener P;Arueodo A Generalizing the wavelet-based multifractal formalism to random vectorfields:application to three-dimensional turbulence velocity and vorticity data[外文期刊] 2004(04)9.苟学强;张义军;董万胜 基于小波的地闪首次回击辐射场的多重分形分析[期刊论文]-地球物理学报 2007(01)10.芩为;杨世峰;薛蓉 基于小波变换模极大法的聚乙烯催化剂表面分形分析[期刊论文]-中国科学B辑 2007(04)11.Chhabra AB;Meneveau C;Jensen RV Direct of the f(a) singularity spectrum and its application tofully developed turbulence 1989(09)
12.National Institutes of Health PhysioNet 2010
本文链接：http://d.wanfangdata.com.cn/Periodical_zgswyxgcxb201005024.aspx