[37] | 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)) |
---|