此程序是应用mex文件来实现潜入维与关联维的计算。.doc

此程序是应用mex 文件来实现潜入维和关联维的计算。 #include mex.h #include math.h ************************************************ * data:the time series * r: the r series * N: the length of the time series * delay : the delay of the time series * rlen: the length of r series * ln_C: the correlation dimention series ************************************************** void G_P(double data[], double r[], int N, int m, int rlen, int delay, double ln_C[],double ln_r[]) { double distance; /*定义相空间中任意两点之间的距离*/ double max=0; /*初始化两点间最大距离为0*/ int cr=0; /*Heaviside 函数计数初始化为0*/ int i,j,k,l; double Cr=0; int M=N-(m-1)*delay; /*相空间点的个数*/ for (l=0; lrlen; l++) /*r 系列的依次取值*/ { for (i=0; iM-1; i++) /*任意两点间距离*/ { for (j=i+1; jM; j++) { for (k=0; km; k++) { distance = fabs(data[i+k*delay]-data[j+k*delay]); if (distancemax) max = distance; } if (r[l]=max) cr++; /*Heaviside 函数取值为1*/ max = 0; } } Cr = cr/(M*(M-1));/*关联积分C2(r,m) cr = 0; ln_C[l] = log(Cr); ln_r[l] = log(r[l]); } } void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]) /*mex 文件入口函数是mexFunction 函数 { double *data; double *r; int m; int delay; double *ln_C; double *ln_r; int datalen,rlen; /*检查正确的参数数目*/ if (nrhs!=4) mexErrMsgTxt(Four inputs required!); else if (nlhs2) mexErrMsgTxt(Too many output arguments!); datalen = mxGetN(prhs[0]); rlen = mxGetN(prhs[1]); m = mxGetScalar(prhs[2]); delay = mxGetScalar(prhs[3]); /*定义两个输出矩阵，以便返回参数*/ plhs[0] = mxCreateDoubleMatrix(1,rlen,mxREAL); plhs[1] = mxCreateDoubleMatrix(1,rlen,mxREAL); /*分配输入输出参数的指针*/ data = mxGetPr(prhs[0]); r = mxGetPr(prhs[1]); ln_C = mxGetPr(plhs[0]); ln_r = mxGetPr(plhs[1]); /*调用G_P 子函数*/ G_P(data,r,datalen,m,rlen,delay,ln_C,ln_r);} 附录3： 此程序实现了时间延迟参数的确定。 #define N 1000 #define MAX 200 #include stdio.h float ave(float *a,int n) { int i; float ave1; long double sum=0; for (i=0;in;i++) sum+=a[i]; printf(%f\n,sum); ave1=sum/n; return ave1; } float autol(float *y,int t,int n,float ave) { float up=0,dowm=0,s=0,c; int i; for(i=0;in;i++) s=(y[i]-ave)*(y[i]-ave); dowm=s/n; for