授業/動的システムの解析と制御2019 の履歴(No.30)

第1回(2019.9.20) 前半（古典制御）概要～第5章周波数応答†

シラバス syllabus.pdf （変更あり）
成績の評価方法と評価項目
- 前半：レポート(36% = 6点×6回=36点)、中間テスト(64%)
- 後半：レポート(36% = 6点×6回=36点)、期末テスト(64%)
- 前半と後半の平均点を最終成績とする。
教科書（前半、古典制御）：「フィードバック制御入門」(杉江俊治、藤田政之著、コロナ社)
- 5. 周波数応答
  - 5.1 周波数応答と伝達関数
  - 5.2 ベクトル軌跡
  - 5.3 ボード線図
  - 5.4 ボード線図の性質
- 6. フィードバック制御系の安定性
- 8. フィードバック制御系の設計法
講義スライド &ref(): File not found: "slide01.pdf" at page "授業/動的システムの解析と制御2019";
演習問題# &ref(): File not found: "exercise01.pdf" at page "授業/動的システムの解析と制御2019"; 第2回の授業で解説する予定です。事前に問題を解いておいてください。
グラフ用紙1 graph01.pdf
グラフ用紙2 graph02.pdf

板書3

Q: わかりやすいですが、ちょっと速いです。最後のはちょっとわかりにくいです。
Q: 速すぎる
Q: スピードがちょっとはやい。
A: カリキュラム上、内容を減らすことができないため、できるだけ無駄を省いて話をしたいと思います。

Q: ゲインや位相などのやり方は分かったが、なぜそうなるのか理解しきれなかった感じがする。
A: ゲイン線図や位相線図は描けるが、なぜそのように描くのか、という質問だとすると、そのように描くと便利だからです。次回、その話（ボード線図の利点）からします。

Q: 伝達関数からの式の展開がやや追いきれなかった。
A: できれば、その場で質問して止めてもらえると助かります。

Q: レポートの \[ m \ddot x = -c\dot x - kx + af \] は \[ \frac{a f}{c\dot x - kx} \] になるのですか
A: なるのかならないのかと言われれば、（前者は等式で後者は式でないので）比較自体ができませんが、質問の意味が分からないため再度聞いてください。

レポート#1 report1_fixed.pdf

第2回(2019.9.27) 第5章周波数応答（つづき）†

講義スライド ... 第1回のつづき（5.4.2 ボード線図の利点）から
演習問題#1（第1回で公開済）の解説をします（解答例 &ref(): File not found: "exercise01A.pdf" at page "授業/動的システムの解析と制御2019"; &ref(): File not found: "exercise01slide.pdf" at page "授業/動的システムの解析と制御2019";）

板書3

Q: G(s)=Ts+1の時の数値が1以外の例えばTs+10の時はそのままボード線図は書けないのでしょうか
A: そのまま書くのは難しいと思うので、時定数 T/10 または折れ点角周波数 10/T を作ってから書いてください。つまり、 \[ Ts+10 = T(s + 10/T) = T (10/T) \frac{s+(10/T)}{10/T} = 10 \frac{s+(10/T)}{10/T} \] とすれば、折れ点角周波数は 10/T、ゲイン10 として書けます。

Q: やはり基礎知識が無い為非常に勤かしい（原文ママ）
A: どこまで簡単（理解できる）で、どこから難しい（理解できない）か、境界を確認しつつ進めると良いと思います。

レポート#2 report2.pdf

第3回(2019.10.4) 第6章フィードバック制御系の安定性†

講義スライド &ref(): File not found: "slide03.pdf" at page "授業/動的システムの解析と制御2019";
演習問題#2 &ref(): File not found: "exercise02.pdf" at page "授業/動的システムの解析と制御2019"; 第4回の授業で解説する予定です（変更の可能性あり）。事前に問題を解いておいてください。

Q: 矢印が＋が分かりにくかった
A: 申し訳ありません。pdfファイルの方を見てください。

Q: むずかしい
A: 予習、復習、教科書の練習問題など、ある程度時間をかける必要はあると思います。

Q: 内部安定性のスライド22の \[ P(s) = \frac{N_P(s)}{D_P(s)}, \quad K(s) = \frac{N_K(s)}{D_K(s)} \] が \[ \phi(s) := D_P(s)D_K(s) + N_P(s)N_K(s) \] になるどうしゅつかていがよくわからない
A: 黒板に r から y までの伝達関数の場合を \[ T_{yr} = \frac{PK}{1+PK} = \frac{\frac{N_P}{D_P}\frac{N_K}{D_K}}{1+\frac{N_P}{D_P}\frac{N_K}{D_K}} \] と書きましたが、ここでさらに分母分子に D_P D_K をかけると、 \[ T_{yr} = \frac{N_P N_K}{D_P D_K + N_P N_K} \] となって、分母多項式がΦ(s)であることがわかります。

