#include?“stdafx.h“

#include?“EigTrain.h“
#include?“math.h“

double?*MatrixMul（double?*Mat1int?iRow1int?iCol1double?*Mat2int?iRow2int?iCol2）
{
	double?*Result;
	int?i=0j=0k=0;
	Result=new?double[iRow1*iCol2];
	if（iCol1==iRow2）
	{
		for（i=0;i		{
			for（j=0;j			{
				double?value=0.0;
				for（k=0;k				{
					value+=Mat1[i*iCol1+k]*Mat2[k*iCol2+j];
				}
				*（Result+i*iCol2+j）=value;
			}
		}
	}
	return?Result;
}
double?*MatrixAdd（double?*Mat1int?iRow1int?iCol1double?*Mat2int?iRow2int?iCol2）
{
	double?*Result;
	int?i=0j=0;
	Result=new?double[iRow1*iCol1];
	if（iRow1==iRow2&&iCol1==iCol2）
	{
		for（i=0;i		{
			for（j=0;j				*（Result+i*iCol1+j）=Mat1[i*iCol1+j]+Mat2[i*iCol1+j];
		}
	}
	return?Result;
}
double?*MatrixSub（double?*Mat1int?iRow1int?iCol1double?*Mat2int?iRow2int?iCol2）
{
	double?*Result;
	int?i=0j=0;
	Result=new?double[iRow1*iCol1];
	if（iRow1==iRow2&&iCol1==iCol2）
	{
		for（i=0;i		{
			for（j=0;j				*（Result+i*iCol1+j）=Mat1[i*iCol1+j]-Mat2[i*iCol1+j];
		}
	}
	return?Result;
}
double?*Transpose（double?*Matint?iRowint?iCol）
{
	double?*Result;
	Result=new?double[iCol*iRow];
	int?i=0j=0;
	for（i=0;i		for（j=0;j		{
			*（Result+j*iRow+i）=Mat[i*iCol+j];
		}
		return?Result;
}
double?*getColv（double?*Matint?iRowint?iColint?n）
{
	double?*Result;
	int?i=0;
	Result=new?double[iRow];
	for（i=0;i		*（Result+i）=Mat[i*iCol+n];
	return?Result;

}
double?*Average（double?*dataint?iRowint?iCol）
{
	int?i=0j=0;
	double?*md;
	md=new?double[iRow];
	for（i=0;i		*（md+i）=0.0;
	for（i=0;i	{
		for（j=0;j		{
			*（md+i）+=data[i*iCol+j]/iCol;
		}
	}
	return?md;
}

double?*FcMatrix（double?*dataint?iRowint?iColdouble?*md）
{
	int?i=0;
	double?*Result;
	double?*temp;
	double?*tempm;
	Result=new?double[iRow*iRow];
	for（i=0;i		*（Result+i）=0.0;
	temp=new?double[iRow];
	tempm=new?double[iRow*iRow];
	for（i=0;i	{
		temp=getColv（dataiRowiColi）;
		temp=MatrixSub（tempiRow1mdiRow1）;
		tempm=MatrixMul（tempiRow1temp1iRow）;
		Result=MatrixAdd（ResultiRowiRowtempmiRowiRow）;
	}
	for（i=0;i		*（Result+i）=*（Result+i）/iCol;

	delete?[]temp;
	delete?[]tempm;
	return?Result;
}

