世纪气象 » 学海笑谈 » 老C模式--仅供娱乐。。。

cwnraul 发表于 2010-6-30 13:19

老C模式--仅供娱乐。。。

最近在学习混沌时间序列预测，jF[;q4? e;^
按照最原始最基本最简单的混沌时间序列相空间重构方法（互信息+虚假邻点+零阶局域），编了个小程序，
n1\&^V)yB0\e6p 然后把SOI（1876年1月至今）、太阳黑子每年5月的数目（1700至今）以及太阳黑子每个月的数目（1750年1月至今）看成混沌时间序列，预测了一下（仅供娱乐）。。。。3Z-{y6E)t[-lIQv

9C~#a5i,IbP.E~ 太阳黑子（年）：根据1700-1999的数据，预测了一下2000-2015的情况，和实况校验结果如下：
z5E:G]5T$O#qCi [attach]42457[/attach]dg^-kJ(vF Nl
由于太阳黑子数的周期性还是挺明显的，看来预测的还算可以。当然这两年黑子数出奇地低，连人家官方的模式都预测不到，何况我这个胡乱编出来的程序。。。BkUQ a&hI$P o6D }

*er6gG#Q@Z([N%d /g}&`$c3Yr
太阳黑子（月），根据1750年1月--2005年12月的情况，预测了一下2006年1月开始的情况，和实况校验结果如下：Jk;~#[%BwAf
[attach]42458[/attach]7otr$Fg(_Rk
月份的数据更多更复杂，但同上，太阳黑子数的周期性还是挺明显的，所以总体的趋势还是能基本预测出来。。。-nd IwkE,ViJ

X(w Z6AmP8? M&qY ?A"nDc.F
SOI（月），根据1876年1月--2005年12月的情况，预测了一下2006年1月开始的情况，和实况校验结果如下："q;_!xNb},WI&hN
[attach]42459[/attach]
!e)DyG'O:[9[ 这个预测结果就非常差了。。。完全预测不到。。。。预测的拉尼娜厄尔尼诺都完全反了。。。。看来SOI数据太过复杂，这个破程序没法搞掂。。。

kitty8123 发表于 2010-6-30 17:47

老C模式可以演算大气环流吗？:y7

cwnraul 发表于 2010-6-30 22:50

[quote]老C模式可以演算大气环流吗？:y7cS,y'{?2c&n
[size=2][color=#999999]kitty8123 发表于 2010-6-30 17:47[/color] [url=http://www.21cma.net/redirect.php?goto=findpost&pid=60201&ptid=7662][img]http://www.21cma.net/images/common/back.gif[/img][/url][/size][/quote]
;} F"@]/P3Q:_7l
w0c;K(c+^} 目前这个简单程序，离演算大气环流还差了很远。。。。

methk 发表于 2010-7-1 12:49

老C！你的程序是用FORTRAN, C++還是其他語言寫的？

cwnraul 发表于 2010-7-1 20:12

用C语言的，图是根据output出来的数据，再用其他软件画的。。。。。

爱我海坛 发表于 2010-7-3 08:00

哈哈，SOI震荡非常厉害，比股市还厉害，所以预测不到是难免呢~~~~~

psc 发表于 2010-7-5 00:23

老C：Q!^,X;l%O
SOI本身是一种人工统计手段，所以不适合混沌模式推导。
D3xnm] P   d-vUJ,E ~$O,b7x,k%}
或许可以将T和D岛气压做混沌延伸，然后再按SOI定义计算结果；或者直接延伸Nino监测区。

cwnraul 发表于 2010-8-13 20:19

