OptimizationForEveryone/cvrt.py at main · Optimization-for-Everyone/OptimizationForEveryone · GitHub

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
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
import numpy as np
from numpy import sin, cos, tan ,cosh, tanh, sinh, abs, exp, mean, pi, prod, sqrt, sum, floor

from sympy import *
import sympy as sym

t=sym.symbols('t')
#t=1

x =[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]
#x = np.random.uniform(0,1,(27,))
print(x)
#k=len(x)/3
I_t=0
U_t=0
V_t=0
I_t_d=0
U_t_d=0
V_t_d=0

#print('I am exactly here')

for i in range(0,int(len(x)/3),3):
    #I_t=I_t+x[i]*(1/(1+pow(math.e,-(x[i+2]*t+x[i+1]))))
    I_t=I_t+x[i]*(1/(1+exp(-(x[i+2]*t+x[i+1]))))
    #print(I_t)

for i in range(int(len(x)/3),int(2*len(x)/3),3):
    #U_t=U_t+x[i]*(1/(1+pow(math.e,-(x[i+2]*t+x[i+1]))))
    U_t=U_t+x[i]*(1/(1+exp(-(x[i+2]*t+x[i+1]))))
    #print(U_t)

for i in range(int(2*len(x)/3),len(x),3):
    #V_t=V_t+x[i]*(1/(1+pow(math.e,-(x[i+2]*t+x[i+1]))))
    V_t=V_t+x[i]*(1/(1+exp(-(x[i+2]*t+x[i+1]))))
    #print(V_t)

for i in range(0,int(len(x)/3),3):
    #I_t_d=I_t_d+x[i]*x[i+2]*((pow(math.e,-(x[i+2]*t+x[i+1])))/(1+pow(math.e,-(x[i+2]*t+x[i+1])))**2)
    I_t_d=I_t_d+x[i]*x[i+2]*((exp(-(x[i+2]*t+x[i+1])))/(1+exp(-(x[i+2]*t+x[i+1])))**2)
    #print(I_t_d)

for i in range(int(len(x)/3),int(2*len(x)/3),3):
    #U_t_d=U_t_d+x[i]*x[i+2]*((pow(math.e,-(x[i+2]*t+x[i+1])))/(1+pow(math.e,-(x[i+2]*t+x[i+1])))**2)
    U_t_d=U_t_d+x[i]*x[i+2]*((exp(-(x[i+2]*t+x[i+1])))/(1+exp(-(x[i+2]*t+x[i+1])))**2)
    #print(U_t_d)

for i in range(int(2*len(x)/3),len(x),3):
    #V_t_d=V_t_d+x[i]*x[i+2]*((pow(math.e,-(x[i+2]*t+x[i+1])))/(1+pow(math.e,-(x[i+2]*t+x[i+1])))**2)
    V_t_d=V_t_d+x[i]*x[i+2]*((exp(-(x[i+2]*t+x[i+1])))/(1+exp(-(x[i+2]*t+x[i+1])))**2)
    #print(V_t_d)

print("I_t: ",I_t)
print("U_t: ",U_t)
print("V_t: ",V_t)
print("Derivatives: ")
print("I_t_d: ",I_t_d)
print("U_t_d: ",U_t_d)
print("V_t_d: ",V_t_d)

# tekrar
I_t=0
U_t=0
V_t=0
I_t_d=0
U_t_d=0
V_t_d=0
t=1
for i in range(0,int(len(x)/3),3):
    #I_t=I_t+x[i]*(1/(1+pow(math.e,-(x[i+2]*t+x[i+1]))))
    I_t=I_t+x[i]*(1/(1+exp(-(x[i+2]*t+x[i+1]))))
    #print(I_t)

for i in range(int(len(x)/3),int(2*len(x)/3),3):
    #U_t=U_t+x[i]*(1/(1+pow(math.e,-(x[i+2]*t+x[i+1]))))
    U_t=U_t+x[i]*(1/(1+exp(-(x[i+2]*t+x[i+1]))))
    #print(U_t)

for i in range(int(2*len(x)/3),len(x),3):
    #V_t=V_t+x[i]*(1/(1+pow(math.e,-(x[i+2]*t+x[i+1]))))
    V_t=V_t+x[i]*(1/(1+exp(-(x[i+2]*t+x[i+1]))))
    #print(V_t)

for i in range(0,int(len(x)/3),3):
    #I_t_d=I_t_d+x[i]*x[i+2]*((pow(math.e,-(x[i+2]*t+x[i+1])))/(1+pow(math.e,-(x[i+2]*t+x[i+1])))**2)
    I_t_d=I_t_d+x[i]*x[i+2]*((exp(-(x[i+2]*t+x[i+1])))/(1+exp(-(x[i+2]*t+x[i+1])))**2)
    #print(I_t_d)

for i in range(int(len(x)/3),int(2*len(x)/3),3):
    #U_t_d=U_t_d+x[i]*x[i+2]*((pow(math.e,-(x[i+2]*t+x[i+1])))/(1+pow(math.e,-(x[i+2]*t+x[i+1])))**2)
    U_t_d=U_t_d+x[i]*x[i+2]*((exp(-(x[i+2]*t+x[i+1])))/(1+exp(-(x[i+2]*t+x[i+1])))**2)
    #print(U_t_d)

for i in range(int(2*len(x)/3),len(x),3):
    #V_t_d=V_t_d+x[i]*x[i+2]*((pow(math.e,-(x[i+2]*t+x[i+1])))/(1+pow(math.e,-(x[i+2]*t+x[i+1])))**2)
    V_t_d=V_t_d+x[i]*x[i+2]*((exp(-(x[i+2]*t+x[i+1])))/(1+exp(-(x[i+2]*t+x[i+1])))**2)
    #print(V_t_d)

#print("I_t: ",I_t)
#print("U_t: ",U_t)
#print("V_t: ",V_t)
#print("Derivatives: ")
#print("I_t_d: ",I_t_d)
#print("U_t_d: ",U_t_d)
#print("V_t_d: ",V_t_d)

# end tekrar

#general parameters
#r=2*pow(10,-2)
#delta=0.7*pow(10,-9)
#alpha=0.00842
#b=50

#f=r/delta
#h=alpha/delta
#g=1/b
# end general parameters

#problem 1
f=0
h=0.15
g=0.02
#end problem 1
N=1

sum_EI=0
for i in range(N):
    sum_EI=sum_EI+(I_t_d-f*I_t+I_t*V_t)**2
EI=sum_EI/N

sum_EU=0
for i in range(N):
    sum_EU=sum_EU+(U_t_d-I_t*V_t+U_t)**2
EU=sum_EU/N

sum_EV=0
for i in range(N):
    sum_EV=sum_EV+(V_t_d-U_t+g*I_t*V_t+h*V_t)**2
EV=sum_EV/N

I0=0.6
U0=0
V0=0.5

EC=(1/3)*((I_t-I0)**2+(U_t-U0)**2+(V_t-V0)**2)

print("EI: ", EI)
print("EU: ", EU)
print("EV: ", EV)

E=EI+EU+EV+EC

print("E: ", E)