模式识别“k-means”C语言程序，就是关于坐标点的划分问题，C语言高手请进!

Question

如果懂模式识别k-means算法的就请直接给我c程序吧，如果不懂看一下介绍：在程序前这样输入15个点int x[15][3] = {{ 3, 1, 0}, { 3, 2, 0}, { 2, 2, 0}, { 3, 3, 0}, { 2, 3, 0}, { 8, 8, 0}，{ 8, 9, 0}, { 9, 8, 0}, { 9, 7, 0}，{ 9, 9, 0}, { 16, 5, 0}, { 16, 4, 0}，{ 15, 5, 0}, { 15, 6, 0}, { 16, 6, 0}} （z坐标为0，相当于二维坐标），选其中任意3个点z1、z2、z3作为基本点，比较其余12点与这3点距离，与哪个点近就和哪个点划分在一起，例如某一点到z2的距离小于到z1、z3距离，它就和z2算是一个类cluster2，这样找出12个点分别属于那个cluster，划分好所有点之后一共分3个cluster，求出各个cluster的所拥有点的中心z1’、z2’、z3’，然后再以这三个新中心点求出原15个点到这3个点距离并划分成新的3个cluster，然后再求新的cluster各点的中心，再求出15个点分别属于哪个新的cluster，然后一次循环下去，直到得到的中心值不变，printf出三个不变的中心值坐标，15点分别属于哪个中心所在的cluster。
由于本人初学C语言不久，请用基本语句，并注明每一步的作用，详细！谢谢！有追加分数！
下面的这个看不懂，呵呵，再给总结下程序步骤：
1.在已知15个点中任意找出3点
2.比较12个点与3点距离并划分成3类
3.求3个类的中心值
4.比较所有15点与3个中心值距离并重新划分成3类
5.再求新中心、比较距离并划分...
6.当新中心值不再变化，即第n与n-1次中心值相同则输出结果。
简单点的改改！

Answer 1

写了一个，参考一下，有什么问题请说。
#include<stdio.h>
#include<math.h>

int x[15][3]={{3,1,0},{3,2,0},{2,2,0},{3,3,0},{2,3,0},{8,8,0},{8,9,0},{9,8,0},{9,7,0},{9,9,0},{16,5,0},{16,4,0},{15,5,0},{15,6,0},{16,6,0}};

double oldcentral[3][3];//旧的中心点的坐标
double newcentral[3][3];//计算距离，分类后求平均值得到的新的中心点的坐标

int clas[15];//15个点各属于哪个类，类的编号从0开始
int clsno;

int minDist(double x,double y,double z) //计算三个数的最小值，返回其序号。
{
if(x<=y&&x<=z)
return 0;
if(y<x&&y<=z)
return 1;
if(z<x&&z<y)
return 2;
}

double distA(int i,int j)//计算两个点的距离，i和j分别是数组x和newcentral中的序号
{
double distx,disty,distz,dist;

distx=(double)x[i][0]-(double)newcentral[j][0];
disty=(double)x[i][1]-(double)newcentral[j][1];
distz=(double)x[i][2]-(double)newcentral[j][2];

dist=sqrt(distx*distx+disty*disty+distz*distz);

return dist;
}

void main()
{
int count[3];//记录每一类的个数；
int i,j;
double dist[3];//求坐标点距离3个中心点距离时用到的变量，
double cenL,cenLx,cenLy,cenLz;//最后求新旧中心点距离的时候用到的变量

printf("The 15 points are:\n");//把所以点的坐标打印一遍，非必要语句
for(i=0;i<15;i++)
printf("%d %d %d \n",x[i][0],x[i][1],x[i][2]);

for(i=0;i<3;i++)//新旧中心点赋初值。
for(j=0;j<3;j++)
{
newcentral[i][j]=(double)x[i][j];
oldcentral[i][j]=-9999.0;
}

for(i=0;i<3;i++)//头3个点作为初始的中心点。
{
clas[i]=i;
cls[i]=i;
}

while(1)//无限循环，退出条件是新旧中心点的距离小于0.005.
{
for(i=0;i<3;i++)//记录每一类的个数的数组赋初值
count[i]=0;

for(i=0;i<3;i++)//新中心点的坐标拷贝到旧中心点数组中，因新中心点需重新计算。
for(j=0;j<3;j++)
oldcentral[i][j]=newcentral[i][j];

for(i=0;i<15;i++)//对15个点分别计算到中心点的距离。
{
for(j=0;j<3;j++)
dist[j]=distA(i,j);

clsno=minDist(dist[0],dist[1],dist[2]);//求距离最小值，返回距离最小的对应中心点坐标。

clas[i]=clsno;//将此点归到距离最小的那一类。
count[clsno]++;

}

for(i=0;i<3;i++)//对新中心点坐标赋初值，进行下面的计算。
for(j=0;j<3;j++)
newcentral[i][j]=0.0;

for(i=0;i<15;i++)//对每一类的坐标点计算其坐标之和。
for(j=0;j<3;j++)
newcentral[clas[i]][j]+=x[i][j];

for(i=0;i<3;i++)//坐标之和除以count数组元素，即得中心点坐标
for(j=0;j<3;j++)
newcentral[i][j]=newcentral[i][j]/count[i];

int flag=0;//标志
for(i=0;i<3;i++)//求新旧中心点的距离
{
cenLx=newcentral[i][0]-oldcentral[i][0];
cenLy=newcentral[i][1]-oldcentral[i][1];
cenLz=newcentral[i][2]-oldcentral[i][2];

cenL=sqrt(cenLx*cenLx+cenLy*cenLy+cenLz*cenLz);

if(cenL>0.005)//只要有一个距离过大，表明未收敛，重新开始循环
flag=1;
}

if(flag!=1)//只有当标志未被设置，退出while循环。
break;
}

for(i=0;i<15;i++) //打印15个点各自所属的类。
{
printf("point %d(%d,%d,%d) belongs to class %d\n",i,x[i][0],x[i][1],x[i][2],clas[i]);
}
}

Answer 2

搬的凳子坐这看