声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1940|回复: 2

[分形与混沌] 混沌工具箱的使用问题

[复制链接]
发表于 2008-11-18 20:47 | 显示全部楼层 |阅读模式

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

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

x
根据混沌工具箱程序: Lorenz=[[0.539589905
0.62942492
0.541100324
0.364385692
0.726098655
0.76503888
0.350448718
0.345463259
0.417826772
0.059067524
0.315777414
0.125562701
0.074241486
0.055517241
0.027329193
0.203525741
0.105401929
0.285198098
0.336062992
0.038396226
0.096384365
0.521884984
0.598346095
0.216837061
0.417821138
0.535605816
0.378103976
0.402092308
0.349152542
0.454915254
0.079970326
0.107425897
0.092779553
0.013129032
0.03265625
0.030870968
0.213841369
0.331534884
0.061875
0.027281399
0.023011094
0.037022654
0.084201954
0.126791277
0.182333333
0.009369085
0.114217252
0.055177994
0.160795107
0.078402458
0.047594554
0.07550702
0.000695103
0.016525822
0.002444444
0.000890688
0.221748808
0.052956381
0.057693523
0.010748092
0.020357143
0.019877489
0.01858679
0
0.466043614
0.215692771
0.421886228
0.458165138
0.29807571
0.039578606
0.151855346
0.39407767
0.333971743
0.463621795
0.381263823
0.279450382
0.460241327
0.538865031
0.020157978
0.102044025
0.105794118
1.145762712
0.002376543
0.361980831
0.584523077
0.703717949
0.715810147
0.5575625
0.549822866
0.219004525
0.499684543
0.320561622
0.71699842
0.409233083
0.379023622
0.645356577
0.498954248
0.371644042
0.141601256
0.335968992
0.415959119
0.194791667
0.528278041
0.619620253
0.478639241
0.469861963
0.736744186
0.37488408
0.313923077
0.36778626
0.765107362
1.338046875
0.779876161
0.532594937
0.375905512
0.75754717
0.707018927
0.518203125
0.604914197
0.098430141
0.073971519
0.689251592
0.42301092
0.748034591
0.46257716
0.195216049
0.730172414
0.388286164
0.787697161
0.51734375
0.480708661
0.61661743
0.487419872
0.622769231
0.627959502
1.079816514
0.525320513
0.653773585
0.863063764
0.456055901
0.500629921
0.271518987
0.484341085
0.190517241
0.534260731
0.733623418
0.32982866
0.832638889
0.424845201
0.426530612
0.61072
0.399919743
0.564593698
0.34546875
0.519968799
0.122705696
0.090677966
0.727639752
0.739921875
0.513276398
0.599761905
0.264285714
0.606807512
0.306877898
0.528463855
0.702580645
0.394937695
0.560610932
0.61908646
0.611399371
0.632242991
0.437114198
0.506418384
0.377315542
0.324247227
0.220518868
0.686862442
0.4503125
0.207098765
0.435669291
0.817840735
0.306340289
0.571187801
0.472203947
0.499516908
0.268103448
0.20485133
0.540546875
0.177586207
0.458623418
0.06637931
0.624882629
0.41271722
0.20312989
0.429508197
0.618223583
0.694332298
0.49984544
0.25851735
0.370077519
0.447196262
0.049292453
0.6
0.514544025
0.509904913
0.530141287
0.787068966
0.248191824
0.502570481
1.364783282
0.69609375
0.60911215
0.605345912
0.451415094
0.485396825
0.653178295
0.727308294
0.190517241
0.504453125
0.592109375
0.349605678
0.582810047
0.491221374
0.681317829
1.587917329
0.566936791
0.537893082
0.670175439
0.776886035
0.407636656
0.520504732
0.130870918
0.479731861
0.15628882
0.086124031
0.546540881
0.295689655
0.629188482
0.271978022
0.12433281
0.53436019
0.838012618
0.538828125
0.745023328
0.555205047
0.480926217
0.503521127
0.037185535
0.622275641
0.717429022
0.351267829
0.701724138
0.664330709
0.5225
0.565492958
0.479769231
0.490945513
0.750319489
0.697645212
0.532972136
0.488309637
0.415789474
0.449249249
1.223172628
0.822413793
0.504675119
0.53351735
0.03514377
0.605259027
0.30399361
0.17700157
0.477406832
0.431338583
0.365207006
0.342679128
0.109289176
0.415789474
0.416428571
0.303834116
0.357032457
0.090944882
0.355477308
0.577760252
0.173865415
0.124643423
0.559572785
0.201811024
0.337380192
0.154213938
0.373828125
0.36607717
0.535298742
0.16576507
0.308241758
0.679462875
0.471178344
0.193961353
0.142879499
0.088932496
0.407120501
0.196056782
0.278426791
];
% 时间序列
% dx/dt = sigma*(y-x)
% dy/dt = r*x - y - x*z
% dz/dt = -b*z + x*y
sigma = 16;             % Lorenz方程参数
r = 45.92;               
b = 4;         
y = [-1;0;1];           % 起始点 (3x1 的列向量)
h = 0.01;               % 积分时间步长
k1 = 30000;             % 前面的迭代点数
k2 = 5000;              % 后面的迭代点数
X = LorenzData(y,h,k1+k2,sigma,r,b);
X = X(k1+1:end,1);      % 时间序列(列向量)
%--------------------------------------------------------------------------
% G-P算法计算关联维
rr = 0.5;
Log2R = -6:rr:0;        % log2(r)
R = 2.^(Log2R);
t = 10;                 % 时延
dd = 1;                 % 嵌入维间隔
D = 2:dd:10;            % 嵌入维   
p = 50;                 % 限制短暂分离,大于序列平均周期(不考虑该因素时 p = 1)   
tic
Log2Cr = log2(CorrelationIntegral(X,t,D,R,p));      %Log2Cr
toc                     % 输出每一行对应一个嵌入维


请问:我用matlab计算过程中老是出现这种情况:“??? Undefined function or variable 'Log2Cr'.
Error in ==> C:\MATLAB6p5p1\work\Main_CorrelationDimension_GP.m
On line 50  ==> plot(Log2R,Log2Cr','k.-'); axis tight; grid on; hold on;”
是怎么回事啊?请各位帮忙分析下啊,万分感谢!
回复
分享到:

使用道具 举报

发表于 2008-11-19 08:05 | 显示全部楼层
log2(CorrelationIntegral(X,t,D,R,p));      这个函数是需要你自己定义的,定义后在能在这个程序中引用,你可以新建一个程序,用function函数定义上面的功能函数,文件名默认成.m格式,这里就可以调用了
至于plot(Log2R,Log2Cr','k.-'); axis tight; grid on; hold on;”我想是因为其中包括了Log2Cr。因为你没能求出这个,所以这里也提示你错误了吧,能力有限,希望帮的上忙
发表于 2008-12-23 09:45 | 显示全部楼层
Chaos toolbox V2.0发我一份行吗??我下载不了……
wujiezhou@163.com
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-12-25 13:36 , Processed in 0.093979 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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