-
Notifications
You must be signed in to change notification settings - Fork 0
/
cec2009_ps.txt
74 lines (69 loc) · 1.58 KB
/
cec2009_ps.txt
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
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
x=(0:0.001:1)';
xx=ones(length(x),3);
xx(:,1)=x;
for i=1:length(x)
for j=2:3
xx(i,j)=sin(6*pi*xx(i,1)+j*pi/30);
end
end
plot3(a(:,1),a(:,2),a(:,3),'go',xx(:,1),xx(:,2),xx(:,3),'-b',a(1,1),a(1,2),a(1,3),'r*')
xlabel('x1');
ylabel('x2');
zlabel('x3');
title('UF1-PS');
x=(0:0.001:1)';
xx=ones(length(x),3);
xx(:,1)=x;
for i=1:length(x)
for j=2:3
if j%2==1
xx(i,j)=(0.3*xx(i,1)*xx(i,1)*cos(24*pi*xx(i,1)+4*j*pi/30)+0.6*xx(i,1))*cos(6*pi*xx(i,1)+j*pi/30);
else
xx(i,j)=(0.3*xx(i,1)*xx(i,1)*cos(24*pi*xx(i,1)+4*j*pi/30)+0.6*xx(i,1))*sin(6*pi*xx(i,1)+j*pi/30);
end
end
end
plot3(a(:,1),a(:,2),a(:,3),'go',xx(:,1),xx(:,2),xx(:,3),'-b',xx(1,1),xx(1,2),xx(1,3),'r*')
xlabel('x1');
ylabel('x2');
zlabel('x3');
title('UF2-PS');
x=(0:0.001:1)';
xx=ones(length(x),3);
xx(:,1)=x;
for i=1:length(x)
for j=2:3
xx(i,j)=xx(i,1)^(0.5*(1+3.0*(j-2)/(30-2)));
end
end
plot3(a(:,1),a(:,2),a(:,3),'go',xx(:,1),xx(:,2),xx(:,3),'-b',xx(1,1),xx(1,2),xx(1,3),'r*')
xlabel('x1');
ylabel('x2');
zlabel('x3');
title('UF3-PS');
x=(0:0.001:1)';
xx=ones(length(x),3);
xx(:,1)=x;
for i=1:length(x)
for j=2:3
xx(i,j)=sin(6*pi*xx(i,1)+j*pi/30);
end
end
plot3(a(:,1),a(:,2),a(:,3),'go',xx(:,1),xx(:,2),xx(:,3),'-b',xx(1,1),xx(1,2),xx(1,3),'r*')
xlabel('x1');
ylabel('x2');
zlabel('x3');
title('UF4-PS');
x=(0:0.001:1)';
xx=ones(length(x),3);
xx(:,1)=x;
for i=1:length(x)
for j=2:3
xx(i,j)=sin(6*pi*xx(i,1)+j*pi/30);
end
end
plot3(a(:,1),a(:,2),a(:,3),'go',xx(:,1),xx(:,2),xx(:,3),'-b',xx(1,1),xx(1,2),xx(1,3),'r*')
xlabel('x1');
ylabel('x2');
zlabel('x3');
title('UF7-PS');