void?JacobiEigenv（double?*aint?ndouble?*vdouble?eps）
{
	int?ijpquwts;
????double?fffmcnsnomegaxyd;
????for?（i=0;?i<=n-1;?i++）
	{?
		v[i*n+i]=1.0;
		for?（j=0;?j<=n-1;?j++）
			if?（i!=j）?v[i*n+j]=0.0;
	}
????ff=0.0;
????for?（i=1;?i<=n-1;?i++）
		for?（j=0;?j<=i-1;?j++）
		{?
			d=a[i*n+j];?ff=ff+d*d;
		}
	ff=sqrt（2.0*ff）;
loop0:
	ff=ff/（1.0*n）;
loop1:
	for?（i=1;?i<=n-1;?i++）
		for?（j=0;?j<=i-1;?j++）
		{?
			d=fabs（a[i*n+j]）;
????????????if?（d>ff）
			{?
				p=i;?
				q=j;
				goto?loop;
			}
		}
		if?（ff			return;
		goto?loop0;

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件????????393??2009-05-15?17:48??GaitRsystem\新建?文本文檔.txt

?????文件???????8647??2005-05-13?17:44??GaitRsystem\GaitRsystem\GaitRsystem.rc

?????文件??????70656??2003-12-13?15:30??GaitRsystem\GaitRsystem\zlib1.dll

?????文件?????153966??2004-01-14?03:16??GaitRsystem\GaitRsystem\jpeg62.dll

?????文件???????1930??2005-05-12?14:50??GaitRsystem\GaitRsystem\Resource.h

?????文件??????38200??2005-05-13?17:44??GaitRsystem\GaitRsystem\GaitRsystem.aps

?????文件??????????0??2005-05-14?15:28??GaitRsystem\GaitRsystem\gaitdata.dat

?????文件???????5118??2005-06-08?15:00??GaitRsystem\GaitRsystem\EigTrain.cpp

?????文件???????3518??2005-06-08?17:47??GaitRsystem\GaitRsystem\GaitRsystem.clw

?????文件???????3669??2005-05-12?13:21??GaitRsystem\GaitRsystem\ReadMe.txt

?????文件???????1379??2005-05-12?13:21??GaitRsystem\GaitRsystem\GaitRsystem.h

?????文件???????2133??2005-05-12?13:21??GaitRsystem\GaitRsystem\GaitRsystem.cpp

?????文件???????4769??2005-05-12?16:28??GaitRsystem\GaitRsystem\GaitRsystem.dsp

?????文件???????1717??2005-06-08?15:14??GaitRsystem\GaitRsystem\GaitRsystemDlg.h

?????文件????????213??2005-05-12?13:21??GaitRsystem\GaitRsystem\StdAfx.cpp

?????文件???????6849??2005-06-05?20:44??GaitRsystem\GaitRsystem\Contour.h

?????文件???????1554??2005-05-12?14:52??GaitRsystem\GaitRsystem\GaitDatabaseDlg.h

?????文件????????967??2005-06-03?15:33??GaitRsystem\GaitRsystem\EigTrain.h

?????文件???????5099??2004-02-11?05:36??GaitRsystem\GaitRsystem\libtiff3.def

?????文件?????441498??2004-02-11?05:36??GaitRsystem\GaitRsystem\libtiff3.dll

?????文件?????381922??2004-02-11?05:34??GaitRsystem\GaitRsystem\libtiff.a

?????文件?????115562??2004-02-11?05:36??GaitRsystem\GaitRsystem\libtiff.dll.a

?????文件????????545??2005-05-12?13:21??GaitRsystem\GaitRsystem\GaitRsystem.dsw

?????文件??????40530??2004-02-11?05:36??GaitRsystem\GaitRsystem\libtiff.lib

?????文件???????2203??2005-06-08?19:12??GaitRsystem\GaitRsystem\mophologic.h

?????文件?????132096??2009-05-16?09:52??GaitRsystem\GaitRsystem\GaitRsystem.ncb

?????文件???????3104??2005-06-05?21:32??GaitRsystem\GaitRsystem\Labelcon.h

?????文件??????21504??2004-02-11?05:36??GaitRsystem\GaitRsystem\libtiff-bcc.lib

?????文件???????2681??2005-04-12?21:30??GaitRsystem\GaitRsystem\pgm.cpp

?????文件????????580??1999-05-01?11:52??GaitRsystem\GaitRsystem\pgm.h

............此處省略38個(gè)文件信息