Q: Π=0のとき、ωを0から∞と変化させたときにベクトル軌跡が"常に"左を見るようなら安定とあるのですが、安定である（註：スライドp.45の右下の図）は右に見る部分があるように見えるのですが、どういうことでしょうか？
A: 進行方向右手に(-1,0)の点が見えるようなωが存在するが、良いのか？という指摘と思います。その意味では、左上と右上の図でも、ωが小さいときにはどちらかと言えば(-1,0)は進行方向右手に存在するので、やはり「安定」と判別できなくなってしまいます。確かにそのように受け取れるわかりにくい表現だと思いますが、ここでは川の流れの場合の右岸か左岸か、という意味だと受け取ってください。そのような説明が不足していてすみません。

Q: 授業のプレゼンを前日までにwebにアップロードしてほしいです。
A: アップロードが遅くなり申し訳ありませんでした。次回から、水曜日のうちにアップロードするようにします。

レポート#3 report3.pdf

第4回(2019.10.11) 第6章フィードバック制御系の安定性（つづき）†

講義スライド ... 第3回のつづき（6.3 ゲイン余裕、位相余裕）から
演習問題#2（第3回で公開済）の解説をします（解答例 &ref(): File not found: "exercise02A.pdf" at page "授業/動的システムの解析と制御2019"; &ref(): File not found: "exercise02slide.pdf" at page "授業/動的システムの解析と制御2019";）

レポート#4 report4.pdf

第5回(2019.10.18) 第8章フィードバック制御系の設計法†

講義スライド &ref(): File not found: "slide05.pdf" at page "授業/動的システムの解析と制御2019"; §8.2 PID補償による制御系設計（p.24）まで
演習問題#3 &ref(): File not found: "exercise03.pdf" at page "授業/動的システムの解析と制御2019"; 次回以降の授業で解説する予定です。授業範囲内の問題は、事前に解いておいてください。

Q: PID制御は制御対象に対してP補償ならPをかければ良いのですか。
A: 良いです。PID制御に限らず、制御をすることは、制御対象に補償器（コントローラ）をかけたものを開ループとすること、なので。制御対象をゲインk倍することは、比例ゲインがkの比例制御と同じです。

レポート#5 report5.pdf

第6回(2019.10.25) 第8章フィードバック制御系の設計法（つづき）†

講義スライド ... 第5回のつづき（8.3 位相進み-遅れ補償による制御系設計）から
演習問題#3（第5回で公開済）の解説をします（解答例 &ref(): File not found: "exercise03A.pdf" at page "授業/動的システムの解析と制御2019"; &ref(): File not found: "exercise03slide.pdf" at page "授業/動的システムの解析と制御2019";）

Q: レポートを期限内に提出しているのに遅れ扱いにされるということがあった。点数の良い悪いにかぎらず、成績に関わることなので、このようなことが起きないよう対策してもらいたい。
A: 締切時刻17:00に不在だったため、TA学生に17:00に目印の紙を投函してもらいました。その上に乗っていたレポートを締切後に投函されたと判断しました。ここからは推測ですが、目印の紙を投函する際にレポートが下に落ちておらず、目印の紙が先に下に落ちた可能性が考えられます。しかし、客観的に判断する材料がなく、大変申し訳ありません。17時に不在にならないようにします。

Q: [例8.3]:どうしてPM>=40°を満たす時、K=1ですか？
A: K=1 とした結果、PM>=40°を満たすことがわかったので、K=1を採用した、ということです。

Q: テストはレポートと演習問題が解ければだいじょうぶですか？
A: 当初説明した通り、レポート課題の中から数値を変えて出題するので、レポートが解ければ大丈夫です。

レポート#6 report6.pdf

第7回(2019.10.25) 第8章フィードバック制御系の設計法（つづき）†

講義スライド ... 第6回のつづき 55枚目から
演習問題#3（第5回で公開済）の解説（つづき、問題2）

Q: グラフが複数重なり合っていて少し見ずらかったです。
A: 時間には余裕があったのでもう少し丁寧に説明すべきでした。pdfファイルで復習してください。

第8回(2019.11.1) 中間テスト†

試験内容: レポート課題の中から、数値を変えて出題する
持ち込み可能なもの: 筆記用具のみ（電卓の持ち込みはできません）
試験時間: 85分

2019.11.6 前半の成績を掲示しています。レポート返却場所にて。確認してください。採点結果に疑義がある場合は申し出ること。

