欢迎您访问程序员文章站本站旨在为大家提供分享程序员计算机编程知识!
您现在的位置是: 首页  >  IT编程

K均值聚类算法的Java版实现代码示例

程序员文章站 2023-12-02 19:41:46
1.简介 k均值聚类算法是先随机选取k个对象作为初始的聚类中心。然后计算每个对象与各个种子聚类中心之间的距离,把每个对象分配给距离它最近的聚类中心。聚类中心以及分配给它们...

1.简介

k均值聚类算法是先随机选取k个对象作为初始的聚类中心。然后计算每个对象与各个种子聚类中心之间的距离,把每个对象分配给距离它最近的聚类中心。聚类中心以及分配给它们的对象就代表一个聚类。一旦全部对象都被分配了,每个聚类的聚类中心会根据聚类中现有的对象被重新计算。这个过程将不断重复直到满足某个终止条件。终止条件可以是没有(或最小数目)对象被重新分配给不同的聚类,没有(或最小数目)聚类中心再发生变化,误差平方和局部最小。

2.什么是聚类

聚类是一个将数据集中在某些方面相似的数据成员进行分类组织的过程,聚类就是一种发现这种内在结构的技术,聚类技术经常被称为无监督学习。

3.什么是k均值聚类

k均值聚类是最著名的划分聚类算法,由于简洁和效率使得他成为所有聚类算法中最广泛使用的。给定一个数据点集合和需要的聚类数目k,k由用户指定,k均值算法根据某个距离函数反复把数据分入k个聚类中。

4.实现

java代码如下:

package org.algorithm;
import java.util.arraylist;
import java.util.random;
/** 
 * k均值聚类算法 
 */
