unit1.cpp

此程序采用K均值的方法对样本数据进行模式分类

//---------------------------------------------------------------------------

#include <vcl.h>
#pragma hdrstop

#include "Unit1.h"
//---------------------------------------------------------------------------
#pragma package(smart_init)
#pragma resource "*.dfm"
TForm1 *Form1;
int IterNum;
//---------------------------------------------------------------------------
__fastcall TForm1::TForm1(TComponent* Owner)
        : TForm(Owner)
{
}
//---------------------------------------------------------------------------

void __fastcall TForm1::FormCreate(TObject *Sender)
{
  SpeedButton2->Enabled=false;
  for(int i=0;i<150;i++)
      for(int j=0;j<4;j++)
         Pattern[i][j]=0;
  //设置载入文件初始目录
  String InitDir=GetCurrentDir();
  OpenDialog1->InitialDir=InitDir;
}
//---------------------------------------------------------------------------

void __fastcall TForm1::SpeedButton1Click(TObject *Sender)
{
  ListView1->Items->Clear();//列表框清空
  ListView2->Items->Clear();
  ListView3->Items->Clear();
  ListView4->Items->Clear();

  Label1->Caption="共计  0  个样本";
  Label2->Caption="共计  0  个样本";
  Label3->Caption="共计  0  个样本";

  StatusBar1->Panels->Items[0]->Text="分类迭代次数：0 次";

  int SampleSN=0;//模式样本序号
  //ListView1->Visible=false;
  int iFileHandle,iFileLength,iBytesRead;
  char *PatternBuffer;
  if(OpenDialog1->Execute())
  {
    try
    {
      iFileHandle=FileOpen(OpenDialog1->FileName,fmOpenRead);
      iFileLength=FileSeek(iFileHandle,0,2);
      FileSeek(iFileHandle,0,0);
      PatternBuffer=new char[iFileLength+1];
      iBytesRead=FileRead(iFileHandle,PatternBuffer,iFileLength);
      FileClose(iFileHandle);
      bool ReadFlag=false;  //数据加入与否标志
      AnsiString EData="";  //读入有效数据
      TListItem *p;
      int DataDim=0;  //样本的维
      for(int i=0;i<iBytesRead;i++)
      {
         if((PatternBuffer[i]>='0'&&PatternBuffer[i]<='9')||PatternBuffer[i]=='.')
            {ReadFlag=true;EData+=PatternBuffer[i];}
         else
            {
              if(ReadFlag)
                 {
                   if(DataDim==0)
                   {
                     p=ListView1->Items->Add();
                     p->Caption=EData;
                   }
                   else{
                        p->SubItems->Add(EData);
                        Pattern[SampleSN][DataDim-1]=atof(EData.c_str());
                        }
                   EData="";
                   ReadFlag=false;
                   DataDim++;
                   if(DataDim==5)
                   {
                     DataDim=0;
                     SampleSN++;
                   }
                   
                 }
            }
      }
      
  delete[] PatternBuffer;
    }
    catch(...)
    {
      Application->MessageBox("不能进行以下文件操作：打开、寻找、读取和关闭。", "File Error", IDOK);
    }
  }

//用加权系数调整样本各维(特征)数据大小
 for(int j=0;j<150;j++)
 {
  Pattern[j][0]=Pattern[j][0]*1.0;
  Pattern[j][1]=Pattern[j][1]*1.0;
  Pattern[j][2]=Pattern[j][2]*10.0;
  Pattern[j][3]=Pattern[j][3]*10.0;
 }

/* //样本各维(特征)数据标准化
 double sum[4];
 for(int i=0;i<4;i++)
 {
   sum[i]=0.0;
   for(int j=0;j<150;j++)sum[i]+=Pattern[j][i];
   sum[i]=sum[i]/150.0;
   for(int k=0;k<150;k++)Pattern[k][i]=Pattern[k][i]/ sum[i];
 }   */

 ListView1->Visible=true;
 SpeedButton2->Enabled=true;
}
//---------------------------------------------------------------------------

void __fastcall TForm1::SpeedButton2Click(TObject *Sender)
{
  Screen->Cursor=crHourGlass;

  //清空上次分类结果
  ListView2->Items->Clear();
  ListView3->Items->Clear();
  ListView4->Items->Clear();

  InitClusters(); //初始化

  //K-Means方法计算
  int Converged; //收敛标志
  int IterNum;  //分类迭代次数
  IterNum=1;
  Converged=0;
  while(Converged==0)
  {
    IterNum++;
    StatusBar1->Panels->Items[0]->Text="分类迭代次数："+IntToStr(IterNum)+"次";
    DistributeSamples();
    Converged=CalcNewClustCenters();
   }

  // 显示分类后各类样本数量
  Label1->Caption="共计"+IntToStr(ClusterNumMembers[0])+"个样本";
  Label2->Caption="共计"+IntToStr(ClusterNumMembers[1])+"个样本";
  Label3->Caption="共计"+IntToStr(ClusterNumMembers[2])+"个样本";

  //显示第一类中心
  TListItem *p;
  p=ListView2->Items->Add();
  p->Caption="中心";
  for(int k=0;k<4;k++)
    p->SubItems->Add(String(int(ClusterCenter[0][k]*10)/10.0));
  //显示第一类各样本数据
  for(int i=0;i<ClusterNumMembers[0];i++)
  {
   p=ListView2->Items->Add();
   p->Caption=String(ClusterMember[0][i]+1);
   for(int k=0;k<4;k++)
      p->SubItems->Add(String(Pattern[(ClusterMember[0][i])][k]));
   }

  // 显示第二类中心
  p=ListView3->Items->Add();
  p->Caption="中心";
  for(int k=0;k<4;k++)
    p->SubItems->Add(String(int(ClusterCenter[1][k]*10)/10.0));
  //  显示第二类各样本数据
  for(int i=0;i<ClusterNumMembers[1];i++)
  {
   p=ListView3->Items->Add();
   p->Caption=String(ClusterMember[1][i]+1);
   for(int k=0;k<4;k++)
      p->SubItems->Add(String(Pattern[(ClusterMember[1][i])][k]));
   }
  //显示第三类中心
  p=ListView4->Items->Add();
  p->Caption="中心";
  for(int k=0;k<4;k++)
    p->SubItems->Add(String(int(ClusterCenter[2][k]*10)/10.0));
  //显示第三类各样本数据
  for(int i=0;i<ClusterNumMembers[2];i++)
  {
   p=ListView4->Items->Add();
   p->Caption=String(ClusterMember[2][i]+1);
   for(int k=0;k<4;k++)
      p->SubItems->Add(String(Pattern[(ClusterMember[2][i])][k]));
   }
   Screen->Cursor=crDefault;        
}
//---------------------------------------------------------------------------

void __fastcall TForm1::SpeedButton3Click(TObject *Sender)
{
  Close();         
}
//---------------------------------------------------------------------------