博客
关于我
聚类分析笔记-K均值matlab算法(一)
阅读量:221 次
发布时间:2019-02-28

本文共 1440 字,大约阅读时间需要 4 分钟。

一:聚类分析的动态聚类算法

聚类分析是一种无监督学习方法,用于将数据点按一定规则分组。常见的动态聚类算法之一是K均值算法(K-means),其思想是通过迭代优化聚类中心,使得每个样本点离其聚类中心的距离平方和最小。

K均值算法的思想

  • 初始聚类中心:随机选取K个样本作为初始聚类中心。
  • 迭代优化
    • 对每个样本点计算其到各个聚类中心的距离。
    • 根据距离将样本点分配到最近的聚类中心。
    • 计算每个聚类中心的新坐标(均值)。
    • 比较当前聚类中心与新聚类中心的差异,若收敛则终止,否则继续迭代。
  • 收敛判断:若聚类中心向量变化小,视为收敛。
  • K均值算法的常用方法

    • K均值算法:适用于已知类别数的情况,效果较好。
    • ISODATA算法(迭代自组织数据分析算法):动态调整类别数,适合未知类别数的情况,常用Matlab实现。

    K均值算法的优点与注意事项

    • 优点
      • 简单易实现。
      • 适用于已知类别数的聚类问题。
      • 计算结果具有几何意义。
    • 注意事项
      • 初始聚类中心的选择会影响结果,需合理选择。
      • 对于高维数据,计算量较大,需注意性能优化。
      • 数据的顺序可能影响聚类结果。

    二:Matlab程序实现示例

    以下是基于K均值算法的Matlab程序实现,用于两类聚类问题:

    function mean = k_means_new% 生成模拟数据num = 100; % 样本总数x1 = rand(1, num) * 5; % 类型1y1 = rand(1, num) * 5;x2 = rand(1, num) * 5 + 3; % 类型2y2 = rand(1, num) * 5 + 3;cities = [x1, y1; x2, y2];% 随机选择初始聚类中心m = round(rand(1, num) * num); % 随机选取一个样本作为聚类中心while true    m2 = round(rand(1, num) * num);    if m ~= m2        break    else        m2    endendu1 = cities(:, m);u2 = cities(:, m2);u_old = [u1, u2];u_new = [u2, u1];while true    diff = u_old ~= u_new    if ~diff        break    end    u_old = u_new;    [c, ~] = min(distances(cities, u_old), 2);    index1 = find(c == 1);    index2 = find(c == 2);    u1 = mean(cities(:, index1), 2);    u2 = mean(cities(:, index2), 2);    u_new = [u1, u2];endmean = u_new;end

    三:实验结果与分析

    通过实验验证,K均值算法在已知类别数的情况下表现良好。随着初始样本数量的增加,聚类精度有所提升。以下是不同初始点数量下的聚类效果对比:

    • 10个初始点:聚类效果较为分散。
    • 20个初始点:聚类效果有所改善,分类准确率提高。
    • 50个初始点:聚类效果更为稳定,分类准确率显著提升。
    • 100个初始点:聚类效果最为理想,分类准确率接近100%。

    实验结果表明,K均值算法在已知类别数的情况下表现出色,适合用于分类问题。

    转载地址:http://sowi.baihongyu.com/

    你可能感兴趣的文章
    Python 初学者需要知道的四条建议
    查看>>
    Python 判断字符串是否包含子字符串
    查看>>
    python 判断当前时间是否为零点
    查看>>
    python 利用pandas读取本地中CSV文件的指定列 列名重命名 并保存回本地
    查看>>
    python 利用pyspark读取HDFS中CSV文件的指定列 列名重命名 并保存回HDFS
    查看>>
    python 利用pyttsx3文字转语音
    查看>>
    python 利用已有Ner模型进行数据清洗合并
    查看>>
    python 到大数据开发工程师_如何成为一个大数据开发工程师?
    查看>>
    python 加密解密(base64, AES)
    查看>>
    Python 包管理器和 Node.js
    查看>>
    Python 单词字母顺序不变且所有倒排
    查看>>
    python 反射机制
    查看>>
    python 启动提示IDLE's subprocess didn't make conne...
    查看>>
    python 命令接口_实现“[命令][操作][参数]”样式的命令行接口?
    查看>>
    Python 和 OpenCV.如何检测图像中的所有(填充)圆形/圆形对象?
    查看>>
    Python 和 RabbitMQ - 聆听来自多个渠道的消费事件的最佳方式?
    查看>>
    python编辑器打不开_关于命令ride.py打不开RF,而是打开pycharm编辑器问题解决思路...
    查看>>
    python编译exe同时支持32位_第三周:同时管理64位和32位版本的Python,并用Pyinstaller打包成exe...
    查看>>
    Python 和Java 哪个更适合做自动化测试?
    查看>>
    Python编程:掌握高级语言程序设计。从零基础到精通,收藏这篇就够了!
    查看>>