swiss roll的LLE降维

clc
close all
clear all
%[Mn1]=size（I）
M=400;
%瑞士卷
t=（3*pi/2）*（1+2*rand（1M））;
h=21*rand（1M）;
X=[t.*cos（t）;h;t.*sin（t）];
x1=t.*cos（t）;
x2=t.*sin（t）;
ColorVector=t‘;
scatter3（X（1:）X（2:）X（3:）3ColorVector‘filled‘）;
%Toroidal------------------------------------------------------------------
%noisesigma=005;
%tt=（1:N）‘/N;
%tt=tt.^2*pi;
%figureplot3（（2+cos（8*tt））.*cos（tt）（2+cos（8*tt））.*sin（tt）sin（8*tt）‘bo‘‘linewidth‘1）;
%--------------------------------------------------------------------------
K=8;
d=2;
[DN]=size（X）;
%1）计算近邻
X2=sum（X.^21）;
distance=repmat（X2N1）+repmat（X2‘1N）-2*X‘*X;
[sortedindex]=sort（distance）;
neighborhood=index（2:（1+K）:）;%K个近邻
%2）
if K>d