-
Notifications
You must be signed in to change notification settings - Fork 0
/
LC_reactive_bound.m
150 lines (150 loc) · 6.46 KB
/
LC_reactive_bound.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
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
function au = LC_reactive_bound(x)
x_norm=[19;5;10;0.1746];
x_min=[-4;-1;-5;-0.0873];
y=(x-x_min)./x_norm(1:4);
m1=1/(1+exp(100*(y(1)-4.969518e-01)));
m2=1-m1;
ub = zeros(1,1);
ub(1)=(m1*(y(1)^0*y(2)^0*y(3)^0*y(4)^0*3.356047e-01+...
y(1)^0*y(2)^0*y(3)^0*y(4)^1*-7.270902e-03+...
y(1)^0*y(2)^0*y(3)^0*y(4)^2*-9.536271e-03+...
y(1)^0*y(2)^0*y(3)^0*y(4)^3*-1.356443e-03+...
y(1)^0*y(2)^0*y(3)^0*y(4)^4*8.594406e-03+...
y(1)^0*y(2)^0*y(3)^1*y(4)^0*8.080632e-02+...
y(1)^0*y(2)^0*y(3)^1*y(4)^1*1.622080e-02+...
y(1)^0*y(2)^0*y(3)^1*y(4)^2*-1.414444e-02+...
y(1)^0*y(2)^0*y(3)^1*y(4)^3*-2.289215e-02+...
y(1)^0*y(2)^0*y(3)^2*y(4)^0*6.292256e-02+...
y(1)^0*y(2)^0*y(3)^2*y(4)^1*2.618491e-02+...
y(1)^0*y(2)^0*y(3)^2*y(4)^2*5.779985e-04+...
y(1)^0*y(2)^0*y(3)^3*y(4)^0*2.920921e-02+...
y(1)^0*y(2)^0*y(3)^3*y(4)^1*1.889219e-02+...
y(1)^0*y(2)^0*y(3)^4*y(4)^0*7.703013e-04+...
y(1)^0*y(2)^1*y(3)^0*y(4)^0*-1.523522e-01+...
y(1)^0*y(2)^1*y(3)^0*y(4)^1*-4.858029e-02+...
y(1)^0*y(2)^1*y(3)^0*y(4)^2*8.658370e-03+...
y(1)^0*y(2)^1*y(3)^0*y(4)^3*3.994045e-02+...
y(1)^0*y(2)^1*y(3)^1*y(4)^0*-4.851016e-02+...
y(1)^0*y(2)^1*y(3)^1*y(4)^1*-2.825565e-02+...
y(1)^0*y(2)^1*y(3)^1*y(4)^2*-1.508591e-02+...
y(1)^0*y(2)^1*y(3)^2*y(4)^0*-8.362546e-03+...
y(1)^0*y(2)^1*y(3)^2*y(4)^1*1.725335e-03+...
y(1)^0*y(2)^1*y(3)^3*y(4)^0*-2.155234e-03+...
y(1)^0*y(2)^2*y(3)^0*y(4)^0*-1.245300e-01+...
y(1)^0*y(2)^2*y(3)^0*y(4)^1*-2.706406e-02+...
y(1)^0*y(2)^2*y(3)^0*y(4)^2*2.842074e-02+...
y(1)^0*y(2)^2*y(3)^1*y(4)^0*-4.782995e-02+...
y(1)^0*y(2)^2*y(3)^1*y(4)^1*-2.271222e-02+...
y(1)^0*y(2)^2*y(3)^2*y(4)^0*-8.902848e-03+...
y(1)^0*y(2)^3*y(3)^0*y(4)^0*-5.939569e-02+...
y(1)^0*y(2)^3*y(3)^0*y(4)^1*6.882775e-03+...
y(1)^0*y(2)^3*y(3)^1*y(4)^0*-1.845722e-02+...
y(1)^0*y(2)^4*y(3)^0*y(4)^0*1.986278e-03+...
y(1)^1*y(2)^0*y(3)^0*y(4)^0*1.585205e-01+...
y(1)^1*y(2)^0*y(3)^0*y(4)^1*5.645107e-02+...
y(1)^1*y(2)^0*y(3)^0*y(4)^2*1.277367e-02+...
y(1)^1*y(2)^0*y(3)^0*y(4)^3*-5.250584e-03+...
y(1)^1*y(2)^0*y(3)^1*y(4)^0*1.162574e-01+...
y(1)^1*y(2)^0*y(3)^1*y(4)^1*4.379941e-02+...
y(1)^1*y(2)^0*y(3)^1*y(4)^2*1.203902e-02+...
y(1)^1*y(2)^0*y(3)^2*y(4)^0*7.140376e-02+...
y(1)^1*y(2)^0*y(3)^2*y(4)^1*2.702971e-02+...
y(1)^1*y(2)^0*y(3)^3*y(4)^0*4.035879e-02+...
y(1)^1*y(2)^1*y(3)^0*y(4)^0*-2.267390e-02+...
y(1)^1*y(2)^1*y(3)^0*y(4)^1*-1.346023e-02+...
y(1)^1*y(2)^1*y(3)^0*y(4)^2*-1.164254e-02+...
y(1)^1*y(2)^1*y(3)^1*y(4)^0*2.347800e-03+...
y(1)^1*y(2)^1*y(3)^1*y(4)^1*-2.751996e-03+...
y(1)^1*y(2)^1*y(3)^2*y(4)^0*1.691771e-03+...
y(1)^1*y(2)^2*y(3)^0*y(4)^0*-5.158867e-02+...
y(1)^1*y(2)^2*y(3)^0*y(4)^1*-2.325297e-02+...
y(1)^1*y(2)^2*y(3)^1*y(4)^0*-1.815789e-02+...
y(1)^1*y(2)^3*y(3)^0*y(4)^0*-4.839067e-02+...
y(1)^2*y(2)^0*y(3)^0*y(4)^0*1.078854e-01+...
y(1)^2*y(2)^0*y(3)^0*y(4)^1*4.173623e-02+...
y(1)^2*y(2)^0*y(3)^0*y(4)^2*1.380640e-02+...
y(1)^2*y(2)^0*y(3)^1*y(4)^0*6.992502e-02+...
y(1)^2*y(2)^0*y(3)^1*y(4)^1*2.792343e-02+...
y(1)^2*y(2)^0*y(3)^2*y(4)^0*4.237808e-02+...
y(1)^2*y(2)^1*y(3)^0*y(4)^0*5.262125e-03+...
y(1)^2*y(2)^1*y(3)^0*y(4)^1*1.044248e-03+...
y(1)^2*y(2)^1*y(3)^1*y(4)^0*7.691514e-03+...
y(1)^2*y(2)^2*y(3)^0*y(4)^0*-1.512666e-02+...
y(1)^3*y(2)^0*y(3)^0*y(4)^0*6.079719e-02+...
y(1)^3*y(2)^0*y(3)^0*y(4)^1*2.378575e-02+...
y(1)^3*y(2)^0*y(3)^1*y(4)^0*3.756700e-02+...
y(1)^3*y(2)^1*y(3)^0*y(4)^0*6.818328e-03+...
y(1)^4*y(2)^0*y(3)^0*y(4)^0*3.249668e-02)...
+m2*(y(1)^0*y(2)^0*y(3)^0*y(4)^0*5.875964e-01+...
y(1)^0*y(2)^0*y(3)^0*y(4)^1*4.916670e-02+...
y(1)^0*y(2)^0*y(3)^0*y(4)^2*-2.950654e-02+...
y(1)^0*y(2)^0*y(3)^0*y(4)^3*-7.120732e-02+...
y(1)^0*y(2)^0*y(3)^0*y(4)^4*-7.776639e-02+...
y(1)^0*y(2)^0*y(3)^1*y(4)^0*-1.287621e-02+...
y(1)^0*y(2)^0*y(3)^1*y(4)^1*2.702638e-02+...
y(1)^0*y(2)^0*y(3)^1*y(4)^2*1.330900e-03+...
y(1)^0*y(2)^0*y(3)^1*y(4)^3*-1.608484e-02+...
y(1)^0*y(2)^0*y(3)^2*y(4)^0*-4.538225e-02+...
y(1)^0*y(2)^0*y(3)^2*y(4)^1*-6.454523e-04+...
y(1)^0*y(2)^0*y(3)^2*y(4)^2*-6.372252e-04+...
y(1)^0*y(2)^0*y(3)^3*y(4)^0*-5.339222e-02+...
y(1)^0*y(2)^0*y(3)^3*y(4)^1*-1.246365e-02+...
y(1)^0*y(2)^0*y(3)^4*y(4)^0*-4.720078e-02+...
y(1)^0*y(2)^1*y(3)^0*y(4)^0*-3.393216e-02+...
y(1)^0*y(2)^1*y(3)^0*y(4)^1*7.707006e-03+...
y(1)^0*y(2)^1*y(3)^0*y(4)^2*-4.989912e-02+...
y(1)^0*y(2)^1*y(3)^0*y(4)^3*-7.661481e-02+...
y(1)^0*y(2)^1*y(3)^1*y(4)^0*-3.227421e-04+...
y(1)^0*y(2)^1*y(3)^1*y(4)^1*1.811629e-02+...
y(1)^0*y(2)^1*y(3)^1*y(4)^2*-5.572832e-03+...
y(1)^0*y(2)^1*y(3)^2*y(4)^0*-2.311986e-02+...
y(1)^0*y(2)^1*y(3)^2*y(4)^1*1.896969e-03+...
y(1)^0*y(2)^1*y(3)^3*y(4)^0*-3.458284e-02+...
y(1)^0*y(2)^2*y(3)^0*y(4)^0*-2.742228e-02+...
y(1)^0*y(2)^2*y(3)^0*y(4)^1*5.136801e-03+...
y(1)^0*y(2)^2*y(3)^0*y(4)^2*-4.751093e-02+...
y(1)^0*y(2)^2*y(3)^1*y(4)^0*1.367947e-02+...
y(1)^0*y(2)^2*y(3)^1*y(4)^1*2.047385e-02+...
y(1)^0*y(2)^2*y(3)^2*y(4)^0*-8.799415e-03+...
y(1)^0*y(2)^3*y(3)^0*y(4)^0*2.158268e-02+...
y(1)^0*y(2)^3*y(3)^0*y(4)^1*2.215696e-02+...
y(1)^0*y(2)^3*y(3)^1*y(4)^0*4.035253e-02+...
y(1)^0*y(2)^4*y(3)^0*y(4)^0*6.390618e-02+...
y(1)^1*y(2)^0*y(3)^0*y(4)^0*2.256822e-02+...
y(1)^1*y(2)^0*y(3)^0*y(4)^1*6.213375e-02+...
y(1)^1*y(2)^0*y(3)^0*y(4)^2*-1.223603e-02+...
y(1)^1*y(2)^0*y(3)^0*y(4)^3*-5.260043e-02+...
y(1)^1*y(2)^0*y(3)^1*y(4)^0*1.898383e-02+...
y(1)^1*y(2)^0*y(3)^1*y(4)^1*3.751403e-02+...
y(1)^1*y(2)^0*y(3)^1*y(4)^2*8.013790e-03+...
y(1)^1*y(2)^0*y(3)^2*y(4)^0*-1.593895e-02+...
y(1)^1*y(2)^0*y(3)^2*y(4)^1*8.754721e-03+...
y(1)^1*y(2)^0*y(3)^3*y(4)^0*-3.092645e-02+...
y(1)^1*y(2)^1*y(3)^0*y(4)^0*-1.235412e-02+...
y(1)^1*y(2)^1*y(3)^0*y(4)^1*2.558354e-02+...
y(1)^1*y(2)^1*y(3)^0*y(4)^2*-2.964552e-02+...
y(1)^1*y(2)^1*y(3)^1*y(4)^0*2.313835e-02+...
y(1)^1*y(2)^1*y(3)^1*y(4)^1*2.806051e-02+...
y(1)^1*y(2)^1*y(3)^2*y(4)^0*-1.770844e-03+...
y(1)^1*y(2)^2*y(3)^0*y(4)^0*-1.002101e-02+...
y(1)^1*y(2)^2*y(3)^0*y(4)^1*2.171069e-02+...
y(1)^1*y(2)^2*y(3)^1*y(4)^0*3.090982e-02+...
y(1)^1*y(2)^3*y(3)^0*y(4)^0*2.604919e-02+...
y(1)^2*y(2)^0*y(3)^0*y(4)^0*1.427637e-02+...
y(1)^2*y(2)^0*y(3)^0*y(4)^1*6.188260e-02+...
y(1)^2*y(2)^0*y(3)^0*y(4)^2*-2.478055e-03+...
y(1)^2*y(2)^0*y(3)^1*y(4)^0*3.000835e-02+...
y(1)^2*y(2)^0*y(3)^1*y(4)^1*4.020574e-02+...
y(1)^2*y(2)^0*y(3)^2*y(4)^0*-1.715281e-04+...
y(1)^2*y(2)^1*y(3)^0*y(4)^0*-1.659904e-02+...
y(1)^2*y(2)^1*y(3)^0*y(4)^1*3.106589e-02+...
y(1)^2*y(2)^1*y(3)^1*y(4)^0*2.999115e-02+...
y(1)^2*y(2)^2*y(3)^0*y(4)^0*-1.463517e-02+...
y(1)^3*y(2)^0*y(3)^0*y(4)^0*-1.142872e-02+...
y(1)^3*y(2)^0*y(3)^0*y(4)^1*5.381773e-02+...
y(1)^3*y(2)^0*y(3)^1*y(4)^0*2.824173e-02+...
y(1)^3*y(2)^1*y(3)^0*y(4)^0*-3.552282e-02+...
y(1)^4*y(2)^0*y(3)^0*y(4)^0*-4.561864e-02)...
)*8/(--9.194929e-01*m1--9.609337e-01*m2)+-4;
au = min(ub);
end