基于有序模式的度量对多变量时间序列进行非线性分析研究(Matlab代码实现)

简介: 基于有序模式的度量对多变量时间序列进行非线性分析研究(Matlab代码实现)

💥1 概述

基于有序模式的度量(Ordinal Pattern-based Measures)是一种用于多变量时间序列非线性分析的方法。它可以通过分析时间序列中的有序模式来获取关于序列动力学和非线性性质的信息。以下是基于有序模式的度量方法的一般研究步骤:


1. 多变量时间序列选择:选择你感兴趣的多变量时间序列数据。这些数据可以代表各种应用领域中的观测变量,例如金融数据、天气数据或生物学数据等。


2. 计算有序模式:将多变量时间序列数据转换为有序模式序列。有序模式是基于变量之间的相对大小关系而形成的序列。这可以通过将原始时间序列中的每个数据点与其他数据点进行比较来实现。


3. 确定有序模式:对于给定的模式长度(pattern length),确定所有可能的有序模式。模式长度是在指定的时间窗口内对变量进行比较的步骤数量。


4. 计算模式频率:计算每个有序模式的频率或出现次数。这可以通过扫描整个有序模式序列并对每个模式进行计数来实现。


5. 模式度量计算:基于模式频率,计算一系列模式度量来描述时间序列的非线性性质。一些常见的模式度量包括熵、复杂度、不规则度等。这些度量可以量化时间序列的复杂性和规律性。


6. 统计分析:对模式度量进行统计分析,以获取关于时间序列非线性性质的信息。这可以包括计算平均值、标准差、相关性等统计指标,以及使用假设检验方法来评估时间序列之间的差异。


7. 结果解释:解释和解读分析结果,以洞察时间序列的非线性动力学和相关特性。这可能需要与领域专家结合,将模式度量结果与问题域的理论知识进行关联。


需要注意的是,基于有序模式的度量方法在时间序列非线性分析中是一种较新的方法,需要对相关理论和算法进行深入了解。此外,选择合适的模式长度、模式度量和统计分析方法等也需要针对具体问题进行权衡和调整。


📚2 运行结果

部分代码:

%% compute permutation entropy in sliding windows
load( 'tonicClonic.mat' );
cfg            = [];
cfg.method     = 'PE'; % compute permutation entropy
cfg.order      = 3;    % ordinal pattens of order 3 (4-points ordinal patterns)
cfg.delay      = 2;    % delay 2 between points in ordinal patterns 
                       % (one point between successive points in ordinal patterns)
cfg.windowSize = 512;  % window size = 512 time steps
cfg.time       = 0:1/102.4:179.999; % OPTIONAL time axis for plotting
cfg.units      = 'seconds';         % OPTIONAL units of time for plotting
outdata        = OPanalysis( cfg, indata );
%% compute permutation entropy and ordinal distributions in sliding windows
load( 'tonicClonic.mat' );
cfg            = [];
cfg.method     = 'opdPE'; % compute permutation entropy
cfg.order      = 3;       % ordinal pattens of order 3 (4-points ordinal patterns)
cfg.orderSeq   = 6;       % ordinal pattens of order 6 for plotting their sequence (7-points ordinal patterns)
cfg.delay      = 1;       % delay 1 between points in ordinal patterns (successive points)
cfg.windowSize = 1024;    % window size = 1024 time steps
cfg.time       = 0:1/102.4:179.999; % OPTIONAL time axis for plotting
cfg.units      = 'seconds';         % OPTIONAL units of time for plotting
outdata        = OPanalysis( cfg, indata );
%% compute all the implemented measures simultaneously for comparison
load( 'tonicClonic.mat' );
cfg                = [];
cfg.method         = 'all';  % compute all implemented ordinal-patterns-based measures
cfg.order          = 4;      % ordinal patterns of order 4 (5-points ordinal patterns)
cfg.delay          = 1;      % delay 1 between points in ordinal patterns
cfg.windowSize     = 512;    % window size = 512 time steps
cfg.lowerThreshold = 0.2;    % the distance considered negligible between points
cfg.upperThreshold = 200;    % the distance between points most probably related to artifact
cfg.time           = 0:1/102.4:179.999; % OPTIONAL time axis for plotting


🎉3 参考文献

部分理论来源于网络,如有侵权请联系删除。

[1] Amigo, J.M., Keller, K. and Unakafova, V.A., 2015. On entropy, entropy-like quantities, and applications. Discrete & Continuous Dynamical Systems-Series B, 20(10).

[2] Bandt C., Pompe B., Permutation entropy: a natural complexity measure for time series. Physical review letters, 2002, APS

[3] Keller, K., and M. Sinn. Ordinal analysis of time series. Physica A: Statistical Mechanics and its Applications 356.1 (2005): 114--120

[4] Keller, K., Unakafov, A.M. and Unakafova, V.A., 2014. Ordinal patterns, entropy, and EEG. Entropy, 16(12), pp.6212-6239.

[5] Zanin, M., Zunino, L., Rosso, O.A. and Papo, D., 2012.

Permutation entropy and its main biomedical and econophysics applications: a review. Entropy, 14(8), pp.1553-1577.

