200532590150.cpp

实现了二维数据的K均值算法以及ISODATA分类算法

// 200532590150.cpp : Defines the class behaviors for the application.
//

#include "stdafx.h"
#include "200532590150.h"

#include "MainFrm.h"
#include "ChildFrm.h"
#include "200532590150Doc.h"
#include "200532590150View.h"
#include "math.h"

#ifdef _DEBUG
#define new DEBUG_NEW
#undef THIS_FILE
static char THIS_FILE[] = __FILE__;
#endif

/////////////////////////////////////////////////////////////////////////////
// CMy200532590150App

BEGIN_MESSAGE_MAP(CMy200532590150App, CWinApp)
	//{{AFX_MSG_MAP(CMy200532590150App)
	ON_COMMAND(ID_APP_ABOUT, OnAppAbout)
		// NOTE - the ClassWizard will add and remove mapping macros here.
		//    DO NOT EDIT what you see in these blocks of generated code!
	//}}AFX_MSG_MAP
	// Standard file based document commands
	ON_COMMAND(ID_FILE_NEW, CWinApp::OnFileNew)
	ON_COMMAND(ID_FILE_OPEN, CWinApp::OnFileOpen)
	// Standard print setup command
	ON_COMMAND(ID_FILE_PRINT_SETUP, CWinApp::OnFilePrintSetup)
END_MESSAGE_MAP()

/////////////////////////////////////////////////////////////////////////////
// CMy200532590150App construction

CMy200532590150App::CMy200532590150App()
{
	// TODO: add construction code here,
	// Place all significant initialization in InitInstance
}

/////////////////////////////////////////////////////////////////////////////
// The one and only CMy200532590150App object

CMy200532590150App theApp;

/////////////////////////////////////////////////////////////////////////////
// CMy200532590150App initialization

BOOL CMy200532590150App::InitInstance()
{
	AfxEnableControlContainer();

	// Standard initialization
	// If you are not using these features and wish to reduce the size
	//  of your final executable, you should remove from the following
	//  the specific initialization routines you do not need.

#ifdef _AFXDLL
	Enable3dControls();			// Call this when using MFC in a shared DLL
#else
	Enable3dControlsStatic();	// Call this when linking to MFC statically
#endif

	// Change the registry key under which our settings are stored.
	// TODO: You should modify this string to be something appropriate
	// such as the name of your company or organization.
	SetRegistryKey(_T("Local AppWizard-Generated Applications"));

	LoadStdProfileSettings();  // Load standard INI file options (including MRU)

	// Register the application's document templates.  Document templates
	//  serve as the connection between documents, frame windows and views.

	CMultiDocTemplate* pDocTemplate;
	pDocTemplate = new CMultiDocTemplate(
		IDR_MY2005TYPE,
		RUNTIME_CLASS(CMy200532590150Doc),
		RUNTIME_CLASS(CChildFrame), // custom MDI child frame
		RUNTIME_CLASS(CMy200532590150View));
	AddDocTemplate(pDocTemplate);

	// create main MDI Frame window
	CMainFrame* pMainFrame = new CMainFrame;
	if (!pMainFrame->LoadFrame(IDR_MAINFRAME))
		return FALSE;
	m_pMainWnd = pMainFrame;

	// Parse command line for standard shell commands, DDE, file open
	CCommandLineInfo cmdInfo;
	ParseCommandLine(cmdInfo);

	// Dispatch commands specified on the command line
	if (!ProcessShellCommand(cmdInfo))
		return FALSE;

	// The main window has been initialized, so show and update it.
	pMainFrame->ShowWindow(m_nCmdShow);
	pMainFrame->UpdateWindow();

	return TRUE;
}


/////////////////////////////////////////////////////////////////////////////
// CAboutDlg dialog used for App About

class CAboutDlg : public CDialog
{
public:
	CAboutDlg();

// Dialog Data
	//{{AFX_DATA(CAboutDlg)
	enum { IDD = IDD_ABOUTBOX };
	//}}AFX_DATA

	// ClassWizard generated virtual function overrides
	//{{AFX_VIRTUAL(CAboutDlg)
	protected:
	virtual void DoDataExchange(CDataExchange* pDX);    // DDX/DDV support
	//}}AFX_VIRTUAL

// Implementation
protected:
	//{{AFX_MSG(CAboutDlg)
		// No message handlers
	//}}AFX_MSG
	DECLARE_MESSAGE_MAP()
};

CAboutDlg::CAboutDlg() : CDialog(CAboutDlg::IDD)
{
	//{{AFX_DATA_INIT(CAboutDlg)
	//}}AFX_DATA_INIT
}

void CAboutDlg::DoDataExchange(CDataExchange* pDX)
{
	CDialog::DoDataExchange(pDX);
	//{{AFX_DATA_MAP(CAboutDlg)
	//}}AFX_DATA_MAP
}

BEGIN_MESSAGE_MAP(CAboutDlg, CDialog)
	//{{AFX_MSG_MAP(CAboutDlg)
		// No message handlers
	//}}AFX_MSG_MAP
END_MESSAGE_MAP()

// App command to run the dialog
void CMy200532590150App::OnAppAbout()
{
	CAboutDlg aboutDlg;
	aboutDlg.DoModal();
}

/////////////////////////////////////////////////////////////////////////////
// CMy200532590150App message handlers

//以最小距离原则分类：计算各样本到各聚类中心的距离并按最小距离分类存放在各聚类数组中
void classify(double x1[],double x2[],double z[][2],int num[],double y1[][2],double y2[][2],double y3[][2],int count)
{
	double temp[100][3],t;
	int a=0,b=0,c=0;
	for(int i=0;i<count;i++)
	{
		//分别计算各样本到各聚类中心的距离
		for(int j=0;j<3;j++)
		{
			t=(x1[i]-z[j][0])*(x1[i]-z[j][0])+(x2[i]-z[j][1])*(x2[i]-z[j][1]);
			temp[i][j]=sqrt(t);
		}
		//寻找样本到三个聚类中心的距离的最小值
		t=min(temp[i][0],temp[i][1]);
		t=min(t,temp[i][2]);
		//将样本归属于距离最近的聚类中心
        if(t==temp[i][0])
		{
			y1[a][0]=x1[i];
			y1[a][1]=x2[i];
			a++;//计算新类群中的样本个数
		}
		else if(t==temp[i][1])
		{
			y2[b][0]=x1[i];
			y2[b][1]=x2[i];
			b++;
		}
		else if(t==temp[i][2])
		{
			y3[c][0]=x1[i];
			y3[c][1]=x2[i];
			c++;
		}
		num[0]=a;
		num[1]=b;
		num[2]=c;
	}
}

//计算聚类中心
void center(double y1[][2],double y2[][2],double y3[][2],double z1[3][2],int num[])
{
	//分别计算各类群的均值向量
	for(int i=0;i<num[0];i++)
	{
		z1[0][0]+=y1[i][0]/num[0];
		z1[0][1]+=y1[i][1]/num[0];
	}
	for(i=0;i<num[1];i++)
	{
		z1[1][0]+=y2[i][0]/num[1];
		z1[1][1]+=y2[i][1]/num[1];
	}
	for(i=0;i<num[2];i++)
	{
		z1[2][0]+=y3[i][0]/num[2];
		z1[2][1]+=y3[i][1]/num[2];
	}
}