2012/4/27 ロバスト制御系設計特論計算機シミュレーション特論シ 川邊武俊 1 モデリング 本講義 2012/04/13 物理モデルの導出( 別資料 ) Black box モデルの導出 ( 別資料 ) ロバスト制御理論 線形ロバスト制御 周波数整形 拡大系 内部安定性 小ゲイン定理 ロバスト
|
|
- まいえ こうい
- 5 years ago
- Views:
Transcription
1 ロバスト制御系設計特論計算機シミュレーション特論シ 川邊武俊 モデリング 本講義 /4/3 物理モデルの導出( 別資料 ) Black bo モデルの導出 ( 別資料 ) ロバスト制御理論 線形ロバスト制御 拡大系 内部安定性 小ゲイン定理 ロバスト安定性 一般化プラント H 制御 H 制御 (McFarelae Glover の方法 ) Sldg Mode 制御 ( 非線形ロバスト制御 )( 別資料 ) 計算機シミュレーションによる制御系の評価 Xah/Sebuld による設計 シミュレーション実習 ( 別資料 ) ( 数値積分法 ) システムの不確かさとロバスト制御 Coroller Noal la はじめに はじめにシステムにはなぜ不確かさが存在するか? モデリング : システムの特性を微分方程式であらわす. 物理モデル Black bo モデル ( 同定モデル ) いずれのモデルにも精度に限界がある : モデル化誤差 (uodeled dac). 制御系の設計に用いた制御対象 ( プラント ) のモデルとプラントとの動特性の差. モデル化誤差はシステムであることに注意すること. モデル化誤差が 大きい と 制御系は不安定になることがある. あるモデル化誤差が存在しても所定の性能を有する制御系をロバスト (Robu) な制御系という. 3 4 モデル同定 ( 伝達関数同定 ) u( k ) ( k ) システム 時間 時間 同定モデル 入力の時系列データと出力の時系列データからパルス伝達関数の係数を推定する. 5 k ( ) = ak ( ) ak ( ) ak ( ) bu( k) b u( k ) ( z) b bz b z u( z) az az = A( q) ( k) = B( q) u( k) q: 進み演算子 q ( k ) = ( k) A( q) = aq aq B( q) = b bq bq 6
2 次の伝達関数であれば ( k ) = a ( k) bu( k) より () = a () bu() (3) = a () bu() () () u() a = Z θ (3) = () u() = b = Zθ 一般的にこう書ける Z - が存在すれば θ = Z = Zθ Z 7 ところが観測雑音が存在する Aq ( ) k ( ) = Bquk ( ) ( ) ek ( ) = Zθˆ Z そこでデータを多くとる () () u() e() (3) () u() e() = θ ( N) ( N ) u( N ) e( N) Z e ( 縦長 ) XZが正則となるX = Zθ e X XZ Xe = θ ( ) ( ) ( ) θˆ ˆ = XZ X θ = XZ X XZ Xe ここが小さくなればよい推定値になる 8 最小 乗法 X = Z ˆ θ = ととる ( Z Z) X = Zθˆ ε ( k ) 正則となる条件 E 条件 u(k) が周波数成分を十分に含む Z J LS = N = ε ε ε = Z ˆ θ N = ( Z ˆ θ) ( Z ˆ θ) N を最小にする θˆ を求める ( ε () ε () ε ( N )) ( θ ) J LS = Z Z ˆ = ˆ θ N より e(k) が白色雑音なら ˆ θ θ N 9 ところが e(k) は一般に白色雑音でないので 結論 Aq ( ) k ( ) = Bquk ( ) ( ) wk ( ) Dq ( ) ek ( ) = wk ( ) Dq ( ) DqAqk ( ) ( ) ( ) = DqBquk ( ) ( ) ( ) wk ( ) 白色雑音 D(q) : 整形フィルタ ( 白色化フィルタ ). 最小 乗法で D( q) A( q) D( q) B( qを推定する ) 次のシステムなら 次以上のシステムとして最小 乗法を用いる DqBq ( ) ( ). DqAq ( ) ( ) を約分する ( モデル低次元化法を使う ). 例 : 航空機のラダー モデル低次元化 Xah wbalace 使用 外挿的 3Hz Hz 同定結果 入力データのパワスペクトラム
3 データに関する注意 = Zθˆ = Zθˆ θ () () u() e() (3) () u() e() = θ ( N) ( N ) u( N ) e( N) を次のように修正すればよい 入力に必要な周波数帯域を確保する後輪舵角からヨーレートの伝達関数同定 ( 実験データ ) 参考文献 [5] 低周波 高周波 θ Z データのオフセットを取り除いておく ( オフセット : センサーのオフセット ) () () u() e() (3) () u () θ e() = θ ( N) ( N) u( N) e( N) この列を加える あるいは オフセットも同定するパラメータとする 3 4 同定結果 外挿的 同定モデル ( 伝達関数同定 ) のまとめ 最小 乗法は大きな次数で そしてモデル低次元化法を用いる. データのオフセットは除いておく あるいはオフセットをパラメータの つとして同定する. 入力のパワスペクトラムに注意 ( 外挿的なのはどこか ) 同定はつ以上の方法で そして比較. 5 6 はじめに 信号の 大きさ? 物理モデル r e Coroller d d 別資料 e () 講義資料 _.df 制御誤差 e() は 小さい 方が良い. 小さい の定義は? 7 8 3
4 システムの 大きさ? はじめに 説明の流れ はじめに Coroller Noal la 目的 :( ロバスト制御 ) モデル化誤差が ある限られた範囲に限定されるのであれば 安定性や制御性能を保証する ( 制御誤差が 小さい ) 制御系を設計する ( ロバスト制御設計特論 / 計算機シミュレーション特論へつづく ) そのためにシステム ( モデル化誤差 ) の 大きさ を定義する : モデル化誤差 (uodeled dac). 制御系の設計に用いた制御対象 ( プラント ) のモデルとプラントとの動特性の差. モデル化誤差はシステムであることに注意すること. モデル化誤差が 大きい と 制御系は不安定になることがある. そのために信号の 大きさ を定義するそのためにベクトルや行列の 大きさ を定義する 説明の順序 参考文献 : 児玉 須田 システム制御のためのマトリクス理論 ( コロナ社 ) 9 行列 行列式 C : 行 列 ( サイズ ) の行列全体の集合. R : 行 列 ( サイズ ) の実数行列全体の集合. 記号 C : 行ベクトル全体の集合. R : 行実数ベクトル全体の集合. A= [ a j ] とする. a j : 行列 A の第 行第 j 列要素. 正方行列 : 行と列との数が等しい行列行列式 : de A または A とかく. A= [ a ] = j = とすると j de A= g ( ) a a a ; g ( ) = ; が遇転置の場合. が奇転置の場合. A : 行列 A の転置行列. A : 行列 A の複素共役行列. A : 行列 A の共役転置行列. A = [ a j ] A = [ a j ] A = a j [ ] ( AB) = B A AB = AB ( A) = A 遇転置 : 数列 { 3 } の要素を偶数回 ( 回 ) 入れ替える操作のことをいう. 奇転置 : 数列 { 3 } の要素を奇数回入れ替える操作のことをいう. {3} {3}: 遇転置 {3} {3}: 奇転置 de A= de A de( AB) = de( BA) = de Ade A 小行列とランク ( 階数 ) 小行列 : 行列 A= [ aj ] C から 第 行 第 行 第 r 行 ( < < r ) かつ第 j 列 第 j 行 第 j r 行 ( j < j < j r ) 以外の要素を取り去った行列をr 次の小行列という. 様々な正方行列 対角行列 : 対角要素以外の要素は である正方行列. a dag( abc ) = b c 例 : の小行列は [] 6 4 など ランク ( 階数 ): r 次以上の小行列の行列式は全てで r 次の小行列には行列式がでないものが存在する. 行列のランク ( 階数 ) は r.rak A = r と書く. 3 ブロック対角行列 : 対角位置が対角行列である行列 A O O dag( ABC ) = O B O O O C O : ゼロ行列. すべての要素が である行列. 4 4
5 様々な正方行列 正則行列 : Aは正方行列とする. AX = XA = I がなりたつ逆行列 X が存在する. Aは正則行列である. X = A とかく. ただし IはAと同じサイズの単位行列を表す. 正則行列の性質 A C 正規行列 : AA = A A 対称行列 : A= A 様々な正規行列 は自然数とする ( 正方行列 ). 正規行列エルミート行列ユニタリ行列 A が正則行列 de A ( AB) = B A ( A ) = ( A ) : = A A A A A AA O A = O A Mar vero lea ( A BD C) = A A B( DCA B) CA 5 直行行列 :( ベクトルの回転を表す ) A = A AA = A A= I エルミート行列 : A ユニタリ行列 : = A A = A AA = A A= I 複素数への拡張 エルミート行列対称行列ユニタリ行列直行行列 6 固有値 固有ベクトル 固有値 固有ベクトル A は正方行列とする. A C C にたいし A = λ λ C が成立する λ をそれぞれ A の固有ベクトル 固有値という. ρ ( A ) = a λ = をスペクトル半径と呼ぶ. de( λi A) = を固有方程式と呼ぶ. 相異なる固有値に対応する固有ベクトルは 互いに線形独立である. 固有値が重複する場合 : 一般化固有ベクトル Jorda 標準形. 7 A はエルミート行列 A の固有値は全て実数である. = λ = A A R A はユニタリ行列 A の固有値は全て絶対値 の複素数である. = λ = = A A 下三角行列 上三角行列の対角要素はその行列の固有値である. 下三角行列 上三角行列 8 = R とする. R に対して > 正定行列 は正定 の固有値は全て正の実数. は正定 λ ( ) λa ( ) ならば は正定行列 > と書く. ならば は準正定行列という. と書く. λ ( ): λ ( ): a の最小固有値. の最大固有値. 9 A C 特異値 とする.A は次のように特異値分解できる. A= UΣV U C U U = I V C V V = I Σ R Σr O Σ= O O Σ r = dag( σ σr) σ > = r r : A のランク σ : Aの特異値. 写像 = Aは ベクトル の 回転 σ ( A): 行列 Aの特異値. 要素別の伸縮( 次元の変化も含む ) 回転 σ a ( A): 行列 Aの最大特異値. と考えられる σ ( A): 行列 Aの最小特異値. 3 5
6 特異値の性質 A C にたいし の固有値は非負になる. AA AA の固有値を λ とすると R σ = λ = とする. σ ( A) A σa ( A) ノルム ノルムの公理 : ノルム は空間の元 に対し次の ) = = ) ( 三角不等式 ) 3) a = a a C を全て満たす. ノルムが定義された空間をノルム空間という. ベクトル 行列 信号に対してノルムを定義することができる. A σ a ( A) a = a A = A σ ( A) = A = 3 3 ノルム = = ベクトルノルム ただし は自然数 R とする. = = = = ( ) ( ユークリッドノルム ) = = a ノルム: = = a = = a ベクトルノルムの等価性 [ ] にたいし 定数 k k が存在し が成立する. k k ある ノルムで評価して が有界なら 別の ノルムで評価しても は有界である. 例 : = = である をノルムで評価すると l = = a 行列のノルム ( 誘導ノルム ) ベクトルノルムから行列のノルムを 誘導 する. ベクトル が A で写像され になったと考える. と の長さ ( ノルム ) の比から A のノルムを 定義する. の次元は同じとは限らない. 行列の ノルム : A A = a C = a = C A ただし A C. A A = a A = a aj = C j 行列の ノルム = a a a a A = a j a j a j a j a a a a を に規格化してもよい. = での最大値を与えるのは A のそれぞれの列の絶対和の最 = A 大値である
7 例 : A = = 誘導ノルムの定義より A A A A = a = λ a ( ) = C 行列の ノルム A = 3. 3 λ ( A) a は A の最大固有値を表す. a A = C を与える は = a A = 4 = C = A とすると 特異値分解により A= U V = A A= V dag( σ σ ) V ここで VV =I 誘導ノルムの定義より A = a A = a a j = C 行列の ノルム j= 例 : A = a = 3. 3 a a a a A = a j a j a j a j a a a a a で A の最大値を与えるのは A のそれぞれの行の絶対和の 最大値である. A = a{ 3. } a A = C を与える は 3 = a A = 4. = C 39 4 定義 : A R とする. r( A) = a = 性質 : ) r( AB) = r( BA) A R B R トレース (race) A= [ A A A ] B = [ B B B ] A R j = Bj R j = j たてながベクトル B B B 横長ベクトル B = B R = (AB と BA とはサイズが異なることに注意 ) A A = = A A A R とおくと とおくと r( AB) = A B = a jbj r( BA) = B A = bjaj j = = = j = = = 4 ) r( ) r( ) AA A A a トレース (race) = = j= = ) よりあきらか r 3) r( AA ) = σ ( A) = 特異値分解より AA = ( UΣV )( VΣ U ) = UΣ U ) より r( AA ) = r( UΣ U ) = r( U UΣ ) = r( Σ ) j 4 7
8 信号のノルム 信号ノルムの定義 信号ノルム : ( ) u( ) R R L ノルム ( 時刻 ) により変化する 次元ベクトルの実数への写像. 次の公理を満たす. u () = u() d = ) u ( ) ) u ( ) = u ( ) 3) au( ) = a u( ) a R 4) u () v () u () v () L ノルム u () = u () ud () L ノルム u () = a u u () ( ) e e e = システムのノルム 線形時不変システムのH ノルム u G システムの入力 u R と出力 R からシステムのノルムを考える. 仮定 : システムは安定かつ線形時不変とする. () = G ( τ ) u ( τ ) d τ とすると () = G ( τ) u( τ) dτ G () = gj () : G ˆ () = LG [ ()] = gˆ j () Gが安定 gj () がすべて絶対可積分 gˆj () が複素左半平面で解析的 ( 極をもたない ) システムGのインパルス応答システムGの伝達関数 ( 行列 ). ただし L[ ] はラプラス変換. 45 定義 : ˆ G( ) = r ( Gˆ( jω) Gˆ( jω) ) dω σ { ( )} G jω dω π = π = = π = = j= = j= g ( jω) dω j パーセバルの定理 () d = ˆ( jω) ˆ( jω) dω より π ˆ G () = r ( G () G ()) d g () d j 入力 出力系について物理的な意味を考察せよ. 46 インパルス応答の状態空間表現 システムの状態方程式を () = A() Bu() () = C () とする. A A ( τ ) () = e () e Bu () τ d τ () = C () より インパルス応答は A G () = Ce B と表される. [ g ( )] j 入力 u から出力 j へのインパルス応答を並べた行列 47 状態空間表現とH ノルムの計算 { } ˆ G () = rblb o ここで L o はリヤプノフ方程式 ALo LA o =CC の正定対称解である. 証明 : ˆ G () = r ( () ()) G G d A A = r Ce B Ce B d = {( ) ( )} r( B e A C Ce B) d A = B r( e C Ce ) db A A 一方 d { A d CCe } = Ae CCe e CCAe A A A A A であるから上式左辺は d { e d C Ce } d = e C Ce = C C A A A A 右辺は Ae = e A を使って次のように書ける. A A Ae A CCe d e CCe Ad A A A そこで Lo = e C Ce A d とおくと A ALo LA o =CC が成り立つ. 48 8
9 状態空間表現と H ノルムの計算 { } ˆ G () = rclc C ここで L 注意 C はリヤプノフ方程式 不安定なシステムのノルムは定義しない. ALC LCA =BB 安定でも強プロパでないシステムのノルの正定対称解である. ムは定義できない. 証明 : 強プロパでないシステム : r ( AA ) =r( A A) であるから () = A() Bu() () = C () Du () D O ˆ G () = r ( G () G ()) d あるいは Gˆ ( ) O = r GG () () d = ( ) A r( Ce BB e C ) d A 49 以下省略 とすると 上式の右辺は : H ノルムの別解釈 H ノルムは 白色雑音に対する定常応答の 乗平均である. E[ u ( )] = E[ uu ( ) ( τ)] = δ( τ) I E[ ( ) ( )] = E { [ g ( α) u( α) ] dα [ g ( β) u( β) ] dβ } E { [ g( α) u( α) ] [ g( β) u( β) ] dαdβ } ところで A = r( A) = r( A ) = r( A) なので E { [ g ( α) u( α) ] [ g ( β) u( β) ] dαdβ } { g α g β u α u β dαdβ } = E r ( ) ( ) ( ) ( ) 5 インパルス応答は確率変数ではないので つづき { gα gβ uα uβ dαdβ } E r ( ) ( ) ( ) ( ) = r g ( α) g( β ) E( u( α ) u( β )) dα dβ = r g ( α) g ( β) δ( β α) Idαdβ まず α で積分すると r g ( α ) g( β ) δ ( β α) Idαdβ = r g ( β) g ( β) dβ 伝達関数の表記法 システムの状態方程式を () = A() Bu() () = C () Du () とするとその伝達関数は G ˆ () = CI ( A) B D この伝達関数を次のように表記する (Dole の表記法 ) ˆ A B G () = C D 因果性 (g = < ) を考慮して変数変換すると r g ( β) g( β) dβ = r g() g() d = Gˆ 5 5 システムGの状態方程式を () = () u() () = () とする. ) インパルス応答 () = e より H ノルムの計算例 ˆ G = e d e = = ) 複素関数の積分 ˆ ˆ ˆ ˆ G = G dω G Gd π j Re ( Gˆ) π = π j = C π j C ˆ ReG = = = 定義 : { } ) Gˆ = u σ ˆ a G( jω) ω あるいは ˆ () ) G = u u u () = u ( ) u = H ノルム Gが 入力 出力系であれば 周波数伝達関数のゲインの上界 あるいは 出力のノルムを最大化する 最悪な 入力にたいする出力のノルムの値である. 3) 可観測グラミアン L ˆ O LO = LO = G = = 53 詳細は後述 54 9
10 線形ロバスト制御に係わる事項 内部モデル原理 ( 古典制御 ) ロバスト安定性ロバスト性能 拡大系 内部安定性 小ゲイン定理 H 制御 (H 制御 ) 線形制御系の設計方針 d () r () e () u () () () () () 以下の伝達関数を考える. () () () () () () = r () d () () ( ) ( ) ( ) ( ) ( ) ( ) () () () () u () = r () d () () ( ) ( ) ( ) ( ) ( ) ( ) ( ) e () = r () d () () ( ) ( ) ( ) ( ) ( ) ( ) 仕様 : 制御誤差 e(jω) は小さい方がよい ( 制御の性能は良くしたい.) ()/r()=/() () r(). 仕様 : 制御入力 u(jω) は小さいほうがよい ( 制御にかかるコストは小さい方が良い.) ( jω) e( jω) = r( jω) d( jω) ( jω) ( jω) ( jω) ( jω) ( jω) ( jω) ( jω) ( jω) ( jω) ( jω) ( jω) u( jω) = r( jω) d( jω) ( jω) ( jω) ( jω) ( jω) ( jω) ( jω) ( jω) (jω) の特性は与えられているから 上式より A) 仕様 を満たすには (jω) あるいは (jω)(jω) は大きい方がよい. B) 仕様 を満たすには (jω) あるいは (jω)(jω) は小さい方がよい. : 同じ周波数で A)B) を同時に満たすことはできないので e(jω) を小さくしなければならない周波数帯域 ( 制御帯域 ) で (jω) あるいは (jω)(jω) を大きくし u(jω) を小さくしなければならない周波数帯域で (jω) あるいは (jω)(jω) が小さくなるようにコントローラ () を設計する. 57 内部モデル原理 d () r () e () u () () 上図のフィードバック制御系において r () = または d () = であるとき l e ( ) = () () () となる ( 定常誤差が に収束する ). 上図のフィードバック制御系は安定であり () は = に少なくとも 個の極を持つ. 58 内部モデル原理 ( 一般化 ) r () e () u () () r () = r () または d () = であるとき a d () a l e ( ) = d () () () () となる ( 定常誤差がに収束する ). 上図のフィードバック制御系は安定であり () はr () d () と同じ極をつ以上持つ. () が無限大のゲインを持つ周波数で外乱の影響は に収束する. 59 一般的な r(jω) (jω) d(jω) の周波数帯域と一巡伝達関数の の基本方針 : () を次のように設計する. a) 制御帯域内の一巡伝達関数のゲインを上げる b) 交差周波数を上げる ( 制御帯域を広げる ). c) 位相余裕を大きくする d) ゲイン余裕を大きくする パワースペクトル ゲイン 位相 -8 r( jω) d( jω ) ( jω) 制御帯域 a) b) c) ω 交差周波数 ( クロスオーバー周波数 ) d) ω ω 6
11 拡大系 d () r () e () () a () u () W() () () () u d r 内部安定性 左の閉ループ系の入出力関係は ( I ) ( I ) r = u ( ) ( I ) d ()W() を制御対象であると仮に考え ()W() を安定化するコントローラ () を設計する. 実装するコントローラはW() () とする. ()W(): 拡大系. 可制御 可観測となるようにW() を選ぶ. W(): 積分器 ( 型サーボの場合 ) ローパスフィルタなど. () の設計には現代制御理論やH 制御理論 H 制御理論を使う. ( I ) ( I ) ( ) ( I ) がすべて安定なとき 閉 ループ系は内部安定であるという. 制御系には内部安定性が必要である. 6 6 内部安定でない制御系 内部安定でない制御系 例 : 極零相殺 () () u() r () = 5( ) u () () u() が発 - 散する 時刻 [ 秒 ] 63 例 : 極零相殺 () r () = ) u() () d( ( ) ( ) d () u ().8 () 時刻 [ 秒 ] 外乱の影響 : 出力が振動する. 64 制御対象の不安定零点をコントローラの極で相殺する場合 u () () () () = () = と分母多項式 分子多項式で表す. d () () d () ここで () = の根は複素右半平面にあるとする. u () () () () () () () d () () d () d () d () u u = = となる. u () () () () d() d() () () = = ( ) ( ) () () d ( ) ( ) d () ( ) d () d () () () d() d() () u () d () u () d () () d ( ) = = = ( ) () () () d () () () () d u () { d ( ) d( ) ( ) ( ) } d ( ) d () は分母の () のために不安定となる (u()/r() が不安定 ). u が安定であったとしても 65 制御対象の不安定極をコントローラの零点で相殺する場合 () du () () () = () d ( ) d ( ) = d ( ) と分母多項式 分子多項式で表す. u ここで d () = の根は複素右半平面にあるとする. u () du () () () () () () = = となる. d ( ) d () d () d () d ( ) u () () () () d () d () () () = = ( ) ( ) () () d () d() () ( ) d () d () () () d () d () () = = = ( ) ( ) d( d ) ( ) ( ) ( d ( ) d ( ) du () d() du ( ) d() () d ( ) () () ) d () u { d () d() () ( ) } は分母の d () のために不安定となる (()/d() が不安定 ). u が安定であったとしても 66
12 小ゲイン定理に基づく制御系設計 小ゲイン (all ga) 定理 u d r M 定義 : { } ) Gˆ = u σ ˆ a G( jω) ω あるいは ˆ () ) G = u u u () = u ( ) u = H ノルム M は安定な有理伝達関数行列であるとする. を満たす全ての に対して 左の閉ループ系が内部安定となる必要十分条件は M < である. Gが 入力 出力系であれば 周波数伝達関数のゲインの上界 あるいは 出力のノルムを最大化する 最悪な 入力にたいする出力のノルムの値である. Nqu の安定判別との関連を考察せよ H ノルム H ノルムの計算 ) と ) の関係 { a } ˆ ˆ () () d = uˆ( jω) G( jω) G( jω) uˆ( jω) dω π ˆ ( ) ˆ( ) ˆ( ) σ π G jω u jω u jω dω = Gˆ u() u() d すなわち () Gˆ = u = u ( ) u u () u = ˆ A B G () = C O とする. である必要十分条件は行列 G ˆ () < A BB H = CC A が虚軸上に固有値をもたないこと. 証明 ( 十分性のみ ) I G G ならば G ˆ () < である. ˆ () ˆ() < 69 7 H ノルムの計算 ( 準備 ) ˆ A B G () = C D とすると { } ˆ A C G () = C( I A) B D = B ( I A ) C D = B D ˆ A BD C BD G () = D C D A B G ˆ () = C D ˆ () ˆ () G = G = とすると C D C D A O B Gˆ () Gˆ () = BC A BD D C C DD 7 以上の公式より H ノルムの計算 A O B A O B ˆ () ˆ() O B O O B I A BB B B ˆ ˆ H I G () G() CC A O O = = O B O O B O I G G = I CC A O = CC A O { } { } すなわち H が虚軸上に固有値をもたない ˆ I G () Gˆ() は虚軸上に極をもたない I Gˆ () Gˆ() = O は虚軸上に根 ( 零点 ) を持たない G ˆ () < である. G( jω) SISO システムの場合 ω 7
13 システムGの状態方程式を () = () u() () = γ () とする. G ˆ () = γ より あきらかに G ˆ () = γ いま H = γ H ノルムの計算 ( 例 ) H の固有値は de( I H) = より G ˆ () のH ノルムを計算するには 係数 γ を用い γ G ˆ () から定義される行列 H の固有値を調べる. ロバスト安定性 a ロバスト安定性制御対象に不確かさ : () < γ a() < γ a が存在しても制御系が安定であること. ただし () : 乗法的不確かさ a () : 加法的不確かさ 不確かさは安定と仮定しているところに注意. =± γ 73 不確かさのゲインは一般に高周波数帯域で大きく 低周波数帯域で小さい 74 ロバスト安定性と小ゲイン定理 ロバスト性能 a モデルに不確かさが存在しても外乱 W Wa () () となるような重み W を用いて (SISO の場合 ) W < となるように を設計する. r r e e W W W d d (d ) から制御誤差 (r) への () ノルムが所定値以下である制御系をロバスト性能な制御系という. e () WS () < d () W (): 制御仕様に関する周波数重み一般に 制御誤差は低周波数帯域 で小さければよい ロバスト性能と制御系設計問題 問題の一般化と一般化プラント r ロバスト性能を有する制御系の設計問題 e W W d () () = ( ) < < () () に対して 制御系が内部安定となる コントローラ を設計せよ. () G u z W w z W G: 一般化プラントと呼ぶ. 制御対象と制御仕様とを表した プラント 拡大系の一般化 以下の状態方程式で表現する. = A B w B u z = C D w D u = C D w D u M () W W 一般化プラントの一例 z: 評価出力 : 測定出力 u: 制御入力 w: 外乱
14 最適制御問題としての一般化 混合感度問題と H 制御 z = F( G ) w l Fl ( G ): G のに関する下方線形分数変換 H 最適制御問題 Fl ( G ) は内部安定 かつ Fl ( G ) を最小化するを求めよ. H 最適制御問題 w u G W z [ ] 混合感度問題 z = z z z w z z () = W () () w () G W W z() = W () S() w () u () () 相補感度 : () = ( ) ( ) 感度 : S () = ( ) ( ) () S() = F ( G ) l F ( G ) l は内部安定 かつを最小化するを求めよ. w から z への伝達関数 [ ] z = W () () () () W S w の Hノルムを最小化するコントローラ を設計せよ.( ロバスト性能の十分条件 ) 79 混合感度問題の問題点を指摘せよ. 8 H 最適制御問題の解 一般化プラント : = A Bw Bu R w R A4) w から への伝達関数は虚軸上に不変零点を持たない. すなわち z = C D w D u z R u R A jωi B = C Dw Du R rak C D = は次の仮定を満たすとする. ω A) (A B ) は可安定 かつ (A C ) は可検出. A) D は縦長で列フルランク かつD は横長で行フルランク. A3) u からz への伝達関数は虚軸上に不変零点を持たない. すなわち A jωi B = ω rak C D さらに次の条件を満たすように変数変換す る. D = O D = O D [ C D ] = [ O I] B O D = D I 8 を満たすコントローラは次のように表せる () = F( M() Q()) Fl ( G ) γ A X XA X ( γ BB B B) X C C = O l A ZL ZB M() = F O I C I O = γ A A BB X BF ZLC F = B X L =Y C ( γ ) Z = I Y X Qは Q () < γ を満たす任意の安定でプロパな伝達関数 X Y は次の Rcca 方程式の正定対称解 ( γ ) YA A Y Y C C C C Y B B = O Q = O と置いた解を中心解と呼ぶ. 8 と混合感度問題 ( jω) ( jω) ( jω) = ( jω) ( jω) ( jω) ( jω) ( jω) ( jω) S( jω) = { ( jω) ( jω) } ( jω) ( jω) ( jω) ( jω) ( jω) ( jω) W S W H 制御の最適解はオールパス関数 ω 83 H 最適制御問題の解 一般化プラント : = A Bw Bu R w R z = C Dw Du z R u R = C Dw Du R は次の仮定を満たすとする. A) (A B ) は可安定 かつ (A C ) は可検出. A) D は縦長で列フルランク かつD は横長で行フルランク. A3) u からz への伝達関数は虚軸上に不変零点を持たない. すなわち A jωi B = ω rak C D A4) w から への伝達関数は虚軸上に不変零点を持たない. すなわち A jωi B rak ω C D = A5) D = O さらに次の条件を満たすように変数変換す る. D = O D [ C D ] = [ O I] B O D D = I 84 4
15 最適なコントローラは A L () = F O A = A B F L C F = B X L =Y C 大型ロケットの姿勢制御系 設計例 X Y は次の Rcca 方程式の正定対称解 A X XA XB B X C C = O YA A Y YC C Y B B = O ( 別資料 )
DVIOUT
最適レギュレータ 松尾研究室資料 第 最適レギュレータ 節時不変型無限時間最適レギュレータ 状態フィードバックの可能な場合の無限時間問題における最適レギュレータについて確定系について説明する. ここで, レギュレータとは状態量をゼロにするようなコントローラのことである. なぜ, 無限時間問題のみを述べるかという理由は以下のとおりである. 有限時間の最適レギュレータ問題の場合の最適フィードバックゲインは微分方程式の解から構成される時間関数として表現される.
More informationMicrosoft PowerPoint - ce07-09b.ppt
6. フィードバック系の内部安定性キーワード : 内部安定性, 特性多項式 6. ナイキストの安定判別法キーワード : ナイキストの安定判別法 復習 G u u u 制御対象コントローラ u T 閉ループ伝達関数フィードバック制御系 T 相補感度関数 S S T L 開ループ伝達関数 L いま考えているのは どの伝達関数,, T, L? フィードバック系の内部安定性 u 内部安定性 T G だけでは不十分
More information航空機の運動方程式
可制御性 可観測性. 可制御性システムの状態を, 適切な操作によって, 有限時間内に, 任意の状態から別の任意の状態に移動させることができるか否かという特性を可制御性という. 可制御性を有するシステムに対し, システムは可制御である, 可制御なシステム という言い方をする. 状態方程式, 出力方程式が以下で表されるn 次元 m 入力 r 出力線形時不変システム x Ax u y x Du () に対し,
More informationMicrosoft PowerPoint - H22制御工学I-10回.ppt
制御工学 I 第 回 安定性 ラウス, フルビッツの安定判別 平成 年 6 月 日 /6/ 授業の予定 制御工学概論 ( 回 ) 制御技術は現在様々な工学分野において重要な基本技術となっている 工学における制御工学の位置づけと歴史について説明する さらに 制御システムの基本構成と種類を紹介する ラプラス変換 ( 回 ) 制御工学 特に古典制御ではラプラス変換が重要な役割を果たしている ラプラス変換と逆ラプラス変換の定義を紹介し
More information航空機の運動方程式
オブザーバ 状態フィードバックにはすべての状態変数の値が必要であった. しかしながら, システムの外部から観測できるのは出力だけであり, すべての状態変数が観測できるとは限らない. そこで, 制御対象システムの状態変数を, システムのモデルに基づいてその入出力信号から推定する方法を考える.. オブザーバとは 次元 m 入力 r 出力線形時不変システム x Ax Bu y Cx () の状態変数ベクトル
More informationMicrosoft PowerPoint - ce07-13b.ppt
制御工学 3 第 8 章 : フィードバック制御系の設計法 8. 設計手順と性能評価キーワード : 設計手順, 性能評価 8. ID 補償による制御系設計キーワード : ( 比例 ),I( 積分 ),D( 微分 ) 8.3 進み 遅れ補償による制御系設計キーワード : 遅れ補償, 進み補償 学習目標 : 一般的な制御系設計における手順と制御系の性能評価について学ぶ. ループ整形の考え方を用いて, 遅れ補償,
More information航空機の運動方程式
過渡応答 定常応答 線形時不変のシステムの入出力関係は伝達関数で表された. システムに対する基本的な 入力に対する過渡応答と定常応答の特性を理解する必要がある.. 伝達関数の応答. 一般的なシステムの応答システムの入力の変化に対する出力の変化の様相を応答 ( 時間応答, 動的応答 ) という. 過渡応答 システムで, 入力がある定常状態から別の定常状態に変化したとき, 出力が変化後の定常状態に達するまでの応答.
More informationMicrosoft PowerPoint - 10.pptx
m u. 固有値とその応用 8/7/( 水 ). 固有値とその応用 固有値と固有ベクトル 行列による写像から固有ベクトルへ m m 行列 によって線形写像 f : R R が表せることを見てきた ここでは 次元平面の行列による写像を調べる とし 写像 f : を考える R R まず 単位ベクトルの像 u y y f : R R u u, u この事から 線形写像の性質を用いると 次の格子上の点全ての写像先が求まる
More informationMicrosoft PowerPoint - chap8.ppt
第 8 章 : フィードバック制御系の設計法 第 8 章 : フィードバック制御系の設計法 8. 設計手順と性能評価 キーワード : 設計手順, 性能評価 8. 補償による制御系設計 キーワード : ( 比例 ),( 積分 ),( 微分 ) 学習目標 : 一般的な制御系設計における手順と制御系の性能評価について学ぶ. 学習目標 : 補償の有効性について理解し, その設計手順を習得する. 第 8 章
More information例 e 指数関数的に減衰する信号を h( a < + a a すると, それらのラプラス変換は, H ( ) { e } e インパルス応答が h( a < ( ただし a >, U( ) { } となるシステムにステップ信号 ( y( のラプラス変換 Y () は, Y ( ) H ( ) X (
第 週ラプラス変換 教科書 p.34~ 目標ラプラス変換の定義と意味を理解する フーリエ変換や Z 変換と並ぶ 信号解析やシステム設計における重要なツール ラプラス変換は波動現象や電気回路など様々な分野で 微分方程式を解くために利用されてきた ラプラス変換を用いることで微分方程式は代数方程式に変換される また 工学上使われる主要な関数のラプラス変換は簡単な形の関数で表されるので これを ラプラス変換表
More informationMicrosoft PowerPoint - 第3回2.ppt
講義内容 講義内容 次元ベクトル 関数の直交性フーリエ級数 次元代表的な対の諸性質コンボリューション たたみこみ積分 サンプリング定理 次元離散 次元空間周波数の概念 次元代表的な 次元対 次元離散 次元ベクトル 関数の直交性フーリエ級数 次元代表的な対の諸性質コンボリューション たたみこみ積分 サンプリング定理 次元離散 次元空間周波数の概念 次元代表的な 次元対 次元離散 ベクトルの直交性 3
More informationPowerPoint Presentation
付録 2 2 次元アフィン変換 直交変換 たたみ込み 1.2 次元のアフィン変換 座標 (x,y ) を (x,y) に移すことを 2 次元での変換. 特に, 変換が と書けるとき, アフィン変換, アフィン変換は, その 1 次の項による変換 と 0 次の項による変換 アフィン変換 0 次の項は平行移動 1 次の項は座標 (x, y ) をベクトルと考えて とすれば このようなもの 2 次元ベクトルの線形写像
More information第 4 週コンボリューションその 2, 正弦波による分解 教科書 p. 16~ 目標コンボリューションの演習. 正弦波による信号の分解の考え方の理解. 正弦波の複素表現を学ぶ. 演習問題 問 1. 以下の図にならって,1 と 2 の δ 関数を図示せよ δ (t) 2
第 4 週コンボリューションその, 正弦波による分解 教科書 p. 6~ 目標コンボリューションの演習. 正弦波による信号の分解の考え方の理解. 正弦波の複素表現を学ぶ. 演習問題 問. 以下の図にならって, と の δ 関数を図示せよ. - - - δ () δ ( ) - - - 図 δ 関数の図示の例 δ ( ) δ ( ) δ ( ) δ ( ) δ ( ) - - - - - - - -
More informationMicrosoft PowerPoint - H22制御工学I-2回.ppt
制御工学 I 第二回ラプラス変換 平成 年 4 月 9 日 /4/9 授業の予定 制御工学概論 ( 回 ) 制御技術は現在様々な工学分野において重要な基本技術となっている 工学における制御工学の位置づけと歴史について説明する さらに 制御システムの基本構成と種類を紹介する ラプラス変換 ( 回 ) 制御工学 特に古典制御ではラプラス変換が重要な役割を果たしている ラプラス変換と逆ラプラス変換の定義を紹介し
More informationスライド タイトルなし
線形代数 演習 (008 年度版 ) 008/5/6 線形代数 演習 Ⅰ コンピュータ グラフィックス, 次曲面と線形代数指南書第七の巻 直交行列, 実対称行列とその対角化, 次曲線池田勉龍谷大学理工学部数理情報学科 実行列, 正方行列, 実対称行列, 直交行列 a a N A am a MN 実行列 : すべての成分 a が実数である行列 ij ji ij 正方行列 : 行の数と列の数が等しい (
More informationMicrosoft Word - 第6章 H∞制御.doc
NAOSIE: Nagaaki Univerity' Ac itle 現代制御理論ノート Author() 辻, 峰男 Citation 現代制御理論ノート ; 7 Iue Date 7 URL http://hdl.handle.net/69/37695 Right hi document i downloaded http://naoite.lb.nagaaki-u.ac.jp 第 6 章 H
More information<4D F736F F D E4F8E9F82C982A882AF82E98D7397F1>
3 三次における行列 要旨高校では ほとんど 2 2 の正方行列しか扱ってなく 三次の正方行列について考えてみたかったため 数 C で学んだ定理を三次の正方行列に応用して 自分たちで仮説を立てて求めていったら 空間における回転移動を表す行列 三次のケーリー ハミルトンの定理 三次における逆行列を求めたり 仮説をたてることができた. 目的 数 C で学んだ定理を三次の正方行列に応用する 2. 概要目的の到達点として
More informationMicrosoft PowerPoint - 10.pptx
0. 固有値とその応用 固有値と固有ベクトル 2 行列による写像から固有ベクトルへ m n A : m n n m 行列によって線形写像 f R R A が表せることを見てきた ここでは 2 次元平面の行列による写像を調べる 2 = 2 A 2 2 とし 写像 まず 単位ベクトルの像を求める u 2 x = v 2 y f : R A R を考える u 2 2 u, 2 2 0 = = v 2 0
More informationMicrosoft PowerPoint - H21生物計算化学2.ppt
演算子の行列表現 > L いま 次元ベクトル空間の基底をケットと書くことにする この基底は完全系を成すとすると 空間内の任意のケットベクトルは > > > これより 一度基底を与えてしまえば 任意のベクトルはその基底についての成分で完全に記述することができる これらの成分を列行列の形に書くと M これをベクトル の基底 { >} による行列表現という ところで 行列 A の共役 dont 行列は A
More informationPowerPoint プレゼンテーション
回転型クレーン / 倒立振子の制御 回転型クレーンの制御 状態方程式 コントローラ設計 ( 極配置法 ) コントローラ設計 ( 最適レギュレータ ) 回転型倒立振子の制御 状態方程式 コントローラ設計 コントローラの形式 : 状態フィードバック P-D コントローラ アームの P-D 振子の P-D 目標値 状態フィードバック制御 回転型クレーン コントローラ で 状態フィードバック制御 回転型クレーン
More informationMicrosoft Word - 補論3.2
補論 3. 多変量 GARC モデル 07//6 新谷元嗣 藪友良 対数尤度関数 3 章 7 節では 変量の対数尤度を求めた ここでは多変量の場合 とくに 変量について対数尤度を求める 誤差項 は平均 0 で 次元の正規分布に従うとする 単純化のため 分散と共分散は時間を通じて一定としよう ( この仮定は後で変更される ) したがって ij から添え字 を除くことができる このとき と の尤度関数は
More informationMicrosoft PowerPoint - 9.pptx
9. 線形写像 ここでは 行列の積によって 写像を定義できることをみていく また 行列の積によって定義される写像の性質を調べていく 行列演算と写像 ( 次変換 3 拡大とスカラー倍 p ' = ( ', ' = ( k, kk p = (, k 倍 k 倍 拡大後 k 倍拡大の関係は スカラー倍を用いて次のように表現できる ' = k ' 拡大前 拡大 4 拡大と行列の積 p ' = ( ', '
More informationディジタル信号処理
ディジタルフィルタの設計法. 逆フィルター. 直線位相 FIR フィルタの設計. 窓関数法による FIR フィルタの設計.5 時間領域での FIR フィルタの設計 3. アナログフィルタを基にしたディジタル IIR フィルタの設計法 I 4. アナログフィルタを基にしたディジタル IIR フィルタの設計法 II 5. 双 次フィルタ LI 離散時間システムの基礎式の証明 [ ] 4. ] [ ]*
More informationMatrix and summation convention Kronecker delta δ ij 1 = 0 ( i = j) ( i j) permutation symbol e ijk = (even permutation) (odd permutation) (othe
Matr ad summato covto Krockr dlta δ ( ) ( ) prmutato symbol k (v prmutato) (odd prmutato) (othrs) gvalu dtrmat dt 6 k rst r s kt opyrght s rsrvd. No part of ths documt may b rproducd for proft. 行列 行 正方行列
More informationパソコンシミュレータの現状
第 2 章微分 偏微分, 写像 豊橋技術科学大学森謙一郎 2. 連続関数と微分 工学において物理現象を支配する方程式は微分方程式で表されていることが多く, 有限要素法も微分方程式を解く数値解析法であり, 定式化においては微分 積分が一般的に用いられており. 数学の基礎知識が必要になる. 図 2. に示すように, 微分は連続な関数 f() の傾きを求めることであり, 微小な に対して傾きを表し, を無限に
More informationMicrosoft PowerPoint - ce07-12c.ppt
制御工学 第 8 章 : フィードバック制御系の設計法 8. 設計手順と性能評価キーワード : 設計手順, 性能評価 8. 補償による制御系設計キーワード : ( 比例 ),( 積分 ),( 微分 ) 8.3 進み 遅れ補償による制御系設計キーワード : 遅れ補償, 進み補償 学習目標 : 一般的な制御系設計における手順と制御系の性能評価について学ぶ. ループ整形の考え方を用いて, 遅れ補償, 進み補償による制御系設計を習得する.
More informationPowerPoint プレゼンテーション
復習 ) 時系列のモデリング ~a. 離散時間モデル ~ y k + a 1 z 1 y k + + a na z n ay k = b 0 u k + b 1 z 1 u k + + b nb z n bu k y k = G z 1 u k = B(z 1 ) A(z 1 u k ) ARMA モデル A z 1 B z 1 = 1 + a 1 z 1 + + a na z n a = b 0
More informationMicrosoft PowerPoint - 9.pptx
9/7/8( 水 9. 線形写像 ここでは 行列の積によって 写像を定義できることをみていく また 行列の積によって定義される写像の性質を調べていく 拡大とスカラー倍 行列演算と写像 ( 次変換 拡大後 k 倍 k 倍 k 倍拡大の関係は スカラー倍を用いて次のように表現できる p = (, ' = k ' 拡大前 p ' = ( ', ' = ( k, k 拡大 4 拡大と行列の積 拡大後 k 倍
More information09.pptx
講義内容 数値解析 第 9 回 5 年 6 月 7 日 水 理学部物理学科情報理学コース. 非線形方程式の数値解法. はじめに. 分法. 補間法.4 ニュートン法.4. 多変数問題への応用.4. ニュートン法の収束性. 連立 次方程式の解法. 序論と行列計算の基礎. ガウスの消去法. 重対角行列の場合の解法項目を変更しました.4 LU 分解法.5 特異値分解法.6 共役勾配法.7 反復法.7. ヤコビ法.7.
More information数学 t t t t t 加法定理 t t t 倍角公式加法定理で α=β と置く. 三角関数
. 三角関数 基本関係 t cot c sc c cot sc t 還元公式 t t t t t t cot t cot t 数学 数学 t t t t t 加法定理 t t t 倍角公式加法定理で α=β と置く. 三角関数 数学. 三角関数 5 積和公式 6 和積公式 数学. 三角関数 7 合成 t V v t V v t V V V V VV V V V t V v v 8 べき乗 5 6 6
More informations とは何か 2011 年 2 月 5 日目次へ戻る 1 正弦波の微分 y=v m sin ωt を時間 t で微分します V m は正弦波の最大値です 合成関数の微分法を用い y=v m sin u u=ωt と置きますと dy dt dy du du dt d du V m sin u d dt
とは何か 0 年 月 5 日目次へ戻る 正弦波の微分 y= in を時間 で微分します は正弦波の最大値です 合成関数の微分法を用い y= in u u= と置きますと y y in u in u (co u co になります in u の は定数なので 微分後も残ります 合成関数の微分法ですので 最後に u を に戻しています 0[ra] の co 値は [ra] の in 値と同じです その先の角
More informationMicrosoft PowerPoint - aep_1.ppt [互換モード]
物理計測法特論 No.1 第 1 章 : 信号と雑音 本講義の主題 雑音の性質を理解することで 信号と雑音の大きさが非常に近い状態での信号の測定技術 : 微小信号計測 について学ぶ 講義の Web http://www.g-munu.t.u-tokyo.ac.jp/mio/note/sig_mes/tokuron.html 物理学の基本は実験事実の積み重ねである そして それは何かを測定することから始まる
More informationMicrosoft PowerPoint - 6.PID制御.pptx
プロセス制御工学 6.PID 制御 京都大学 加納学 Division of Process Control & Process Systems Engineering Department of Chemical Engineering, Kyoto University manabu@cheme.kyoto-u.ac.jp http://www-pse.cheme.kyoto-u.ac.jp/~kano/
More informationPowerPoint プレゼンテーション
データ解析 第 7 回 : 時系列分析 渡辺澄夫 過去から未来を予測する 観測データ 回帰 判別分析 解析方法 主成分 因子 クラスタ分析 時系列予測 時系列を予測する 無限個の確率変数 ( 確率変数が作る無限数列 ){X(t) ; t は整数 } を生成する情報源を考える {X(t)} を確率過程という 確率過程に ついて過去の値から未来を予測するにはどうしたらよいだろうか X(t-K),X(t-K+1),,X(t-1)
More informationPowerPoint プレゼンテーション
PID 制御の基礎 ON/OFF 制御 PID 制御 P 制御 過渡特性を改善しよう PD 制御と P-D 制御 定常特性を改善しよう PI-D 制御 4.2 節 I-PD 制御 角度制御実験装置 0 [deg] 30 [deg] 角度制御実験装置 目標値 コントローラ ( マイコン ) アクチュエータ (DC モータ ) 制御対象 ( アーム ) 角度 センサ ( ロータリエンコーダ ) ON/OFF
More information1. 線形シフト不変システムと z 変換 ここで言う システム とは? 入力数列 T[ ] 出力数列 一意変換 ( 演算子 ) 概念的には,, x 2, x 1, x 0, x 1, x 2, を入力すると, y 2, y 1, y 0, y 1, y 2, が出力される. 線形システム : 線形シ
1. 線形シフト不変システムと z 変換 ここで言う システム とは? 入力数列 T[ ] 出力数列 一意変換 ( 演算子 ) 概念的には,, x 2, x 1, x, x1, x2, を入力すると, y 2, y 1, y, y1, y2, が出力される. 線形システム : 線形システムの例 x nxn 1 yn= 2 線形でないシステムの例 xn yn={ 2 xn xn othewise なぜ線形システム?
More information2018年度 東京大・理系数学
08 東京大学 ( 理系 ) 前期日程問題 解答解説のページへ関数 f ( ) = + cos (0 < < ) の増減表をつくり, + 0, 0 のと sin きの極限を調べよ 08 東京大学 ( 理系 ) 前期日程問題 解答解説のページへ n+ 数列 a, a, を, Cn a n = ( n =,, ) で定める n! an qn () n とする を既約分数 an p として表したときの分母
More information景気指標の新しい動向
内閣府経済社会総合研究所 経済分析 22 年第 166 号 4 時系列因子分析モデル 4.1 時系列因子分析モデル (Stock-Watson モデル の理論的解説 4.1.1 景気循環の状態空間表現 Stock and Watson (1989,1991 は観測される景気指標を状態空間表現と呼ば れるモデルで表し, 景気の状態を示す指標を開発した. 状態空間表現とは, わ れわれの目に見える実際に観測される変数は,
More informationClipboard
機械システム工学実験 Ⅲ 現代制御実験 授業 ( 時間程度 シミュレーション ( 時間 実験 ( 時間 - 時間半 課題 ( 時間 - 時間半 レポート提出に関して日時 : 翌週の月曜 時 3 分場所 : 9 号館 553 室 質問があれば鎌田研究室 (9 号館 35 室 まで 制御とは? 対象とする物 ( またはシステム を自分の思うように操る 制御するためには何が必要か? コントローラ ( 制御器
More informationMicrosoft Word - Chap17
第 7 章化学反応に対する磁場効果における三重項機構 その 7.. 節の訂正 年 7 月 日. 節 章の9ページ の赤枠に記載した説明は間違いであった事に気付いた 以下に訂正する しかし.. 式は 結果的には正しいので安心して下さい 磁場 の存在下でのT 状態のハミルトニアン は ゼーマン項 と時間に依存するスピン-スピン相互作用の項 との和となる..=7.. g S = g S z = S z g
More information多次元レーザー分光で探る凝縮分子系の超高速動力学
波動方程式と量子力学 谷村吉隆 京都大学理学研究科化学専攻 http:theochem.kuchem.kyoto-u.ac.jp TA: 岩元佑樹 iwamoto.y@kuchem.kyoto-u.ac.jp ベクトルと行列の作法 A 列ベクトル c = c c 行ベクトル A = [ c c c ] 転置ベクトル T A = [ c c c ] AA 内積 c AA = [ c c c ] c =
More informationMicrosoft Word - thesis.doc
剛体の基礎理論 -. 剛体の基礎理論初めに本論文で大域的に使用する記号を定義する. 使用する記号トルク撃力力角運動量角速度姿勢対角化された慣性テンソル慣性テンソル運動量速度位置質量時間 J W f F P p .. 質点の並進運動 質点は位置 と速度 P を用いる. ニュートンの運動方程式 という状態を持つ. 但し ここでは速度ではなく運動量 F P F.... より質点の運動は既に明らかであり 質点の状態ベクトル
More information数学 ⅡB < 公理 > 公理を論拠に定義を用いて定理を証明する 1 大小関係の公理 順序 (a > b, a = b, a > b 1 つ成立 a > b, b > c a > c 成立 ) 順序と演算 (a > b a + c > b + c (a > b, c > 0 ac > bc) 2 図
数学 Ⅱ < 公理 > 公理を論拠に定義を用いて定理を証明する 大小関係の公理 順序 >, =, > つ成立 >, > > 成立 順序と演算 > + > + >, > > 図形の公理 平行線の性質 錯角 同位角 三角形の合同条件 三角形の合同相似 量の公理 角の大きさ 線分の長さ < 空間における座漂とベクトル > ベクトルの演算 和 差 実数倍については 文字の計算と同様 ベクトルの成分表示 平面ベクトル
More informationMicrosoft PowerPoint - ip02_01.ppt [互換モード]
空間周波数 周波数領域での処理 空間周波数 (spatial frquncy) とは 単位長さ当たりの正弦波状の濃淡変化の繰り返し回数を表したもの 正弦波 : y sin( t) 周期 : 周波数 : T f / T 角周波数 : f 画像処理 空間周波数 周波数領域での処理 波形が違うと 周波数も違う 画像処理 空間周波数 周波数領域での処理 画像処理 3 周波数領域での処理 周波数は一つしかない?-
More informationMicrosoft Word - 第9章 PID制御.doc
NAOSITE: Nagaak Unry' Ac Tl 自動制御の理論と応用 Auhr() 辻, 峰男 Can 自動制御の理論と応用 ; 5 Iu Da 5 URL h://hdl.handl.n/69/35886 Rgh Th dcumn dwnladd h://na.lb.nagaak-u.ac.j 第 9 章 PID 制御 これまで, どのような制御器を用いるかということはあまり触れなかったが,
More information横浜市環境科学研究所
周期時系列の統計解析 単回帰分析 io 8 年 3 日 周期時系列に季節調整を行わないで単回帰分析を適用すると, 回帰係数には周期成分の影響が加わる. ここでは, 周期時系列をコサイン関数モデルで近似し単回帰分析によりモデルの回帰係数を求め, 周期成分の影響を検討した. また, その結果を気温時系列に当てはめ, 課題等について考察した. 気温時系列とコサイン関数モデル第 報の結果を利用するので, その一部を再掲する.
More information<4D F736F F F696E74202D2091E6824F82538FCD8CEB82E88C9F8F6F814592F990B382CC8CB4979D82BB82CC82505F D E95848D8682CC90B69
第 章 誤り検出 訂正の原理 その ブロック符号とその復号 安達文幸 目次 誤り訂正符号化を用いる伝送系誤り検出符号誤り検出 訂正符号 7, ハミング符号, ハミング符号生成行列, パリティ検査行列の一般形符号の生成行列符号の生成行列とパリティ検査行列の関係符号の訂正能力符号多項式 安達 : コミュニケーション符号理論 安達 : コミュニケーション符号理論 誤り訂正符号化を用いる伝送系 伝送システム
More informationMicrosoft PowerPoint - 資料04 重回帰分析.ppt
04. 重回帰分析 京都大学 加納学 Division of Process Control & Process Sstems Engineering Department of Chemical Engineering, Koto Universit manabu@cheme.koto-u.ac.jp http://www-pse.cheme.koto-u.ac.jp/~kano/ Outline
More informationMicrosoft PowerPoint - CSA_B3_EX2.pptx
Computer Science A Hardware Design Excise 2 Handout V2.01 May 27 th.,2019 CSAHW Computer Science A, Meiji University CSA_B3_EX2.pptx 32 Slides Renji Mikami 1 CSAHW2 ハード演習内容 2.1 二次元空間でのベクトルの直交 2.2 Reserved
More informationMicrosoft Word - 第2章 ブロック線図.doc
NAOSIE: Nagaaki Univriy' Ac il ディジタル制御システム Auhor() 辻, 峰男 Ciaion ディジタル制御システム ; 06 Iu Da 06 URL hp://hdl.handl.n/0069/3686 Righ hi documn i downloadd hp://naoi.lb.nagaaki-u.ac.jp 第 章ブロック線図. インパルス列を用いた z
More informationDVIOUT-SS_Ma
第 章 微分方程式 ニュートンはリンゴが落ちるのを見て万有引力を発見した という有名な逸話があります 無重力の宇宙船の中ではリンゴは落ちないで静止していることを考えると 重力が働くと始め静止しているものが動き出して そのスピードはどんどん大きくなる つまり速度の変化が現れることがわかります 速度は一般に時間と共に変化します 速度の瞬間的変化の割合を加速度といい で定義しましょう 速度が変化する, つまり加速度がでなくなるためにはその原因があり
More informationMicrosoft PowerPoint - spe1_handout11.ppt
目次 信号処理工学 Ⅰ 第 回 : ディジタルフィルタ 電気通信大学電子工学専攻電子知能システム学講座 ディジタルフィルタ ( 復習 ) FIR フィルタの補足 IIR フィルタの設計 IIR フィルタの実現 FIR フィルタと IIR フィルタの比較 最後の課題 長井隆行 ディジタルフィルタ ( 復習 ) 線形位相 FIR フィルタの補足 FIR フィルタ フィードバックがない インパルス応答が有限
More informationMicrosoft Word ã‡»ã…«ã‡ªã…¼ã…‹ã…žã…‹ã…³ã†¨åłºæœ›å•¤(佒芤喋çfl�)
Cellulr uo nd heir eigenlues 東洋大学総合情報学部 佐藤忠一 Tdzu So Depren o Inorion Siene nd rs Toyo Uniersiy. まえがき 一次元セルオ-トマトンは数学的には記号列上の行列の固有値問題である 固有値問題の行列はふつう複素数体上の行列である 量子力学における固有値問題も無限次元ではあるが関数環上の行列でその成分は可換環である
More informations と Z(s) の関係 2019 年 3 月 22 日目次へ戻る s が虚軸を含む複素平面右半面の値の時 X(s) も虚軸を含む複素平面右半面の値でなけれ ばなりません その訳を探ります 本章では 受動回路をインピーダンス Z(s) にしていま す リアクタンス回路の駆動点リアクタンス X(s)
と Z の関係 9 年 3 月 日目次へ戻る が虚軸を含む複素平面右半面の値の時 X も虚軸を含む複素平面右半面の値でなけれ ばなりません その訳を探ります 本章では 受動回路をインピーダンス Z にしていま す リアクタンス回路の駆動点リアクタンス X も Z に含まれます Z に正弦波電流を入れた時最大値 抵抗 コイル コンデンサーで作られた受動回路の ラプラスの世界でのインピーダンスを Z とします
More information画像解析論(2) 講義内容
画像解析論 画像解析論 東京工業大学長橋宏 主な講義内容 信号処理と画像処理 二次元システムとその表現 二次元システムの特性解析 各種の画像フィルタ 信号処理と画像処理 画像解析論 処理の応答 記憶域 入出力の流れ 信号処理系 実時間性が求められる メモリ容量に対する制限が厳しい オンラインでの対応が厳しく求められる 画像処理系 ある程度の処理時間が許容される 大容量のメモリ使用が容認され易い オフラインでの対応が容認され易い
More informationベイズ統計入門
ベイズ統計入門 条件付確率 事象 F が起こったことが既知であるという条件の下で E が起こる確率を条件付確率 (codtoal probablt) という P ( E F ) P ( E F ) P( F ) 定義式を変形すると 確率の乗法公式となる ( E F ) P( F ) P( E F ) P( E) P( F E) P 事象の独立 ある事象の生起する確率が 他のある事象が生起するかどうかによって変化しないとき
More information数学の世界
東京女子大学文理学部数学の世界 (2002 年度 ) 永島孝 17 6 行列式の基本法則と効率的な計算法 基本法則 三次以上の行列式についても, 二次の場合と同様な法則がなりたつ ここには三次の場合を例示するが, 四次以上でも同様である 1 単位行列の行列式の値は 1 である すなわち 1 0 0 0 1 0 1 0 0 1 2 二つの列を入れ替えると行列式の値は 1 倍になる 例えば a 13 a
More information数値計算で学ぶ物理学 4 放物運動と惑星運動 地上のように下向きに重力がはたらいているような場においては 物体を投げると放物運動をする 一方 中心星のまわりの重力場中では 惑星は 円 だ円 放物線または双曲線を描きながら運動する ここでは 放物運動と惑星運動を 運動方程式を導出したうえで 数値シミュ
数値計算で学ぶ物理学 4 放物運動と惑星運動 地上のように下向きに重力がはたらいているような場においては 物体を投げると放物運動をする 一方 中心星のまわりの重力場中では 惑星は 円 だ円 放物線または双曲線を描きながら運動する ここでは 放物運動と惑星運動を 運動方程式を導出したうえで 数値シミュレーションによって計算してみる 4.1 放物運動一様な重力場における放物運動を考える 一般に質量の物体に作用する力をとすると運動方程式は
More information第6章 実験モード解析
第 6 章実験モード解析 6. 実験モード解析とは 6. 有限自由度系の実験モード解析 6.3 連続体の実験モード解析 6. 実験モード解析とは 実験モード解析とは加振実験によって測定された外力と応答を用いてモードパラメータ ( 固有振動数, モード減衰比, 正規固有モードなど ) を求める ( 同定する ) 方法である. 力計 試験体 変位計 / 加速度計 実験モード解析の概念 時間領域データを利用する方法
More informationChap2
逆三角関数の微分 Arcsin の導関数を計算する Arcsin I. 初等関数の微積分 sin [, ], [π/, π/] cos sin / (Arcsin ) 計算力の体力をつけよう π/ π/ E. II- 次の関数の導関数を計算せよ () Arccos () Arctan E. I- の解答 不定積分あれこれ () Arccos n log C C (n ) n e e C log (log
More informationmemo
数理情報工学特論第一 機械学習とデータマイニング 4 章 : 教師なし学習 3 かしまひさし 鹿島久嗣 ( 数理 6 研 ) kashima@mist.i.~ DEPARTMENT OF MATHEMATICAL INFORMATICS 1 グラフィカルモデルについて学びます グラフィカルモデル グラフィカルラッソ グラフィカルラッソの推定アルゴリズム 2 グラフィカルモデル 3 教師なし学習の主要タスクは
More information学習指導要領
(1) 数と式 ア数と集合 ( ア ) 実数数を実数まで拡張する意義を理解し 簡単な無理数の四則計算をすること 絶対値の意味を理解し適切な処理することができる 例題 1-3 の絶対値をはずせ 展開公式 ( a + b ) ( a - b ) = a 2 - b 2 を利用して根号を含む分数の分母を有理化することができる 例題 5 5 + 2 の分母を有理化せよ 実数の整数部分と小数部分の表し方を理解している
More information1/30 平成 29 年 3 月 24 日 ( 金 ) 午前 11 時 25 分第三章フェルミ量子場 : スピノール場 ( 次元あり ) 第三章フェルミ量子場 : スピノール場 フェルミ型 ボーズ量子場のエネルギーは 第二章ボーズ量子場 : スカラー場 の (2.18) より ˆ dp 1 1 =
/ 平成 9 年 月 日 ( 金 午前 時 5 分第三章フェルミ量子場 : スピノール場 ( 次元あり 第三章フェルミ量子場 : スピノール場 フェルミ型 ボーズ量子場のエネルギーは 第二章ボーズ量子場 : スカラー場 の (.8 より ˆ ( ( ( q -, ( ( c ( H c c ë é ù û - Ü + c ( ( - に限る (. である 一方 フェルミ型は 成分をもち その成分を,,,,
More information4.6: 3 sin 5 sin θ θ t θ 2t θ 4t : sin ωt ω sin θ θ ωt sin ωt 1 ω ω [rad/sec] 1 [sec] ω[rad] [rad/sec] 5.3 ω [rad/sec] 5.7: 2t 4t sin 2t sin 4t
1 1.1 sin 2π [rad] 3 ft 3 sin 2t π 4 3.1 2 1.1: sin θ 2.2 sin θ ft t t [sec] t sin 2t π 4 [rad] sin 3.1 3 sin θ θ t θ 2t π 4 3.2 3.1 3.4 3.4: 2.2: sin θ θ θ [rad] 2.3 0 [rad] 4 sin θ sin 2t π 4 sin 1 1
More information(Microsoft Word - 10ta320a_\220U\223\256\212w\223\301\230__6\217\315\221O\224\274\203\214\203W\203\201.docx)
6 章スペクトルの平滑化 スペクトルの平滑化とはフーリエスペクトルやパワ スペクトルのギザギザを取り除き 滑らかにする操作のことをいう ただし 波のもっている本質的なものをゆがめてはいけない 図 6-7 パワ スペクトルの平滑化 6. 合積のフーリエ変換スペクトルの平滑化を学ぶ前に 合積とそのフーリエ変換について説明する 6. データ ウィンドウデータ ウィンドウの定義と特徴について説明する 6.3
More informationファイナンスのための数学基礎 第1回 オリエンテーション、ベクトル
時系列分析 変量時系列モデルとその性質 担当 : 長倉大輔 ( ながくらだいすけ 時系列モデル 時系列モデルとは時系列データを生み出すメカニズムとなるものである これは実際には未知である 私たちにできるのは観測された時系列データからその背後にある時系列モデルを推測 推定するだけである 以下ではいくつかの代表的な時系列モデルを考察する 自己回帰モデル (Auoregressive Model もっとも頻繁に使われる時系列モデルは自己回帰モデル
More informationアルゴリズムとデータ構造
講義 アルゴリズムとデータ構造 第 2 回アルゴリズムと計算量 大学院情報科学研究科情報理工学専攻情報知識ネットワーク研究室喜田拓也 講義資料 2018/5/23 今日の内容 アルゴリズムの計算量とは? 漸近的計算量オーダーの計算の方法最悪計算量と平均計算量 ポイント オーダー記法 ビッグオー (O), ビッグオメガ (Ω), ビッグシータ (Θ) 2 お風呂スケジューリング問題 お風呂に入る順番を決めよう!
More informationMicrosoft PowerPoint - 2.ppt [互換モード]
0 章数学基礎 1 大学では 高校より厳密に議論を行う そのために 議論の議論の対象を明確にする必要がある 集合 ( 定義 ) 集合 物の集まりである集合 X に対して X を構成している物を X の要素または元という 集合については 3 セメスタ開講の 離散数学 で詳しく扱う 2 集合の表現 1. 要素を明示する表現 ( 外延的表現 ) 中括弧で 囲う X = {0,1, 2,3} 慣用的に 英大文字を用いる
More information学習指導要領
(1) 数と式 学習指導要領 数と式 (1) 式の計算二次の乗法公式及び因数分解の公式の理解を深め 式を多面的にみたり目的に応じて式を適切に変形したりすること 東京都立町田高等学校学力スタンダード 整式の加法 減法 乗法展開の公式を利用できる 式を1 つの文字におき換えることによって, 式の計算を簡略化することができる 式の形の特徴に着目して変形し, 展開の公式が適用できるようにすることができる 因数分解因数分解の公式を利用できる
More information2015-2017年度 2次数学セレクション(複素数)解答解説
05 次数学セレクション解答解説 [ 筑波大 ] ( + より, 0 となり, + から, ( (,, よって, の描く図形 C は, 点 を中心とし半径が の円である すなわち, 原 点を通る円となる ( は虚数, は正の実数より, である さて, w ( ( とおくと, ( ( ( w ( ( ( ここで, w は純虚数より, は純虚数となる すると, の描く図形 L は, 点 を通り, 点 と点
More information<4D F736F F D20824F B CC92E8979D814696CA90CF95AA82C691CC90CF95AA2E646F63>
1/1 平成 23 年 3 月 24 日午後 6 時 52 分 6 ガウスの定理 : 面積分と体積分 6 ガウスの定理 : 面積分と体積分 Ⅰ. 直交座標系 ガウスの定理は 微分して すぐに積分すると元に戻るというルールを 3 次元積分に適用した定理になります よく知っているのは 簡単化のため 変数が1つの場合は dj ( d ( ににします全微分 = 偏微分 d = d = J ( + C d です
More informationデータ解析
データ解析 ( 前期 ) 最小二乗法 向井厚志 005 年度テキスト 0 データ解析 - 最小二乗法 - 目次 第 回 Σ の計算 第 回ヒストグラム 第 3 回平均と標準偏差 6 第 回誤差の伝播 8 第 5 回正規分布 0 第 6 回最尤性原理 第 7 回正規分布の 分布の幅 第 8 回最小二乗法 6 第 9 回最小二乗法の練習 8 第 0 回最小二乗法の推定誤差 0 第 回推定誤差の計算 第
More informationMicrosoft Word docx
有限図形の代数的表現について 三角形や星型を式で表現したいという思いから以下のことを 考察をしまし た 有限個の点と辺で 構成される図形を 関数で表現する そのため 基礎 体として 素数の有限体を考える 但し 扱うのは 点の数と辺の数が等しい 特別場合である 先ず P5 のときから 始めることにします. グラフと写像と関数について ( 特別な場合 ) 集合 F {,,,, } について 写像 f :
More informationMicrosoft Word - ModernControl_Scilab.doc
Scilab による現代制御の学習 2014 年 4 月 7 日 目次 1. はじめに 2. 行列の演算 3. 時間応答の数値計算 4. システムの特性解析 5. 状態フィードバック制御の設計 6. オブザーバの設計と動的制御器 7. サーボ系の設計 8. 最適制御系の設計 9. LQG 理論によるサーボ系の設計 1. はじめに現代制御理論では, 状態方程式に基づいて制御系設計が行われる. 現代制御理論で用いられる計算としては,
More informationに対して 例 2: に対して 逆行列は常に存在するとは限らない 逆行列が存在する行列を正則行列 (regular matrix) という 正則である 逆行列が存在する 一般に 正則行列 A の逆行列 A -1 も正則であり (A -1 ) -1 =A が成り立つ また 2 つの正則行列 A B の積
2 逆行列 逆行列の計算は 連立一次方程式を数値的に解くために利用される 気象学の分野では線形系の応答問題を数値的に解くときに用いられることも多い ここでは計算機を用いて逆行列を求める方法を学ぶ 2.1 はじめにたとえば 次のような連立一次方程式を解くことを考える このような 2 元連立一次方程式は 代入法や消去法によって容易に解くことができる 解法をプログラミング言語によって記述することも困難ではない
More information重要例題113
04_ 高校 数学 Ⅱ 必須基本公式 定理集 数学 Ⅱ 第 章式の計算と方程式 0 商と余り についての整式 A をについての整式 B で割ったときの商を Q, 余りを R とすると, ABQ+R (R の次数 ) > 0
More information以下 変数の上のドットは時間に関する微分を表わしている (ex. 2 dx d x x, x 2 dt dt ) 付録 E 非線形微分方程式の平衡点の安定性解析 E-1) 非線形方程式の線形近似特に言及してこなかったが これまでは線形微分方程式 ( x や x, x などがすべて 1 次で なおかつ
以下 変数の上のドットは時間に関する微分を表わしている (e. d d, dt dt ) 付録 E 非線形微分方程式の平衡点の安定性解析 E-) 非線形方程式の線形近似特に言及してこなかったが これまでは線形微分方程式 ( や, などがすべて 次で なおかつそれらの係数が定数であるような微分方程式 ) に対して安定性の解析を行ってきた しかしながら 実際には非線形の微分方程式で記述される現象も多く存在する
More information( 最初の等号は,N =0, 番目は,j= のとき j =0 による ) j>r のときは p =0 から和の上限は r で十分 定義 命題 3 ⑵ 実数 ( 0) に対して, ⑴ =[] []=( 0 または ) =[6]+[] [4] [3] [] =( 0 または ) 実数 に対して, π()
伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊 伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊 数研通信 70 号を読んで チェビシェフの定理の精密化 と.5 の間に素数がある 伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊 さい才 の 野 せ瀬 いちろう 一郎 伊伊伊伊伊伊伊伊伊伊伊伊伊伊伊 0. はじめに このたび,
More informationMicrosoft Word - Chap11
第 章 次元回転群とそのリー代数. SO のリー代数. 節でリー代数を定義したが 以下にその定義を再録する なお 多くの教科書に従って本章以降は ep t A の代わりに ep t と書くこととする 定義.. G を 次の線型リー群とすると 任意の実数 t に対して ep t G となる gl C の全体をGのリー代数 またはリー環 という 例えば ep t が 次の特殊直交群 SO の元であれば
More information学習指導要領
(1) いろいろな式 学習指導要領紅葉川高校学力スタンダードア式と証明展開の公式を用いて 3 乗に関わる式を展開すること ( ア ) 整式の乗法 除法 分数式の計算ができるようにする 三次の乗法公式及び因数分解の公式を理解し そ 3 次の因数分解の公式を理解し それらを用いて因数れらを用いて式の展開や因数分解をすること また 分解することができるようにする 整式の除法や分数式の四則計算について理解し
More informationMicrosoft PowerPoint - H17-5時限(パターン認識).ppt
パターン認識早稲田大学講義 平成 7 年度 独 産業技術総合研究所栗田多喜夫 赤穂昭太郎 統計的特徴抽出 パターン認識過程 特徴抽出 認識対象から何らかの特徴量を計測 抽出 する必要がある 認識に有効な情報 特徴 を抽出し 次元を縮小した効率の良い空間を構成する過程 文字認識 : スキャナ等で取り込んだ画像から文字の識別に必要な本質的な特徴のみを抽出 例 文字線の傾き 曲率 面積など 識別 与えられた未知の対象を
More information学力スタンダード(様式1)
(1) 数と式 学習指導要領ア数と集合 ( ア ) 実数数を実数まで拡張する意義を理解し 簡単な無理数の四則計算をすること 稔ヶ丘高校学力スタンダード 有理数 無理数の定義や実数の分類について理解し ている 絶対値の意味と記号表示を理解している 実数と直線上の点が一対一対応であることを理解 し 実数を数直線上に示すことができる 例 実数 (1) -.5 () π (3) 数直線上の点はどれか答えよ
More informationProbit , Mixed logit
Probit, Mixed logit 2016/5/16 スタートアップゼミ #5 B4 後藤祥孝 1 0. 目次 Probit モデルについて 1. モデル概要 2. 定式化と理解 3. 推定 Mixed logit モデルについて 4. モデル概要 5. 定式化と理解 6. 推定 2 1.Probit 概要 プロビットモデルとは. 効用関数の誤差項に多変量正規分布を仮定したもの. 誤差項には様々な要因が存在するため,
More informationMicrosoft Word - 微分入門.doc
基本公式 例題 0 定義式 f( ) 数 Ⅲ 微分入門 = の導関数を定義式にもとづいて計算しなさい 基本事項 ( f( ), g( ) が微分可能ならば ) y= f( ) g( ) のとき, y = y= f( ) g( ) h( ) のとき, y = ( f( ), g( ) が微分可能で, g( ) 0 ならば ) f( ) y = のとき, y = g ( ) とくに, y = のとき,
More informationMicrosoft Word - 実験テキスト2005.doc
7. プロセスの動特性 [Ⅰ] 目的液レベル制御実験および同シミュレーションを通して ステップ応答に基づくプロセス伝達関数の同定方法 ステップ応答法による PI 制御パラメータの調整方法 および PI 制御パラメータが制御性能へ与える影響について習熟する さらに 制御シミュレーションを通して むだ時間を有するプロセスに対するスミス補償型制御の有効性を確認する [Ⅱ] 理論 2.1 ステップ応答実験による伝達関数の同定
More information学習指導要領
(1) 数と式 ア数と集合 ( ア ) 実数数を実数まで拡張する意義を理解し 簡単な無理数の四則計算をすること 自然数 整数 有理数 無理数の包含関係など 実数 の構成を理解する ( 例 ) 次の空欄に適当な言葉をいれて, 数の集合を表しなさい ア イ 無理数 整数 ウ 無理数の加法及び減法 乗法公式などを利用した計 算ができる また 分母だけが二項である無理数の 分母の有理化ができる ( 例 1)
More information学習指導要領
(1) 数と式 学習指導要領ア数と集合 ( ア ) 実数数を実数まで拡張する意義を理解し 簡単な無理数の四則計算をすること 千早高校学力スタンダード 自然数 整数 有理数 無理数の用語の意味を理解す る ( 例 ) 次の数の中から自然数 整数 有理 数 無理数に分類せよ 3 3,, 0.7, 3,,-, 4 (1) 自然数 () 整数 (3) 有理数 (4) 無理数 自然数 整数 有理数 無理数の包含関係など
More informationPowerPoint Presentation
応用数学 Ⅱ (7) 7 連立微分方程式の立て方と解法. 高階微分方程式による解法. ベクトル微分方程式による解法 3. 演算子による解法 連立微分方程式 未知数が複数個あり, 未知数の数だけ微分方程式が与えられている場合, これらを連立微分方程式という. d d 解法 () 高階微分方程式化による解法 つの方程式から つの未知数を消去して, 未知数が つの方程式に変換 のみの方程式にするために,
More informationIntroduction to System Identification
y(t) モデルベースデザイン 制御系設計のためのシステム同定入門 s 2 Teja Muppirala t s 2 3s 4 2012 The MathWorks, Inc. 1 モデルベースデザイン 正確なモデルがあることが大前提 実行可能な仕様書 シミュレーションによる設計 モデル 連続したテスト 検証 コード生成による実装 2 動的システムのモデリング モデリング手法 第一原理モデリング データドリブンモデリング
More information学習指導要領
(1) 数と式 ア数と集合 ( ア ) 実数数を実数まで拡張する意義を理解し 簡単な無理数の四則計算をすること 自然数 整数 有理数 無理数 実数のそれぞれの集 合について 四則演算の可能性について判断できる ( 例 ) 下の表において それぞれの数の範囲で四則計算を考えるとき 計算がその範囲で常にできる場合には を 常にできるとは限らない場合には を付けよ ただし 除法では 0 で割ることは考えない
More informationMicrosoft Word - 町田・全 H30学力スタ 別紙1 1年 数学Ⅰ.doc
(1) 数と式 学習指導要領 都立町田高校 学力スタンダード ア 数と集合 ( ア ) 実数 根号を含む式の計算 数を実数まで拡張する意義を理解し 簡単な 循環小数を表す記号を用いて, 分数を循環小数で表 無理数の四則計算をすること すことができる 今まで学習してきた数の体系について整理し, 考察 しようとする 絶対値の意味と記号表示を理解している 根号を含む式の加法, 減法, 乗法の計算ができる
More informationMicrosoft PowerPoint - 物情数学C(2012)(フーリエ前半)_up
年度物理情報工学科 年生秋学期 物理情報数学 C フーリエ解析 (Fourier lysis) 年 月 5 日 フーリエ ( フランス ) (768~83: ナポレオンの時代 ) 歳で Ecole Polyechique ( フランス国立理工科大学 ) の教授 ナポレオンのエジプト遠征に従軍 (798) 87: 任意の関数は三角関数によって級数展開できる という フーリエ級数 の概念を提唱 ( 論文を提出
More information学習指導要領
(1 ) 数と式 ア数と集合 ( ア ) 実数数を実数まで拡張する意義を理解し 簡単な無理数の四則計算をすること 自然数 整数 有理数 無理数の包含関係など 実 数の構成を理解する ( 例 ) 次の空欄に適当な言葉をいれて, 数の集合を表しなさい 実数の絶対値が実数と対応する点と原点との距離で あることを理解する ( 例 ) 次の値を求めよ (1) () 6 置き換えなどを利用して 三項の無理数の乗法の計
More informationフィードバック ~ 様々な電子回路の性質 ~ 実験 (1) 目的実験 (1) では 非反転増幅器の増幅率や位相差が 回路を構成する抵抗値や入力信号の周波数によってどのように変わるのかを調べる 実験方法 図 1 のような自由振動回路を組み オペアンプの + 入力端子を接地したときの出力電圧 が 0 と
フィードバック ~ 様々な電子回路の性質 ~ 実験 (1) 目的実験 (1) では 非反転増幅器の増幅率や位相差が 回路を構成する抵抗値や入力信号の周波数によってどのように変わるのかを調べる 実験方法 図 1 のような自由振動回路を組み オペアンプの + 入力端子を接地したときの出力電圧 が 0 となるように半固定抵抗器を調整する ( ゼロ点調整のため ) 図 1 非反転増幅器 2010 年度版物理工学実験法
More informationスライド 1
作成 : 群馬大学電気電子教員 電子回路設計 OP アンプ (2) 小林春夫 桑名杏奈 Email: koba@gunma-u.ac.jp Tel: 277-3-788 オフィスアワー : AM9:~AM:( 平日 ) 電気電子棟 (3 号館 )4F 44 室 電子回路設計 授業の内容 第 回講義内容の説明と電子回路設計の基礎知識 第 2 回キルヒホッフ則を用いた回路解析と演習 第 3 回集積回路のデバイス
More information線型代数試験前最後の 3 日間 できるようになっておきたい計算問題 ( 特に注意 まぁ注意 ) シュミットの直交化とその行列表示 (P5) ユニタリ行列による行列の対角化 (P8) 数列, 微分方程式の解法 対角可能な条件もおさえておきたい とりあえず次の問題を ( まだやっていない人は ) やって
線型代数試験前最後の 日間 できるようになっておきたい計算問題 特に注意 まぁ注意 シュミットの直交化とその行列表示 P ユニタリ行列による行列の対角化 P8 数列 微分方程式の解法 対角可能な条件もおさえておきたい とりあえず次の問題を まだやっていない人は やってください 8 年 月 日 理二三 組の線型代数担当志甫先生の過去問から持ってきました 結構計算が大変だったと思います これが難なくできる人は以下の総復習編はさらっと目を通すだけで
More informationシステム工学実験 パラメータ推定手順
システム工学実験パラメータ推定手順 大木健太郎 2014/11/14 2014 年度システム工学実験 : フレキシブルリンク 1 アウトライン 1. 線形システムと周波数情報 2. パラメータ推定 3. 実際の手順 2014/11/14 2014 年度システム工学実験 : フレキシブルリンク 2 線形時不変システムと伝達関数 入力と出力の関係が線形な定係数微分方程式で与えられるとき, この方程式を線形時不変システムという
More information学習指導要領
(1) 数と式 学習指導要領ア数と集合 ( ア ) 実数数を実数まで拡張する意義を理解し 簡単な無理数の四則計算をすること 第 1 章第 節実数 東高校学力スタンダード 4 実数 (P.3~7) 自然数 整数 有理数 無理数 実数のそれぞれの集 合について 四則演算の可能性について判断できる ( 例 ) 下の表において, それぞれの数の範囲で四則計算を考えるとき, 計算がその範囲で常にできる場合には
More informationMicrosoft Word - Time Series Basic - Modeling.doc
時系列解析入門 モデリング. 確率分布と統計的モデル が確率変数 (radom varable のとき すべての実数 R に対して となる確 率 Prob( が定められる これを の関数とみなして G( Prob ( とあらわすとき G( を確率変数 の分布関数 (probablt dstrbuto ucto と呼 ぶ 時系列解析で用いられる確率変数は通常連続型と呼ばれるもので その分布関数は (
More information学習指導要領
(1) 数と式 ア整式 ( ア ) 式の展開と因数分解二次の乗法公式及び因数分解の公式の理解を深め 式を多面的にみたり目的に応じて式を適切に変形したりすること (ax b)(cx d) acx (ad bc)x bd などの基本的な公式を活用して 二次式の展開や因数分解ができる また 式の置き換えや一文字に着目するなどして 展開 因数分解ができる ( 例 ) 次の問に答えよ (1) (3x a)(4x
More information