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
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
|
#!/usr/bin/env python3
from random import seed, random
from time import time
import math
class X:
def __init__(self, x: float, y: float, z: float):
self.x_now = 0
self.y_now = 0
self.z_now = 0
self.x_k = x
self.y_k = y
self.z_k = z
self.x_k1 = 0
self.y_k1 = 0
self.z_k1 = 0
self.k = 0
def neighbour(self):
self.X_k_to_now()
self.x_now = self.x_now + (
random() / 100 if random() > 0.5 else -random() / 100
)
self.y_now = self.y_now + (
random() / 100 if random() > 0.5 else -random() / 100
)
self.z_now = self.z_now + (
random() / 100 if random() > 0.5 else -random() / 100
)
self.k = self.k + 1
def function(self, x: float, y: float, z: float):
return (
(x ** 5) - ((x ** 2) * y * z) + (z * x) + (y ** 2) - (z ** 3) - 10
)
def X_now(self):
return self.function(self.x_now, self.y_now, self.z_now)
def X_k(self):
return self.function(self.x_k, self.y_k, self.z_k)
def X_k_to_now(self):
self.x_now = self.x_k
self.y_now = self.y_k
self.z_now = self.z_k
def X_now_to_k1(self):
self.x_k1 = self.x_now
self.y_k1 = self.y_now
self.z_k1 = self.z_now
def X_k_to_k1(self):
self.x_k1 = self.x_k
self.y_k1 = self.y_k
self.z_k1 = self.z_k
def X_k1_to_k(self):
self.x_k = self.x_k1
self.y_k = self.y_k1
self.z_k = self.z_k1
def SA():
K_max = 600
T_zero = 30
alpha = 0.90
alpha_step = 20
x = X(1, 0, -1)
T = T_zero
seed(time())
for k in range(1, K_max):
x.neighbour()
if x.X_now() <= x.X_k():
x.X_now_to_k1()
else:
p = math.e ** ((-(x.X_now() - x.X_k())) / T)
seed(time())
r = random()
if p >= r:
x.X_now_to_k1()
else:
x.X_k_to_k1()
if k % alpha_step == 0:
T = T * alpha
x.X_k1_to_k()
print(x.x_k, x.y_k, x.z_k)
print("k=", x.X_k())
def main():
SA()
if __name__ == "__main__":
main()
|
