同时计算关联维和Kolmogorov熵

同时计算时间序列的关联维数及Kolmogorov熵

// 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];
	//关联维输出矩阵
	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  同时计算关联维和Kolmogorov熵\GPKolmogorov\GPKolmogorov.cpp

     文件       4933  2004-11-12 10:24  同时计算关联维和Kolmogorov熵\GPKolmogorov\GPKolmogorov.dsp

     文件        549  2004-11-12 08:59  同时计算关联维和Kolmogorov熵\GPKolmogorov\GPKolmogorov.dsw

     文件       1116  2004-11-15 20:02  同时计算关联维和Kolmogorov熵\GPKolmogorov\GPKolmogorov.h

     文件      50176  2004-11-15 20:03  同时计算关联维和Kolmogorov熵\GPKolmogorov\GPKolmogorov.ncb

     文件      96256  2004-11-15 20:03  同时计算关联维和Kolmogorov熵\GPKolmogorov\GPKolmogorov.opt

     文件       1422  2004-11-15 08:23  同时计算关联维和Kolmogorov熵\GPKolmogorov\GPKolmogorov.plg

     文件        703  2004-11-12 10:01  同时计算关联维和Kolmogorov熵\GPKolmogorov\Henon.h

     文件       2883  2004-11-12 09:36  同时计算关联维和Kolmogorov熵\GPKolmogorov\HyperRossler.h

     文件       2557  2004-11-12 10:14  同时计算关联维和Kolmogorov熵\GPKolmogorov\Lorenz.h

     文件       1244  2004-11-12 08:59  同时计算关联维和Kolmogorov熵\GPKolmogorov\ReadMe.txt

     文件        299  2004-11-12 08:59  同时计算关联维和Kolmogorov熵\GPKolmogorov\StdAfx.cpp

     文件        667  2004-11-12 08:59  同时计算关联维和Kolmogorov熵\GPKolmogorov\StdAfx.h

     目录          0  2011-03-14 09:35  同时计算关联维和Kolmogorov熵\GPKolmogorov

     目录          0  2011-03-14 09:35  同时计算关联维和Kolmogorov熵

----------- ---------  ---------- -----  ----

               166182                    15