[quote]老C： VV:rH@
SOI本身是一种人工统计手段，所以不适合混沌模式推导。
0@3Y1i1A7A.Z5f G({F   
D,H7~R2Bor] 或许可以将T和D岛气压做混沌延伸，然后再按SOI定义计算结果；或者直接延伸Nino监测区。
@s3I x'P0GF!L L)~ [size=2][color=#999999]psc 发表于 2010-7-5 00:23[/color] [url=http://www.21cma.net/redirect.php?goto=findpost&pid=60377&ptid=7662][img]http://www.21cma.net/images/common/back.gif[/img][/url][/size][/quote]6}CF1P!k&[0Ox

UT r8vQMm5g? 非常有道理。。。
*R~5AAP? 于是我在更新了算法之余，换成用NINO3.4-Index来进行预测&D _E3a)D`N
.c].E#pFN
sunspot1这个图，基于1700-1960年的太阳黑子数据，预测1961-2010年的"`\[a4kc
可见已经预测得比较准确了。%]9^lq!VUj'B?
[attach]44715[/attach]
$E4qy,w}.{_ o5g$u F)\E|!KF4b
sunspot2这个图，基于1700-2009年，预测后面的。。。。看看仅供娱乐
,PC7q4HG&BwI [attach]44716[/attach]8mNt:F!\J&J`D
:K+Z&xxw^8EM
参考LS的意见，采用NINO3.4数据来预测
o6oH l5w1vPks 基于1865-1995的数据（每月的数据都有），预测1996-2010的C(|)hC.HDZ)N^hP
虽然不能说很准确，但至少趋势上有点似了吧。
eb3NK2m
y [attach]44717[/attach]

通灵草 发表于 2010-8-13 20:46

严重同意sunspot2这个图，这是我一直怀疑的。,~-beIT%Vh'l
  
)wn&[X9N$M.? 至于NINO3.4结果，1997/1998年居然跟现有模型一样意外，看来1997/1998年有明显外因干扰。
1e-q$hQ ^%gf2u)o   6[A-SH ~ _C
另外，人类观测数据也是有限的。

cwnraul 发表于 2010-8-13 20:58

[quote]严重同意sunspot2这个图，这是我一直怀疑的。
7v5A!b}o2Sk N,v ^   l P]&GB
至于NINO3.4结果，1997/1998年居然跟现有模型一样意外，看来1997/1998年有明显外因干扰。j `Af$i
~
  
)w;r M#?t&N4uT 另外，人类观测数据也是有限的。0z`b%k*n"`
[size=2][color=#999999]通灵草 发表于 2010-8-13 20:46[/color] [url=http://www.21cma.net/redirect.php?goto=findpost&pid=63016&ptid=7662][img]http://www.21cma.net/images/common/back.gif[/img][/url][/size][/quote]
2Ay
y+n6q5r*Q.s )O2K(GwZm_R
1997-1998那个强厄，从时间序列这个角度，我换了很多方法都没有能够套出来。。lJi,]W\0b
其实我选从1996开始预测窗口，就是为了看看能不能套出这个强厄的
)J\0@$ti ~%TB 很遗憾不能。。。

通灵草 发表于 2010-8-13 21:03

[b] [url=http://www.21cma.net/redirect.php?goto=findpost&pid=63020&ptid=7662]10#[/url] [i]cwnraul[/i] [/b](a.cElY y?V
然后对本次ENSO过程，是否怎么调整模式也是反应微弱？

cwnraul 发表于 2010-8-13 21:05

[quote] 10# cwnraul
%B&QAyp6B E 然后对本次ENSO过程，是否也怎么调整模式也是反应微弱？Y_.J:hP_M
[size=2][color=#999999]通灵草 发表于 2010-8-13 21:03[/color] [url=http://www.21cma.net/redirect.php?goto=findpost&pid=63021&ptid=7662][img]http://www.21cma.net/images/common/back.gif[/img][/url][/size][/quote]
t4b$QK"Kskiw?3? S }(JS8P8\s)m|q
未试过，如果要测本次ENSO过程，应该将最新的数据也用作训练集了
*|r-Vy0n(f&zD0Q%P 前面那个是基于到1995的数据的，预测2010之后不靠谱
v
U~DOg7O 以后有空运行一下1856-2009作训练集的看看

0001 发表于 2010-8-14 02:00

老C模式能否将黑子模型数据集关联到NINO3.4模型数据集中，然后运算NINO3.4，假设黑子也受到跟ENSO一样的未知外因子干扰，估计97/98、09/10 ENSO就能体现。w;Tb_4q9lLz}"u$e
  &Uy T2R W/rT&E I
不过，那不是黑子的关系，只是“黑子也受到跟ENSO一样的未知外因子干扰”，所以可以从复合混沌模式中提取到共同因子。尽管不知那是什么因子，但可以找出这个因子的存在，然后再追踪来源。
uS+Z9Z$WG)A Y9C   
t8cAYaX6neZ 防腐、调香同理。
2m.kYgc+a
`N   
6J:S VKR1i 知可知之不可知，知不可知之可知——独裁9式：总决式——仅供娱乐。。。

liuxinmin 发表于 2010-8-14 11:54

我很怀疑98年ENSO过大的正偏离,可能会在2010年后被纠正(过大的负偏离),以至达到总体的平衡,目前看海温的演变正向此靠拢

雪中红 发表于 2011-1-13 01:39

按目前情况看，第8楼的“看看仅供娱乐”将成功运算出至少相当于道尔顿极小期的到来，即是相当于慈禧时代{:soso_e120:}

cwnraul 发表于 2011-6-15 21:29

最后时刻，顶顶自己的帖子最后留个名字吧。。。挖坟。。
)r1lpu'U#p
:q3M#a E\-]^X 其实相关算法的论文已经投稿了，不过至今没有回音。
(z!Tgz!\9aH 那边没有回音，这边却已经要结业了。。。。。

查看完整版本: 老C模式--仅供娱乐。。。