K-Means 聚类的目标函数:簇内误差平方和

作者:F_D_Z日期:2026/6/16

1. 什么是 K-Means?

K-Means 是一种无监督迭代式的聚类算法:
给定数据集 {x₁, x₂, …, xₙ} 与预设簇数 K,算法把样本划分为 K 个不相交的簇 C₁, C₂, …, Cₖ,使得同一簇内样本尽可能相似,不同簇间样本尽可能远离

核心思想:
> “让簇内‘抱团’,让簇间‘疏远’。”


2. 目标函数 J:簇内误差平方和(WCSS)

K-Means 用几何距离衡量相似性,目标函数 J 定义为:
J=∑k=1K∑x∈Ck∥x−μk∥2 J = \sum_{k=1}^{K} \sum_{x \in C_k} \|x - \mu_k\|^2J=k=1∑K​x∈Ck​∑​∥x−μk​∥2

  • μₖ:第 k 个簇的质心(centroid)
  • ‖x − μₖ‖²:样本到所属质心的欧氏距离平方
  • J 的物理意义:Within-Cluster Sum of Squares (WCSS),即“簇内误差平方和”

> 算法目标:找到使 J 最小的簇划分 {C₁,…,Cₖ} 与质心 {μ₁,…,μₖ}


3. 迭代两步:坐标下降求 J

K-Means 采用坐标下降策略,交替更新两个变量:

步骤固定量优化量公式
E步 (Assignment)质心 μₖ样本归属 CₖCₖ = {x : ‖x − μₖ‖² ≤ ‖x − μⱼ‖², ∀j}
M步 (Update)簇 Cₖ质心 μₖμₖ = (1/Cₖ) ∑_{x∈Cₖ} x

示例

1def kmeans(X, K, max_iter=100):
2    n, d = X.shape
3    mu = X[torch.randperm(n)[:K]]          # 随机初始化 K 个质心
4    for _ in range(max_iter):
5        # E步:计算距离并分配样本
6        dist = torch.cdist(X, mu)           # (n, K)
7        labels = torch.argmin(dist, dim=1)  # (n,)
8        # M步:重新计算质心
9        for k in range(K):
10            mask = labels == k
11            if mask.sum() > 0:
12                mu[k] = X[mask].mean(dim=0)
13    return labels, mu
14

K-Means 聚类的目标函数:簇内误差平方和》 是转载文章,点击查看原文


相关推荐


STM32F4 TIM定时器HAL库源码完全解析
冉卓电子2026/6/9

一、前言    在STM32嵌入式开发中,TIM定时器是使用频率最高、功能最丰富的外设之一,几乎所有嵌入式项目都会用到定时中断、PWM输出、输入捕获、编码器解码等功能,而这些功能的底层支撑都是TIM外设。   很多开发者只会调用HAL库现成函数,却不懂底层实现逻辑,遇到定时器卡死、PWM波形异常、捕获数据不准、DMA传输失效等问题时无从排查。   本文基于STM32F4xx HAL库 tim.c 完整源码,保姆级拆解TIM定时器全套底层架构、功能模块、核心函数、回调机制与使用规范,带你从源


【GaussDB】会话里出现大量idle in transaction状态的问题排查
DarkAthena2026/6/1

【GaussDB】会话里出现大量idle in transaction状态的问题排查 背景 客户DBA给应用开发人员宣讲GaussDB数据库的idle in transaction会话的监控告警后,开发人员自查,发现在开发环境中存在大量的idle in transaction会话,一时恐慌,于是来找DBA帮忙排查是不是应用代码哪里写得有问题。 数据库内核版本是 506.0 SPC0500 集中式。 检查应用 应用软件使用的spring框架接mybaitis,连接使用基于DBCP修改的连


Milvus 实战:当 RAG 遇上向量数据库,从"玩具 Demo"到"生产可用的"那一步
Lee川2026/5/25

Milvus 实战:当 RAG 遇上向量数据库,从"玩具 Demo"到"生产可用的"那一步 前面的文章讲过,RAG 的第四步是"向量存储"。在原型阶段,一个 MemoryVectorStore 就能跑通全链路。但问题是:程序一关,所有向量灰飞烟灭;数据一多,内存直接爆掉;更别提多用户隔离、权限控制、分布式检索这些生产环境的刚需。 MemoryVectorStore 是 Demo 的终点,却是向量数据库的起点。 本文通过一段与 Milvus(Zilliz Cloud)交互的真实代码,拆解向量数据库


RAG 系列(五):Embedding 模型——语义理解的核心
冬奇Lab2026/5/3

