[[授業]]
//[[最新回へ>#zce48c3f]]
担当：小林、TA:M1岡村、竹村
[[スケジュール2020]]（変更の可能性あり。特に試験の実施について）
*第1回(2020.9.4) &color(red){前半（古典制御）};概要〜第5章 周波数応答 [#l52e32ae]
- シラバス &ref(syllabus.pdf); （変更あり）
- 成績の評価方法と評価項目（変更の可能性あり。試験が実施できなくなった場合、レポートのみで最終成績をつけます。）
-- 前半：レポート(36% = 6点×6回=36点)、中間テスト(64%)
-- 後半：レポート(36% = 6点×6回=36点)、期末テスト(64%)
-- 前半と後半の平均点を最終成績とする。
- 教科書&color(red){（前半、古典制御）};：[[「フィードバック制御入門」(杉江 俊治、藤田 政之 著、コロナ社)>http://www.coronasha.co.jp/np/isbn/9784339033038/]]
-- 5. 周波数応答
--- 5.1 周波数応答と伝達関数
--- 5.2 ベクトル軌跡
--- 5.3 ボード線図
--- 5.4 ボード線図の性質
-- 6. フィードバック制御系の安定性
-- 8. フィードバック制御系の設計法
- 講義スライド &ref(slide01.pdf);
- 演習問題# &ref(exercise01.pdf); 第2回の授業で解説する予定です。事前に問題を解いておいてください。
- グラフ用紙1 &ref(graph01.pdf);
- グラフ用紙2 &ref(graph02.pdf);
レポート#1 &ref(report1.pdf);
#ref(2020.09.04-1.jpg,left,noimg,板書1)
-Q: この講義は、計測制御工学の内容を理解していることを前提としている科目ですか？
-A: はい、前提としています。シラバスに「計測制御工学 (2 年 2 学期) あるいは同等の内容を履修済みであること.」と記載されている通りです。ただし、該当する内容についても、質問があれば回答します。
*第2回(2020.9.11) 第5章 周波数応答（つづき） [#ybd74b9e]
- 講義スライド ... 第1回のつづき（5.4.2 ボード線図の利点）から
- 演習問題#1（第1回で公開済）の解説をします（&color(red){解答例}; &ref(exercise01A.pdf); &ref(exercise01slide.pdf);）
レポート#2 &ref(report2.pdf);
#ref(2020.09.11-1.jpg,left,noimg,板書1)
-Q: 時定数はゲイン曲線、位相曲線を書くにあたってどのように配慮すればよいですか。(時定数があるとどう変化しますか。)
-A: 時定数によって折れ点角周波数が変化します。答えになっていない場合は、また聞いてください。
-Q: 演習課題などが載っているURLを教えてください&br;
よろしくお願いします
-A: 授業のホームページのことでしょうか。ならば、このページなので、この回答自体が見れないですね。。。
*第3回(2020.9.18) 第6章 フィードバック制御系の安定性 [#oa2ae761]
- 講義スライド &ref(slide03.pdf);
- 演習問題#2 &ref(exercise02.pdf); 第4回の授業で解説する予定です（変更の可能性あり）。事前に問題を解いておいてください。
レポート#3 &ref(report3.pdf);
#ref(2020.09.18-1.jpg,left,noimg,板書1)
-Q: 特にありません
-Q: 毎回のレポートの答えを教えてもらいたいです。復習やあってたかどうかの確認をしたいです。
-A: もっともな指摘です。例年は採点後のレポートをすぐ返却しているのですが、今年度はそれができず、中間試験も実施可能か方針が決まっておらず、現在検討中です。申し訳ありません。
*第4回(2020.9.25) 第6章 フィードバック制御系の安定性（つづき） [#a35e0c54]
- 講義スライド ... 第3回のつづき（6.3 ゲイン余裕、位相余裕）から
- 演習問題#2（第3回で公開済）の解説をします
（&color(red){解答例}; &ref(exercise02A.pdf); &ref(exercise02slide.pdf);）
レポート#4 &ref(report4.pdf);
#ref(2020.09.25-1.jpg,left,noimg,板書1)
*第5回(2020.10.2) 第8章 フィードバック制御系の設計法 [#t468364a]
- 講義スライド &ref(slide05.pdf); §8.2 PID補償による制御系設計（p.24）まで
- 演習問題#3 &ref(exercise03.pdf); 次回以降の授業で解説する予定です。授業範囲内の問題は、事前に解いておいてください。
レポート#5 &ref(report5.pdf);
#ref(2020.10.02-1.jpg,left,noimg,板書1)
#ref(2020.10.02-2.jpg,left,noimg,板書2)
*第6回(2020.10.9) 第8章 フィードバック制御系の設計法（つづき） [#f676c50c]
- 講義スライド ... 第5回のつづき（8.3 位相進み-遅れ補償による制御系設計）から
- 演習問題#3（第5回で公開済）の解説をします（&color(red){解答例}; &ref(exercise03A.pdf); &ref(exercise03slide.pdf);）
//■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■
&color(black,red){&size(20){!!!!!!!!!!!!!!!!!!!!!!!!!!!!!以下は過去の情報です!!!!!!!!!!!!!!!!!!!!!!!!!!};};
#ref(20191025_1.jpg,left,noimg,板書1)
#ref(20191025_2.jpg,left,noimg,板書2)
#ref(20191025_3.jpg,left,noimg,板書3)
#ref(20191025_4.jpg,left,noimg,板書4)
-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 &ref(report6.pdf);
*第7回(2019.10.25) 第8章 フィードバック制御系の設計法（つづき） [#w47bf472]
- 講義スライド ... 第6回のつづき 55枚目から
- 演習問題#3（第5回で公開済）の解説（つづき、問題2）
#ref(20191025_5.jpg,left,noimg,板書1)
#ref(20191025_6.jpg,left,noimg,板書2)
-Q: グラフが複数重なり合っていて少し見ずらかったです。
-A: 時間には余裕があったのでもう少し丁寧に説明すべきでした。pdfファイルで復習してください。
*第8回(2019.11.1) 中間テスト [#s3fbea9d]
- 試験内容: レポート課題の中から、数値を変えて出題する
- 持ち込み可能なもの: 筆記用具のみ（電卓の持ち込みはできません）
- 試験時間: 85分
&color(red){&size(25){2019.11.6 前半の成績を掲示しています。レポート返却場所にて。確認してください。採点結果に疑義がある場合は申し出ること。};};
*第1回(2019.11.8) &color(red){後半（現代制御）};概要〜第１章 システムを状態方程式で記述する [#z139db6e]
&color(blue,orange){この授業の目的：与えられた物理系に対して、最適制御系を設計できるようになる};
- シラバス &ref(syllabus.pdf);
- 成績の評価方法と評価項目：レポート(36% = 6点×6回=36点)、期末テスト(64%)
-教科書：[[「演習で学ぶ現代制御理論」(森 泰親 著、森北出版)>https://www.morikita.co.jp/books/book/2368/]]
- 古典制御と現代制御の違い（長所と短所）
- 伝達関数と状態空間表現の関係（簡単な運動方程式を例に）
#ref(20191108_1.jpg,left,noimg,板書1)
#ref(20191108_2.jpg,left,noimg,板書2)
#ref(20191108_3.jpg,left,noimg,板書3)
#ref(20191108_4.jpg,left,noimg,板書4)
-Q: レポートの(2)では出力が
\[
\dot z
\] なので sZ(s) = GU(s) となる G を求めればよいのか。
-A: よいです。ただし、U(s)が力fのラプラス変換ならば。
レポート#7 &ref(report7.pdf);
*第2回(2019.11.15)第２章 システムの応答と安定性 [#xf82b2a3]
- 解析問題と設計問題: 解析が基本、今日は解析を扱う。特に安定性
- 古典制御における安定性：インパルス応答がt→∞で0に収束する
- ⇔ 現代制御における安定性：任意の初期値に対する応答がt→∞で0に収束する
- ⇔ 伝達関数の極の実部がすべて負
- ⇔ A行列の固有値の実部がすべて負
- 復習：たたみこみのラプラス変換、行列の固有値と固有ベクトル、逆行列
- 状態遷移行列の定義、状態遷移行列のラプラス変換
#ref(20191115_1.jpg,left,noimg,板書1)
#ref(20191115_2.jpg,left,noimg,板書2)
#ref(20191115_3.jpg,left,noimg,板書3)
#ref(20191115_4.jpg,left,noimg,板書4)
#ref(20191115_5.jpg,left,noimg,板書5)
-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 &ref(report8.pdf);
*第3回(2019.11.22) 第3章 可制御性(3.3可制御性とその条件), 第5章 極配置法(5.1 フィードバック係数ベクトルを直接計算する) [#m441440a]
-解析から設計へ、出力フィードバック（難)→状態フィードバック（簡単、基本）
-閉ループ系のA行列 = A - BF
-A-BFの固有値をFによって任意に指定できる=可制御性
-例1(Fによる固有値の指定不可)、例1'(指定可)
-可制御性の定義
-可制御性の判定方法(可制御性行列の正則性)
-レポート対策
#ref(20191122_1.jpg,left,noimg,板書1)
#ref(20191122_2.jpg,left,noimg,板書2)
#ref(20191122_3.jpg,left,noimg,板書3)
#ref(20191122_4.jpg,left,noimg,板書4)
#ref(20191122_5.jpg,left,noimg,板書5)
#ref(20191122_6.jpg,left,noimg,板書6)
&color(red){&size(25){... 最後の[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 &ref(report9.pdf);
*第4回(2019.11.29) 第6章 最適レギュレータ §6.1 評価関数と最適制御 [#e9c4aa61]
- 簡単な例題（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(20191129_1.jpg,left,noimg,板書1)
#ref(20191129_2.jpg,left,noimg,板書2)
#ref(20191129_3.jpg,left,noimg,板書3)
#ref(20191129_4.jpg,left,noimg,板書4)
#ref(20191129_5.jpg,left,noimg,板書5)
-Q: 式が多い
-A: 基本的には数学なので、ある程度は仕方ないように思います。
レポート#10 &ref(report10.pdf);
*第5回(2019.12.6) 第6章つづき〜 §6.2 重み行列と正定・半正定、§6.3最適制御系の安定性 [#xd1595ff]
- 行列へ拡張: q→Q≧0, r→R>0
- (半)正定行列の定義
- 対称行列の固有値と正定性の関係(演習6.6), 対称行列の固有値はすべて実数(演習6.8) &ref(授業/動的システムの解析と制御2015/proof3.pdf);
- 最適制御⇒安定かつJが最小 &ref(授業/動的システムの解析と制御2015/proof4.pdf,left,証明);
- 最適レギュレータの設計(演習6.3,6.4)：リカッチ方程式の解P、P > 0 ⇒ 対応する F が閉ループ系を安定化
#ref(20191211_1.jpg,left,noimg,板書1)
#ref(20191211_2.jpg,left,noimg,板書2)
#ref(20191211_3.jpg,left,noimg,板書3)
#ref(20191211_4.jpg,left,noimg,板書4)
-Q: §6.2（例３）で
\[
\left[\begin{array}{cc} x_1 & x_2 \end{array}\right]
\left[\begin{array}{cc} 2 & 1 \\ 1 & 2 \end{array}\right]
\left[\begin{array}{c} x_1 \\ x_2 \end{array}\right]
= (x_1 + x_2)^2 + x_1^2 + x_2^2 > 0
\]
としてもよいか。
-A: 良いです。
-Q: レポート対策の
\[
(\mbox{リカッチの左辺}) = \left[\begin{array}{cc} 1 & p_1 - p_2 \\ \ast & 2(p_2 - p_3) + 1 \end{array}\right] - \left[\begin{array}{cc} p_2^2 & p_2 p_3 \\ \ast & p_3^2 \end{array} \right] = 0
\]
この式の形になるのがわからない ←（註：＊二箇所と、2(p_2 - p_3) + 1 への矢印）どうなったらこの形になるのか
-A: 口頭で説明しましたが、対称行列なので、左下は記載を省略（＊）しました。また、2(p_2 - p_3) + 1 については、PA の (2,2) 要素 p_2 - p_3 と、A^T P の (2,2)要素 p_2 - p_3 と、Q の (2,2)要素 1 の和です。
レポート#11 &ref(report11.pdf);
*第6回(2019.12.13) §9.1状態観測器の構造 [#mabe8aa8]
- 状態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(20191213_1.JPG,left,noimg,板書1)
#ref(20191213_2.JPG,left,noimg,板書2)
#ref(20191213_3.JPG,left,noimg,板書3)
#ref(20181218.jpg,left,noimg,板書4)
#ref(20191213_4.JPG,left,noimg,板書5)
- 板書□7を撮り忘れました、すみません。&br;
代わりに去年の板書（板書4の□8）を掲載します。板書内容は同じものとなっています。
レポート#12 &ref(report12.pdf);
*第7回(2019.12.20) §9.3併合系の固有値 [#qd296dbc]
- 方法3で安定化できる理由：閉ループ系の固有値 = A - BF と A - LC の固有値
- 方法3が評価関数Jの最小値に与える影響？
- A-LC の固有値の実部が-∞に近づくように L を設定するとどうなるか？
#ref(20191220_1.jpg,left,noimg,板書1)
#ref(20191220_2.jpg,left,noimg,板書2)
#ref(20191220_3.jpg,left,noimg,板書3)
#ref(20191220_4.jpg,left,noimg,板書4)
//■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■
//&color(black,red){&size(20){!!!!!!!!!!!!!!!!!!!!!!!!!!!!!以下は過去の情報です!!!!!!!!!!!!!!!!!!!!!!!!!!};};
-Q: 大学院の制御特論とこの授業はどのようにつながるのか。
-A: この授業では、制御対象が伝達関数または状態空間表現として一つに（一点で）与えられる状況を考えました。制御工学特論では、制御対象が一点ではなく集合で与えられる場合を考えます。その上で、集合の制御対象に対して制御系が安定性で、かつある性能仕様を満たすコントローラの設計問題を扱います。
//■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■
//&color(black,red){&size(20){!!!!!!!!!!!!!!!!!!!!!!!!!!!!!以下は過去の情報です!!!!!!!!!!!!!!!!!!!!!!!!!!};};
*第8回(2019.12.27) 期末テスト, 授業アンケート実施 [#z7ad6cb6]
- 試験内容: レポート課題の中から、数値を変えて出題する
- 持ち込み可能なもの: 筆記用具のみ
- 試験時間: 85分
- 授業アンケート（本科目の前半・後半をまとめて）
&color(red){&size(25){2019.12.27 前半・後半を含めた総合成績を掲示しています。レポート返却場所にて。確認してください。採点結果に疑義がある場合は1月6日（月）までに申し出ること。};};
![[PukiWiki]](image/logo_crop.png "[PukiWiki]"){kind=logo}
