世纪气象云

 找回密码
 注册

QQ登录

只需一步,快速开始

查看: 8405|回复: 20
收起左侧

[学习交流] 老C模式--仅供娱乐。。。

[复制链接]
发表于 2010-6-30 13:19 | 显示全部楼层 |阅读模式
最近在学习混沌时间序列预测,
按照最原始最基本最简单的混沌时间序列相空间重构方法(互信息+虚假邻点+零阶局域),编了个小程序,
然后把SOI(1876年1月至今)、太阳黑子每年5月的数目(1700至今)以及太阳黑子每个月的数目(1750年1月至今)看成混沌时间序列,预测了一下(仅供娱乐)。。。。

太阳黑子(年):根据1700-1999的数据,预测了一下2000-2015的情况,和实况校验结果如下:
sunspot year.jpg
由于太阳黑子数的周期性还是挺明显的,看来预测的还算可以。当然这两年黑子数出奇地低,连人家官方的模式都预测不到,何况我这个胡乱编出来的程序。。。


太阳黑子(月),根据1750年1月--2005年12月的情况,预测了一下2006年1月开始的情况,和实况校验结果如下:
sunspot month.jpg
月份的数据更多更复杂,但同上,太阳黑子数的周期性还是挺明显的,所以总体的趋势还是能基本预测出来。。。


SOI(月),根据1876年1月--2005年12月的情况,预测了一下2006年1月开始的情况,和实况校验结果如下:
soi.jpg
这个预测结果就非常差了。。。完全预测不到。。。。预测的拉尼娜厄尔尼诺都完全反了。。。。看来SOI数据太过复杂,这个破程序没法搞掂。。。

评分

1

查看全部评分

本帖被以下淘专辑推荐:

发表于 2010-6-30 17:47 | 显示全部楼层
老C模式可以演算大气环流吗?:y7
 楼主| 发表于 2010-6-30 22:50 | 显示全部楼层
发表于 2010-7-1 12:49 | 显示全部楼层
老C!你的程序是用FORTRAN, C++還是其他語言寫的?
 楼主| 发表于 2010-7-1 20:12 | 显示全部楼层
用C语言的,图是根据output出来的数据,再用其他软件画的。。。。。
发表于 2010-7-3 08:00 | 显示全部楼层
哈哈,SOI震荡非常厉害,比股市还厉害,所以预测不到是难免呢~~~~~
发表于 2010-7-5 00:23 | 显示全部楼层
老C:
SOI本身是一种人工统计手段,所以不适合混沌模式推导。
  
或许可以将T和D岛气压做混沌延伸,然后再按SOI定义计算结果;或者直接延伸Nino监测区。

评分

1

查看全部评分

 楼主| 发表于 2010-8-13 20:19 | 显示全部楼层
老C:
SOI本身是一种人工统计手段,所以不适合混沌模式推导。
  
或许可以将T和D岛气压做混沌延伸,然后再按SOI定义计算结果;或者直接延伸Nino监测区。
psc 发表于 2010-7-5 00:23


非常有道理。。。
于是我在更新了算法之余,换成用NINO3.4-Index来进行预测

sunspot1这个图,基于1700-1960年的太阳黑子数据,预测1961-2010年的
可见已经预测得比较准确了。
sunspot1.jpg

sunspot2这个图,基于1700-2009年,预测后面的。。。。看看仅供娱乐
sunspot2.jpg

参考LS的意见,采用NINO3.4数据来预测
基于1865-1995的数据(每月的数据都有),预测1996-2010的
虽然不能说很准确,但至少趋势上有点似了吧。
NINO34.jpg
发表于 2010-8-13 20:46 | 显示全部楼层
严重同意sunspot2这个图,这是我一直怀疑的。
  
至于NINO3.4结果,1997/1998年居然跟现有模型一样意外,看来1997/1998年有明显外因干扰。
  
另外,人类观测数据也是有限的。
 楼主| 发表于 2010-8-13 20:58 | 显示全部楼层
严重同意sunspot2这个图,这是我一直怀疑的。
  
至于NINO3.4结果,1997/1998年居然跟现有模型一样意外,看来1997/1998年有明显外因干扰。
  
另外,人类观测数据也是有限的。
通灵草 发表于 2010-8-13 20:46


1997-1998那个强厄,从时间序列这个角度,我换了很多方法都没有能够套出来。。
其实我选从1996开始预测窗口,就是为了看看能不能套出这个强厄的
很遗憾不能。。。
您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|手机版|世纪气象云 | | | 关于我们

GMT+8, 2019-9-19 18:25 , Processed in 0.776240 second(s), 22 queries .

Powered by Discuz! X3.4

© 2001-2017 Comsenz Inc.

快速回复 返回顶部 返回列表