mvudemo.m

流形学习中的重要方法MVU的源代码

function mvudemo;
%mvuDemo tests mvu with a simple spiral consisting of 100 points.
%please accustom the code to your SDP solver (CSDP/SEDUMI)
%
%function mvuDemo;
%
%tt=linspace(0,4*pi,100);
%X=[tt.*sin(tt);tt.*cos(tt)];
%figure;
%subplot(2,1,1);
%scatter(X(1,:),X(2,:),'o',tt,'filled'); axis equal;
%title('Original');
%drawnow; 
%try
% [Y,D]=mvuCSDP(X);       % CSDP 
% fprintf('Seems like mvuCSDP is working!\n');
%catch
%  [Y,D]=mvuNT(X);           % SEDUMI 
%  fprintf('Seems like mvuSeDuMi is working!\n');
%end;
%subplot(2,1,2);
%scatter(Y(1,:),Y(2,:),'o',tt,'filled'); axis equal
%title('Reduced Dimensionality');

tt=linspace(0,4*pi,100);
X=[tt.*sin(tt);tt.*cos(tt)];
%X=X-repmat(mean(X,2),1,length(X));
figure;
pars.slack=1;

try
 fprintf('Computing distances...');
 Dis=distance(X);  
 fprintf('done\n');
catch
 error('ERROR! Are you sure distance.m is in the path?');
  
end;

subplot(2,1,1);
scatter(X(1,:),X(2,:),60,tt,'filled'); axis equal;
title('Original');
drawnow; 
try
 pars.solver=2;
 [Y,D]=mvu(Dis,3,pars);       % CSDP 
 fprintf('\n\nCSDP is working!\n');
catch
  pars.solver=1;
  [Y,D]=mvu(Dis,3,pars);           % SEDUMI 
  fprintf('\n\nCSDP does not seem to be installed correctly.\n');
  fprintf('SeDuMi is working!\n');
end;
subplot(2,1,2);
scatter(Y(1,:),Y(2,:),60,tt,'filled'); axis equal;
title('Reduced Dimensionality');