public class kmeans {
	private int k;
	// 分成多少簇 
	private int m;
	// 迭代次数 
	private int datasetlength;
	// 数据集元素个数,即数据集的长度 
	private arraylist<float[]> dataset;
	// 数据集链表 
	private arraylist<float[]> center;
	// 中心链表 
	private arraylist<arraylist<float[]>> cluster;
	// 簇 
	private arraylist<float> jc;
	// 误差平方和,k越接近datasetlength,误差越小 
	private random random;
	/** 
   * 设置需分组的原始数据集 
   * 
   * @param dataset 
   */
	public void setdataset(arraylist<float[]> dataset) {
		this.dataset = dataset;
	}
	/** 
   * 获取结果分组 
   * 
   * @return 结果集 
   */
	public arraylist<arraylist<float[]>> getcluster() {
		return cluster;
	}
	/** 
   * 构造函数,传入需要分成的簇数量 
   * 
   * @param k 
   *      簇数量,若k<=0时,设置为1,若k大于数据源的长度时,置为数据源的长度 
   */
	public kmeans(int k) {
		if (k <= 0) {
			k = 1;
		}
		this.k = k;
	}
	/** 
   * 初始化 
   */
	private void init() {
		m = 0;
		random = new random();
		if (dataset == null || dataset.size() == 0) {
			initdataset();
		}
		datasetlength = dataset.size();
		if (k > datasetlength) {
			k = datasetlength;
		}
		center = initcenters();
		cluster = initcluster();
		jc = new arraylist<float>();
	}
	/** 
   * 如果调用者未初始化数据集,则采用内部测试数据集 
   */
	private void initdataset() {
		dataset = new arraylist<float[]>();
		// 其中{6,3}是一样的,所以长度为15的数据集分成14簇和15簇的误差都为0 
		float[][] datasetarray = new float[][] { { 8, 2 }, { 3, 4 }, { 2, 5 }, 
		        { 4, 2 }, { 7, 3 }, { 6, 2 }, { 4, 7 }, { 6, 3 }, { 5, 3 }, 
		        { 6, 3 }, { 6, 9 }, { 1, 6 }, { 3, 9 }, { 4, 1 }, { 8, 6 } };
		for (int i = 0; i < datasetarray.length; i++) {
			dataset.add(datasetarray[i]);
		}
	}
	/** 
   * 初始化中心数据链表,分成多少簇就有多少个中心点 
   * 
   * @return 中心点集 
   */
	private arraylist<float[]> initcenters() {
		arraylist<float[]> center = new arraylist<float[]>();
		int[] randoms = new int[k];
		boolean flag;
		int temp = random.nextint(datasetlength);
		randoms[0] = temp;
		for (int i = 1; i < k; i++) {
			flag = true;
			while (flag) {
				temp = random.nextint(datasetlength);
				int j = 0;
				// 不清楚for循环导致j无法加1 
				// for(j=0;j<i;++j) 
				// { 
				// if(temp==randoms[j]); 
				// { 
				// break; 
				// } 
				// } 
				while (j < i) {
					if (temp == randoms[j]) {
						break;
					}
					j++;
				}
				if (j == i) {
					flag = false;
				}
			}
			randoms[i] = temp;
		}
		// 测试随机数生成情况 
		// for(int i=0;i<k;i++) 
		// { 
		// system.out.println("test1:randoms["+i+"]="+randoms[i]); 
		// } 
		// system.out.println(); 
		for (int i = 0; i < k; i++) {
			center.add(dataset.get(randoms[i]));
			// 生成初始化中心链表
		}
		return center;
	}
	/** 
   * 初始化簇集合 
   * 
   * @return 一个分为k簇的空数据的簇集合 
   */
	private arraylist<arraylist<float[]>> initcluster() {
		arraylist<arraylist<float[]>> cluster = new arraylist<arraylist<float[]>>();
		for (int i = 0; i < k; i++) {
			cluster.add(new arraylist<float[]>());
		}
		return cluster;
	}
	/** 
   * 计算两个点之间的距离 
   * 
   * @param element 
   *      点1 
   * @param center 
   *      点2 
   * @return 距离 
   */
	private float distance(float[] element, float[] center) {
		float distance = 0.0f;
		float x = element[0] - center[0];
		float y = element[1] - center[1];
		float z = x * x + y * y;
		distance = (float) math.sqrt(z);
		return distance;
	}
	/** 
   * 获取距离集合中最小距离的位置 
   * 
   * @param distance 
   *      距离数组 
   * @return 最小距离在距离数组中的位置 
   */
	private int mindistance(float[] distance) {
		float mindistance = distance[0];
		int minlocation = 0;
		for (int i = 1; i < distance.length; i++) {
			if (distance[i] < mindistance) {
				mindistance = distance[i];
				minlocation = i;
			} else if (distance[i] == mindistance) // 如果相等,随机返回一个位置 
			{
				if (random.nextint(10) < 5) {
					minlocation = i;
				}
			}
		}
		return minlocation;
	}
	/** 
   * 核心,将当前元素放到最小距离中心相关的簇中 
   */
	private void clusterset() {
		float[] distance = new float[k];
		for (int i = 0; i < datasetlength; i++) {
			for (int j = 0; j < k; j++) {
				distance[j] = distance(dataset.get(i), center.get(j));
				// system.out.println("test2:"+"dataset["+i+"],center["+j+"],distance="+distance[j]);
			}
			int minlocation = mindistance(distance);
			// system.out.println("test3:"+"dataset["+i+"],minlocation="+minlocation); 
			// system.out.println(); 
			cluster.get(minlocation).add(dataset.get(i));
			// 核心,将当前元素放到最小距离中心相关的簇中
		}
	}
	/** 
   * 求两点误差平方的方法 
   * 
   * @param element 
   *      点1 
   * @param center 
   *      点2 
   * @return 误差平方 
   */
	private float errorsquare(float[] element, float[] center) {
		float x = element[0] - center[0];
		float y = element[1] - center[1];
		float errsquare = x * x + y * y;
		return errsquare;
	}
	/** 
   * 计算误差平方和准则函数方法 
   */
	private void countrule() {
		float jcf = 0;
		for (int i = 0; i < cluster.size(); i++) {
			for (int j = 0; j < cluster.get(i).size(); j++) {
				jcf += errorsquare(cluster.get(i).get(j), center.get(i));
			}
		}
		jc.add(jcf);
	}
	/** 
   * 设置新的簇中心方法 
   */
	private void setnewcenter() {
		for (int i = 0; i < k; i++) {
			int n = cluster.get(i).size();
			if (n != 0) {
				float[] newcenter = { 0, 0 };
				for (int j = 0; j < n; j++) {
					newcenter[0] += cluster.get(i).get(j)[0];
					newcenter[1] += cluster.get(i).get(j)[1];
				}
				// 设置一个平均值 
				newcenter[0] = newcenter[0] / n;
				newcenter[1] = newcenter[1] / n;
				center.set(i, newcenter);
			}
		}
	}
	/** 
   * 打印数据,测试用 
   * 
   * @param dataarray 
   *      数据集 
   * @param dataarrayname 
   *      数据集名称 
   */
	public void printdataarray(arraylist<float[]> dataarray, 
	      string dataarrayname) {
		for (int i = 0; i < dataarray.size(); i++) {
			system.out.println("print:" + dataarrayname + "[" + i + "]={" 
			          + dataarray.get(i)[0] + "," + dataarray.get(i)[1] + "}");
		}
		system.out.println("===================================");
	}
	/** 
   * kmeans算法核心过程方法 
   */
	private void kmeans() {
		init();
		// printdataarray(dataset,"initdataset"); 
		// printdataarray(center,"initcenter"); 
		// 循环分组,直到误差不变为止 
		while (true) {
			clusterset();
			// for(int i=0;i<cluster.size();i++) 
			// { 
			// printdataarray(cluster.get(i),"cluster["+i+"]"); 
			// } 
			countrule();
			// system.out.println("count:"+"jc["+m+"]="+jc.get(m)); 
			// system.out.println(); 
			// 误差不变了,分组完成 
			if (m != 0) {
				if (jc.get(m) - jc.get(m - 1) == 0) {
					break;
				}
			}
			setnewcenter();
			// printdataarray(center,"newcenter"); 
			m++;
			cluster.clear();
			cluster = initcluster();
		}
		// system.out.println("note:the times of repeat:m="+m);//输出迭代次数
	}
	/** 
   * 执行算法 
   */
	public void execute() {
		long starttime = system.currenttimemillis();
		system.out.println("kmeans begins");
		kmeans();
		long endtime = system.currenttimemillis();
		system.out.println("kmeans running time=" + (endtime - starttime) 
		        + "ms");
		system.out.println("kmeans ends");
		system.out.println();
	}
}

