改进C0+复杂度及其应用.pdfVIP

第47卷第6期 复旦学报(自然科学版) Vbl．47No．6 2008年12月 ofFudan Dec．2008 Joumal UniverSi够(Natl】ralSci印ce) 文章编号：0427—7104(2008)06一0791一06 改进的Co复杂度及其应用 蔡志杰，孙洁 (复旦大学数学科学学院，上海200433) 摘要：引入参数r，改进了Co复杂度，给出了改进的一维及二维c0复杂度的有关性质及若干实例，并应用于 脑电复杂性的分析中．实验证明改迸的co复杂度比原来的C0复杂度更具鲁棒性，实际应用时适应性更强． 关键词：复杂度；Co复杂度；二维复杂度；logistic映射；脑电信号 911．73 中图分类号：TN9儿．72；TN 文献标识码：A 许多实际信号(例如EEG信号、心电信号等)都具有强烈的非线性性，传统的线性分析方法无法反映 它们的本质特性．而大多数非线性分析方法需要长程数据才能得到鲁棒的结果，这又与这类信号所具有的 强烈的非平稳性相矛盾．因此非常需要寻找一种适合于分析同时具有非线性性和非平稳性的时间序列的 J．但是，大多数实际信号都是连 非线性分析方法．kmpel．Ziv复杂度(CLZ)是一种合适的非线性方法u-3 续的，丽G2只能定义在有限字符集上．因此，在计算Qz值之前需要对原始时间序列进行粗粒化处理，例 如将原始信号转换成二进制序列．对原始信号作这种过分的粗粒化可能会改变其动力学性质，从而有可能 得到错误的结论【4．5】．所以，在对这类信号进行分析时应遵循两条标准：1)对短程数据也能得到鲁棒的结 果；2)对连续时间序列不需要作过分的粗粒化处理，而是直接定义． C。复杂度【6-9J满足这两条标准．Co复杂度的定义如下：计算序列功率谱的均方值，高于均方值的部分 保持不变，而将低于均方值的部分置为零．计算新的信号频谱的Fourier逆变换，得到新的时间序列．新时 间序列与原时间序列之差的平方和与原始信号的平方和之比定义为原始信号的c0复杂度．文献[9]证明 了C0复杂度的重要性质，说明了C0复杂度可以作为描述时间序列随机程度的复杂性指标．但在实际应用 时，也发现了G复杂度的一些缺陷，如对零均值的时间序列，其计算结果与理论值不相符合．究其原因，是 因为时间序列频谱值的平方值与其期望值非常接近，实际计算时，无法预计它们之间的大小关系，从而难 以估计最终的Co复杂度值． 本文引入参数，一，改进了C。复杂度，使其在均值为零时也能得到鲁棒的、与理论值相吻合的结果．我 们给出了改进的一维和二维Co复杂度的主要性质，列举了若干实例验证改进的C。复杂度的合理性，并将 其运用于EEG信号的分析中，得到了满意的结果． 1 改进的Co复杂度 记{厂(忌)，惫=O，1，2，…，N—l}是一个长度为N的时间序列， 歹=o，1，2，…，N—l FN(歹)=考∑，(是)骶茸 表示其F0urier变换序列，其中wⅣ=eXp(等)，i=厅表示虚数单位． 设{FN(』)，歹=0，1，2，…，N—l}的均方值为 收稿日期：2008一03—14 基金项目：上海市自然科学基金资助项目(08ZRl401100) 作者简介：蔡志杰(1968一)，男，博士，副教授． 万方数据 复旦学报(自然科学版) 第47卷 GN=志姜J氏(川2， 引人参数，．，保留超过均方值r倍的频谱，而将其余部分置为零，即 瓦∽={≯o’耋；2暑f：三墨：筹： 其中r(1)是一个给定的正常数．对{FN(歹)，歹=O，1，2，…，N一1}作F。u