Kmea的ns C语言代码.doc

#include stdio.h #include math.h #include time.h #include stdlib.h #define TRUE 1 #define FALSE 0 int N;//数据个数 int K;//集合个数 int * CenterIndex;//初始化质心数组的索引 double * Center;//质心集合 double * CenterCopy;//质心集合副本 double * AllData;//数据集合 double ** Cluster;//簇的集合 int * Top;//集合中元素的个数，也会用作栈处理 //随机生成k个数x(0=x=n-1)作为起始的质心集合 void CreateRandomArray(int n, int k,int * center) { int i=0; int j=0; srand((unsigned)time(NULL)); for( i=0;ik;++i)//随机生成k个数 { int a=rand()%n; //判重 for(j=0;ji;j++) { if(center[j]==a)//重复 { break; } } if(j=i)//如果不重复，加入 { center[i]=a; } else { i--; //如果重复，本次重新随机生成 } } } //返回距离最小的质心的序号 int GetIndex(double value,double * center) { int i=0; int index=i;//最小的质心序号 double min=fabs(value-center[i]);//距质心最小距离 for(i=0;iK;i++) { if(fabs(value-center[i])min)//如果比当前距离还小，更新最小的质心序号和距离值 { index=i; min=fabs(value-center[i]); } } return index; } //拷贝质心数组到副本 void CopyCenter() { int i=0; for(i=0;iK;i++) { CenterCopy[i]=Center[i]; } } //初始化质心，随机生成法 void InitCenter() { int i=0; CreateRandomArray(N,K,CenterIndex);//产生随机的K个N的不同的序列 for(i=0;iK;i++) { Center[i]=AllData[CenterIndex[i]];//将对应数据赋值给质心数组 } CopyCenter();//拷贝到质心副本 } //加入一个数据到一个Cluster[index]集合 void AddToCluster(int index,double value) { Cluster[index][Top[index]++]=value;//这里同进栈操作 } //重新计算簇集合 void UpdateCluster() { int i=0; int tindex; //将所有的集合清空，即将TOP置0 for(i=0;iK;i++) { Top[i]=0; } for(i=0;iN;i++) { tindex=GetIndex(AllData[i],Center);//得到与当前数据最小的质心索引 AddToCluster(tindex,AllData[i]); //加入到相应的集合中 } } //重新计算质心集合，对每一簇集合中的元素加总求平均即可 void UpdateCenter() { int i=0; int j=0;