技能一覧 › 大学院数学 変分法・最適制御

制御に制約のある最適制御

制御が のような範囲に制限されるとき、各時刻でハミルトニアンを について区間上で最小化し(2 次なら区間への切り詰め、1 次ならバンバン制御)、切り替え時刻を求めて評価関数の最小値を求められる

前提

この技能までの道筋(72 技能)
算数中学高校大学大学院

読み込み中…

右端(太枠)がこの技能です。左へドラッグすると、さらに前の前提をさかのぼれます。技能を押すと、そこまでの道筋が光ります。

一覧で見る(前提が先)
  1. 分数のたし算・ひき算 (算数 計算と数のきまり)
  2. 小数の四則計算 (算数 計算と数のきまり)
  3. 正負の数の加法・減法 (中学数学 数と式)
  4. 分数のかけ算・わり算 (算数 計算と数のきまり)
  5. 四則の混じった計算の順序 (算数 計算と数のきまり)
  6. 正負の数の乗除・累乗と四則混合 (中学数学 数と式)
  7. 分配法則を使った計算のくふう (算数 計算と数のきまり)
  8. 文字式の計算(同類項・分配法則) (中学数学 数と式)
  9. 乗法公式による展開 (中学数学 数と式)
  10. 因数分解 (中学数学 数と式)
  11. 3 次式までの展開と因数分解 (数学I 数と式)
  12. 微分係数と導関数の定義 (数学II 微分法と積分法)
  13. 多項式の微分 (数学II 微分法と積分法)
  14. 等間隔にならんだ数の和 (算数 計算と数のきまり)
  15. 等差・等比数列の和と Σ (数学B 数列)
  16. 数列の極限 (数学III 極限)
  17. 関数の極限 (数学III 極限)
  18. 積と商の微分 (数学III 微分法)
  19. 合成関数の微分 (数学III 微分法)
  20. 偏微分と全微分 (大学数学 微分積分)
  21. 多項式の不定積分 (数学II 微分法と積分法)
  22. 円の周の長さと面積 (算数 平面図形)
  23. おうぎ形の弧の長さと面積 (算数 平面図形)
  24. 弧の長さからおうぎ形の中心角と面積を求める (中学数学 平面図形)
  25. 三角形・四角形の面積 (算数 平面図形)
  26. 約数の見つけ方 (算数 整数の性質)
  27. 素因数分解と約数の個数 (算数 整数の性質)
  28. 平方根の意味と根号の簡単化 (中学数学 数と式)
  29. 三平方の定理 (中学数学 相似と三平方の定理)
  30. 特別な直角三角形の辺の比 (中学数学 相似と三平方の定理)
  31. 弧度法と単位円による三角関数 (数学II 三角関数)
  32. の極限 (数学III 極限)
  33. 加法定理と 2 倍角・半角・積和の公式 (数学II 三角関数)
  34. 三角関数の導関数 (数学III 微分法)
  35. 根号の乗除と加減 (中学数学 数と式)
  36. 指数の拡張と指数関数 (数学II 指数関数と対数関数)
  37. 対数の定義と性質 (数学II 指数関数と対数関数)
  38. 自然対数の底 (数学III 極限)
  39. 指数関数・対数関数の導関数 (数学III 微分法)
  40. 基本的な関数の不定積分 (数学III 積分法)
  41. 置換積分 (数学III 積分法)
  42. 剰余の定理・因数定理と高次方程式 (数学II 複素数と方程式)
  43. 分数関数の積分(部分分数分解) (数学III 積分法)
  44. 変数分離形 (大学数学 微分方程式)
  45. 小数と分数のまじった計算 (算数 計算と数のきまり)
  46. 逆算(□ を求める) (算数 計算と数のきまり)
  47. 1 次方程式の解き方 (中学数学 方程式)
  48. 比例・反比例の式とグラフ (中学数学 関数)
  49. 関数 の式と変域 (中学数学 関数)
  50. 平方完成と 2 次関数のグラフ (数学I 2次関数)
  51. 2 次不等式 (数学I 2次関数)
  52. 増減表と極値 (数学II 微分法と積分法)
  53. 関数の連続性と中間値の定理 (数学III 極限)
  54. 平均値の定理 (数学III 微分法)
  55. 第 2 次導関数とグラフの凹凸 (数学III 微分法)
  56. 増減・凹凸・漸近線とグラフ (数学III 微分法)
  57. テイラー展開 (大学数学 微分積分)
  58. 無限等比級数と無限級数 (数学III 極限)
  59. べき級数と収束半径 (大学数学 微分積分)
  60. 複素数と 2 次方程式の虚数解 (数学II 複素数と方程式)
  61. オイラーの公式 (大学数学 微分積分)
  62. 三角関数の合成 (数学II 三角関数)
  63. 定数係数 2 階線形(斉次) (大学数学 微分方程式)
  64. 部分積分 (数学III 積分法)
  65. 1 階線形微分方程式 (大学数学 微分方程式)
  66. 定数係数 2 階線形(非斉次) (大学数学 微分方程式)
  67. オイラー・ラグランジュ方程式 (大学院数学 変分法・最適制御)
  68. オイラー型(コーシー・オイラー型)方程式 (大学数学 微分方程式)
  69. 自由端と自然境界条件 (大学院数学 変分法・最適制御)
  70. 非斉次の連立線形微分方程式 (大学数学 微分方程式)
  71. ポントリャーギンの最大原理 (大学院数学 変分法・最適制御)
  72. 定義域つきの2次関数の最大・最小 (数学I 2次関数)

