kmeanlbg.m

patten regnization source从1-14章能运行

function [x,esq,j] = kmeanlbg(d,k)
%KMEANLBG Linde-Buzo-Gray 舅绊府硫阑 捞侩茄 氦磐 剧磊拳 [X,ESQ,J]=(D,K)
%涝仿 牢磊:
% D 单捞磐 氦磐(one per row)
% K 吝缴 肮荐
%
%免仿 牢磊:
% X 免仿 凯 氦磐(K 俺狼 凯)
% ESQ 乞闭 磊铰 坷瞒(mean square error)
% J 阿 单捞磐 氦磐啊 绢蠢 吝缴俊 加窍绰啊狼 祸牢 
%
%
%  Implements LBG K-means algorithm:
%  Linde, Y., A. Buzo, and R. M. Gray,
%  "An Algorithm for vector quantiser design,"
%  IEEE Trans Communications, vol. 28, pp.84-95, Jan 1980.
%  VOICEBOX 傍俺 家胶 荐沥
%  VOICEBOX home page: http://www.ee.ic.ac.uk/hp/staff/dmb/voicebox/voicebox.html
%
nc=size(d,2);
[x,esq,j]=kmeans(d,1);
m=1;
iter = 0;
while m<k
   iter = iter + 1;
   n=min(m,k-m);
   m=m+n;
   e=1e-4*sqrt(esq)*rand(1,nc);
   [x,esq,j]=kmeans(d,m,[x(1:n,:)+e(ones(n,1),:); x(1:n,:)-e(ones(n,1),:); x(n+1:m-n,:)]);

   % 唱传 单捞磐 焊咯林扁
   color = 'rgbkcmy';
   for i = 1:m
      r = (j == i);
      z = d(r,:);
      scatter(z(:,1),z(:,2),10,color(i));
   end
   scatter(x(:,1),x(:,2),'b','filled'); % 盎脚等 努矾胶磐 吝缴
   title([' LBG : 馆汗荐:',num2str(iter)]);
   pause
   drawnow;
end