[6] Unakafova, V.A., Keller, K., 2013. Efficiently measuring complexity on the basis of real-world Data. Entropy, 15(10), 4392-4415.

[7] Unakafova, V.A., 2015. Investigating measures of complexity for dynamical systems and for time series (Doctoral dissertation, University of Luebeck).

[8] Bian, C., Qin, C., Ma, Q.D. and Shen, Q., 2012. Modified permutation-entropy analysis of heartbeat dynamics. Physical Review E, 85(2), p.021906.

[9] Amigo, J.M., Zambrano, S. and Sanjuan, M.A., 2008. Combinatorial detection of determinism in noisy time series. EPL (Europhysics Letters), 83(6), p.60005.

[10] Cao, Y., Tung, W.W., Gao, J.B. et al., 2004. Detecting dynamical changes in time series using the permutation entropy. Physical Review E, 70(4), p.046217.

[11] Riedl, M., Muller, A. and Wessel, N., 2013. Practical considerations of permutation entropy. The European Physical Journal Special Topics, 222(2), pp.249-262.


🌈4 Matlab代码实现

相关文章
|
1天前
|
机器学习/深度学习 算法 数据挖掘
基于PSO优化的CNN-GRU-Attention的时间序列回归预测matlab仿真
摘要: 本文介绍了运用粒子群优化(PSO)调整深度学习模型超参数以提升时间序列预测性能的方法。在比较了优化前后的效果(Ttttttttttt12 vs Ttttttttttt34)后,阐述了使用matlab2022a软件的算法。文章详细讨论了CNN、GRU网络和注意力机制在时间序列预测中的作用,以及PSO如何优化这些模型的超参数。核心程序展示了PSO的迭代过程,通过限制和调整粒子的位置(x1)和速度(v1),寻找最佳解决方案(gbest1)。最终,结果保存在R2.mat文件中。
基于PSO优化的CNN-GRU-Attention的时间序列回归预测matlab仿真
|
3天前
|
算法 数据安全/隐私保护 数据格式
基于混沌序列的图像加解密算法matlab仿真,并输出加解密之后的直方图
该内容是一个关于混沌系统理论及其在图像加解密算法中的应用摘要。介绍了使用matlab2022a运行的算法,重点阐述了混沌系统的特性,如确定性、非线性、初值敏感性等,并以Logistic映射为例展示混沌序列生成。图像加解密流程包括预处理、混沌序列生成、数据混淆和扩散,以及密钥管理。提供了部分核心程序,涉及混沌序列用于图像像素的混淆和扩散过程,通过位操作实现加密。
|
3天前
|
数据采集 算法 数据可视化
MATLAB、R用改进Fuzzy C-means模糊C均值聚类算法的微博用户特征调研数据聚类研究
MATLAB、R用改进Fuzzy C-means模糊C均值聚类算法的微博用户特征调研数据聚类研究
10 1
|
5天前
|
算法 数据建模
MATLAB随机波动率SV、GARCH用MCMC马尔可夫链蒙特卡罗方法分析汇率时间序列
MATLAB随机波动率SV、GARCH用MCMC马尔可夫链蒙特卡罗方法分析汇率时间序列
15 6
|
5天前
Matlab用向量误差修正VECM模型蒙特卡洛Monte Carlo预测债券利率时间序列和MMSE 预测
Matlab用向量误差修正VECM模型蒙特卡洛Monte Carlo预测债券利率时间序列和MMSE 预测
23 14
|
7天前
|
机器学习/深度学习 传感器 数据可视化
MATLAB用深度学习长短期记忆 (LSTM) 神经网络对智能手机传感器时间序列数据进行分类
MATLAB用深度学习长短期记忆 (LSTM) 神经网络对智能手机传感器时间序列数据进行分类
24 1
MATLAB用深度学习长短期记忆 (LSTM) 神经网络对智能手机传感器时间序列数据进行分类
|
8天前
|
机器学习/深度学习 算法 数据挖掘
基于PSO优化的CNN-LSTM-Attention的时间序列回归预测matlab仿真
该文档介绍了使用MATLAB2022A中PSO优化算法提升时间序列预测模型性能的过程。PSO优化前后对比显示了优化效果。算法基于CNN、LSTM和Attention机制构建CNN-LSTM-Attention模型,利用PSO调整模型超参数。代码示例展示了PSO的迭代优化过程及训练、预测和误差分析环节。最终,模型的预测结果以图形形式展示,并保存了相关数据。
|
12天前
|
数据可视化 Python
Matlab正态分布、历史模拟法、加权移动平均线 EWMA估计风险价值VaR和回测标准普尔指数 S&P500时间序列
Matlab正态分布、历史模拟法、加权移动平均线 EWMA估计风险价值VaR和回测标准普尔指数 S&P500时间序列
36 11
|
3月前
|
Perl
【MFAC】基于全格式动态线性化的无模型自适应控制(Matlab代码)
【MFAC】基于全格式动态线性化的无模型自适应控制(Matlab代码)
|
3月前
【数值分析】迭代法求方程的根(附matlab代码)
【数值分析】迭代法求方程的根(附matlab代码)

热门文章

最新文章