基于有限穿越可视图的时间序列网络模型

提出了一种改进的时间序列有限穿越可视图建网方法,并对三种可视图(可视图、水平可视图、有限穿越可视图)网络度分布进行了评价.结果表明:水平可视图网络均无法有效识别各类时间序列信号(周期、分形、混沌);对分形信号,可视图及有限穿越可视图网络均具有无标度幂律形式,但抗噪能力较差;对周期信号及混沌信号,有限穿越可视图网络比可视图具有更强的抗噪性.在此基础上,采用有限穿越可视图网络从油气水三相流电导波动信号中提取了度分布特征参数,通过其特征参数组合实现了对三种典型三相流流型(水包油泡状流、水包油泡状-段塞过渡流型及水包油段塞流)较好的辨识效果.     

一种新的心率变异性度量方法

心率变异性的复杂波动反映了心脏的自主调节功能.本文提出了一种新的心率变异性度量方法——ICBN方法,该方法通过改进的自适应噪声完备集合经验模态分解方法对心率变异性信号进行分解,得到多个模态分量,计算每个模态分量的bubble熵得到熵值向量,把该向量映射成复杂网络,通过计算网络的特征参数,对心率变异性在不同时频尺度状态下的非线性特征之间的耦合关系进行度量.首先,采用时域、频域和ICBN分析方法对29名充血性心力衰竭病人和29名正常窦性心律对象的心率变异性进行分析,结果表明:时域指标三角指数HRVTi,频域指标LF/HF,网络层级加权值WB,平均点权值PW,特征路径长度CL具有统计学差异;基于网络层级加权值WB,特征路径长度CL,频域指标LF/HF和Fisher判别方法的识别模型对充血性心力衰竭病人的识别正确率达到89.66%.然后,又对43名房颤心律失常患者和43名正常窦性心律对象的心率变异性进行分析,结果表明:时域指标SDNN,pNN50,RMSSD,频域指标LF/HF,网络层级加权值WB,平均点权值PW具有统计学差异;时域指标pNN50,RMSSD,频域指标LF/HF和网络层级加权值WB,平均点权值PW作为特征向量,Fisher判别方法作为分类器,对房颤心律失常患者的识别正确率达到91.86%.综合以上实验结果可知,本文为心率变异性的度量研究提供了一种新的思路.     

两相流有限穿越可视图演化动力学研究

针对小管径两相流流动特性, 全新优化设计弧形对壁式电导传感器. 通过动态实验在获取传感器测量信号的基础上, 采用有限穿越可视图理论构建对应于不同流型的两相流复杂网络. 通过分析发现, 有限穿越可视图网络异速生长指数和网络平均度值的联合分布可实现对小管径两相流的流型辨识; 有限穿越可视图度分布曲线峰值可有效刻画与泡径大小分布相关的流动物理结构细节特征; 网络平均度值可表征流动结构的宏观特性; 网络异速生长指数对流体动力学复杂性十分敏感, 可揭示不同流型演化过程中的细节演化动力学特性. 两相流测量信号的有限穿越可视图分析为揭示两相流流型的形成及演化动力学机理提供了新途径. 关键词： 两相流 复杂网络 有限穿越可视图 网络异速生长指数     

昼夜节律下心率变异性信号的熵信息和谱特征

生理系统产生的复杂波动信号能够反映其潜在的动力学特征.采用基本尺度熵和功率谱的方法分析24 h心率变异性信号.结果表明,心脏系统昼夜节律下生理和病理的变化伴随着变化的基本尺度熵和功率谱分布,但是对于近似熵,其变化却不明显;同时发现,基本尺度熵的变化能够反映相应的自主神经调控的变化,由于充血性心力衰竭患者迷走神经的调控被抑制,交感神经的调控占优势,所以数据中会出现更多变化的矢量模式组合,因此心力衰竭患者心率变异性信号的熵值较高;在夜间睡眠状态时,由于迷走神经的调控增强,交感神经的调控减少,所以健康人和心力衰竭患者的基本尺度熵都比白天清醒状态时产生了下降趋势.     

两种熵测度在量化射击运动员短时心率变异性信号复杂度上的一致性

