用C++实现的CURE算法的源码

/****************************************************************************
*                                                                           *
*  CURE                                                                     *
*                                                                           *
*****************************************************************************/

#include 
#include 
#include 
#include 
#include 
#include 

#define         SUCCESS         1
#define         FAILURE         0
#define         TRUE            1
#define         FALSE           0
#define         NUMPATTERNS     10   // 数据点的个数
#define         SIZEVECTOR      2     // 数据集的维数
#define         NUMCLUSTERS     2     // 类的个数
#define         NUMPRE          3    // 代表点个数
#define         SHRINK          0.5   // 收缩率

// ***** 定义结构和类 *****
struct ClustNode {
	int          Member[NUMPATTERNS];  //类中的数据项
	int          NumMembers;
	double       Means[SIZEVECTOR];  //均值点
	double       Pre[NUMPRE][SIZEVECTOR];//代表点
	int          PreMember[NUMPRE];//作为代表点的数据项
	int          closet;  //最近的簇
	double       MinDist;//与最近的簇之间的距离
};

class System {
private:
	double      Pattern[NUMPATTERNS][SIZEVECTOR];
	ClustNode *p[NUMPATTERNS];
    void        BuildClustList（）; //建立初始簇链表
	ClustNode * Merge（ClustNode *ClustNode *）; //合并两个簇
	int         MinClust（）;    //找到最近的两个簇
	double      dist（ClustNode*ClustNode *）;//两个簇之间的最短距离
	double      EucNorm（double *double *）;//代表点之间的距离
	void        ShrinkPre（ClustNode *）;  //收缩代表点

public:
	System（）;
	~System（）;
	int LoadPatterns（char *fname）;      // 处理数据集
	void RunCure（）;                     // 聚类的过程
	void ShowClusters（）;                // 显示聚类结果
};
System::System（）
{
}
System::~System（）
{
	int i;
	for（i=0;i		if（p[i]==NULL）
			continue;
		else
			delete p[i];
	}
}
int System::LoadPatterns（char *fname）{
	int    ijlen;
	char buff[300];
	char *ptr*ptr1;
	ifstream infile（fnameios::in）;    //将fname打开作为输入
	if（!infile）{
		cout<<“不能打开输入文件:\n“;
		return FAILURE;
	}
	for（i=0; i		infile.getline（buff300）;     //将300个字符传入buff中，buff是指针
		len=0;
		ptr=buff;
		for （j=0; j			ptr1=strchr（ptr‘‘）;   //找出sptr指向的字符串中第一次出现““的位置
			*ptr1=‘\0‘;
			Pattern[i][j]=atof（ptr）;    //将ptr指向的字符串转化为浮点数
			len=strlen（ptr）;
			ptr=ptr+len+1;
		} /* endfor */
	} /* endfor */
	infile.close（）;
	return SUCCESS;
}
void  System::BuildClustList（）{
	int ij;
	double MinDistd;
	for（i=0;i		p[i]=new ClustNode;
		p[i]->NumMembers=1;
		p[i]->Member[0]=i;
		p[i]->PreMember[0]=i;
		for（j=0;j			p[i]->Means[j]=Pattern[i][j];//均值
		}
		for（j=0;j			p[i]->Pre[0][j]=Pattern[i][j]+（SHRINK）*（p[i]->Means[j]-Pattern[i][j]）;//计算代表点
		}
	}
	for（i=0;i		p[i]->MinDist=9.9e+99;
		for（j=0;j