授業
担当：小林、TA:M1永井
第1回(2016.11.4) 概要〜第１章 システムを状態方程式で記述する†
この授業の目的：与えられた物理系に対して、最適制御系を設計できるようになる
- Q: □5の一般の場合にでてくるIとは何ですか？
- A: 単位行列です。
- Q: □6の例の導出がよくわからなかった。
- A: この前まではわかっているとすると、伝達関数行列の式
\[ C(sI - A)^{-1} B \]
に、わかっている行列の中身を代入するのは良いと思います。
あとは、(sI-A) の逆行列の計算でしょうか。2×2 の逆行列は良くでてくるので、
覚えておくと良いと思います。
つまり、
\[ \left[ \begin{array}{cc} a & b \\ c & d \end{array} \right]^{-1} = \frac{\left[\begin{array}{cc} d & -b \\ -c & a \end{array}\right]}{ad - bc} \]
です（実際にもとの行列とかければ単位行列となることから逆行列であることが確認できます）。
- Q:
\[ \left\{\begin{array}{l} \dot x = A x + B u \\ y = C x \end{array} \right. \]
の A, B, C の求め方がいまいち理解できなかった。
- A: 与えられた伝達関数に対応する行列 A, B, C が予めわかっているので（伝達関数の係数 ai, bi で表されているので）、伝達関数の係数を読み取って、それを A, B, C に代入するだけです。
なぜ行列 A, B, C が予めそのようになるのか？という質問の場合は、その A, B, C を用いた状態空間表現に対応する伝達関数を求めれば確認できます。
- Q: よくわからなかった
- A: できれば、わからない箇所がでてきた時点で指摘してもらうか、または、わからなかった部分を教えてもらえれば回答します。
- Q:
\[ \dot{\tilde x} T A T^{-1} \mbox{ or } T^T \]
どちらでしたか。座った位置からはっきりしませんでした。
- A: 前者（Tの逆行列）です。字が少し小さかったかもしれません。次回もう少し大きくしますが、席は十分に空いているので、前の方に着席をお願いします。
第2回(2016.11.8)（11/11休講分の補講）第２章 システムの応答と安定性†
- 解析問題と設計問題: 解析が基本、今日は解析を扱う。特に安定性
- 古典制御における安定性：インパルス応答がt→∞で0に収束する
- ⇔ 現代制御における安定性：任意の初期値に対する応答がt→∞で0に収束する
- ⇔ 伝達関数の極の実部がすべて負
- ⇔ A行列の固有値の実部がすべて負
- 復習：たたみこみのラプラス変換、行列の固有値と固有ベクトル、逆行列
- 状態遷移行列の定義、状態遷移行列のラプラス変換
- Q:
\[ \left\{ \begin{array}{l} \dot x = Ax + Bu \\ y = C x \end{array} \right. \]
をラプラス変換・逆変換して初期値応答が0に収束するということはA行列の固有値の実部が負ということが要点ですか？すべて負であれば、制御系が安定ということですか？
- A: ラプラス変換は、初期値応答を求めるために使用しただけなので、安定性には直接関係ありません。初期値応答が0に収束するということとA行列の固有値の実部がすべて負ということが等価です。要点は、小さな黒板に書いた4つがすべて等価ということです。
答になっていない場合はまた聞いてください。
- Q: なぜラプラス変換＆逆ラプラス変換したのでしょうか？
- A: 微分方程式を解くために用いました。他の手法でも良いです。
- Q: いきなり出てきた a < 0 が分かりません。aとは何ですか？
- A: 行列Aがスカラの場合を A = a（実数）としました。a < 0 は安定条件です。
- Q: 矢印でいったりきたりするのが少し分かりにくかった。
- A: 黒板の番号が先に進んだにもかかわらず、書く手間の都合で、前の番号の箇所に戻った所が一箇所ありました。その部分のことだと思います。確かに、板書を写している人にはわかりにくいだろうな、と思いつつ行いました。以後、気を付けます。
- Q:
\[ {}^\forall x(0) \]
は何ですか？
- A: 全ての x(0) という意味です。
第3回(2016.11.18) 第3章 可制御性(3.3可制御性とその条件), 第5章 極配置法(5.1 フィードバック係数ベクトルを直接計算する)†
- 解析から設計へ、出力フィードバック（難)→状態フィードバック（簡単、基本）
- 閉ループ系のA行列 = A - BF
- A-BFの固有値をFによって任意に指定できる=可制御性
- 例1(Fによる固有値の指定不可)、例1'(指定可)
- 可制御性の定義
- 可制御性の判定方法(可制御性行列の正則性)
- レポート対策
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!以下は過去の情報です!!!!!!!!!!!!!!!!!!!!!!!!!!
- Q: \[ z = \left[ \begin{array}{cc} 0 & 1 \end{array} \right] \left[ \begin{array}{c} \dot z \\ z \end{array} \right] \] で [0 1] でない場合、コントローラを作る方法がないのか。また、
\[ z = \left[ \begin{array}{cc} 0 & 1 \end{array} \right] \left[ \begin{array}{c} \dot z \\ z \end{array} \right] \] とならなかった場合に、
\[ \dot z = \left[ \begin{array}{cc} 1 & 0 \end{array} \right] \left[ \begin{array}{c} \dot z \\ z \end{array} \right] \] ならば成り立つことが確認できれば、コントローラを作ることができると考えてよいか。
- A: y = C x の式は、状態 x に C 行列をかけたものを出力 y とする、という y の定義式です。C がどう与えられるかは問題によります。1 の部分が（性能の劣化などで）0.9 の場合も有り得る、などと説明しましたが、その場合は、0.9 z が出力されます。
- Q: 高専のとき、可制御かの判別に、可制御性行列がフルランクかどうかを使用していたような気がする
- A: フルランクで正しいです。一般には行列Ucは正方でなくフルランクかどうかをみます（正方の場合、正則とフルランクは等価です）。授業では話を簡単にするためにUcが正方の場合のみ話しましたが、一般にはフルランクかどうかで判別するが、という説明を一言すべきでした。すみません。
第4回(2015.11.20) 第6章 最適レギュレータ §6.1 評価関数と最適制御†
- 簡単な例題（a, x, b, u, f がすべてスカラ、b > 0）... 最適制御を直感的に理解できる
- (i) x(t)→0の収束する速さ、(ii) u(t)を小さく抑えること、にトレードオフがある
- 極配置法では、(i),(ii)のバランスを客観的に取りにくい
- 最適制御では、(i),(ii)のバランスを客観的に取れる
- 最適レギュレータ問題と、その解（リカッチ方程式、P>0）
- 最適制御の結果を使わずに J を最小化する f を求める ... f の二次方程式、閉ループが安定となる解を選ぶ（行列の場合：&ref(): File not found: "proof4.pdf" at page "授業/動的システムの解析と制御2016";）
- ＊1：f の二次方程式とリカッチ方程式の関係
- ＊2：閉ループ系の安定性と P > 0 の関係
- ＊3：Jの最小値を求める
- q を大きくする（(i) を重視する）と、|a-bf| が大きくなる
#ref(): File not found: "2015-11-20-1.jpg" at page "授業/動的システムの解析と制御2016"
#ref(): File not found: "2015-11-20-2.jpg" at page "授業/動的システムの解析と制御2016"
#ref(): File not found: "2015-11-20-3.jpg" at page "授業/動的システムの解析と制御2016"
- Q: 私の持っている教科書では、評価関数を「時間区間[0,T]にわたる生産率u(t)を制御して、在庫量x(t)を希望状態xd(t)に制御する場合、
\[ J(u) = \int_0^T [q(t)(x(t)-x_d(t))^2 + r(t) u^2(t)] dt \]
q(t), r(t) は正値関数（q(t)は非負も可）」としていたのですがこれは黒板の評価関数と同義なのでしょうか？それとも評価関数は複数あるのでしょうか？
- A: 積分区間が有限か無限か、重みが時不変か時変か、などの違いで評価関数は様々に定義されます。上記のように重み関数が時変の場合はリカッチ方程式の解も時間依存となります。この授業では取り扱いが最も簡単な場合＝積分区間は無限、重みは定数行列、リカッチの解も定数行列、を扱います。
第5回(2015.11.27) 第6章つづき〜 §6.2 重み行列と正定・半正定、§6.3最適制御系の安定性†
#ref(): File not found: "2015.11.27-1.jpg" at page "授業/動的システムの解析と制御2016"
#ref(): File not found: "2015.11.27-2.jpg" at page "授業/動的システムの解析と制御2016"
#ref(): File not found: "2015.11.27-3.jpg" at page "授業/動的システムの解析と制御2016"
#ref(): File not found: "2015.11.27-4.jpg" at page "授業/動的システムの解析と制御2016"
#ref(): File not found: "2015.11.27-5.jpg" at page "授業/動的システムの解析と制御2016"
- Q: M=0 となってしまった場合は、半正定、半負定のどちらにも当てはまるという考え方で正しいか？
- A: 正しいです。
- Q: テストは何日ですか？26日ですか？
- A: 12月25日（金）に実施します。26日は昨年の情報です。まぎらわしいので、以下の情報を今年度のものに更新しました。
第6回(2015.12.4) §9.1状態観測器の構造†
- 状態xが使えない場合
- (方法1) 状態の代わりに出力yを使う = 静的出力フィードバック ⇒ ダメ
- 別の方法：状態を推定して、それをxの代わりに使う
- 状態観測器の定義：t→∞で誤差 x(t) の推定誤差が0となる
- (方法2) 状態観測器? (演習9.1) + 状態フィードバック ⇒ ダメ
- (方法3) 状態観測器((9.3)式, 演習9.2) + 状態フィードバック
- 状態観測器を作る（（9.3)式の導出、方法２は状態観測器にならない）
- 可観測性（可制御性との関係）
- 演習9.3'：A - L C を安定（固有値の実部がすべて負）とする L の求め方
#ref(): File not found: "2015.12.04-1.jpg" at page "授業/動的システムの解析と制御2016"
#ref(): File not found: "2015.12.04-2.jpg" at page "授業/動的システムの解析と制御2016"
#ref(): File not found: "2015.12.04-3.jpg" at page "授業/動的システムの解析と制御2016"
#ref(): File not found: "2015.12.04-4.jpg" at page "授業/動的システムの解析と制御2016"
#ref(): File not found: "2015.12.04-5.jpg" at page "授業/動的システムの解析と制御2016"
- Q: 受け取り忘れた返却レポートを受け取れるか？
- A: 受け取りに来てください。
第7回(2015.12.11) §9.3併合系の固有値, 授業アンケート実施†
- 方法3で安定化できる理由：閉ループ系の固有値 = A - BF と A - LC の固有値
- 方法3が評価関数Jの最小値に与える影響？
- A-LC の固有値の実部が-∞に近づくように L を設定すると何が起こるか？
授業アンケート（本科目の前半・後半をまとめて）
#ref(): File not found: "2015.12.11-1.jpg" at page "授業/動的システムの解析と制御2016"
#ref(): File not found: "2015.12.11-2.jpg" at page "授業/動的システムの解析と制御2016"
#ref(): File not found: "2015.12.11-3.jpg" at page "授業/動的システムの解析と制御2016"
... 下記QAの誤りあり、すみません。
- Q: K(s) = Y(s) / U(s) どうしてですか？
- A: K(s) = U(s) / Y(s) の誤りです。ごめんなさい。
第8回(2015.12.25) 期末テスト, 授業アンケート実施†
- 試験内容: レポート課題の中から、数値を変えて出題する
- 持ち込み可能なもの: 筆記用具のみ
- 試験時間: 90分
- 授業アンケート（本科目の前半・後半をまとめて）
2015.12.28 前半・後半を含めた総合成績を掲示しています。レポート返却場所にて。確認してください。
![[PukiWiki]](image/logo_crop.png "[PukiWiki]"){kind=logo}
syllabus.pdf
板書1
板書2
板書3
板書4
板書1
板書2
板書3
板書4
板書5
板書6
板書1
板書2
板書3
板書4