Changes between Version 10 and Version 11 of misc/setup_for_ML/example_beamsize

Timestamp:

09/13/23 23:26:33 (10 months ago)

Author:

obina

Comment:

--

misc/setup_for_ML/example_beamsize

@@ -74 +74 @@
     beta_y = [p.beta_y for p in tws]
 
+    titlestr = "QM1.k1=%.2f, QM2.k1=%.2f, sx=%.3f, sy=%.3f"%(x0, x1, sx, sy)
     plt.plot(s, beta_x, 'b-');
     plt.plot(s, beta_y, 'r-');
@@ -79 +80 @@
     plt.xlabel('s [m]')
     plt.ylabel('beta [m]')
+    plt.title(titlestr)
     plt.grid(True)
     plt.show()
@@ -132 +134 @@
 
 == 次のステップ ==
-適当に色々と変更して試すと良い
- - ビームの初期パラメータを変える：簡単のためには円形ビーム（βx = βy ）にして、αx=αy=0 にすると、もっと単純
+色々と変更して試すと良い
+ - ビームの初期パラメータを変える：簡単のためには円形ビーム（βx = βy ）にして、αx=αy=0 にすると、もっと単純な形になる
  - 解析解との比較
  - 評価関数を色々変える：sx+sy, sx*sy, log(sx+sy+), -1/(sx+sy), ....ゼロ割に気をつけて適当なオフセットを入れると-1/((sx+10)*(sy+10)) とか....
     * それぞれで最小値を見つける速度が変わる
     * 単なるビームサイズだけではなく、他の条件（個々のQM k値を出来るだけ小さくするとか、途中のビームサイズを小さくするとか(betatron関数を小さくする））などを入れてみる
- - Optics的には Bendを入れてdispersionが入った状態にするとか、空間電荷効果を入れるとか、より実践的な形に
+    * 最小にするのではなく、設計値に設定する
+ - Optics的には Bendを入れてdispersionが入った状態にするとか、QMの数を増やすなど
+ - Ocelotで空間電荷効果を入れるとか、より実践的な形に→計算時間との兼ね合い
  - 測定したビームサイズに意図的に「ノイズ」を入れて、GPの応答がどうなるかを確認する（重要！）
  - GPのカーネルを変える
  - 他のGPツール との比較
  - 他の最適化手法（Nelder–Mead(downhill simplex)とか、最急降下法とか、遺伝的アルゴリズムとか）との比較、あるいは組み合わせを考える
+ - Grid Search との比較, Random Search との比較