利用两种基于熵的非线性复杂度测度：近似熵和样本熵,研究了专业射击运动员两种不同状态下(休息和练习赛)心率变异性信号的复杂度.计算结果表明：射击运动员休息时其心率变异性信号的熵值大于射击比赛时信号的熵值,这意味着运动员一旦进行射击比赛时,其心率变异性信号复杂度降低了,心跳变得更为规则了.为了更好地应用这两种基于熵的方法,进一步分析了算法中的两个重要影响因素：矢量匹配容差r和序列长度N对算法性能的影响.分析结果表明：只要参数选择在合适的范围内,近似熵和样本熵都能够正确地区分出两种不 关键词： 近似熵 样本熵 复杂度 射击     

应用多尺度化的基本尺度熵分析心率变异性

本文采用多尺度化的基本尺度熵方法,针对心率变异性信号进行了分析,研究发现多尺度化的基本尺度熵可以区分不同生理病理信号,包括健康人、充血性心力衰竭患者和房颤心律失常患者的心率变异性信号,以及健康人白天黑夜的心率变异性信号.通过对健康人代理数据的分析,发现房颤心律失常患者与代理数据的熵值趋势相似,研究结果表明房颤心律失常患者的心率变异性信号更多的是反映生理信号的线性特征,而对环境变化不能很好的进行自我调节. 关键词： 多尺度化的基本尺度熵 心率变异性 充血性心力衰竭 房颤心律失常     

基于多路可视图的健康与心梗患者心电图信号复杂网络识别

可视图(visibility graph, VG)算法已被证明是将时间序列转换为复杂网络的简单且高效的方法,其构成的复杂网络在拓扑结构中继承了原始时间序列的动力学特性.目前,单维时间序列的可视图分析已趋于成熟,但应用于复杂系统时,单变量往往无法描述系统的全局特征.本文提出一种新的多元时间序列分析方法,将心梗和健康人的12导联心电图(electrocardiograph, ECG)信号转换为多路可视图,以每个导联为一个节点,两个导联构成可视图的层间互信息为连边权重,将其映射到复杂网络.由于不同人群的全连通网络表现为完全相同的拓扑结构,无法唯一表征不同个体的动力学特征,根据层间互信息大小重构网络,提取权重度和加权聚类系数,实现对不同人群12导联ECG信号的识别.为判断序列长度对识别效果的影响,引入多尺度权重度分布熵.由于健康受试者拥有更高的平均权重度和平均加权聚类系数,其映射网络表现为更加规则的结构、更高的复杂性和连接性,可以与心梗患者进行区分,两个参数的识别准确率均达到93.3%.     

颠倒睡眠状态调制心率变异性信号的功率谱和基本尺度熵分析

心率变异性(HRV)信号能够提供心脏活动状态的重要信息.通过建立颠倒睡眠模型,联合功率谱和基本尺度熵方法分析颠倒睡眠状态下24 h的HRV信号,研究颠倒睡眠对自主神经相互作用以及HRV信号混沌强度的调制.结果表明,颠倒睡眠导致自主神经在昼夜间的活动节律发生颠倒,基本尺度熵在昼夜的变化趋势也随之发生逆转,因此HRV信号的混沌强度与自主神经的相互作用密切相关,进一步研究两者之间的关系发现:HRV信号的混沌强度与交感神经的调制强度正相关,与迷走神经的调制强度负相关.     

基于Poincare差值散点图的心率变异性分析方法研究

在心率变异性的非线性分析中,Poincare散点图分析是一种重要手段.本文基于Poincare差值散点图(modified Poincare plot)提出了两个参数—区域分布熵和区域分布系数,用于定量描述所考察区域内散点的分布趋势,并提出对散点在4个象限中的分布进行分别计算.通过对MIT-BIH数据库中健康年轻人、健康老年人和充血性心力衰竭患者样本数据的分析,发现两参数值均呈现显著的组间差异;同时,不同象限的分析结果显示了四个象限具有不同的区分敏感性,而其中尤以第一象限的区分度为最高,反映出充血性心力衰竭患者相对健康人迷走神经调控功能的改变最为显著,与以往的生理学研究结论相符.经验证,该方法可用于短时数据,更易于扩展至临床应用.     

基于Poincaré差值散点图的心率变异性分析方法研究

在心率变异性的非线性分析中, Poincaré散点图分析是一种重要手段.本文基于Poincaré差值散点图(modified Poincaré plot)提出了两个参数——区域分布熵和区域分布系数, 用于定量描述所考察区域内散点的分布趋势, 并提出对散点在4个象限中的分布进行分别计算.通过对MIT-BIH数据库中健康年轻人、健康老年人和充血性心力衰竭患者样本数据的分析, 发现两参数值均呈现显著的组间差异;同时, 不同象限的分析结果显示了四个象限具有不同的区分敏感性, 而其中尤以第一象限的区分度为最高, 反映出充血性心力衰竭患者相对健康人迷走神经调控功能的改变最为显著, 与以往的生理学研究结论相符.经验证, 该方法可用于短时数据, 更易于扩展至临床应用.     

