2008年10月26日日曜日

交通エレクトロニクス2007年10月19日マショレック博士特別講義のスライド

もう大分前のことになりますが、昨年度冬学期に開講の電子情報学専攻大学院講義「交通エレクトロニクス」において、2007年10月19日に、ベルリン工科大学のマショレック博士に、ヨーロッパの鉄道信号保安技術の標準化について「Interoperability of Train Protection Systems in Europe」 と言う題目の講演をしていただきました。マショレック氏が再来日し、その際のスライドのファイルを公開用にまとめて提供してくださりましたので、ここから入手していただけるようにしました。(ファイルの大きさは約5.8MBあります。)

2008年10月25日土曜日

2008/10/23 AGSセミナー 「電気鉄道と環境」御聴講の皆様へ

 拙い講演をお聴きくださりありがとうございました。

 その際、スライド中にて、補足資料はPW付きでこのページからリンクを張ると申し上げましたが、具体的に資料の追加を要するようなご質問、ご議論がございませんでしたので、このページに特段の追加の記載はいたしません。講演資料については、堀教授のものもあわせてAGSクラブの会員向けページに掲載されると伺っております。

 せっかく、このページを訪れていただいた方、追加の内容が無くて申し訳ございません
 なお、せっかく来ていただきましたので、お時間のございますときに、このページの右側にリンクがございます、古関研から独逸に留学の機会を得た3名の学生(うち2名は既に社会人として活躍しております)の滞欧記をご覧になってみてください。三者三様の若い感性で、それぞれに興味深い内容が記されていると思います。

2008年10月8日水曜日

電子・情報系進学内定者へのガイダンスを終えて(10/22リンクづけして改訂)

 本日電子・情報系進学内定の皆さんには1日にわたるいろんな先生の話を聞いていただきました。我々B2コースの話は、そのスケジュールの中でも最後だったので、聞く側としても大変疲れていたと思いますが、最後まで真面目にお付き合いいただきありがとうございました。中島先生の2030年ビジョンの話を含む、我々のコースのエネルギー・環境・グローバルシステムの技術に関わる教育・研究内容の説明に使用した資料は、今回配布しませんでしたが、近日中に学科の公式ページからダウンロードできるようになるはずです。

 当座の措置として、上記の中島先生、古関、熊田先生の講演に関わる資料はここから入手可能です。一人でも多くの前途有望な学生さんが、この分野に関心を持っていただけることを願っております。

2008年10月1日水曜日

制御CAD演習 実施要領(10/01初稿 10/8改訂)

3年生冬学期、月曜日午後には演習が行われます。堀教授と古関が担当するEA3 制御CAD演習ではPC上の制御CADソフトを用いて実践的な制御系の応答計算や設計の演習を扱います。

このブログの昨年度の記載でその様子がお分かりになると思います。本年度もほぼ同じ内容です。
(この形式での演習は本年度が最後になります。)

実際に演習に使うテキストは、ここから入手してください。(10/20に改訂版へのリンクをはり直しました。)
また、図の様式に一部不備がありますが、5, 6日目の堀教授の発展課題のテキストは、必要に応じてここから入手してください。演習当日の時間を有効に活用するため、事前にテキストに目を通し予習をして演習に臨むことを強く推奨します。

(10/02に配布された全体の実験計画の中で確認した結果、実施日程を以下の通りとします。)

場所: 工学部2号館12階古関研(123D1室)
日時: 月曜13:00-16:00

10/2(木) 午後 実験/演習全体のガイダンス

前半グループ (*はレポート課題出題日)
基礎演習: (1) 10/06 (2) 10/20* (3) 10/27 (4) 11/10*
応用演習: (5) 11/17 (6) 11/21*
レポート締切日: (2) 11/10(13:00@演習) (4) 11/21(13:00@演習) (6) 01/13(17:00@事務室)

後半グループ (*はレポート課題出題日)
基礎演習: (1) 12/08 (2) 12/15* (3) 12/22 (4) 12/24*
応用演習: (5) 01/08 (6) 01/19*
レポート締切日: (2) 12/24(13:00@演習) (4) 01/19(13:00@演習) (6) 02/16(17:00@事務室)

担当者: 堀教授 古関准教授 TA 福正博之 (博士課程1年生)

