声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2530|回复: 5

[分形与混沌] 混沌时间序列的预测误差

[复制链接]
发表于 2014-3-17 10:59 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
用多步预测方法预测混沌时间序列,但是预测值是哪个数据啊?好心人帮个忙吧,以下是代码:
  1. clear all
  2. data=load('shuju(359x96).txt');   
  3. %A=data(:,[4,8,12,16,20,24,28,32,36,40,44,48,52,56,60,64,68,72,76,80,84,88,92,96]);
  4. A=data(1:337,[4,8,12,16,20,24,28,32,36,40,44,48,52,56,60,64,68,72,76,80,84,88,92,96]);
  5. m1=336;
  6. m2=m1+1;
  7. for i=1:m2
  8.   for j=1:24
  9.     X((i-1)*24+j,1)=A(i,j);
  10.    end
  11. end
  12. roll=X;%取X数据
  13. len=m1*24;
  14. oldroll=roll;
  15. AverRoll=roll-mean(roll);
  16. roll=AverRoll;
  17. mtbp=4;
  18. taow=mtbp;
  19. deltaT=1;
  20. tao=2;
  21. taowdivtao=round(taow/tao);
  22. d=3;%taowdivtao+1;
  23. prestep=24;
  24. N = length(roll)-prestep;
  25. MaxStep = 24;
  26. Step = 1 : 1 : MaxStep;
  27. [PredictedData] = FunctionChaosPredict(roll(1:N),N,mtbp,deltaT,tao,d,MaxStep);
  28. subplot(2,1,1);plot(roll(1:len+MaxStep)+mean(oldroll));
  29. hold on;
  30. subplot(2,1,1);
  31. plot(Step + N,PredictedData+mean(oldroll),'r:');
  32. hold on;
  33. subplot(2,1,2);
  34. plot(Step + N,PredictedData+mean(oldroll),'r:');
  35. hold on;
  36. subplot(2,1,2);
  37. plot(Step + N,oldroll(Step + N),'b-');




复制代码


回复
分享到:

使用道具 举报

发表于 2014-3-17 14:45 | 显示全部楼层
从代码上看 PredictedData+mean(oldroll) 应该是程序给出的预测值
 楼主| 发表于 2014-3-17 16:39 | 显示全部楼层
gghhjj 发表于 2014-3-17 14:45
从代码上看 PredictedData+mean(oldroll) 应该是程序给出的预测值

谢谢,您讲的是正确的。还想问您下,风速时间序列一般用什么方法重构相空间啊
发表于 2014-3-18 08:56 | 显示全部楼层
原来是小飞 发表于 2014-3-17 16:39
谢谢,您讲的是正确的。还想问您下,风速时间序列一般用什么方法重构相空间啊

相空间重构方法很多,常用的有伪最近邻点法、奇异值分解法、自相关和互信息法等。但都存在这样或者那样的一些问题,个人对风速时间序列的特征不太了解,不知道具体什么方法更好,这需要你结合时间序列的特征以及你所掌握的的数据情况进行选择尝试,这方面的文献应该也是有的,你可以参考参考
发表于 2014-4-6 13:41 | 显示全部楼层
风功率呢  我是初学者 不太懂  能耽误你宝贵的时间给我讲解一下吗
发表于 2014-4-8 15:04 | 显示全部楼层
ABC2014 发表于 2014-4-6 13:41
风功率呢  我是初学者 不太懂  能耽误你宝贵的时间给我讲解一下吗

两则应该是类似的。上面讲的是相空间重构
如果是预测的话,现在一般用神经网络,ARMA模型,或者统计模型较多
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-12-23 01:41 , Processed in 0.052996 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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