説明

実際の制御問題では、操作量 に「推力は 以上 以下」「バルブの開度は から まで」のような上下限があるのが普通です。このとき の解が許される範囲からはみ出すことがあり、そのままでは使えません。この技能では、各時刻にハミルトニアン を の関数とみて区間 上で最小化することで最適制御を決めます。 が の 2 次式なら区間への切り詰め(クリップ)、1 次式なら端点だけを使うバンバン制御になり、制御の形が切り替わる時刻を求めて評価関数の最小値を計算します。

なぜそうしてよいのか

評価関数 を最小にする問題で、 とおきます。前提のポントリャーギンの最大原理で使った条件 は、本来は「最適な は各時刻で を について最小にする」という条件の、内点での形にすぎません。直観的には、ある短い時間 だけ制御を別の値 に取り替えると、評価関数はおよそ だけ変化します(随伴変数 は、状態の小さなずれが評価関数に与える影響を表す量だからです)。最適ならどの に取り替えても は減らないので、 がすべての許される で成り立つ、つまり です。最小点が区間の内部にあれば微分が 0 となり に戻りますが、端にあるときは微分は 0 とは限りません。随伴方程式 と横断条件 はそのまま使えます。

あとは「定義域つきの 2 次関数の最小」と同じ考え方です。( は によらない項)で なら、下に凸な放物線の頂点 が区間内ならそこ、区間より左なら 、右なら が最小点です。すなわち (頂点を区間に切り詰めた値)になります。

( は によらない項)のように について 1 次なら、最小は必ず端点でとられ、 なら 、 なら です(バンバン制御)。 を切り替え関数といい、 の符号が変わる時刻で制御が切り替わります。 が区間上で続く場合(特異制御)は、この条件だけでは が決まらないので別の考察が必要です。

手順

  1. ハミルトニアン を作り、 について何次式か(2 次か 1 次か)を確認する

  2. 随伴方程式 と横断条件 から を求める( が を含まなければ は単独で求まる)

  3. 各時刻で を 上で最小化する。2 次( の係数が正)なら頂点を区間に切り詰めて ( は の解)、1 次なら切り替え関数 の符号で または

  4. , (2 次の場合)または (1 次の場合)を解いて切り替え時刻を求め、 を区間に分けて を区分的に書く

  5. 状態方程式を初期条件から区間ごとに解き、切り替え時刻で が連続につながるようにする

  6. 評価関数を区間ごとに積分して足し合わせ、最小値を求める

例題

状態方程式 , で、制御は に制限されているとする。評価関数 の最小値を求めよ。

例題の解答

手順 1: , より について 2 次で、 の係数 は正である。

手順 2: なので 。終端コストはないので で、横断条件は 。よって 。

手順 3: の解は 。これを区間 に切り詰めて

手順 4: で だから下限には当たらない。 となるのは 。よって切り替え時刻は で

手順 5: では , より 、特に 。 では とおくと , ()より 。

手順 6:
-
-

よって の最小値は 。

(確認: と書き直せるので となり、被積分関数を各 で 上最小化すれば上と同じ が得られる。この問題では最小原理の条件を満たす制御が実際に最小を与えている。)

答え:

よくある誤り

この説明は AI が書き、別の AI が例題を独立に解いて答えを確かめ、説明を一歩ずつ査読しました。

入口の問題

この技能を使う問題