demo1.m

来自「流形学习中的重要方法MVU的源代码」· M 代码 · 共 32 行

tt=linspace(0,4*pi,100);
X=[tt.*sin(tt);tt.*cos(tt)];
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=0;
[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');