Specific Emitter Identification Based on Visibility Graph Entropy

The specific emitter identification(SEI) technique determines the unique emitter of a given signal by using some external feature measurements of the signal. It has recently attracted a great deal of attention because many applications can benefit from it. This work addresses the SEI problem using two methods, namely, the normalized visibility graph entropy(NVGE) and the normalized horizontal visibility graph entropy(NHVGE)based on treating emitters as nonlinear dynamical systems. Firstly, the visibility graph(VG) and the horizontal visibility graph(HVG) are used to convert the instantaneous amplitude, phase and frequency of received signals into graphs. Then, based on the information captured by the VG and the HVG, the normalized Shannon entropy(NSE) calculated from the corresponding degree distributions are utilized as the rf fingerprint. Finally, four emitters from the same manufacturer are utilized to evaluate the performance of the two methods. Experimental results demonstrate that both the NHVGE-based method and NVGE-based method are quite effective and they perform much better than the method based on the normalized permutation entropy(NPE) in the case of a small amount of data. The NVGE-based method performs better than the NHVGE-based method since the VG can extract more information than the HVG does. Moreover, our methods do not distinguish between the transient signal and the steady-state signal, making it practical.     

Improved visibility graph fractality with application for the diagnosis of Autism Spectrum Disorder

Recently, the visibility graph (VG) algorithm was proposed for mapping a time series to a graph to study complexity and fractality of the time series through investigation of the complexity of its graph. The visibility graph algorithm converts a fractal time series to a scale-free graph. VG has been used for the investigation of fractality in the dynamic behavior of both artificial and natural complex systems. However, robustness and performance of the power of scale-freeness of VG (PSVG) as an effective method for measuring fractality has not been investigated. Since noise is unavoidable in real life time series, the robustness of a fractality measure is of paramount importance. To improve the accuracy and robustness of PSVG to noise for measurement of fractality of time series in biological time-series, an improved PSVG is presented in this paper. The proposed method is evaluated using two examples: a synthetic benchmark time series and a complicated real life Electroencephalograms (EEG)-based diagnostic problem, that is distinguishing autistic children from non-autistic children. It is shown that the proposed improved PSVG is less sensitive to noise and therefore more robust compared with PSVG. Further, it is shown that using improved PSVG in the wavelet-chaos neural network model of Adeli and c-workers in place of the Katz fractality dimension results in a more accurate diagnosis of autism, a complicated neurological and psychiatric disorder.     

Entropy Analysis of Heart Rate Variability in Different Sleep Stages

How the complexity or irregularity of heart rate variability (HRV) changes across different sleep stages and the importance of these features in sleep staging are not fully understood. This study aimed to investigate the complexity or irregularity of the RR interval time series in different sleep stages and explore their values in sleep staging. We performed approximate entropy (ApEn), sample entropy (SampEn), fuzzy entropy (FuzzyEn), distribution entropy (DistEn), conditional entropy (CE), and permutation entropy (PermEn) analyses on RR interval time series extracted from epochs that were constructed based on two methods: (1) 270-s epoch length and (2) 300-s epoch length. To test whether adding the entropy measures can improve the accuracy of sleep staging using linear HRV indices, XGBoost was used to examine the abilities to differentiate among: (i) 5 classes [Wake (W), non-rapid-eye-movement (NREM), which can be divide into 3 sub-stages: stage N1, stage N2, and stage N3, and rapid-eye-movement (REM)]; (ii) 4 classes [W, light sleep (combined N1 and N2), deep sleep (N3), and REM]; and (iii) 3 classes: (W, NREM, and REM). SampEn, FuzzyEn, and CE significantly increased from W to N3 and decreased in REM. DistEn increased from W to N1, decreased in N2, and further decreased in N3; it increased in REM. The average accuracy of the three tasks using linear and entropy features were 42.1%, 59.1%, and 60.8%, respectively, based on 270-s epoch length; all were significantly lower than the performance based on 300-s epoch length (i.e., 54.3%, 63.1%, and 67.5%, respectively). Adding entropy measures to the XGBoost model of linear parameters did not significantly improve the classification performance. However, entropy measures, especially PermEn, DistEn, and FuzzyEn, demonstrated greater importance than most of the linear parameters in the XGBoost model.300-s270-s.     

