2電力系統の数理モデル(a) ∆ωi(b) 実線：δi ，破線：∠V i(c) 実線：Ei ，破線：|V i |(d) 実線：Pi ，破線：Qi

図 2.10 初期値に摂動を加えた場合の時間応答（青：母線 1，黒：母線 2，赤：母線 3）© 公益社団法人 計測自動制御学会 2022，「電力系統のシステム制御工学」（2022年，コロナ社発行書籍関連資料）

2. 電力系統の数理モデル(a) ∆ωi(b) 実線：δi ，破線：∠V i(c) 実線：Ei ，破線：|V i |(d) 実線：Pi ，破線：Qi3

図 2.11 機械入力に摂動を加えた場合の時間応答（青：母線 1，黒：母線 2，赤：母線 3）© 公益社団法人 計測自動制御学会 2022，「電力系統のシステム制御工学」（2022年，コロナ社発行書籍関連資料）

4

2. 電力系統の数理モデル(a) ∆ωi(b) 実線：δi ，破線：∠V i(c) 実線：Ei ，破線：|V i |(d) 実線：Pi ，破線：Qi

図 2.12 界磁入力に摂動を加えた場合の時間応答（青：母線 1，黒：母線 2，赤：母線 3）© 公益社団法人 計測自動制御学会 2022，「電力系統のシステム制御工学」（2022年，コロナ社発行書籍関連資料）

3
電力系統モデルの数値シミュレーション(a) 表 3.1 の定常潮流状態
図 3.2(b) 表 3.2 の定常潮流状態

初期値変動に対する角周波数偏差の時間応答（青実線：∆ω1 ，赤破線：∆ω3 ）© 公益社団法人 計測自動制御学会 2022，「電力系統のシステム制御工学」（2022年，コロナ社発行書籍関連資料）

6

3. 電力系統モデルの数値シミュレーション(a) 抵抗が 1%増加した場合(b) 抵抗が 1%減少した場合

図 3.3 負荷の変化に対する角周波数偏差の時間応答（表 3.1 の定常潮流状態，線種は図 3.2 と同様）

(a) 抵抗が 1%増加した場合(b) 抵抗が 1%減少した場合

