//?GPKolmogorov.cpp?:?Defines?the?entry?point?for?the?console?application.
//

#include?“stdafx.h“
#include?
#include?“Henon.h“
#include?“HyperRossler.h“
#include?“Lorenz.h“
#include?“GPKolmogorov.h“

#define?MAXLEN?10000
#define?MLEN?12
#define?RLEN?50


int?main（int?argc?char*?argv[]）
{
	double?Seq[MAXLEN]={0};
	//int?Henon_len=3000;
	//int?Lorenz_len=3000;
	int?HyperRossler_len=3000;
	//double?*Henon_X=new?double[Henon_len];
	//double?*Henon_Y=new?double[Henon_len];
	//Henon映射初值
	//double?Henon_x0=1;
	//double?Henon_y0=0;
	//lorenz映射初值
	//double?Lorenz_init[3]?=?{151337};
	//double?*Lorenz_X=new?double[Lorenz_len];
	//double?*Lorenz_Y=new?double[Lorenz_len];
	//double?*Lorenz_Z=new?double[Lorenz_len];
	//超Rossler映射初值
	double?HyperRossler_init[4]?=?{-200015};
????double?*HyperRossler_X=new?double[HyperRossler_len];
	double?*HyperRossler_Y=new?double[HyperRossler_len];
	double?*HyperRossler_Z=new?double[HyperRossler_len];
	double?*HyperRossler_W=new?double[HyperRossler_len];
	//嵌入維
	int?m[MLEN];
	//延時
	int?delay=15;
	double?r[RLEN];
	//關(guān)聯(lián)維輸出矩陣
	double?ln_C[MLEN][RLEN];
	double?*lncc[MLEN];
	double?ln_r[RLEN];
	double?mean_x[MLEN]={0};
	double?mean_y[MLEN]={0};
	double?lxx[MLEN]={0};
	double?lxy[MLEN]={0};
	double?sum_lxx=0;
	double?sum_lxy=0;
	double?a_D;
	double?b_D[MLEN];
	double?Kolmogorov[MLEN-1];?//得到的Kolmogorov矩陣
????int?ij;

	for?（i=0;?i	{
		lncc[i]?=?ln_C[i];
	}
????//定義嵌入維矩陣
	for?（i=0;?i	{
		m[i]?=?8+2*i;
	}
	for?（i=0;?i	{
		r[i]?=?3+2*i;
	}

	//Henon2（Henon_XHenon_YHenon_x0Henon_y0Henon_len）;
	//lorenz（Lorenz_XLorenz_YLorenz_ZLorenz_initLorenz_len）;
	hyperrossler（HyperRossler_XHyperRossler_Y
		HyperRossler_ZHyperRossler_WHyperRossler_initHyperRossler_len）;
/*	for?（i=0;?i	{
		cout<		if?（i%7==6）
			cout<	}*/
?/*???for?（i=0;?i	{
		cout<		if?（i%7==6）
			cout<	}*/
????for?（i=0;?i	{
		cout<		if?（i%7==6）
			cout<	}
	GP_Kol（HyperRossler_XrHyperRossler_lenmRLENMLENdelaylnccln_r）;
????for?（i=0;?i	{
		for?（j=0;?j			cout<		if?（j%7==6）
			cout<	}
	//計算xy的均值
	for?（i=0;?i	{
		for?（j=0;?j		{
			mean_y[i]?+=?ln_C[i][j];
			mean_x[i]?+=?ln_r[j];
		}
		mean_y[i]?/=?RLEN;
????????mean_x[i]?/=?RLEN;
	}
	//計算lxylxx的值及其和．
	for?（i=0;?i	{
		for?（j=0;?j		{
			lxy[i]?+=?（ln_C[i][j]-mean_y[i]）*（ln_r[j]-mean_x[i]）;
			lxx[i]?+=?（ln_r[j]-mean_x[i]）*（ln_r[j]-mean_x[i]）;
		}
		sum_lxy?+=?lxy[i];
		sum_lxx?+=?lxx[i];
	}
	a_D?=?sum_lxy/sum_lxx;
	cout<<“嵌入維為：“<	for?（i=0;?i	{
		b_D[i]?=?mean_y[i]-a_D*mean_x[i];
	}
	//得到Kolmogorov熵
	cout<<“Kolmogorov熵為：“<	for?（i=0;?i	{
		Kolmogorov[i

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????3377??2004-11-15?20:01??同時計算關(guān)聯(lián)維和Kolmogorov熵\GPKolmogorov\GPKolmogorov.cpp

?????文件???????4933??2004-11-12?10:24??同時計算關(guān)聯(lián)維和Kolmogorov熵\GPKolmogorov\GPKolmogorov.dsp

?????文件????????549??2004-11-12?08:59??同時計算關(guān)聯(lián)維和Kolmogorov熵\GPKolmogorov\GPKolmogorov.dsw

?????文件???????1116??2004-11-15?20:02??同時計算關(guān)聯(lián)維和Kolmogorov熵\GPKolmogorov\GPKolmogorov.h

?????文件??????50176??2004-11-15?20:03??同時計算關(guān)聯(lián)維和Kolmogorov熵\GPKolmogorov\GPKolmogorov.ncb

?????文件??????96256??2004-11-15?20:03??同時計算關(guān)聯(lián)維和Kolmogorov熵\GPKolmogorov\GPKolmogorov.opt

?????文件???????1422??2004-11-15?08:23??同時計算關(guān)聯(lián)維和Kolmogorov熵\GPKolmogorov\GPKolmogorov.plg

?????文件????????703??2004-11-12?10:01??同時計算關(guān)聯(lián)維和Kolmogorov熵\GPKolmogorov\Henon.h

?????文件???????2883??2004-11-12?09:36??同時計算關(guān)聯(lián)維和Kolmogorov熵\GPKolmogorov\HyperRossler.h

?????文件???????2557??2004-11-12?10:14??同時計算關(guān)聯(lián)維和Kolmogorov熵\GPKolmogorov\Lorenz.h

?????文件???????1244??2004-11-12?08:59??同時計算關(guān)聯(lián)維和Kolmogorov熵\GPKolmogorov\ReadMe.txt

?????文件????????299??2004-11-12?08:59??同時計算關(guān)聯(lián)維和Kolmogorov熵\GPKolmogorov\StdAfx.cpp

?????文件????????667??2004-11-12?08:59??同時計算關(guān)聯(lián)維和Kolmogorov熵\GPKolmogorov\StdAfx.h

?????目錄??????????0??2011-03-14?09:35??同時計算關(guān)聯(lián)維和Kolmogorov熵\GPKolmogorov

?????目錄??????????0??2011-03-14?09:35??同時計算關(guān)聯(lián)維和Kolmogorov熵

-----------?---------??----------?-----??----

???????????????166182????????????????????15