|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
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;”是怎么回事啊?请各位帮忙分析下啊,万分感谢! |
|