MATLAB1770太阳黑子,基于MATLAB的太阳黑子时间序列与仿真
2012.3
26 基于MATLAB 的太阳黑子时间序列
分析与仿真
周园 肖洪祥 董俊飞
桂林理工大学信息科学与工程学院 广西 541004
摘要:本文研究了时间序列的分析方法,具体分析了基于最大Lyapunov 指数的方法在太阳黑子时间序列分析中的应用。介绍利用MATLAB 对太阳黑子时间序列进行分析与仿真的方法,并给出相关的流程、程序和相应的仿真结果。最终证明太阳黑子时间序列是一个混沌时间序列。
关键词:混沌时间序列;最大Lyapunov 指数;太阳黑子数;仿真
0 引言
在非线性系统中,初始条件的微小变化,往往会导致结果以指数级的大小发生分离,这时我们称这个系统存在混沌。时间序列是非线性动力系统的一种模型。如果时间序列对初始条件敏感,采用传统线性时间序列分析方法将很难予以分析,因此传统时间序列预测模型对混沌时间序列的拟合和预测准确度都很差。经过混沌学的发展,可以使用序列本身的规律对其进行预测。Lyapunov 指数法即是其中之一。通过最大Lyapunov 指数的数值,可以判断一个时间序列是否是混沌时间序列,亦即该非线性系统中是否存在着混沌。本文对太阳黑子序列进行分析,证明其是一个混沌时间序列。
1 基于Lyapunov 指数的时间序列分析方法
对时间序列进行分析,首先必须进行相空间重构。根据有限的数据重构吸引子以研究系统动力行为的方法即是相空间重构。主要思想为:系统中每个分量的演化皆是由与之联系的其他分量所决定的,相关分量的信息隐含在任意其他分量的变化过程中,即是运用系统的任何一个观察量可以重构出整个系统的模型。
设时间序列为{}t x ,其中1,2,...,t N =。重构相空间m
R 的元素组为:
(1)(,,)(,,...,),T=1,2,3,...,T T T T m X m N X X X p τττ++-= (1)
其中,N 为重构相空间维数;τ为延迟时间间隔数,且为正整数;(1)p N m τ=--为时间序列嵌入相空间的向量数,N 为时间序列的数据点数。
由Tokens 定理,在理论条件下可任选τ。但在现实条件下时间序列都是有限长且有噪声的。因而在重构相空间时,
τ的选取至关重要。目前所采用的方法大多是通过经验来选
择τ, 从而使得T X 和T X τ+相互独立并不完全相关。
Lyapunov 指数是描述奇异吸引子性质的数据量。在m 维离散系统中存在m 个Lyapunov 指数,即Lyapunov 指数族。正的Lyapunov 指数意为在此维度方向,系统以指数级速度分离。1983年,G.Grebogi 证明了若最大Lyapunov 指数
max 0λ>,则系统一定存在着混沌。因此要判断一个时间序
列是否为混沌时间序列,必须求出其最大Lyapunov 指数。为了保证领域点沿着不同的轨道运动,最近邻域点间必须有分离间隔。此处取分离间隔为/w T t =?,其中T 为用FFT 计算出的序列平均周期;t ?为序列的采样周期。
2 计算机仿真步骤
仿真步骤如图1所示。输入太阳黑子年平均序列,通过
FFT 算法计算得到其平均周期T 。计算分离间隔作为时间窗
ωτ。由公式得到嵌入维数m 。运用所得的参数使用Wolf 法
算出最大Lyapunov 指数。进而判断该序列是否是混沌序列。
MATLAB1770太阳黑子,基于MATLAB的太阳黑子时间序列与仿真相关推荐
- matlab频分复用,基于MATLAB的频分复用系统的仿真_.doc
基于MATLAB的频分复用系统的仿真_ 基于MATLAB的频分复用系统的仿真_毕业论文(设计) Abstract With the development of communication techn ...
- matlab故障识别,基于Matlab的电力系统故障分析与仿真(V2.1)最新版
<基于Matlab的电力系统故障分析与仿真.doc>由会员分享,可免费在线阅读全文,更多与<基于Matlab的电力系统故障分析与仿真(V2.1)>相关文档资源请在帮帮文库(ww ...
- matlab非线性系统频域标识,基于MATLAB的最小二乘法系统辨识与仿真
第 29卷 第 2期 许昌学院学报 Vol. 29. No. 2 2010年 3月 JOURNAL OF XUCHANG UN IVERSITY Mar. 2010 收稿日期: 2008 - 10 - ...
- tcsc工作原理matlab仿真,基于Matlab的TCSC建模与仿真研究.doc
基于Matlab的TCSC建模与仿真研究 基于Matlab的TCSC建模与仿真研究 第17卷第5期 2006年1O月 巾原T学院 JOURNALOFZHONGYUANINSTIT[ITEOFTECHN ...
- 【基于MATLAB的火灾疏散模拟仿真】——安全隐患提前发现,疏散方案优化
[基于MATLAB的火灾疏散模拟仿真]--安全隐患提前发现,疏散方案优化 随着城市化进程的不断加速,人口密度越来越大,特别是在高层建筑中.万一发生火灾,往往会对人的生命和财产造成严重损失.因此,火灾疏 ...
- matlab模拟Fraunhofer衍射,基于Matlab的夫琅禾费衍射光学仿真.doc
基于Matlab的夫琅禾费衍射光学仿真 基于Matlab的夫琅禾费衍射光学仿真 摘要 计算机仿真技术是以多种学科和理论为基础,以计算机及其相应的软件为工具,通过虚拟试验的方法来分析和解决问题的一门综合 ...
- matlab link offset,基于MATLAB教学型机器人空间轨迹仿真
基于MATLAB教学型机器人空间轨迹仿真 robotic toolbox for matlab工具箱 1. PUMA560的MATLAB仿真 要建立PUMA560的机器人对象,首先我们要了解PUMA5 ...
- 【机械仿真】基于matlab GUI曲柄摇杆机构运动仿真【含Matlab源码 1608期】
一.获取代码方式 获取代码方式1: 完整代码已上传我的资源:[机械仿真]基于matlab GUI曲柄摇杆机构运动仿真[含Matlab源码 1608期] 点击上面蓝色字体,直接付费下载,即可. 获取代码 ...
- 低通采样定理 matlab,基于matlab的低通抽样定理仿真
基于matlab的低通抽样定理仿真 DSP 课 程 设 计 专业: 电子信息技术工程 年级: 2011 级 不 姓名: 陈兰兰 学号: 20113015 指导教师: 刘 德 春 阿坝师专电子信息工程系 ...
最新文章
- Adam那么棒,为什么还对SGD念念不忘 (1) —— 一个框架看懂优化算法
- 高性能存储之--快速理解redis(简版)
- 如何查询当前手机的cpu架构,so库导入工程又出异常了?
- https跨域到http问题解决
- PHP:第三章——PHP中控制函数的函数
- 160523、Oracle建立表空间和用户
- C语言 strcat_s 函数 - C语言零基础入门教程
- Duplicate Symbol链接错误的原因总结和解决方法[转]
- golang---文件读写
- Android 系统性能优化(77)---电量优化
- Hive - Create TableDrop Table ALTER Table(中)
- STM32 ADC没有输入电压时,采集结果不为0
- EFM32片内外设--TFT 汉字字库
- 【算法导论】有向图的可达矩阵
- python 基于模板生成ppt_【Python3】通过模板实现PPT的自动生成
- jetson nano 3(Archiconda安装)
- iOS 给三方日历加上农历
- web生成pdf后的直接打印,隐藏iframe
- 期货止损篇2:ATR止损,高手都在悄悄用的止损利器
- 【MATLAB】遗传算法(GA)求解TSP问题