-
Notifications
You must be signed in to change notification settings - Fork 8
/
Dollarsign.m
55 lines (41 loc) · 1.2 KB
/
Dollarsign.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
function [X] = Dollarsign(~)
% X: generated data set
% m : label
n=16000;
noise=.03;
% choose S part over | part with prob based on area of each
% S is 3*pi, | is 6
p = (3*pi)/(6+3*pi);
X = zeros(n,3);
Y = zeros(n,1);
m = zeros(n,1);
for i = 1:n,
if rand() > p,
m(i,1) = 0;
% surface is along x=0, y=[0,1], z=[-3,3]
X(i,:) = [0 rand() 6*rand()-3];
Y(i,1) = 13 + X(i,3) + randn(); % will be from 10 to 16 +/- 2
else
m(i,1) = 1;
angle = 1.5*rand()*pi;
x = cos(angle);
y = rand();
z = sin(angle);
% choose between top and bottom of S
if rand() > 0.5,
X(i,:) = [x y z+1];
Y(i,1) = angle; % Y is from 0 to 3*pi
else,
X(i,:) = [-x y -z-1];
Y(i,1) = 3*pi - angle;
end
end
end
X(find(m),3) = X(find(m),3) + randn(length(find(m)),1)*noise;
Y = Y + randn(n,1)*0.1;
% randomly shuffle the points, so we can easily subsample later
randidx = randperm(n);
X = X(randidx,:);
Y = Y(randidx,:);
m = m(randidx,:);
% [group,path]=Path_Based_Cluster_LandMarks( D ,10, 70 , 15);