使用matlab进行降维

clc;
for i=1:168
    X(i,1:20)=totalPDSdat(i,5:24);
    
    labels(i,1)=0;
end

for i=169:386
    X(i,1:20)=totalNPdata(i-168,5:24);
    labels(i,1)=1;
end

item=10;  %迭代次数
sigma=1; %高斯函数方差

filename = ['-sigma-' num2str(sigma) '--' num2str(item) '.gif']; %必须预先建立gif文件  
 
% 估计本质维数，即最终降到的维数
no_dims = round(intrinsic_dim(X, 'MLE'));
disp(['MLE estimate of intrinsic dimensionality: ' num2str(no_dims)]);


for i = 1:item 
    [mappedX, mapping] = compute_mapping(X, 'Laplacian', no_dims,12,i);
    tname  = ['-sigma-' num2str(sigma) 'Result of Laplacian Eigenmaps'];  
    figure(i),
    scatter3(mappedX(:,1), mappedX(:,2), mappedX(:,3),5, labels(mapping.conn_comp),'filled');
    if isequal(labels(mapping.conn_comp),labels)
           sprintf('%s','sucess')
    else sprintf('%s','error')
   end

    title(tname);
    sigma=sigma+1;
    
    % 自动保存为gif图像  
    frame = getframe(i);  
    im = frame2im(frame);  
    [I,map] = rgb2ind(im,256);  
    if i==1  
        imwrite(I,map,filename,'gif','Loopcount',inf, 'DelayTime',0.8);  
    else  
        imwrite(I,map,filename,'gif','WriteMode','append','DelayTime',0.8);  
    end  
    
end

% Isomap降维
[mappedX, mapping] = compute_mapping(X, 'Isomap', no_dims);
if isequal(labels(mapping.conn_comp),labels)
        sprintf('%s','sucess')
else sprintf('%s','error')
end

figure
scatter3(mappedX(:,1), mappedX(:,2), mappedX(:,3),5, labels(mapping.conn_comp),'filled')
title('Result of Isomap')
drawnow