无监督学习k-means简单实现

%随机获取150个点%X = [randn(50,2)+ones(50,2);randn(50,2)-ones(50,2);randn(50,2)+[ones(50,1),-ones(50,1)]];X = load('test.txt') %二维高斯拟合函数opts = statset('Display','final'); %调用Kmeans函数%X N*P的数据矩阵%Idx N*1的向量,存储的是每个点的聚类标号%Ctrs K*P的矩阵,存储的是K个聚类质心位置%SumD 1*K的和向量,存储的是类间所有点与该类质心点距离之和%D N*K的矩阵，存储的是每个点与所有质心的距离; [Idx,Ctrs,SumD,D] = kmeans(X,3,'Replicates',3,'Options',opts); %画出聚类为1的点。X(Idx==1,1),为第一类的样本的第一个坐标；X(Idx==1,2)为第二类的样本的第二个坐标plot(X(Idx==1,1),X(Idx==1,2),'r.','MarkerSize',14)hold onplot(X(Idx==2,1),X(Idx==2,2),'b.','MarkerSize',14)hold onplot(X(Idx==3,1),X(Idx==3,2),'g.','MarkerSize',14) %绘出聚类中心点,kx表示是圆形plot(Ctrs(:,1),Ctrs(:,2),'kx','MarkerSize',14,'LineWidth',4)plot(Ctrs(:,1),Ctrs(:,2),'kx','MarkerSize',14,'LineWidth',4)plot(Ctrs(:,1),Ctrs(:,2),'kx','MarkerSize',14,'LineWidth',4) legend('Cluster 1','Cluster 2','Cluster 3','Centroids','Location','NW') 

数据集test.txt

