模式识别k均值算法解析.docVIP

模式识别导论大作业 （2015-2016第一学期） 姓名 学 号 161310224 任课教师 南京航空航天大学 9 月 19 日 K均值聚类 功能描述： 利用K-均值算法将150个模式样本分成3类别。分别计算最后算法所用的迭代次数，最终 带注释的源代码 #include stdio.h #include math.h typedef struct mydata{ float data[4]; //保存原始数据 int index; }MyData; MyData myData[150]; float classdata[3][150][4]; int saveindex[3][150]; //保存下表 float center[3][4]; //保存三个聚类中心 float newcenter[3][4];//保存三个新的聚类中心 int findMinCenter(float *temp){//找到当前数据的类别 float dis[3]={0,0,0}; float min = 100000; int returnk=0; for(int k=0;k3;k++){ float mytemp = 0; for(int m=0;m4;m++) mytemp += (temp[m]-center[k][m])*(temp[m]-center[k][m]); dis[k] = sqrt(mytemp); //printf(%f\n,dis[k]); if(dis[k]min){ min = dis[k]; returnk = k; }} return returnk; } void counterNewCenter(float tempdata[][4],int len,int num){ float temp[4] = {0,0,0,0}; //计算新的聚类中心方法 for(int j=0;jlen;j++){ for(int i=0;i4;i++){ temp[i] += tempdata[j][i]; } } for(int k=0;k4;k++){ newcenter[num][k] = temp[k]/len; printf(k = %f\n,newcenter[num][k]); }} int main(){ int k=0; int kt[3] = {0,0,0}; int i=0,u=0; FILE *fp; if ((fp = fopen(D:\\Iris.txt, r)) == NULL){ //打开文件 printf(打开文件失败\n); return 0; } float numtemp; while (!feof(fp)){ //读取数据 fscanf(fp, %f, numtemp); myData[u].data[i++] = numtemp; if(i==4){ u++; i=0;} } fclose(fp);//关闭文件 for(k=0;k4;k++){ center[0][k] = myData[0].data[k]; center[1][k] = myData[1].data[k]; center[2][k] = myData[2].data[k]; }//设置初始聚类中心 int s=0;int num = 0;i=0;//聚类所用的迭代次数 while(i150){ i++; k = findMinCenter(myData[i].data); //printf(k = %d\n,k); saveindex[k][kt[k]] = i; for(s=0;s4;s++){ classdata[k][kt[k]][s] = myData[i].data[s];}//每次将数据分到不同的类中 kt[k]++;if (i == 150){ printf(******************); for (int m = 0; m3; m++)//计算新的聚类中心 counterNewCenter(classdata[m],kt[m],m); bool flag = true;//标示前后聚类中心是否相等的布尔类型 for(int n=0;n3;n++){ for(int t=0;t4;t++){ if(newcenter[n][t] != center[n][t]){ flag = false; break