混沌时间序列的预测误差

原来是小飞

用多步预测方法预测混沌时间序列，但是预测值是哪个数据啊？好心人帮个忙吧，以下是代码：

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-');
    
38. 
    
39. 
    
40. 
    
41. 
    

复制代码

gghhjj

从代码上看 PredictedData+mean(oldroll) 应该是程序给出的预测值

原来是小飞

gghhjj 发表于 2014-3-17 14:45
从代码上看 PredictedData+mean(oldroll) 应该是程序给出的预测值

谢谢，您讲的是正确的。还想问您下，风速时间序列一般用什么方法重构相空间啊

gghhjj

原来是小飞 发表于 2014-3-17 16:39
谢谢，您讲的是正确的。还想问您下，风速时间序列一般用什么方法重构相空间啊

相空间重构方法很多，常用的有伪最近邻点法、奇异值分解法、自相关和互信息法等。但都存在这样或者那样的一些问题，个人对风速时间序列的特征不太了解，不知道具体什么方法更好，这需要你结合时间序列的特征以及你所掌握的的数据情况进行选择尝试，这方面的文献应该也是有的，你可以参考参考

ABC2014

风功率呢  我是初学者 不太懂  能耽误你宝贵的时间给我讲解一下吗

gghhjj

ABC2014 发表于 2014-4-6 13:41
风功率呢  我是初学者 不太懂  能耽误你宝贵的时间给我讲解一下吗

两则应该是类似的。上面讲的是相空间重构
如果是预测的话，现在一般用神经网络，ARMA模型，或者统计模型较多