为什么换个 Embedding 模型,检索效果天差地别? 前面四篇文章,我们搞定了 Pipeline 搭建、参数调优和分块策略。但有一个问题一直没细说: 你的文档被切成 Chunk 之后,是怎么变成向量的? 这个过程叫 Embedding(嵌入),它把人类可读的文本变成计算机可算的向量。Embedding 模型的选择,直接决定了: "苹果"和"iPhone"能不能被识别为相关 "数据库连接池耗尽"和"Too many connections"能不能被匹配到一起 中文成语、专业术语、缩写能不


Python全栈项目实战:自建高效多媒体处理工具
天天进步20152026/4/24

在数字化时代,视频剪辑、格式转换、音频提取等需求已成为日常。虽然市面上有很多成熟的工具,但作为开发者,**亲手构建一个属于自己的“全栈多媒体处理平台”**不仅能深度掌握 Python 生态,还能解决隐私安全和批量化定制的痛点。 本博文将带你梳理一个 Python 全栈多媒体处理工具的核心设计与实现方案。 一、 项目核心功能 一个实用的多媒体工具至少应具备以下“硬核”功能: 视频处理:格式转换(MP4/WebM/AVI)、视频抽帧、添加水印、调整分辨率。 音频处理:音频提取、格式压


Frida 源码编译全流程:自己动手编译 frida-server
CYRUS_STUDIO2026/4/15

版权归作者所有,如有转发,请注明文章出处:cyrus-studio.github.io/blog/ 下载 Frida 源码 Frida 源码:github.com/frida/frida 官方文档:frida.re/docs/buildi… 下载源码 git clone https://github.com/frida/frida.git 安装相关依赖: sudo apt-get install build-essential git lib32stdc++-9-dev \ libc


用 3100 个数字造一台计算机
jump_jump2026/4/7

你有没有想过,一台计算机最少需要什么? 不是说你桌上那台——那个有几十亿个晶体管、跑着操作系统和浏览器的庞然大物。我说的是最本质的那个东西:能算数、能画画、能放音乐、能响应你的键盘和鼠标。 答案可能会让你意外:一个数组就够了。 Little Virtual Computer 是一台用 TypeScript 写的虚拟计算机,原作者是 jsdf。我在他的基础上做了不少重构和优化——把代码拆分成了清晰的模块结构,加了音频系统、断点调试、内存追踪、中英文切换等功能。3100 个内存槽位,23 条指令,你


别再被误导!try...catch性能大揭秘
小码哥_常2026/3/30

别再被误导!try...catch性能大揭秘 开头:抛出问题,引发好奇 家人们,最近我在代码审查的时候,被狠狠质疑了一把。我在代码里用了好些try...catch,结果就收到了这样的意见:“try...catch用太多会影响性能,得优化一下”。当时我就在想,try...catch真有这么大罪过吗?平常开发的时候,我们为了处理各种可能出现的异常,try...catch可没少用,它真的会严重影响性能吗 ?今天咱就来好好唠唠这个话题,一起把这层迷雾给拨开! 历史担忧:曾经的性能痛点 (一)早期 Jav


Claude Skills 新手笔记
言萧凡_CookieBoty2026/3/22

开篇:这份笔记能帮你解决什么问题 如果你正在用 Claude(或任何大模型)做写作、总结、审阅、编码、整理资料等重复性工作,你很快会遇到两个痛点: 同一类任务,每次都要把背景、规则、格式重新讲一遍; 输出质量容易波动:有时很稳,有时跑偏,难以复用。 Claude Skills 可以把“可重复任务的知识 + 流程 +资源”打包成可复用的能力单元,让模型在需要时按需加载,从而更稳定地完成特定工作。这篇文章会用最小规范、加载机制、以及一个完整的中文示例,带你快速建立可落地的理解。 1. Clau


破解企业安全软件网络拦截实战记录
来点vc2026/3/14

一、背景与现象 问题描述 个人电脑访问某技术网站时,被企业安全软件(零信任安全架构)拦截,提示: "您访问的站点或应用因不合规,企业安全软件 已自动拦截此次访问" 关键特征: 这是个人电脑,未连接公司 VPN,未访问公司内网 浏览器和系统设置中看不到任何代理配置 重启电脑后问题依然存在 二、排查与分析过程 第一阶段:DNS 层面排查 初始假设:DNS 被劫持到企业服务器 $ scutil --dns | grep nameserver nameserver[0] : 10.x.x.x

首页编辑器站点地图

本站内容在 CC BY-SA 4.0 协议下发布

Copyright © 2026 聚合阅读