level 2
此处放上Maximize函数
findopt[m0_, b0_, g0_, pb1_, pb2_, pb3_, pb4_, pb5_, pg1_, pg2_, pg3_,
pg4_, pg5_] :=
Maximize[{m0 - pb1 x1 - pg1 y1 - 0.02*pb1 Abs[x1] -
0.01*pg1 Abs[y1] - pb2 x2 - pg2 y2 - 0.02*pb2 Abs[x2] -
0.01*pg2 Abs[y2] - pb3 x3 - pg3 y3 - 0.02*pb3 Abs[x3] -
0.01*pg3 Abs[y3] - pb4 x4 - pg4 y4 - 0.02*pb4 Abs[x4] -
0.01*pg4 Abs[y4] - pb5 x5 - pg5 y5 - 0.02*pb5 Abs[x5] -
0.01*pg5 Abs[y5] + pb5 (x1 + x2 + x3 + x4 + x5 + b0) +
pg5 (y1 + y2 + y3 + y4 + y5 + g0),
If[t[[i + 9]] == 0, y1 == 0, Element[y1, Reals]],
If[t[[i + 10]] == 0, y2 == 0, Element[y2, Reals]],
If[t[[i + 11]] == 0, y3 == 0, Element[y3, Reals]],
If[t[[i + 12]] == 0, y4 == 1, Element[y4, Reals]],
If[t[[i + 13]] == 0, y5 == 0, Element[y5, Reals]],
m0 - pb1 x1 - pg1 y1 - 0.02*Abs[x1] - 0.01*Abs[y1] >= 0,
x1 + b0 >= 0, y1 + g0 >= 0,
m0 - pb1 x1 - pg1 y1 - 0.02*pb1 Abs[x1] - 0.01*pg1 Abs[y1] -
pb2 x2 - pg2 y2 - 0.02*pb2 Abs[x2] - 0.01*pg2 Abs[y2] >= 0,
x1 + x2 + b0 >= 0, y1 + y2 + g0 >= 0,
m0 - pb1 x1 - pg1 y1 - 0.02*pb1 Abs[x1] - 0.01*pg1 Abs[y1] -
pb2 x2 - pg2 y2 - 0.02*pb2 Abs[x2] - 0.01*pg2 Abs[y2] - pb3 x3 -
pg3 y3 - 0.02*pb3 Abs[x3] - 0.01*pg3 Abs[y3] >= 0,
x1 + x2 + x3 + b0 >= 0, y1 + y2 + y3 + g0 >= 0,
m0 - pb1 x1 - pg1 y1 - 0.02*pb1 Abs[x1] - 0.01*pg1 Abs[y1] -
pb2 x2 - pg2 y2 - 0.02*pb2 Abs[x2] - 0.01*pg2 Abs[y2] - pb3 x3 -
pg3 y3 - 0.02*pb3 Abs[x3] - 0.01*pg3 Abs[y3] - pb4 x4 - pg4 y4 -
0.02*pb4 Abs[x4] - 0.01*pg4 Abs[y4] >= 0,
x1 + x2 + x3 + x4 + b0 >= 0, y1 + y2 + y3 + y4 + g0 >= 0,
m0 - pb1 x1 - pg1 y1 - 0.02*pb1 Abs[x1] - 0.01*pg1 Abs[y1] -
pb2 x2 - pg2 y2 - 0.02*pb2 Abs[x2] - 0.01*pg2 Abs[y2] - pb3 x3 -
pg3 y3 - 0.02*pb3 Abs[x3] - 0.01*pg3 Abs[y3] - pb4 x4 - pg4 y4 -
0.02*pb4 Abs[x4] - 0.01*pg4 Abs[y4] - pb5 x5 - pg5 y5 -
0.02*pb5 Abs[x5] - 0.01*pg5 Abs[y5] >= 0,
x1 + x2 + x3 + x4 + x5 + b0 >= 0,
y1 + y2 + y3 + y4 + y5 + g0 >= 0}, {x1, x2, x3, x4, x5, y1, y2, y3,
y4, y5}, WorkingPrecision -> 100]
2022年02月19日 08点02分
