pca_test.m

patten regnization source从1-14章能运行

% 葛电 函荐甸苞 弊覆甸阑 瘤款促. 
  clear all;
  close all;

% 鸥盔葛剧狼 盒器甫 啊柳蠢 单捞磐 痢甸阑 积己矫挪促. 
  x(1,:) = randn(1,100);
  x(2,:) = randn(1,100)*3;

% 焊咯林扁 困秦 鸥盔 葛剧 单捞磐 盒器甫 雀傈茄促. 
  [p(1,:),p(2,:)] = cart2pol(x(1,:),x(2,:));
  p(1,:) = p(1,:)-pi/3;
  [x(1,:),x(2,:)] = pol2cart(p(1,:),p(2,:));

% 单捞磐甫 敲吩茄促. 
  scatter(x(1,:),x(2,:));
  axis equal;
  drawnow;
  pause;

% 林己盒(PC)阑 拌魂窍促. 
  [pc, latent, explained] = pcacov(cov(x'));

% 单捞磐 惑俊 林己盒阑 弊赴促. 
  hold on;
  plot([-4 4]*pc(1,1),[-4 4]*pc(2,1),'r-');
  plot([-2 2]*pc(1,2),[-2 2]*pc(2,2),'g-');
  pause;

% 林己盒阑 绵栏肺 单捞磐甫 雀傈茄促.
  y = (x'*pc)';

% 单捞磐甫 敲吩茄促. 
  figure;
  scatter(y(1,:),y(2,:));
  axis equal;
  drawnow;
  pause;

% 绵惑栏肺 林己盒捞 初咯乐绰啊甫 犬牢窍扁 困秦辑 促矫 林己盒(PC)阑 拌魂茄促. 
  [pc2, latent, explained] = pcacov(cov(y'));

% 单捞磐 惑俊 林己盒阑 弊赴促.
  hold on;
  plot([-4 4]*pc2(1,1),[-4 4]*pc2(2,1),'r-');
  plot([-2 2]*pc2(1,2),[-2 2]*pc2(2,2),'g-');
  pause;

% 老瞒盔栏肺 绵家窍扁 滴锅掳 己盒阑 0栏肺 汲沥茄促. 
  y(2,:) = 0;

% 盔 单捞磐甫 开函券茄促. 
  x = (y'*inv(pc))';

% 单捞磐甫 敲吩茄促. 
  figure;
  scatter(x(1,:),x(2,:));
  axis equal;
  drawnow;
  pause;

% 肯己
  close all;