以下に、後期実験テキストに掲載予定の、演習の紹介(の一部を最新情報に基づき微修正したもの)を転載しておきます。

講義だけでは何となく納得のいかなかった部分が、自分でシミュレーションをし、時間応答波形や周波数特性を見る中で、直感的に理解できるようになることを目指していますので、エネルギーコースの人のみならず、情報系、電子系の科目を中心に履修している皆さんにも積極的にご参加いただきたいと思います。

------------------

MATLAB, Simulinkを用いた制御系CAD演習(2008年度)
(本稿、古関の不手際から、実験書には昨年度のページがそのまま掲載されておりますことをお詫びします。本年度版の記述もほぼ変わりませんが、こちらをご覧ください。)
担当: 堀洋一, 古関隆章


1. はじめに

 制御工学は電気系のみならずあらゆる工学分野で用いられている基礎的な技術であり、この演習の受講者の多くはすでに夏学期の「制御工学I」を履修していることであろう。また、駒場4学期でも制御工学と非常に関連の深いラプラス変換を用いた微分方程式の取扱を必修科目として学んでいるはずである。しかし、そこでは、例題があくまでも紙の上の計算で扱える範囲の比較的単純なものに限定されざるを得ず、面白さと言う点でも、現実的なより複雑な問題を解決するためのトレーニングを行うという点でも物足りなさを感じたに違いない。

 堀や古関の学生時代は、現在のように能力のすぐれた計算機・便利なソフトウェアが簡単に使える環境ではなく、時間応答のディジタルシミュレーションを行うことは、それ自体で数日を要する作業であった。たとえば、分母多項式を導出し、ニュートンラフソン法で多項式の解を求めるプログラムを作成、留数の定理などを用いて各時間におけるラプラス逆変換を数値的に求め、さらに、それをグラフで表現するためのプログラムを作成するという作業を経てようやく一つの応答波形を可視化することができるという具合であった。さらに制御工学Iの授業で習ったような、いわゆる「古典制御理論」が盛んに研究された頃は、解析式の導出とアナログのアンプを組み合わせて積分器や近似的な微分器の回路をつくり、その出力波形を観測しながら研究が進められていたと聞く。(これをアナログコンピュータと呼んでいた。)

 今日では、幸いにして、高速の計算機と、プログラムやモデル作成の容易なグラフィック機能も含むソフトウェアが比較的簡単に手に入るようになっている。その代表的なものがこの演習で用いるMATLABで、もともとは行列の演算を扱う数値ライブラリ(サブルーチン集)を使いよくするために、変数をわたす部分のインターフェースを良くするツールとして生まれた。その後、それをベースにグラフィクスの機能などが強化され、商用の技術計算用言語として頒布されるようになった。さらに、フィルタ設計や最適化問題あるいは制御のシミュレーションなどに便利な、Toolboxと呼ばれるアプリケーションごとのルーチン集が別売りされるようになって、CADツールとして研究機関や産業界で広く用いられるようになった。MATLAB本体はあくまでもキャラクタベースでコマンドをスクリプトとして記述しプログラムを行う形態のツールであるが、制御の分野でMATLABが広く用いられる様になったのは、ブロック線図をカットアンドペーストでグラフィカルに作成することで非常に簡単に過渡応答の計算ができてしまうSimulinkというツールがMATLAB上で動くようになってからであろう。

 本演習でもMATLABの基本機能と、Simulinkを必要に応じて使い分けながら作業を進めて行くが、その便利な機能を堪能すると同時に、安直にその便利さに溺れることなく、問題そのもののもつ物理的本質やディジタルシミュレーションに伴う様々な問題点に注意しながら、頭を使って結果を考察するようにしよう。

 与えられた問題の他にも、自分自身で発展例題を作成しあるいは参考書を参照して興味のある問題を探してそれを解いてみることは大いに推奨される。この演習で用いているMATLAB/Simulinkは、以前書籍の付録として配布されていたStudent Editionであるため非常に複雑な問題を考えると、変数の数やブロックの上限が問題となる可能性もあることに注意してほしい。現在では、最新商用ソフトの学生版 Student Version として大学生協の店舗などで、学生個人がよりよいものを購入することができる。

