| 查看: 1155 | 回复: 0 | ||
[求助]
matlab的聚类分析
|
|
n = 300 ; c = 3; t = randperm(n) ; x = [randn(1,n/3)-2 randn(1,n/3) randn(1,n/3)+2 ; randn(1,n/3) randn(1,n/3)+4 randn(1,n/3)]' ; m = x(t(1:c),  ;x2 = sum(x.^2,2) ; s0(1:c,1)=inf; for o = 1 : 100 m2 = sum(m.^2,2) ; [d,y] = min(repmat(m2,1,n)+repmat(x2',c,1)-2*m*x') ; for j = 1 : c m(j, = mean(x(y==j,) ;s(j,1) = mean(d(y==j)) ; end % norm（A） 返回矩阵A的二范数 % 矩阵的二范数 ： A' * A 的最大特征值开平方 % A = [0 1 2 ; 3 4 5 ; 6 7 8] ; % norm（A） = sqrt(max(eig(A'*A))) = 14.2267 if norm(s-s0) < 0.001 break ; end s0 = s ; end figure(1) ;%创建图形窗口 clf ;%清除当前图像窗口 hold on ; plot(x(y==1,1),x(y==1,2),'bo') ; plot(x(y==2,1),x(y==2,2),'rx') ; plot(x(y==3,1),x(y==3,2),'gv') ; 本人初学者，有点看不懂，尤其是从for开始，不知道聚类为什么这样写，请赐教 |
回复此楼» 猜你喜欢
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有7人回复
-
申请2026年博士
已经有5人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有5人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有6人回复
-
2025冷门绝学什么时候出结果
已经有7人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有6人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有7人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