5.说明:

具体代码是从网上找的,根据自己的理解加了注释和进行部分修改,若注释有误还望指正

6.测试

package org.test;
import java.util.arraylist;
import org.algorithm.kmeans;
public class kmeanstest {
	public static void main(string[] args) 
	  {
		//初始化一个kmean对象,将k置为10 
		kmeans k=new kmeans(10);
		arraylist<float[]> dataset=new arraylist<float[]>();
		dataset.add(new float[]{1,2});
		dataset.add(new float[]{3,3});
		dataset.add(new float[]{3,4});
		dataset.add(new float[]{5,6});
		dataset.add(new float[]{8,9});
		dataset.add(new float[]{4,5});
		dataset.add(new float[]{6,4});
		dataset.add(new float[]{3,9});
		dataset.add(new float[]{5,9});
		dataset.add(new float[]{4,2});
		dataset.add(new float[]{1,9});
		dataset.add(new float[]{7,8});
		//设置原始数据集 
		k.setdataset(dataset);
		//执行算法 
		k.execute();
		//得到聚类结果 
		arraylist<arraylist<float[]>> cluster=k.getcluster();
		//查看结果 
		for (int i=0;i<cluster.size();i++) 
		    {
			k.printdataarray(cluster.get(i), "cluster["+i+"]");
		}
	}
}

总结:测试代码已经通过。并对聚类的结果进行了查看,结果基本上符合要求。至于有没有更精确的算法有待发现。具体的实践还有待挖掘

总结

以上就是本文关于k均值聚类算法的java版实现代码示例的全部内容,希望对大家有所帮助。感兴趣的朋友可以继续参阅本站其他相关专题。如有不足之处,欢迎留言指出。感谢朋友们对本站的支持!