C++实现K-means聚类算法 _ 质心迭代与欧式距离计算【源码】

风萱君_2184

风萱君_2184

2026-04-08

710人浏览

原创

kmeans()循环不收敛主因是质心更新时误用未同步的新标签,导致分组依据失效;应严格分两阶段:①全部点重新分配标签,②统一用新标签重算质心,并辅以标签稳定性与位移双条件判断收敛。

c++实现k-means聚类算法 _ 质心迭代与欧式距离计算【源码】

为什么 kmeans() 循环不收敛?检查质心更新是否用了旧标签

常见现象是聚类结果在几轮后停滞,或质心坐标疯狂跳变。根本原因常在于:计算新质心时,仍用上一轮未更新的 labels 数组去分组点——但该数组可能已在本轮开头被重写,而你误把它当作“本轮真实归属”来用。

正确做法是:先完整遍历所有点,用当前质心算出每个点的新标签,存入临时数组;确认全部点都完成归属判定后,再统一用这个新标签数组重新计算所有质心。否则会出现“边分组边改标签,分组依据已失效”的竞态逻辑。

  • 每次迭代必须严格分两阶段:① 所有点重新分配 → ② 所有质心一次性重算
  • 避免在循环中直接修改 labels[i] 后立刻用它参与质心累加(这会让同一轮里部分点按新标签、部分按旧标签归组)
  • 若用 std::vector<:vector>> 存数据,注意内层 vector 大小要一致(比如每点都是 2 维或 3 维),否则 sqrt(sum(pow(x_i - c_i, 2))) 会越界

欧式距离计算别手写循环:用 std::inner_product 更安全

手动写 for 循环算平方和容易漏掉维度检查,尤其当点向量长度和质心向量长度不一致时,运行时崩溃或静默错误都很常见。C++ 标准库的 std::inner_product 能天然对齐两个范围,还能复用已有迭代器逻辑。

示例:计算点 p 到质心 c 的平方欧氏距离(避免开方节省性能):

