1 | function ex14_1_6 |
---|
2 | |
---|
3 | yalmip('clear') |
---|
4 | |
---|
5 | sdpvar x1 x2 x3 x4 x5 x6 x7 x8 x9 x10 objvar; |
---|
6 | |
---|
7 | F = set([]); |
---|
8 | |
---|
9 | F = F + set( - x10 + objvar == 0); |
---|
10 | |
---|
11 | F = F + set( (1 - x1*x2)*x3*(exp(x5*(0.485 + (-0.0052095*x7) - 0.0285132*x8)) - 1) + 23.3037*x2 - x10 <= 28.5132); |
---|
12 | |
---|
13 | F = F + set( (1 - x1*x2)*x3*(exp(x5*(0.752 + (-0.0100677*x7) - 0.1118467*x8)) - 1) + 101.779*x2 - x10 <= 111.8467); |
---|
14 | |
---|
15 | F = F + set( (1 - x1*x2)*x3*(exp(x5*(0.869 + (-0.0229274*x7) - 0.1343884*x8)) - 1) + 111.461*x2 - x10 <= 134.3884); |
---|
16 | |
---|
17 | F = F + set( (1 - x1*x2)*x3*(exp(x5*(0.982 + (-0.0202153*x7) - 0.2114823*x8)) - 1) + 191.267*x2 - x10 <= 211.4823); |
---|
18 | |
---|
19 | F = F + set( (-(1 - x1*x2)*x3*(exp(x5*(0.485 + (-0.0052095*x7) - 0.0285132*x8)) - 1)) - 23.3037*x2 - x10 <= -28.5132); |
---|
20 | |
---|
21 | F = F + set( (-(1 - x1*x2)*x3*(exp(x5*(0.752 + (-0.0100677*x7) - 0.1118467*x8)) - 1)) - 101.779*x2 - x10 <= -111.8467); |
---|
22 | |
---|
23 | F = F + set( (-(1 - x1*x2)*x3*(exp(x5*(0.869 + (-0.0229274*x7) - 0.1343884*x8)) - 1)) - 111.461*x2 - x10 <= -134.3884); |
---|
24 | |
---|
25 | F = F + set( (-(1 - x1*x2)*x3*(exp(x5*(0.982 + (-0.0202153*x7) - 0.2114823*x8)) - 1)) - 191.267*x2 - x10 <= -211.4823); |
---|
26 | |
---|
27 | F = F + set( (1 - x1*x2)*x4*(exp(x6*(0.116 + 0.0233037*x9 - 0.0052095*x7)) - 1) - 28.5132*x1 - x10 <= -23.3037); |
---|
28 | |
---|
29 | F = F + set( (1 - x1*x2)*x4*(exp(x6*(0.101779*x9 - 0.0100677*x7 - 0.502)) - 1) - 111.8467*x1 - x10 <= -101.779); |
---|
30 | |
---|
31 | F = F + set( (1 - x1*x2)*x4*(exp(x6*(0.166 + 0.111461*x9 - 0.0229274*x7)) - 1) - 134.3884*x1 - x10 <= -111.461); |
---|
32 | |
---|
33 | F = F + set( (1 - x1*x2)*x4*(exp(x6*(0.191267*x9 - 0.0202153*x7 - 0.473)) - 1) - 211.4823*x1 - x10 <= -191.267); |
---|
34 | |
---|
35 | F = F + set( 28.5132*x1 - (1 - x1*x2)*x4*(exp(x6*(0.116 + 0.0233037*x9 - 0.0052095*x7) ) - 1) - x10 <= 23.3037); |
---|
36 | |
---|
37 | F = F + set( 111.8467*x1 - (1 - x1*x2)*x4*(exp(x6*(0.101779*x9 - 0.0100677*x7 - 0.502) ) - 1) - x10 <= 101.779); |
---|
38 | |
---|
39 | F = F + set( 134.3884*x1 - (1 - x1*x2)*x4*(exp(x6*(0.166 + 0.111461*x9 - 0.0229274*x7) ) - 1) - x10 <= 111.461); |
---|
40 | |
---|
41 | F = F + set( 211.4823*x1 - (1 - x1*x2)*x4*(exp(x6*(0.191267*x9 - 0.0202153*x7 - 0.473) ) - 1) - x10 <= 191.267); |
---|
42 | |
---|
43 | F = F + set( x1*x3 - x2*x4 == 0); |
---|
44 | |
---|
45 | |
---|
46 | F = F + set(0 < [x1 x2 x3 x4 x5 x6 x7 x8 x9] < 10); |
---|
47 | |
---|
48 | solvesdp(F,objvar,sdpsettings('solver','bmibnb','bmibnb.upper','none','allownon',1,'bmibnb.maxiter',10)) |
---|