Visibility graph analysis on quarterly macroeconomic series of China based on complex network theory

The visibility graph approach and complex network theory provide a new insight into time series analysis. The inheritance of the visibility graph from the original time series was further explored in the paper. We found that degree distributions of visibility graphs extracted from Pseudo Brownian Motion series obtained by the Frequency Domain algorithm exhibit exponential behaviors, in which the exponential exponent is a binomial function of the Hurst index inherited in the time series. Our simulations presented that the quantitative relations between the Hurst indexes and the exponents of degree distribution function are different for different series and the visibility graph inherits some important features of the original time series. Further, we convert some quarterly macroeconomic series including the growth rates of value-added of three industry series and the growth rates of Gross Domestic Product series of China to graphs by the visibility algorithm and explore the topological properties of graphs associated from the four macroeconomic series, namely, the degree distribution and correlations, the clustering coefficient, the average path length, and community structure. Based on complex network analysis we find degree distributions of associated networks from the growth rates of value-added of three industry series are almost exponential and the degree distributions of associated networks from the growth rates of GDP series are scale free. We also discussed the assortativity and disassortativity of the four associated networks as they are related to the evolutionary process of the original macroeconomic series. All the constructed networks have “small-world” features. The community structures of associated networks suggest dynamic changes of the original macroeconomic series. We also detected the relationship among government policy changes, community structures of associated networks and macroeconomic dynamics. We find great influences of government policies in China on the changes of dynamics of GDP and the three industries adjustment. The work in our paper provides a new way to understand the dynamics of economic development.     

结合可视图的多状态交通流时间序列特性分析

交通流时间序列的研究主要采用数据挖掘和机器学习的方法,这些"黑箱"挖掘方法很难直观反映序列特性.为增强交通流时间序列及其特征分析的可视化性,结合可视图理论来构建交通流时间序列的关联网络,从复杂网络角度实现交通流时间序列的特性分析.在网络构建的过程中,考虑到不同交通状态下交通流表征具有的差异性,首先利用交通流参量的相关性对交通流状态进行分类,然后构建不同交通状态下的时间序列复杂网络,并对这些网络的特征属性给出统计分析,如度分布、聚类系数、网络直径、模块化等.研究表明,可视图法可为交通流时间序列映射到网络提供有效途径,并且不同状态下交通流时间序列构建的复杂网络的模块化、聚类系数和度分布等统计特征呈现一定的变化规律,为交通流运行态势的研究提供了可视化的分析角度.     

基于时间序列符号化模式表征的有向加权复杂网络

时间序列复杂网络分析近些年已发展成为非线性信号分析领域的一个国际热点课题.为了能更有效地挖掘时间序列(特别是非线性时间序列)中的结构特征,同时简化时间序列分析的复杂度,提出了一种新的基于时间序列符号化结合滑窗技术模式表征的有向加权复杂网络建网方法.该方法首先按照等概率区段划分的方式将时间序列做符号化处理,结合滑窗技术确定不同时刻的符号化模式作为网络的节点;然后将待分析时间序列符号化模式的转换频次和方向作为网络连边的权重和方向,从而建立时间序列有向加权复杂网络.通过对Logistic系统不同参数设置对应的时间序列复杂网络建网测试结果表明,相比经典的可视图建网方法,本文方法的网络拓扑能更简洁、直观地展示时间序列的结构特征.进而,将本文方法应用于规则排列采集的自然风场信号分析,其网络特性指标能较准确地预测采集信号的排布规律,而可视图建网方法的网络特性指标没有任何规律性的结果.     

Extraction of Nuclear Matter Properties from Nuclear Masses by a Model of Equation of State

The extraction of nuclear matter properties from measured nuclear masses is investigated in the energy density functional formalism of nuclei.It is shown that the volume energy a1 and the nuclear incompressibility Ko depend essentially on μ_nN + μ_pZ - 2E_N,whereas the symmetry energy J and the density symmetry coefficient L as well as symmetry incompressibility Ks depend essentially on μ_n - μ_p,where μ_p ＝μp - ∂E_c/∂Z,μ_n and μ_p are the neutron and proton chemical potentials respectively,E_N the nuclear energy,and E_c the Coulomb energy.The obtained symmetry energy is J ＝ 28.5 MeV,while other coefficients are uncertain within ranges depending on the model of nuclear equation of state.