datagen.m

patten regnization source从1-14章能运行

function X=datagen(Nvec,mean_var)
% 2D 啊快矫救 盒器 单捞磐 积己
% Nvec: (1xclass) c俺狼 啊快矫救 盒器狼 阿阿俊 措茄 单捞磐 肮荐
% mean_var: (3 x class): 阿 努贰胶俊 措窍咯 贸澜 2 青篮 mean 唱赣瘤 1 青篮 variance

[m,c]=size(mean_var);
if m ~= 3 | c ~=length(Nvec),
   disp(' dimension not match, break ')
   break
end
X=[];
for i=1:c,   
   randn('seed',sum(100*clock));
   tmp=sqrt(mean_var(3,i))*randn(Nvec(i),2); % scaled by variance
   mean=mean_var(1:2,i);  % mean is a 2 by 1 vector
   X=[X;tmp+ones(Nvec(i),2)*diag(mean)];
end