14.9,27.5514.45,28.3514.15,28.813.75,28.0513.35,28.4513,29.1513.45,27.513.6,26.512.8,27.3512.4,27.8512.3,28.412.2,28.6513.4,25.112.95,25.9512.9,26.511.85,2711.35,2811.15,28.711.25,27.410.75,27.710.5,28.359.65,28.4510.25,27.2510.75,26.5511.7,26.3511.6,25.911.9,25.0512.6,24.0511.9,24.511.1,25.210.55,25.1510.05,25.959.35,26.69.3,27.259.2,27.87.5,28.258.55,27.458.5,27.058.05,27.27.85,26.87.3,27.46.8,26.857,26.57.55,26.38.55,26.39,25.858.6,25.659.4,25.558.45,25.058.85,24.69.65,24.710.55,24.3511.05,23.910.55,23.559.45,23.359.2,23.98.35,23.97.35,24.757.4,25.456.6,25.756.1,265.8,26.955.65,25.85.3,26.16.4,25.45.4,25.255.35,24.74.8,25.054.2,25.556.4,24.86.55,24.37.4,24.255.45,24.24.3,244,24.253.35,23.34.85,23.054.3,22.755.85,23.45.9,23.557.55,23.76.85,23.257.65,23.16.95,22.556.1,22.65.5,22.64.7,22.13.8,21.854.65,21.24.15,20.355.3,20.45.6,20.755.8,21.956.4,21.956.55,21.157.45,21.957.4,21.557.75,21.27.65,20.656.95,19.86.6,20.16.05,20.25.4,19.655.35,19.055.8,18.256.3,19.17,18.97.15,17.97.35,18.28.2,20.058.3,19.458.3,18.58.75,18.89.05,18.29.35,17.78.9,17.658.45,17.210.05,17.210.4,16.758.6,20.98.65,21.38.65,21.98.65,22.58.95,22.89.95,22.658.95,22.29.65,21.910.55,22.310.9,22.8511.35,23.4512.05,23.412.3,22.7511.7,22.1511.15,22.0510.85,21.510.85,21.059.6,21.39.85,20.79.35,20.69.25,19.659.95,19.810.7,20.3511.3,20.712.35,21.613.1,21.312.85,20.7512,2011,19.8510.35,199.9,18.6510.6,18.1511.4,18.311.4,19.2512.35,18.812.8,19.7512.15,18.111.05,17.511.95,17.2512.25,17.513.05,17.413.75,18.1513.5,18.6513.65,19.2514,19.915.2,18.215.5,17.1513.9,17.113.75,16.612.15,16.47.8,13.78.85,13.359,12.79.7,12.18.05,12.97.7,13.256.8,13.26.6,13.456.2,12.555.4,12.855.7,12.255.2,11.95.15,11.355.85,11.26.1,11.757,12.357.05,12.457.9,12.58.55,12.17.85,11.857.1,11.956.9,11.56.85,10.96.4,10.75.9,10.36.4,10.257.05,10.057.35,10.57.65,11.18.1,11.28.8,11.48.3,10.559,10.99.35,10.510.15,1110.4,10.5510.9,1011.55,10.211.75,10.8510.1,8.6511.05,9.111.85,9.812.85,10.6512.9,11.713.6,11.114.05,11.7514.5,11.814.3,12.4517,12.915.8,12.615.85,1216.7,12.216.25,11.715.55,11.1514.8,11.3514.45,10.7513.75,10.4512.8,10.113.15,9.812.45,9.311.8,8.9511.1,8.4510.35,7.710.1,6.7511.3,7.9512.35,8.4513.1,8.9513.2,9.3514.1,10.0511.5,7.511.35,6.911.95,6.7512.4,7.112.25,7.612.95,7.613.45,7.9513.35,8.2513.75,914.3,9.314.85,9.5515.1,10.2515.45,10.5516.35,10.8516.75,11.516.25,10.215.4,10.115.45,9.715.15,9.315.25,8.6515.55,8.214.25,8.714.25,8.2515.05,7.814.3,7.513.55,7.4514.3,6.9513.95,6.713.05,6.9513.05,6.211.55,6.310.8,5.8510.6,5.0511.35,5.5512.15,5.412.4,5.812.8,5.713.65,5.913.9,5.313.1,5.112.55,4.911.5,4.7511.35,4.0512.4,4.3511.75,3.4512.65,3.713.4,4.3513.9,4.9512.75,313.55,3.1513.7,3.6514.1,4.114.65,5.0514.35,5.7514.5,6.5515.15,7.113.6,2.5514.45,2.414.6,3.0515,3.415.25,3.514.7,4.114.7,4.515.25,2.715.65,2.0515.95,2.816.1,3.5515.9,415.6,4.7515.55,5.0515.35,5.515.15,5.9515.5,6.7515.7,6.3516.2,5.916.35,5.3516.2,4.5516.55,4.216.95,4.7517.05,5.117.3,4.817.3,4.1517.6,4.317.05,3.717.25,3.0516.65,2.816.55,2.1517.2,2.0518.15,1.9518.05,2.4518.15,3.0518.6,3.4518.4,3.618.85,3.219.1,2.6519.45,2.6519,2.119.9,2.0520.45,2.819.8,3.2519.45,3.918.65,4.218.4,4.618.65,4.7518.75,5.1519.1,4.5517.9,5.417.65,5.717.05,6.0517.4,6.516.6,6.8515.7,7.1515.75,7.7516.6,7.9520.4,3.420.7,3.4521.15,2.8521.75,2.6522,3.2522.2,3.521.45,3.7521.1,4.0520.15,4.320.8,4.720.7,5.1519.75,5.0519.85,5.520.4,5.6520.55,5.7518.7,5.7519.25,5.9518.4,618.45,6.617.65,7.0516.7,7.418.65,7.318.05,7.3517.85,7.7517.5,8.2517.15,8.617.05,916.4,8.716.05,8.9516.05,9.616.5,9.7517.25,9.617.6,9.917.8,9.318,8.5518.8,8.118.8,8.3519.4,7.619.25,6.620.05,6.9519.8,7.520.05,6.3521.15,5.721.65,4.8522.15,4.3523.05,3.3523.05,3.823.15,4.422.5,4.7522.15,5.224.15,4.5523.5,5.0523.1,5.323,5.7522.2,5.7521.85,6.220.75,6.5521,7.1520.75,7.6520,8.219.5,8.6518.85,9.0518.75,9.5518.6,1016.95,10.3517.35,10.8518,10.6518.5,10.5518.1,11.117.55,11.317.95,11.918.3,1218,12.519,11.6519.5,11.0519.45,10.5519.4,9.6520.1,9.420.05,9.9520.05,10.219.35,12.219.2,12.2520.05,11.620.6,11.1520.7,10.6521.3,11.6521.8,11.1521.85,10.721.65,10.0520.95,10.220.9,9.721.65,9.4521.2,9.2520.75,8.7520.55,8.7521.1,821.65,8.6521.75,8.221.95,7.5522,6.7522.8,6.4522.65,6.6522.75,7.0523,7.3522.55,7.922.2,8.722.9,8.4522.35,9.222.75,9.3522.4,10.0523.05,10.923.3,9.8523.95,9.823.65,9.123.7,8.8524.25,8.2524.85,7.9523.5,7.8523.85,7.3523.95,6.923.65,6.523.6,5.724.3,5.6524.8,6.434.05,3.533.05,3.8532,3.831.9,4.431.05,4.7530.4,5.6530.75,6.130,6.730.1,7.429.5,8.1530.75,830.85,7.3531.5,6.7531.75,5.9532.35,6.4532.8,632.05,5.132.8,4.832.65,4.433.65,4.633.05,5.1533.6,5.4534.5,5.0534.9,4.6535.45,4.134.6,4.0534.2,4.236.3,5.235.55,5.3535.95,6.0534.8,5.8533.7,6.1533.95,6.633.7,7.0532.75,7.132.3,7.6533,7.931.95,8.1531.15,8.6530.35,8.8529.85,930.7,9.1529.7,9.930.45,9.9530.95,9.8531.8,9.4532.45,8.833.55,8.634.35,7.734.7,834.6,7.2535,6.835.5,7.3536.1,7.536.55,736,8.235.35,8.0536.55,8.6536.4,9.135.5,9.134.55,8.8535.25,9.434.4,9.533.5,9.333.85,9.832.5,9.6532.3,10.2533.3,10.331.6,10.530.6,10.530.4,11.130.9,11.4530.7,11.6530.4,12.0531.2,1231.95,11.3531.65,11.0532.95,11.1532.65,11.732.25,12.2532.05,12.2531.3,12.731.95,12.9532.75,13.133.15,13.233.1,12.7533.15,12.134.3,11.7534,10.8534.65,1134.8,10.135.65,9.8536.35,1035.55,10.7535.8,11.5535.2,11.7534.7,11.7534.95,12.7534.05,12.5534.05,13.0533.25,13.733.2,14.1533.25,14.733,15.1532.95,15.6532.6,16.1532.45,16.7532.65,17.0532.75,17.331.75,17.231.7,17.6531,17.531.15,17.930.45,18.0530.05,18.830.55,18.830.5,19.330.25,19.429.6,19.8529.15,20.5530.25,20.4530.7,20.0531,19.931.55,19.6531.5,18.5532.05,18.631.95,19.132.6,18.3532.85,17.9533.45,17.4533.7,1734.25,17.3534.3,18.0533.85,18.433.05,18.8533.25,18.9532.8,19.1532.3,19.8532.8,2031.75,20.2531.75,20.7532.1,21.1531.55,21.630.65,21.329.95,21.629.5,21.630.35,22.0531.05,2231.55,22.230.95,22.6530.3,23.129.6,23.1529.35,22.5529.2,23.8530.75,2430.95,24.1531.45,23.731.95,23.1532.55,22.0532.6,22.5533.25,22.2533.65,21.933.5,21.333.1,20.7533.8,20.433.85,2034.15,19.334.85,18.8535.3,18.5535.4,19.3534.55,19.7535.05,2035.95,19.8536.35,20.635.5,20.5534.45,20.6534.4,21.2535,21.0535.75,21.335.05,21.534.6,22.0534.2,21.7536.25,21.9535.7,22.335.5,22.935.85,23.2536.3,23.835.45,24.134.9,23.534.2,22.933.85,23.333.25,23.3532.45,23.733.6,23.934.25,23.9534.25,24.135.4,24.735.15,25.334.4,24.933.7,24.8532.25,24.4532.5,24.731.45,24.4531.55,25.231.05,2530.25,24.329.8,24.829.6,25.529.7,26.0530.5,25.530.65,2631.25,26.0531.45,26.9530.75,26.930.65,27.1531.25,27.8531.85,27.7532.7,28.233.25,27.5532.4,27.132.15,26.6532.35,25.9532.95,25.533.85,26.0533.05,26.533.65,2734.1,27.3534.2,27.9534.65,26.8535.25,2635.7,26.1534.4,25.621.3,20.820.15,20.919.2,21.3519.1,21.8518.45,22.818.75,22.9519.4,2319.55,22.2519.8,21.8520.5,21.8521.45,21.4521.7,21.921.4,22.321,22.621.15,22.9520.5,22.8519.75,23.6519.2,23.718.45,24.3520.65,23.8520.65,24.319.7,24.620.15,25.0522.15,25.121.6,24.6521.7,23.821.9,23.6522.55,23.522.55,24.323.3,24.4524.25,24.3523.8,25.2523.4,23.822.9,23.222.3,22.822.2,22.423.1,21.722.85,21.922.65,21.123.15,22.624.1,21.924.7,22.224.3,22.624.15,23.323.9,23.455.2,2.156.35,1.956.75,2.35.9,2.45.4,2.74.85,2.94.85,3.355.15,3.455.7,3.456.2,36.2,3.27.65,2.157.2,2.756.75,3.26.75,3.556.65,3.85.8,44.95,4.055.1,4.355.7,4.455.45,4.856.7,4.86.55,5.057.2,4.96.2,4.257.1,4.37.85,4.57.6,4.157.25,3.557.8,3.358.05,2.758.5,3.258.1,3.558.15,4