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

授業

第1回(2023.9.8) ガイダンス、前半（古典制御）概要～第5章周波数応答†

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

レポート#1 report1.pdf

↑

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

レポート#1を返却しています
講義スライド ... 第1回のつづき（5.4.2 ボード線図の利点）から
演習問題#1（第1回で公開済）の解説をします（解答例 exercise01A.pdf exercise01slide.pdf）

レポート#2 report2.pdf

↑

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

講義スライド slide03.pdf
演習問題#2 exercise02.pdf 第4回の授業で解説する予定です。事前に問題を解いておいてください。

レポート#3 report3.pdf

↑

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

講義スライド ... 第3回のつづき（6.3 ゲイン余裕、位相余裕）から
演習問題#2（第3回で公開済）の解説をします（解答例 exercise02A.pdf exercise02slide.pdf）

レポート#4 report4.pdf

↑

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

講義スライド slide05.pdf §8.2 PID補償による制御系設計（p.24）まで
演習問題#3 exercise03.pdf 次回以降の授業で解説する予定です。授業範囲内の問題は、事前に解いておいてください。

レポート#5 report5.pdf

↑

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

講義スライド ... 第5回のつづき（8.3 位相進み-遅れ補償による制御系設計）から

レポート#6 report6.pdf

↑

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

講義スライド ... 第6回のつづき 55枚目から
演習問題#3（第5回で公開済）の解説 exercise03.pdf 解答例 &ref(): File not found: "exercise03A.pdf" at page "授業/動的システムの解析と制御2023"; exercise03slide.pdf

↑

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

学生証を提示の上、指定された席で受験すること。
筆記用具（定規含む）のみ持ち込み可、電卓などは不可。
85分間

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

↑

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

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

教科書：「演習で学ぶ現代制御理論」(森泰親著、森北出版)
古典制御と現代制御の違い（長所と短所）
伝達関数と状態空間表現の関係（簡単な運動方程式を例に）

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

レポート#7 report7.pdf

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

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

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

Q: やったことのある内容だったがほとんど覚えていなかったので復習したいと思った
A: やったことがある、と覚えていれば話は早いです。

Q: 古典制御と比べて現代制御の欠点は何ですか。
A: 明確に回答することはできませんが、一つの例を挙げるとすれば、PID制御は現在でも産業界などで幅広く利用されています。パラメータが三つだけでコントローラの動作が直感的にイメージしやすい（現代制御ではしにくい）ということはあるかもしれません。

Q: なぜ伝達関数と状態空間表現の変換を行う必要があるのですか。
A: 古典制御だけ、または、現代制御だけ、で満足するなら「変換」の必要はありません。一方、両者を習得する必要がある場合、両者は完全な別物ではなく関係があるので、「変換」を行うことが「習得」の助けになると思います。

Q: ・SSR→TFは可制御正準系などの特定の形を利用しなければ変換することが出来ないのですか。
A: （可制御正準形で与えられているかどうかによらず）SSRの A,B,C 行列が与えられれば、対応する TF は一意に定まります。授業で説明した通りです。

Q: ・可制御正準形以外の変換の形はありますか。また、それがある場合どの形を利用するかの判断方法はありますか。
A: 例えば、可制御正準形の双対の、可観測正準形があります。今日示した通り、伝達関数の係数との対応関係が明らかです。一方、他に教科書で説明されているものに対角正準形があり、伝達関数の極（固有値）がA行列中に対角に並んだ構造をしていることから、変換後の状態ベクトルをバラバラに扱える（状態ベクトルの各成分が、それぞれの固有値に依存した指数関数で振る舞う）利点があります。つまり、目的次第で使い分けることになります。

↑

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

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

レポート#8 report8.pdf

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

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

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

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

↑

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

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

レポート#9 report9.pdf

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

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

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

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

↑

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

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

レポート#10 report10.pdf

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

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

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

↑

第13回(2022.12.2) 第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 が閉ループ系を安定化

レポート#11 report11.pdf

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

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

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

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

↑

第14回(2022.12.9) §9.1状態観測器の構造†

状態xが使えない場合
(方法1) 状態の代わりに出力yを使う = 静的出力フィードバック ⇒ NG
別の方法：状態を推定して、それをxの代わりに使う
状態観測器の定義：t→∞で x(t) の推定誤差が0となる
(方法2) 状態観測器? (演習9.1) + 状態フィードバック ⇒ NG
(方法3) 状態観測器((9.3)式, 演習9.2) + 状態フィードバック
状態観測器を作る（（9.3)式の導出、方法２は状態観測器にならない）
可観測性（可制御性との関係）
演習9.3'：A - L C を安定（固有値の実部がすべて負）とする L の求め方

レポート#12 report12.pdf

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

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

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

↑

第15回(2022.12.16) §9.3併合系の固有値†

方法3で安定化できる理由：閉ループ系の固有値 = A - BF と A - LC の固有値（分離定理）
方法3が評価関数Jの最小値に与える影響
古典制御との関係：スカラの場合を例に
LQRからLQG、ロバスト制御へ

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

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

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

↑

第8回(2022.12.23) 期末テスト, 授業アンケート実施来年に延期†

週末の寒波の影響や大学へのバスの復旧目途が立っていないため、年内は対面で授業・試験等を実施しないよう教務係より指示がありました。よって、期末試験の実施を来年に延期します。新しい試験日程については後日連絡します。(2022.12.22 )

試験内容: レポート課題の中から、数値を変えて出題する
持ち込み可能なもの: 筆記用具のみ
学生証を提示の上、指定された席で受験すること。

↑

(2022.12.23) 期末試験は実施しません†

来年に延期した期末試験について、日程を検討した結果、実施可能な期間が短く、受講生全員を集めて対面で試験を実施することが困難であることがわかりました。そこで、期末試験の実施は無しとし、後半の成績は既に提出済の6回分のレポート36点満点を100点満点に換算して評価することとします。つまり、成績の評価方法について、初回のガイダンスで説明した内容を次のように変更します： ●変更前

       前半：レポート(36% = 6点×6回=36点)、中間テスト(64%)
       後半：レポート(36% = 6点×6回=36点)、期末テスト(64%)
       前半と後半の平均点を最終成績とする。

↓ ●変更後

       前半：レポート(36% = 6点×6回=36点)、中間テスト(64%)
       後半：レポート(100% = (6点×6回=36点)÷36×100)
       前半と後半の平均点を最終成績とする。

話が二転三転して大変申し訳ありません。レポートのみでの成績評価は、新型コロナウィルス対応のため昨年度及び一昨年度にも実施しており、今回の変更もそれに倣った措置です。

なお、この変更が原因で不合格となる学生はいません（もし期末試験が100点満点だったら合格となったのに今回の変更のために不合格となる学生はいない）。

前半・後半を含めた総合成績を掲示しました。レポート返却場所にて。大学に登校できるようになったら確認の上、採点結果に疑義がある場合は1月6 日（月）午前中までに申し出てください。ただ、前半の成績は10月に掲示して周知済であること、後半のレポートはすべて返却していることから、総合成績は自分で計算できるはずです。

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

第1回(2023.9.8) ガイダンス、前半（古典制御） 概要～第5章 周波数応答†

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

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

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

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

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

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