図 3.4 負荷の変化に対する角周波数偏差の時間応答(表 3.2 の定常潮流状態，線種は図 3.2 と同様）

© 公益社団法人 計測自動制御学会 2022，「電力系統のシステム制御工学」（2022年，コロナ社発行書籍関連資料）

3. 電力系統モデルの数値シミュレーション(a) 50 [ms] の地絡

7(b) 100 [ms] の地絡

図 3.5 地絡に対する角周波数偏差の時間応答（表 3.1 の定常潮流状態，線種は図 3.2 と同様）

(a) 50 [ms] の地絡(b) 100 [ms] の地絡図 3.6 地絡に対する角周波数偏差の時間応答（表 3.2 の定常潮流状態，線種は図 3.2 と同様）

© 公益社団法人 計測自動制御学会 2022，「電力系統のシステム制御工学」（2022年，コロナ社発行書籍関連資料）

83. 電力系統モデルの数値シミュレーション(a) iL の時間応答(b) vC の時間応答

図 3.12 LC 並列回路の時間応答（青実線：ode45，赤破線：ode15s）© 公益社団法人 計測自動制御学会 2022，「電力系統のシステム制御工学」（2022年，コロナ社発行書籍関連資料）

4
電力系統モデルの定態安定性解析(a) ∆ω lin(b) δ lin(c) E lin(d) P lin

図 4.2 近似線形モデルの初期値応答（青：発電機 1，黒：発電機 2，赤：発電機 3）© 公益社団法人 計測自動制御学会 2022，「電力系統のシステム制御工学」（2022年，コロナ社発行書籍関連資料）

10

4. 電力系統モデルの定態安定性解析(a) D = (10, 10, 10)，Y = Y 0(c) D = (10, 10, 10)，Y =
図 4.3

Y0100(b) D = (1 1 1, )，Y = Y 010 10 10111(d) D = ( 10, 10, 10)，Y =Y0100

近似線形モデルが安定となるパラメータの領域(a) 受動送電条件 (ii) を満たす場合(b) 受動送電条件 (ii) を満たさない場合

図 4.5 例 4.2 に対する蓄積関数の時間変化（青：WF ，赤：WG ，黒：WF + WG ）© 公益社団法人 計測自動制御学会 2022，「電力系統のシステム制御工学」（2022年，コロナ社発行書籍関連資料）

4. 電力系統モデルの定態安定性解析(a)G(jω) + GT (−jω)の固有値2(b)11H(jω) − H T (−jω)の固有値
2j

図 4.9 G(s) の正実性と H(s) の負虚性
（青：受動送電条件 (ii) を満たす，赤：満たさない）© 公益社団法人 計測自動制御学会 2022，「電力系統のシステム制御工学」（2022年，コロナ社発行書籍関連資料）

124. 電力系統モデルの定態安定性解析(a) 実線：∆ω lin ， 破線：∆ω̂ lin(b) 実線：δ lin ，破線：δ̂ lin(c) 実線：E lin ，破線：Ê lin(d) 実線：P lin ，破線：P̂ lin

図 4.10 低次元近似を適用した場合の時間応答（青：発電機 1，黒：発電機 2，赤：発電機 3）© 公益社団法人 計測自動制御学会 2022，「電力系統のシステム制御工学」（2022年，コロナ社発行書籍関連資料）

4. 電力系統モデルの定態安定性解析(a) D = (10, 10, 10)，Y = Y 0(c) D = (10, 10, 10)，Y =
図 4.11

Y0100

13(b) D = (0.1, 0.1, 0.1)，Y = Y 0(d) D = (0.1, 0.1, 0.1)，Y =Y0100

近似線形モデルが安定となるパラメータの領域© 公益社団法人 計測自動制御学会 2022，「電力系統のシステム制御工学」（2022年，コロナ社発行書籍関連資料）

5電力系統モデルの安定化制御(a) 自動発電制御なし
図 5.2(a) 青実線：P1 ，赤破線：P3
図 5.3(b) 自動発電制御あり

消費電力増加に対する角周波数偏差の時間応答（青実線：∆ω1 ，赤破線：∆ω3 ）(b) 送電損失 P1 + P2 + P3

寄与係数の変化に対する有効電力の時間応答© 公益社団法人 計測自動制御学会 2022，「電力系統のシステム制御工学」（2022年，コロナ社発行書籍関連資料）

5. 電力系統モデルの安定化制御(a) 定常潮流状態 1 に対応する初期値

15(b) 定常潮流状態 2 に対応する初期値

図 5.5 初期値応答に対する蓄積関数の時間変化（青：Wx⋆ ，赤：Wx⋆ ，緑：Wξ⋆ ，黒：総和）F(a) ∥∆ω∥L2
図 5.8G(b) |V | − |V ⋆ | L

2

回転子偏角差の初期値に対する過渡安定度評価（青実線：(i)，黒破線：(ii)，赤鎖線：(iii)）© 公益社団法人 計測自動制御学会 2022，「電力系統のシステム制御工学」（2022年，コロナ社発行書籍関連資料）

16

5. 電力系統モデルの安定化制御(a) 自動電圧調整器なし(b) 自動電圧調整器あり

図 5.9 角周波数偏差の初期値応答（青実線：∆ω1 ，赤破線：∆ω3 ）(a) 自動電圧調整器なし(b) 自動電圧調整器あり

図 5.10 角周波数偏差の初期値応答（青実線：∆ω1 ，赤破線：∆ω3 ）© 公益社団法人 計測自動制御学会 2022，「電力系統のシステム制御工学」（2022年，コロナ社発行書籍関連資料）

5. 電力系統モデルの安定化制御(a) δ1 (0) − δ3 (0) = −117(b) δ1 (0) − δ3 (0) = −1.57

図 5.12 角周波数偏差の初期値応答（系統安定化装置あり，線種は図 5.9 と同様）

(a) δ3 と δ̂3 の時系列データ(b) E3 と Ê3 の時系列データ

図 5.15 無作為な励振入力に対する時間応答（青実線：δ3 , E3 ，赤破線：δ̂3 , Ê3 ）(a) Q(δ3 , δ̂3 )
図 5.16(b) Q(E3 , Ê3 )同定パラメータに対する内部状態の近似誤差© 公益社団法人 計測自動制御学会 2022，「電力系統のシステム制御工学」（2022年，コロナ社発行書籍関連資料）

6
大規模モデルの数値シミュレーション例(a) 機械入力と界磁入力が定数の場合(b) 機械入力が定数の場合(c) 機械入力を制御する場合(d) 機械入力を制御する場合（拡大）図 6.2

負荷変動に対する角周波数偏差の時間応答© 公益社団法人 計測自動制御学会 2022，「電力系統のシステム制御工学」（2022年，コロナ社発行書籍関連資料）

6. 大規模モデルの数値シミュレーション例(a) Ei ejδi [pu](b) V i [pu](c) I i [pu](d) P + jQ [pu]

図 6.3

19

負荷変動に対する発電機変数と母線変数の変化（丸：初期時刻，四角：終端時刻）© 公益社団法人 計測自動制御学会 2022，「電力系統のシステム制御工学」（2022年，コロナ社発行書籍関連資料）