聚类小练习1

    xiaoxiao2025-10-29  8

    数据来源为华为杯全国研究生数学建模比赛2015年B题第一问,它提供的数据格式为.mat(maltab默认保存数据类型)。

    数据一共有200组,每组数据有100个特征。

    编写运用maltab。

    clear clc load 1.mat; sud=0; jig=ones(200,200); d=ones(200,1); for i=1:200;     for j=1:200;         sd=sum((data(:,i)-data(:,j)).^2);         sud=sd+sud;         jig(i,j)=sqrt(sd)/2;     end end for i=1:200;     ac=sort(jig(i,:));     d(i)=ac(200); end md=sort(d); maxd=md(200); [c,d]=find(jig==maxd); err=0.01; dd=ones([1,2]); u=zeros([2,200]); errma=1; E=0; n=1; while errma>err;     n=n+1;     if n==1000000;         break;     end     for i=1:200;         for j=1:2;         dd(1,j)=jig(c(j),i);         if j==2 && dd(1,j)>dd(1,1);              u(1,i)=1;         end         if j==2 && dd(1,j)<=dd(1,1);           u(2,i)=1;         end         end     end     adds=0;     for j=1:2;        adds=adds+jig(c(j),:)*u(j,:)';     end     errma=abs(adds-E);     E=adds;     for j=1:2;         sum1=[];      dd1=find(u(j,:)==1);       for k=1:length(dd1);           sum2=0;           for i=1:length(dd1);              sum2=sum2+jig(dd1(k),dd1(i));           end           sum1=[sum1 sum2];       end       mb=sort(sum1);       dd2=find(sum1==mb(1));       c(i)=dd1(dd2);     end   end

    最终结果为

    第一类(组号)

    41    42    43    44    45    46    47    48    49    50    51    52    53    54    55    56    57    58    59    60    61    62    63 64    65    66    67    68    69    70    71    72    73    74    75    76    77    78    79    80    81    82    83    84    85    86 87    88    89    90    91    92    93    94    95    96    97    98    99   100   101   102   103   104   105   106   107   108   109 110   111   112   113   114   115   116   117   118   119   120   121   122   123   124   125   126   127   128   129   130   131   132 133   134   135   136   137   138   139   140 第二类(组号)   1     2     3     4     5     6     7     8     9    10    11    12    13    14    15    16    17    18    19    20    21    22    23 24    25    26    27    28    29    30    31    32    33    34    35    36    37    38    39    40   141   142   143   144   145   14 147   148   149   150   151   152   153   154   155   156   157   158   159   160   161   162   163   164   165   166   167   168   169 170   171   172   173   174   175   176   177   178   179   180   181   182   183   184   185   186   187   188   189   190   191   192 193   194   195   196   197   198   199   200

    转载请注明原文地址: https://ju.6miu.com/read-1303628.html
    最新回复(0)