(詳細はソフトウェアのマニュアルを参照のこと。例:http://www.mathworks.com/access/helpdesk/help/toolbox/control/

 MATLAB/Simulinkの使い方に関しては、演習の現場でも指導を行うが、現場で閲覧できるStudent Edition のマニュアル(英語版)のほか、生協にもいくつかの成書があるのでそれを各自で適宜選択して購入するとよい。たとえばこの演習の内容に近いものとしては

西村, 野波:「MATLABによる制御理論の基礎」東京電気大学出版局1998年

などが最近出版されている。またMATLABの概要を知りいくつかのデモなどがみたければ、このソフトを日本で販売しているサイバネット社のホームページ http://www.cybernet.co.jp/products/matlabを覗いてみると良いであろう。

 演習で用いているMATLABは、オンラインヘルプもマニュアルもすべて英語で書かれているが、少し慣れれば難しいことはない。このようなマニュアルは、翻訳されたものよりも原文を読んだほうがずっとわかりやすいことを実際の研究の中でしばしば経験する。また、プログラム中のコメントやグラフの軸などもすべて英語で書かねばならない。もちろんローマ字で日本語を書くことは可能である。技術的な英語を用いることにアレルギーを起こさず、良い訓練の機会と思って積極的に活用し、必要なtechnical termを覚えてしまうと良い。


2. 日程と学習内容

 本演習は、電子情報学専攻貸与のノートPCを用いて、工学部2号館12階の古関研究室にて、教員およびTAの指導の下で行う。週1回のペースで6週で完結する。その概要を以下の表に示す。主として制御工学Iの内容に準拠している。制御工学Iの教科書の該当箇所を読み事前に準備することが可能である。具体的演習問題説明などは配布プリントとして、この古関隆章ブログ からテキストをダウンロードをしていただきたい。 制御工学IIの内容として演習の課題となりうるものとしては

2.1. ディジタル制御:Z変換
微分と疑似微分:Bode線図と波形の応答/Tustin 変換と厳密なZ変換の応答の比較
MatLabを用いたディジタルフィルタの設計

2.2. 状態方程式と伝達関数:相互の書き換え 極配置に基づく状態フィードバック
時不変線形システムの最適レギュレータ問題と結果として与えられる極配置の関係

2.3. 状態観測器の設計

2.4. 定常カルマンフィルタ

などが挙げられる。これらを本演習の基礎編で本格的に扱うことは時間の制約上残念ながらできないが、今は、制御工学第一の範囲で状態空間法入門を扱っているので、これに関する演習も原理的に可能である。4回までの演習でMATLABを用いた作業に習熟したら、その後の選択課題の中で、これら冬学期の講義内容にも関係する課題に積極的に挑戦し、授業で習った内容の理解に役立て欲しい。



第1回: MATLAB/Simulink入門、微分方程式とラプラス変換

第2回: システムの応答, システムの周波数特性とボーデ線図, 安定性の解析: ラウスの判定法, フルビッツの判定法, ナイキストの判定法

第3回: 二次系の性質と極の位置, フィードバック制御の特性, フィードバック制御の設計

第4回: PI(D)制御とI-P(D)制御と分子多項式の影響/二自由度制御, フィードバックとフィードフォワード, 状態空間法入門

選択課題: 第5, 6回 発展課題:磁気浮上, 電気自動車の制御, カオス現象の解析と制御など

3. 参考書

なお、演習に役立つであろう制御理論に関する参考書としては、講義で指定している堀の教科書「制御工学の基礎」の他、

[1] 茅陽一:「制御工学第一」

[2] 細江繁幸:「システムと制御」オーム社

[3] 伊藤正美:「制御理論演習」 昭晃堂

[4] 平井一正, 羽根田博正, 北村新三:「システム制御工学」森北出版

などがあるので適宜自分の気に入ったものを参照して欲しい。


4. 有用なフリーソフトに関して

 このガイダンスの文書は、Sun Microsystemsが本学の情報基盤センターを通じて学内関係者にフリーで使用権を認めてくれているStarSuiteというソフトを用いて作成している。本演習ではレポート用のワープロ、作図、数式エディタを用いてレポートをまとめたいという履修者の希望に対応するため、フリーソフトであるopenoffice.orgを2002年度から導入した。その後、本来商用ソフトであるStarSuiteがGoogleパックスの形で実質的に自由に使えるようになったため、StarSuiteあるいはLaTexを使用してレポーティングを行うことを推奨する。

-----------------------------------------------
 また、本演習で用いているMATLABは、比較的高価な商用ソフトであり、自宅学習用に配布することはできない。(演習で用いているソフト自体が、機能限定のついた学生版の古いバージョンのものである。)しかし、MATLABに対しても俗称MATLABクローンと呼ばれる優秀なフリーのソフトウェアが存在する。自習用にはフリーソフトSciLabをhttp://www-rocq.inria.fr/scilab/からダウンロードして、自分のPCにインストールし使用してみることも奨める。これらフリーソフトに関しては


早稲田大学の大石先生による数値計算の基礎教育のページ
http://www.oishi.info.waseda.ac.jp/~oishi/lec2001/l-1.htm

メディアラボのホームページにあるフリーの科学技術計算用のプログラム
http://www.mlb.co.jp/linux/science/

特に、シミュリンク相当のツールもあるプログラムの紹介
http://www.geocities.jp/rui_hirokawa/scilab/

などのページが参考になろう。(2008年10月01日リンク有効確認済。)

インターネットのキーワード検索で「Matlabクローン」などの用語で検索すると、様々な関連プログラムの情報が得られる。なお、これらのフリーソフトウェアに関する情報は、古関の個人HP (http://www.geocities.jp/takafumikoseki/)にも掲載している。Linux環境上で理工学系の仕事をすることは、最近のユーザフレンドリなOSの普及によりかつてなく容易になっている。ここで述べた計算やレポート作成のための環境は

Ubuntu: http://www.ubuntulinux.jp/

Fedora9: http://fedoraproject.org/ja/

などを用いることで、今日では非常に安価に個人でも構成することが可能である。

3年生冬学期講義「制御工学第二」の進め方に関して(10/1暫定版)

(変更は、適宜修正・付記します。)

本講義の指定教科書はありません。
 「新: 制御理論の基礎 昭晃堂」は、特に制御工学全般を短時間で見渡すのに優れた参考書であると思います。
 下に推奨参考書として挙げている、堀洋一著「応用制御工学」(丸善出版)は現在入手困難と思われます。図書室などで必要に応じてご覧ください。

 (講義にご参加いただいている中でより詳しく学びたい皆様には、
著者である堀先生のご厚意により堀先生の書かれた教科書用原稿電子ファイルを上巻相当の内容下巻相当の内容)から特別に御覧いただくことができます。なお、左記のリンク情報は2週間ほどをめどに削除の予定です。)

10/7の講義で説明した、夏学期 制御工学第一の状態フィードバックに関係する期末テストの解説(デモ)に用いたMatlabのソースを以下に示しておきます。

% ex081007.m: CTRL example  program coded by KOSEKI, OCT 2008
 
clear all
%-----Definition of physical constants
 %%%---Fundamental state feedback---
%
A0=[0,1; 10,0 ]
B0=[0;1] 
C0=[1,0]
D=0

%---Pole placement---
Te=1/sqrt(5) % Equivalent settling time
% --3rd order Kessler canonical form---
H=[Te^2/2, Te, 1];
P=roots(H)

F=acker(A0, B0, P)

sys0=ss(A0, B0, C0, D)
sys1=ss((A0-B0*F), B0, C0, D)


figure(1); pzmap(sys0); sgrid
figure(2); pzmap(sys1);sgrid

figure(3); step(sys1)


--以下は第1回目の講義の導入資料--
制御工学第II (火曜日第一限 242教室)
古関隆章・堀 洋一
Ext. 26676, takafumikoseki @ ieee.org

制御工学第二講義予定
(01) 10/07 状態空間におけるシステムの取扱(1):夏学期の復習、状態変数と状態方程式、状態遷移行列
(02) 10/14 状態空間におけるシステムの取扱(2): 実現問題、時間応答シミュレーション
(03) 10/21 状態空間におけるシステムの取扱(3):状態フィードバック、最適制御、演習出題
(04) 10/28 堀教授特別講義 (電気自動車とモーションコントロール)
(05) 11/04 状態空間におけるシステムの取扱(4): 状態推定、システム同定
(06) 11/11 状態空間法、最適制御演習
(07) 11/18 ディジタル制御(1): 空間量子化と時間量子化
(08) 11/25 ディジタル制御(2): サンプリング動作の数学的表現、サンプリングを含む系の取扱とZ変換
(09) 12/02 呉助教特別講義 (福祉制御工学)
(10) 12/09 ディジタル制御(3): ディジタル系安定性判別、ディジタル系のアナログ近似、演習出題
(11) 12/16 ディジタル制御(4): 有限整定とデッドビート応答
(12) 01/13 状態推定法、ディジタル制御演習
(13) 01/20 非線形制御入門 演習出題
(14) 01/27 予備日
(##) 02/03 or 02/10 期末テスト

制御工学第一の内容(復習)
1.序論
 1.1  制御工学とは何か
 1.2  制御システムの例,分類
 1.3  閉ループ制御と開ループ制御
 1.4  制御工学の歴史

2.システム動特性の表現
 2.1  信号伝達と状態遷移
 2.2  動作点まわりの線形化
 2.3  線形システムの表現
(2.4  ラプラス変換法の基礎)
 2.5  ブロック図とその合成
 2.6  特性の計測法

3.制御システムの安定性
 3.1  線形システムの安定性
 3.2  ラウスの安定判別法
 3.3  ナイキストの安定判別法
 3.4  フィードバック系の安定度指標

4.フィードバック制御系の基本特性
 4.1  入力追従特性と外乱抑圧特性
 4.2  定常誤差と誤差係数
 4.3  2次系の過渡応答
 4.4  2次系の周波数応答
 4.5  高次系の代表根

5.線形フィードバック系の補償
 5.1  フィードバック制御系設計の基本指針
 5.2  直列補償とフィードバック補償
 5.3  ニコルズ線図を用いた制御系設計
 5.4  根軌跡を用いた制御系設計
 5.5  直列補償とPID調節計
 5.6  フィードバック補償

制御工学第二から制御工学第一への移行
 古典的制御の復習と状態空間法との関係(フィードバック制御と多項式法など)
 状態空間におけるシステム序論:連続量と離散量

---

2. 教科書・参考書
新: 制御理論の基礎 昭晃堂
堀・大西: 応用制御工学  丸善
ただし、これは教科書指定ではなく、推薦参考書なので、授業では該当ページを示すのみで板書の内容と一対一の対応にはなっていない。主として、自宅学習や演習問題を解く際の参考として活用してほしい。具体的にはいろいろと探してみて自分に合うと思う本を買って学ぶのが良い。

平井・羽根田・北村: システム制御工学  森北出版
金原・黒須: ディジタル制御入門  日刊工業新聞社
小郷・美多: システム制御理論入門 実教出版社
正田: 制御工学  培風館 
前田・杉江: アドバンスト制御のためのシステム制御理論 
システム制御情報学会編 朝倉書店

MATLAB: SimuLink, Control tool boxなど
計測制御学会誌、電気学会雑誌、論文誌など

3. 授業の受け方

式を自分で追ってみる-------演習問題。
月曜午後の演習のようにSCILABなどを用いたシミュレーションもできると良い。
演習問題のレポートをすべて期限内に提出することを、期末試験を受ける条件とする。

A4の方眼紙をノートとせよ: 授業は基本的に板書を中心に行う。
色鉛筆を用意する-------作図の理解に便利。
昨今は振替休日のため、月曜日の講義の回数が少なく、最後に時間がなくなって非線形制御の解説は手薄になりがち。冬休みなどに「基礎制御工学」の相当部分を自主的に読んでもらえるとありがたい。

本資料の更新版をはじめ、講義中の資料は必要に応じて、http://takafumikoseki.blogspot.com/ から入手可能なようにアップロードします。

2008年度冬学期のスタート(10/01初稿)

しばらく、本ブログの更新が止まっておりましたが、本日冬学期の初日を迎えました。本学期は、古関は制御工学第二と制御CAD演習を主として担当いたします。

 制御CAD演習は、昨年度も述べたとおり、必ずしも演習の履修者でなくても、このページからテキストを入手し、フリーの制御CADソフトを活用することで自習ができます。このページから、テキストをダウンロードしご活用ください。

 制御CAD演習の履修者は、時間を有効に使えるよう予習をして演習に臨むよう努めてください。良い学習成果を挙げられるよう、皆様御自身の御努力と講義・演習へのご協力ををよろしくお願いします。