第1回(2019.11.8) 後半（現代制御）概要～第１章　システムを状態方程式で記述する†

この授業の目的：与えられた物理系に対して、最適制御系を設計できるようになる

シラバス syllabus.pdf
成績の評価方法と評価項目：レポート(36% = 6点×6回=36点)、期末テスト(64%)
教科書：「演習で学ぶ現代制御理論」(森泰親著、森北出版)
古典制御と現代制御の違い（長所と短所）
伝達関数と状態空間表現の関係（簡単な運動方程式を例に）

Q: レポートの(2)では出力が \[ \dot z \] なので sZ(s) = GU(s) となる G を求めればよいのか。
A: よいです。ただし、U(s)が力fのラプラス変換ならば。

レポート#7 report7.pdf

第2回(2019.11.15)第２章システムの応答と安定性†

解析問題と設計問題: 解析が基本、今日は解析を扱う。特に安定性
古典制御における安定性：インパルス応答がt→∞で0に収束する
⇔ 現代制御における安定性：任意の初期値に対する応答がt→∞で0に収束する
⇔ 伝達関数の極の実部がすべて負
⇔ A行列の固有値の実部がすべて負
復習：たたみこみのラプラス変換、行列の固有値と固有ベクトル、逆行列
状態遷移行列の定義、状態遷移行列のラプラス変換

Q: レポート対策の \[ A = \left[ \begin{array}{ccc} 0 & 1 & 0 \\ 0 & 0 & 1 \\ -8 & 0 & 0 \end{array} \right] \] 固有値の一つは -2 の問題は \[|\lambda I - A | = \left[ \begin{array}{ccc} \lambda & 0 & 0 \\ 0 & \lambda & 0 \\ 0 & 0 & \lambda \end{array} \right] - \left[ \begin{array}{ccc} 0 & 1 & 0 \\ 0 & 0 & 1 \\ -8 & 0 & 0 \end{array} \right] = \left[ \begin{array}{ccc} \lambda & -1 & 0 \\ 0 & \lambda & -1 \\ 8 & 0 & \lambda \end{array} \right] \] になるのであっていますか。

A: 行列式と行列をしっかり区別してください。つまり最後は3×3の行列ではなく行列式なので \[\left| \begin{array}{ccc} \lambda & -1 & 0 \\ 0 & \lambda & -1 \\ 8 & 0 & \lambda \end{array} \right| \] と書く必要があります（途中の行列も| |で挟む必要がある）。それ以外は合っています。

レポート#8 report8.pdf

第3回(2019.11.22) 第3章可制御性(3.3可制御性とその条件), 第5章極配置法(5.1 フィードバック係数ベクトルを直接計算する)†

解析から設計へ、出力フィードバック（難)→状態フィードバック（簡単、基本）
閉ループ系のA行列 = A - BF
A-BFの固有値をFによって任意に指定できる=可制御性
例1(Fによる固有値の指定不可)、例1'(指定可)
可制御性の定義
可制御性の判定方法(可制御性行列の正則性)
レポート対策

... 最後の[2 4 9] は [2 9 4] の誤記です。単純ミスで申し訳ありません。

