demo_isomap.m

基于测地距离不变性的非线性降维算法源码

load swiss_roll_data
D = L2_distance(X_data(:,1:1000), X_data(:,1:1000), 1); 

%To run Isomap with K = 7 neighbors/point, and produce embeddings
%in dimensions 1, 2, ..., 10, type these commands: 

options.dims = 1:10;
[Y, R, E] = Isomap(D, 'k', 7, options);