#?include
#?include
#?include?
#?include?
#?include?
#?include
using?namespace?Eigen;
using?namespace?std;

int?main（）
{


	int?m?=?3;
	int?n?=?5;
	double?dl[3][5]?=?{?{?5??0?3?7?4}?{?1?2?3?4?5?}?{?1?2?4?6?2?}?};
	double?data[3][5]?=?{?0?};
	double?datam[3]?=?{?0?};
	double?sb[3][3]?=?{?0?};
	double?sw[5][5]?=?{?0?};
	double?me[3]?=?{0};
	double?mu=0;
	double?meanb?=?0;
	double?meanl?=?0;
	MatrixXf?SW（5?5）;
	MatrixXf?SW_1（5?5）;
	double?*aMat?=?new?double[25];



	double?covzo（double?data_b[3][5]?double?data_a[3][5]double?sw[5][5]）;
	int?i?jjj;

	cout?<	for?（i?=?0;?i?	{
		for?（j?=?0;?j?		{
			meanl?=?meanl?+?dl[i][j];?	
		}
		me[i]?=?meanl?/?n;//sw?mean
		cout?<		meanb?=?meanb?+?meanl;
		meanl?=?0;
	}
	cout?<
	mu?=?meanb?/?（n*m）;//sb?mean
	cout?<	for?（i?=?0;?i?	{
		for?（j?=?0;?j?		{
			data[i][j]?=?dl[i][j]-me[i];//substract?mean
			cout?<		}
		datam[i]?=?me[i]?-?mu;
		cout?<	}
	cout?<	for?（j?=?0;?j?	{
		for?（jj?=?0;?jj?		{
			sb[j][jj]?=?datam[j]?*?datam[jj];
			cout?<
		}
		cout?<
	}
	cout?<	covzo（data?datasw）;
	
	for?（j?=?0;?j?	{
		for?（jj?=?0;?jj?		{
			cout?<			aMat[j*n+jj]?=?sw[j][jj];
		}
		cout?<
	}
	
	Map?dymMat（aMat?5?5）;
	cout?<	SW_1=SW.reverse（）;
	cout?<	


	
	system（“pause“）;
	return?0;
}

double?covzo（double?data_b[3][5]?double?data_a[3][5]double?sw[5][5]）
{//m類n維數據
	
	double?s[5][5]?=?{?0?}?data_bt[5][3];
	int?i?j?jj?k;
	int?n?=?5;
	int?m?=?3;
//	double（*p）[5];

	/*
	p?=?new?double*[n];
	for?（i?=?0;?i?	{
		p[i]?=?new?double[n];
	}
	

	if?（!*p）
	{
		cout?<	}
	*/


	cout?<	for?（k?=?0;?k?	{
		for?（i?=?0;?i?
	{
		data_bt[k][i]?=?data_b[i][k];
		cout?<
	}?
		cout?<	}
		
	cout?<		for?（j?=?0;?j?	{
		for?（jj?=?0;?jj?		{
			for?（i?=?0;?i?			{
				sw[j][jj]?+=?data_bt[j][i]?*?data_a[i][jj];
			}
			cout?<
		}
		cout?<
	}
	//	p?=?s1;
		return?（1）;

	}
	
	/*
	for?（i?=?0;?i?	{
		for?（j?=?0;?j?		{
			for?（jj?=?0;?jj?			{
				s1[i][j]?+=?data_bt[i][jj]?*?data_a[jj][j];

			}
			cout?<		}
		cout?<	}
		for?（k?=?0;?k?		{
			for?（i?=?0;?i?
			{
				s[k][i]?=?s[k][i]+s1[k][i];
				cout?<			}
			cout?<		}
		*/

	/*free（s1）;
		free（s）;
	free（data_bt）;
	free（i）;
	free（j）;
		free（jj）;
		free（k）;
	fre

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

?????文件?????396456??2017-03-26?17:53??PCAtest\Debug\ldaz.obj

?????文件????6524333??2017-03-31?11:11??PCAtest\Debug\LDA_Z.obj

?????文件??????31668??2017-03-27?18:25??PCAtest\Debug\pca2.obj

?????文件???????1667??2017-03-31?12:28??PCAtest\Debug\PCAtest.log

?????文件???????4742??2017-03-31?12:28??PCAtest\Debug\PCAtest.tlog\cl.command.1.tlog

?????文件?????189112??2017-03-31?12:28??PCAtest\Debug\PCAtest.tlog\CL.read.1.tlog

?????文件???????7180??2017-03-31?12:28??PCAtest\Debug\PCAtest.tlog\CL.write.1.tlog

?????文件???????3776??2017-03-31?12:28??PCAtest\Debug\PCAtest.tlog\link.command.1.tlog

?????文件???????2852??2017-03-31?12:28??PCAtest\Debug\PCAtest.tlog\link.read.1.tlog

?????文件????????606??2017-03-31?12:28??PCAtest\Debug\PCAtest.tlog\link.write.1.tlog

?????文件????????188??2017-03-31?12:28??PCAtest\Debug\PCAtest.tlog\PCAtest.lastbuildstate

?????文件????6976660??2017-03-31?12:28??PCAtest\Debug\pcaz.obj

?????文件?????356514??2017-03-30?21:15??PCAtest\Debug\PCA_Z.obj

?????文件????4402176??2017-03-31?12:28??PCAtest\Debug\vc120.idb

?????文件???26669056??2017-03-31?12:28??PCAtest\Debug\vc120.pdb

?????文件???????3031??2017-03-26?17:52??PCAtest\ldaz.cpp

?????文件???????5421??2017-03-27?22:31??PCAtest\LDA_Z.cpp

?????文件???????3491??2017-03-31?12:02??PCAtest\PCAtest.vcxproj

?????文件????????945??2017-03-31?12:02??PCAtest\PCAtest.vcxproj.filters

?????文件????????165??2016-12-29?20:12??PCAtest\PCAtest.vcxproj.user

?????文件???????4150??2017-03-31?12:28??PCAtest\pcaz.cpp

?????文件???????1976??2017-03-30?21:15??PCAtest\PCA_Z.cpp

?????目錄??????????0??2017-03-31?15:18??PCAtest\Debug\PCAtest.tlog

?????目錄??????????0??2017-03-31?15:18??PCAtest\Debug

?????目錄??????????0??2017-03-31?15:20??PCAtest

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

?????????????45586165????????????????????25