Q: 可制御性行列の時に例えば□9の(例1')の \[ U_c = \left| \begin{array}{cc} B & AB \end{array} \right| \] の AB は \[ \left[ \begin{array}{cc} 1 & 0 \\ 0 & -1 \end{array} \right] \] になって 0 を消して \[ \left[ \begin{array}{c} 1 \\ -1 \end{array} \right] \] になるということでしょうか？
A: 「0を消す」の意味がわかりませんが、AB は通常の意味（A と B の積）です。質問の答になっていない場合はまた聞いてください。

レポート#9 report9.pdf

第4回(2019.11.29) 第6章最適レギュレータ §6.1 評価関数と最適制御†

簡単な例題（a, x, b, u, f がすべてスカラ）... 最適制御を直感的に理解できる
(i) x(t)→0の収束する速さ、(ii) u(t)を小さく抑えること、にトレードオフがある
極配置法では、(i),(ii)のバランスを客観的に取りにくい
最適制御では、(i),(ii)のバランスを客観的に取れる
最適レギュレータ問題と、その解（リカッチ方程式、P>0）
最適制御の結果を使わずに J を最小化する f を求める ... f の二次方程式、閉ループが安定となる解を選ぶ
＊1：f の二次方程式とリカッチ方程式の関係
＊2：閉ループ系の安定性と P > 0 の関係
＊3：Jの最小値を求める
r を大きくする（(ii) を重視する）と、|f| と |a-bf| が小さくなる

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!以下は過去の情報です!!!!!!!!!!!!!!!!!!!!!!!!!!

#ref(): File not found: "2018.12.07_1.jpg" at page "授業/動的システムの解析と制御2019"

#ref(): File not found: "2018.12.07_2.jpg" at page "授業/動的システムの解析と制御2019"

#ref(): File not found: "2018.12.07_3.jpg" at page "授業/動的システムの解析と制御2019"

#ref(): File not found: "2018.12.07_4.jpg" at page "授業/動的システムの解析と制御2019"

Q: リカッチ代数方程式のスカラーに対応する式が分かりにくかった。
A: f=bp/r の下で、どちらも同じ二次方程式（一方はfについて、もう一方はpについて）である、ということが分かってもらえれば十分です。

レポート#10 &ref(): File not found: "report10_fixed.pdf" at page "授業/動的システムの解析と制御2019";

第5回(2018.12.14) 第6章つづき～ §6.2 重み行列と正定・半正定、§6.3最適制御系の安定性†

行列へ拡張: q→Q≧0, r→R>0
(半)正定行列の定義
対称行列の固有値と正定性の関係(演習6.6), 対称行列の固有値はすべて実数(演習6.8) proof3.pdf
最適制御⇒安定かつJが最小証明
最適レギュレータの設計(演習6.3,6.4)：リカッチ方程式の解P、P > 0 ⇒ 対応する F が閉ループ系を安定化

#ref(): File not found: "2018.12.14_1.jpg" at page "授業/動的システムの解析と制御2019"

#ref(): File not found: "2018.12.14_2.jpg" at page "授業/動的システムの解析と制御2019"

#ref(): File not found: "2018.12.14_3.jpg" at page "授業/動的システムの解析と制御2019"

#ref(): File not found: "2018.12.14_4.jpg" at page "授業/動的システムの解析と制御2019"

#ref(): File not found: "2018.12.14_5.jpg" at page "授業/動的システムの解析と制御2019"

Q: 朝起きられません
A: 夜早めに寝てみてください。

Q: レポート対策で、Pは正定だから対角成分が正と仮定しているのに、次の問でPが正定であることを示すのに違和感があった。
A: はっきり説明すべきでしたが、リカッチ方程式の（正定とは限らない）解は複数存在します（今日説明した例でも、p2 = ±1、p3 = 1, -3 などの可能性があった通りです）。よって、リカッチ方程式の解が一つ求まったとしても、それが正定であるとは限りません。p3 = -3 と選んでしまうとその時点で正定でないことが確定してしまうので、p3 = 1 と選ぶ必要がありますが、これはPが正定であるために十分ではありません。そのため、(2)で正定性を調べています。とても良い質問で、もし授業中に質問してくれていたら、不十分な説明の良い補足になったと思います。

レポート#11 report11.pdf

第6回(2018.12.18) §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: "2018.12.18_1.jpg" at page "授業/動的システムの解析と制御2019"

#ref(): File not found: "2018.12.18_2.jpg" at page "授業/動的システムの解析と制御2019"

#ref(): File not found: "2018.12.18_3.jpg" at page "授業/動的システムの解析と制御2019"

#ref(): File not found: "2018.12.18_4.jpg" at page "授業/動的システムの解析と制御2019"

#ref(): File not found: "2018.12.18_5.jpg" at page "授業/動的システムの解析と制御2019"

Q: □7（方法2）... L=0とした（方法3）
聞き逃しかもしれませんが、（方法2）..._____（方法3）
と方法が2通り混在して書かれている意味が見返すとわからなくなってしまいました。
A: （方法2）は（方法3）の特別な場合、つまり L=0 とした場合、y を使用しない場合です。答えになっていない場合はまた聞いてください。

Q: レポート#12はいつ返却予定か？
A: 21日（木）10:00に返却予定です。

レポート#12 report12.pdf

第7回(2018.12.21) §9.3併合系の固有値†

方法3で安定化できる理由：閉ループ系の固有値 = A - BF と A - LC の固有値
方法3が評価関数Jの最小値に与える影響？
A-LC の固有値の実部が-∞に近づくように L を設定するとどうなるか？

#ref(): File not found: "2018.12.21_1.jpg" at page "授業/動的システムの解析と制御2019"

#ref(): File not found: "2018.12.21_2.jpg" at page "授業/動的システムの解析と制御2019"

#ref(): File not found: "2018.12.21_3.jpg" at page "授業/動的システムの解析と制御2019"

#ref(): File not found: "2018.12.21_4.jpg" at page "授業/動的システムの解析と制御2019"

Q: 制御を学ぶ楽しさが少し分かった。
A: よくがんばりました。これを聞いて私も少し報われた気がします。