Lyapunov 指数求解及预测程序

一个朋友做博士论文，论文里面要用到李雅普诺夫指数做工程预测。我查了很多资料最后才帮他解决了问题，其中花费了很多心血，解决问题的过程让人很痛苦。因为这不单单是编程设计的问题，其中的算法很费解（如果学过混沌的朋友应该会比较容易）。这个程序包里面包含了主要计算的调用示例，供有相似问题的朋友方便使用。

function P=ave_period（data）

year=length（data）;
wolfer=data;
% plot（yearwolfer）
% title（‘ signal‘）

Y = fft（wolfer）;
Y（1）=[];

% plot（Y‘ro‘）
% title（‘Fourier Coefficients in the Complex Plane‘）;
% xlabel（‘Real Axis‘）;
% ylabel（‘Imaginary Axis‘）;

n=length（Y）;
power = abs（Y（1:floor（n/2）））.^2;
nyquist = 1/2;
freq = （1:n/2）/（n/2）*nyquist;
% plot（freqpower）
% xlabel（‘cycles/year‘）
% title（‘Periodogram‘）

period=1./freq;
% plot（periodpower）;
% ylabel（‘Power‘）;
% xlabel（‘Period （Years/Cycle）‘）;

% hold on;
index=find（power==max（power））;
mainPeriodStr=num2str（period（index））;
% plot（period（index）power（index）‘r.‘ ‘MarkerSize‘25）;
% text（period（index）+2power（index）[‘Period = ‘mainPeriodStr]）;
% hold off;

P=round（period（index））

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件     325279  2011-03-07 02:20  使用图示.png

     文件        775  2006-03-02 21:08  Lyapunov预测程序\ave_period.m

     文件      29184  2005-12-12 16:57  Lyapunov预测程序\buffer.dll

     文件        800  2011-03-08 02:09  Lyapunov预测程序\callLargest_lyapunov_exponent.asv

     文件        932  2011-03-08 02:20  Lyapunov预测程序\callLargest_lyapunov_exponent.m

     文件       1076  2011-03-08 01:54  Lyapunov预测程序\CallPre.asv

     文件       1076  2011-03-08 02:28  Lyapunov预测程序\CallPre.m

     文件        829  2011-03-08 02:40  Lyapunov预测程序\callReconstitution.m

     文件      20480  2006-03-01 20:41  Lyapunov预测程序\cao1.dll

     文件        961  2006-03-10 11:24  Lyapunov预测程序\cao_m.m

     文件       3063  2006-03-01 17:38  Lyapunov预测程序\cao_old.m

     文件      20480  2006-01-05 12:14  Lyapunov预测程序\ccFunction0.dll

     文件       1603  2006-04-08 14:30  Lyapunov预测程序\CC_Method.m

     文件      20480  2005-11-25 17:33  Lyapunov预测程序\CorrelationIntegral.dll

     文件       2315  2006-03-06 22:43  Lyapunov预测程序\FunctionChaosPredict.m

     文件       1510  2006-02-16 10:54  Lyapunov预测程序\GP_Algorithm.m

     文件        659  2006-04-02 20:48  Lyapunov预测程序\LargestLyapunov.m

     文件       1744  2011-03-07 01:34  Lyapunov预测程序\largest_lyapunov_exponent.m

     文件       1699  2006-03-08 16:19  Lyapunov预测程序\largest_lyapunov_exponent_revised.m

     文件        228  2006-04-02 20:55  Lyapunov预测程序\lmd1.m

     文件     122368  2005-12-16 22:39  Lyapunov预测程序\lyapunov_buffer.dll

     文件      20480  2006-03-02 20:10  Lyapunov预测程序\Lyapunov_luzhenbo.dll

     文件       2431  2006-02-14 11:25  Lyapunov预测程序\lyapunov_rosenstein.m

     文件      20480  2006-03-09 20:00  Lyapunov预测程序\lyapunov_small.dll

     文件       5574  2011-03-07 01:36  Lyapunov预测程序\lyapunov_wolf.m

     文件      20480  2006-03-09 17:47  Lyapunov预测程序\lyapunov_wolf_dll.dll

     文件        653  2006-03-10 11:45  Lyapunov预测程序\main.m

     文件       1924  2006-03-12 19:29  Lyapunov预测程序\main_process.m

     文件       5391  2006-03-06 22:15  Lyapunov预测程序\mutual.c

     文件      20480  2006-03-10 21:23  Lyapunov预测程序\mutual.dll

............此处省略45个文件信息