double sq_dist = std::inner_product(
    p.begin(), p.end(), c.begin(), 0.0,
    std::plus{},
    [](double a, double b) { return (a - b) * (a - b); }
);
  • 不用自己写 for (size_t i = 0; i ,省去下标越界风险
  • 如果 p.size() != c.size(),std::inner_product 在短序列结束时自动停止,不会崩溃(但结果错误——所以仍需在调用前断言尺寸一致)
  • 返回的是平方距离,K-means 中比较大小完全够用,比每次调用 sqrt() 快 2–3 倍

质心初始化不能全设为零:用 std::uniform_real_distribution 在数据范围内采样

所有质心初始值设成 {0,0} 或固定值,会导致大量点距离相等,std::min_element 返回首个索引,造成严重偏差——第一簇永远吃掉最多点,后续迭代难以修正。

C++14
C++14

C++14 对 C++11 的修正与增强版本,适合旧系统维护和较老工具链兼容。

下载

正确做法是从原始数据点中随机选 k 个作为初始质心(k-means++ 改进版更优,但随机采样已解决 90% 的发散问题):

  • 先遍历一遍数据,记录每维的最小/最大值:min_val[i], max_val[i]
  • 用 std::random_device + std::mt19937 初始化引擎,对每维用 std::uniform_real_distribution(min_val[i], max_val[i]) 生成坐标
  • 或者更简单:用 std::shuffle 打乱数据索引,取前 k 个点的值直接赋给 centroids
  • 切忌用 rand() % 100 这类整数随机数填充浮点质心——范围错、精度丢、分布不均

判断收敛别只看质心位移:必须同时检查标签稳定性

仅监控质心坐标的 L2 变化量(如 max_norm_diff )不够。某些数据集会出现质心微调但点的大规模重分配,此时算法实际未稳定,却提前退出。

健壮的收敛条件应是“标签不再变化”或“变化点比例低于阈值”,因为 K-means 的目标是最小化簇内平方和,而标签决定该目标函数值:

  • 每次迭代后用 std::equal(labels_old.begin(), labels_old.end(), labels_new.begin()) 判断是否全等
  • 若不等,统计不同位置数量:std::inner_product(labels_old.begin(), labels_old.end(), labels_new.begin(), 0, std::plus{}, std::not_equal_to{})
  • 建议设置双条件:位移 且 标签变化率

维度高、样本少时,质心可能在几个位置间震荡,标签反复切换——这时需要加最大迭代轮数兜底,否则死循环。

C++免费学习笔记(深入):立即使用
在学习笔记中,你将探索 C++ 的入门与实战技巧!

相关文章

PHP速学视频免费教程(入门到精通)
PHP速学视频免费教程(入门到精通)

PHP怎么学习?PHP怎么入门?PHP在哪学?PHP怎么学才快?不用担心,这里为大家提供了PHP速学教程(入门到精通),有需要的小伙伴保存下载就能学习啦!

下载

相关标签:

c++

本站声明:本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系admin@php.cn

相关专题

更多
页面置换算法
页面置换算法

页面置换算法是操作系统中用来决定在内存中哪些页面应该被换出以便为新的页面提供空间的算法。本专题为大家提供页面置换算法的相关文章,大家可以免费体验。

2023.08.14

6776

4

Postman身份认证怎么配置
Postman身份认证怎么配置

本专题从实际API接口调试场景出发,讲解Postman身份认证的完整配置流程,包括常见认证协议介绍、授权参数设置、访问令牌填写、OAuth 2.0授权及认证失败排查,帮助用户正确完成接口授权与访问验证。

2026.11.06

0

12

Postman请求头Headers怎么设置
Postman请求头Headers怎么设置

围绕Postman HTTP请求头配置展开,讲解Headers选项的具体使用方式,包括Key和Value参数填写、请求头信息添加、常见Header字段设置及请求头修改方法,帮助用户解决接口调试过程中遇到的请求头配置问题。

2026.11.06

0

13

Postman接口请求发送教程
Postman接口请求发送教程

详细介绍Postman接口请求的创建与发送方法,涵盖GET和POST请求、URL配置、请求参数、Headers请求头、Body请求体设置以及响应结果查看,帮助用户掌握Postman接口请求的完整操作流程。

2026.11.06

0

10

DBeaver怎么导入数据
DBeaver怎么导入数据

DBeaver数据导入专题,整理Excel、CSV、SQL、TXT等常见文件的数据导入方法,涵盖目标数据库选择、数据表创建、字段映射、数据类型设置及导入参数配置等操作,帮助用户将本地文件中的数据写入数据库,并了解批量导入、重复数据处理和导入失败等常见问题的解决方法。

2026.11.05

20

12

DBeaver安全防护教程大全
DBeaver安全防护教程大全

DBeaver安全防护专题,围绕数据库连接安全、账号密码保护和访问权限管理展开,整理SSL/TLS加密连接、SSH隧道、数据库凭据保存及证书验证等配置方法,同时介绍敏感信息保护、远程连接风险和高权限账号使用等注意事项,帮助用户提高数据库日常管理的安全性。

2026.11.05

20

12

DBeaver连接创建与配置教程
DBeaver连接创建与配置教程

本专题聚焦DBeaver数据库连接配置与管理,覆盖从新建连接、选择数据库类型到填写连接参数和验证连接的完整流程。针对MySQL、PostgreSQL、SQL Server等常见数据库,整理相应的连接设置方法,并提供驱动缺失、认证失败和连接异常等问题的排查方向。

2026.11.04

20

16

XORM框架结构体映射与表结构同步
XORM框架结构体映射与表结构同步

本专题围绕XORM框架的结构体映射展开,说明表名与字段的命名策略、xorm标签常用属性、主键自增与时间字段自动填充,并演示Sync2自动建表、增量字段升级、索引声明与反向生成结构体的完整流程。

2026.11.04

20

15

XORM框架安装配置与数据库连接
XORM框架安装配置与数据库连接

本专题聚焦XORM框架的配置与连接环节,拆解MySQL、PostgreSQL、SQLite、SQL Server的DSN写法与驱动选择,讲清连接池参数、多数据源拆分、读写分离与超时重连的处理方式,并补充容器与兼容数据库的适配经验。

2026.11.04

20

15

热门下载

更多
网站特效
/
网站源码
/
网站素材
/
前端模板

精品课程

更多
相关推荐
/
热门推荐
/
最新课程
C++ Core Guidelines
C++ Core Guidelines

共0课时 | 0人学习

C++ Reference
C++ Reference

共0课时 | 0人学习

C++ 官方标准说明
C++ 官方标准说明

共0课时 | 0人学习