- 最適化(極小値)探索して見るとか。
問
題の出典はこ
このを使う。まず、制約条件尽きの例
x^2 + ( y - 0.5)^2, y>=0 && y >= x + 1
CCTLの計算結果 |
Mathematicaの計算結果 |
arg[0]=-0.259332
arg[1]=0.741514
minval=0.1255820560818098
|
-0.25
0.75
0.125 |
次の例は凹凸がたくさんあってこの種の問題においては意地悪なケースである。
CCTLの
計算結果 |
Mathematicaの
計算結果 |
arg[0]=-0.024279
arg[1]=0.210291
minval=-3.3063125936525632
確かにこの点は最も深い箇所である!!
|
-0.343687
0.44071
-3.00383<==失敗例
-0.0244031
-0.210612
-3.30687
|
ShekelのFoxholesは全部の深さが異なり,一般に最適化するのが難しい例である.
1/(0.002 +
1/(25 + pow(-32 + x,6) + pow(-32 + y,6)) + 1/(24 + pow(-16
+ x,6) + pow(-32 + y,6)) +
1/(23 + pow(x,6) + pow(-32 + y,6))
+ 1/(22 + pow(16 + x,6) + pow(-32 + y,6)) +
1/(21 + pow(32 + x,6) + pow(-32 + y,6)) +
1/(20 + pow(-32 + x,6) + pow(-16 + y,6)) +
1/(19 + pow(-16 + x,6) + pow(-16 + y,6)) +
1/(18 + pow(x,6) + pow(-16 + y,6)) +
1/(17 + pow(16 + x,6) + pow(-16 + y,6)) +
1/(16 + pow(32 + x,6) + pow(-16 + y,6)) +
1/(15 + pow(-32 + x,6) + pow(y,6)) + 1/(14 +
pow(-16 + x,6) + pow(y,6)) +
1/(13 + pow(x,6) + pow(y,6)) + 1/(12 + pow(16
+ x,6) + pow(y,6)) +
1/(11 + pow(32 + x,6) + pow(y,6)) + 1/(10 +
pow(-32 + x,6) + pow(16 + y,6)) +
1/(9 + pow(-16 + x,6) + pow(16 + y,6)) + 1/(8
+ pow(x,6) + pow(16 + y,6)) +
1/(7 + pow(16 + x,6) + pow(16 + y,6)) + 1/(6 +
pow(32 + x,6) + pow(16 + y,6)) +
1/(5 + pow(-32 + x,6) + pow(32 + y,6)) + 1/(4
+ pow(-16 + x,6) + pow(32 + y,6)) +
1/(3 + pow(x,6) + pow(32 + y,6)) + 1/(2 +
pow(16 + x,6) + pow(32 + y,6)) + 1/(1 + pow(32 + x,6) + pow(32 + y,6)))
|
CCTLの計算結果
|
Mathematicaの計算結果 |
arg[0]=-31.958227
arg[1]=-32.012368
minval=0.9980038432748928
確かにこの点は最も深い箇所である!!
|
21.9884 <==
失敗例
11.7187 <== 失敗例
17.3744 <== 失敗例
13.6186 <== 失敗例
17.3744 <== 失敗例
0.998004
(-31.9783,-31.9783)
|
画像はクリックで大きく表示できます。
※これら計算結果はWCCTLによって得られたものです
|