有人用过陆振波博士求时间序列关联维数和熵的程序么？

warden

RT,我最近做时间序列的分析，在网上找的陆振波博士的计算时间序列关联维数和熵的程序，对时间序列进行HHT滤波处理之后，导入陆博士的程序，有些得到的是nan和0.00000，这些说明什么情况呢？有人做过吗？求教我一下

下面的是主程序：

1. clc
   
2. clear all
   
3. close all
   
4. x=importdata('D:\Documents\Desktop\新的要求\10-30%-5.mat'); % 时间序列(列向量)
   
5. x=x';
   
6. %-----------------------------------  ---------------------------------------
   
7. % G-P算法计算关联维
   
8. k=length(x);
   
9. rr = 0.5;
   
10. Log2R = -6:rr:0;        % log2(r)
    
11. R = 2.^(Log2R);
    
12. fs =100;               % 信号采样频率
    
13. t = 10;                 % 时延
    
14. dd = 4;                 % 嵌入维间隔
    
15. D = 4:dd:72;            % 嵌入维   
    
16. p = 50;                 % 限制短暂分离，大于序列平均周期(不考虑该因素时 p = 1)  
    
17. 
    
18. tic
    
19. Log2Cr = log2(CorrelationIntegral(x,t,D,R,p));   % 输出每一行对应一个嵌入维
    
20. toc
    
21. %--------------------------------------------------------------------------
    
22. % 结果作图
    
23. figure
    
24. plot(Log2R,Log2Cr','k.-'); axis tight; hold on; grid on;
    
25. xlabel('log2(r)');
    
26. ylabel('log2(C(r))');
    
27. title(['signal, length = ',num2str(k)]);
    
28. %--------------------------------------------------------------------------
    
29. % 最小二乘拟合
    
30. Linear = [6:10];                            % 线性似合区域
    
31. [a,B] = LM2(Log2R,Log2Cr,Linear);           % 最小二乘求斜率和截距
    
32. 
    
33. disp(sprintf('Correlation Dimension = %.4f',a));
    
34. for i = 1:length(D)
    
35.     Y = polyval([a,B(i)],Log2R(Linear),1);
    
36.     plot(Log2R(Linear),Y,'r');
    
37. end
    
38. hold off;
    
39. %--------------------------------------------------------------------------
    
40. % 求梯度
    
41. 
    
42. Slope = diff(Log2Cr,1,2)/rr;                % 求梯度
    
43. xSlope = Log2R(1:end-1);                    % 梯度所对应的log2(r)
    
44. 
    
45. figure;
    
46. plot(xSlope,Slope','k.-'); axis tight; grid on;
    
47. xlabel('log2(r)');
    
48. ylabel('Slope');
    
49. title(['signal, length = ',num2str(k)]);
    
50. 
    
51. %--------------------------------------------------------------------------
    
52. % 差分求K熵
    
53. 
    
54. KE = -diff(B)/(dd*t)*fs*log(2);             % 用采样频率 fs 和公式 log(x) = log2(x)*log(2) 将单位转化成 nats/s
    
55. D_KE = D(1:end-1);                          % K熵所对应的嵌入维
    
56. 
    
57. figure;
    
58. plot(D_KE,KE,'k.-'); grid on; hold on;
    
59. xlabel('m');
    
60. ylabel('Kolmogorov Entropy (nats/s)');
    
61. title(['signal, length = ',num2str(k)]);
    
62. 
    
63. %--------------------------------------------------------------------------
    
64. % 输出显示
    
65. 
    
66. disp(sprintf('Kolmogorov Entropy = %.4f',min(KE)));

复制代码

处理一组数据后得到的图形如下：

happy

没有数据，没有代码，没有说明输出是什么，这个问题无从回答
建议看看http://forum.vibunion.com/thread-36219-1-1.html
把问题说明清楚

章紫风

我试过之后感觉出来的Kolmogorov Entropy图不对，结果并不收敛！！！

happy

章紫风 发表于 2012-12-3 18:54

                               
登录/注册后可看大图

我试过之后感觉出来的Kolmogorov Entropy图不对，结果并不收敛！！！

你的时间序列是仿真结果还是实验结果？

章紫风

happy 发表于 2012-12-4 09:40

                               
登录/注册后可看大图

你的时间序列是仿真结果还是实验结果？

仿真结果，就是Lorenz模型，Kolmogorov Entropy随着m值变化，结果不收敛，图像是先下降然后居然又上升。还有那个程序有使用次数限制吗？

happy

章紫风 发表于 2012-12-5 20:51

                               
登录/注册后可看大图

仿真结果，就是Lorenz模型，Kolmogorov Entropy随着m值变化，结果不收敛，图像是先下降然后居然又上升。还 ...

1. 检查一下仿真得到的时间序列是否已经提出了瞬态响应数据
2. 没听说过他的工具包有这方面的限制

章紫风

happy 发表于 2012-12-6 10:58

                               
登录/注册后可看大图

1. 检查一下仿真得到的时间序列是否已经提出了瞬态响应数据
2. 没听说过他的工具包有这方面的限制

非常感谢，我的问题差不多解决了！！！

happy

章紫风 发表于 2012-12-6 18:32

                               
登录/注册后可看大图

非常感谢，我的问题差不多解决了！！！

建议解决问题后把相关的情况和大家分享一下

warden

章紫风 发表于 2012-12-6 18:32

                               
登录/注册后可看大图

非常感谢，我的问题差不多解决了！！！

同求解决，不知道能否交流一下？

高永强

我想请问相关的问题

牛小贱

章紫风 发表于 2012-12-6 18:32
非常感谢，我的问题差不多解决了！！！

希望你分享一下你的经验，方便大家学习！谢了

liuxiaoxu2011

你好，我是中国矿业大学的研究生，我在网上看到你发的一篇帖子，我对里面的内容很感兴趣。不知道你现在还有没有涉及柯尔莫戈洛夫熵值的领域。我想和你就我们共同感兴趣的地方交流一下。我的QQ是652247469.如果你方便的话请加我QQ好友。谢谢！

沐雨柠檬

31行程序的LM2是什么？子函数么？matlab没这个命令。如果是子函数，那是求什么的？子函数具体什么样？求解答！