(57)【特許請求の範囲】
【請求項１】
太陽光発電用パワーコンディショナ（以下、「ＰＣＳ」という。）の単独運転検出に用いた周波数フィードバック方式によるハンチング現象を解析するために、周波数フィードバック方式の伝達関数と、ＰＣＳが接続される電力系統の伝達関数とに基づく下記の数式（１）の現状解析用の開ループ伝達関数（G1(s)）を、その変数（ｓ）に基づいて解析装置に演算させてグラフの中にナイキスト軌跡を描かせるものであって、
周波数フィードバック方式の伝達関数は、基準周波数の偏差（Δf0）に対する電力系統の周波数偏差（Δfc）の和である第１の周波数偏差の検出に要するむだ時間（τ）の特性を示すむだ時間要素と、単独運転系統の周波数を検出した２つの時期の時間差の特性をその時定数（Tr）に基づいて示すハイパスフィルタ要素と、２つの時期のうち最新時期における周波数の検出時間の特性をその時定数（Td）に基づいて示すローパスフィルタ要素と、２つの時期の周波数の差である第２の周波数偏差に応じてＰＣＳが出力する無効電力のゲイン（Ｇ）の特性を示すゲイン要素との積で表現し、
電力系統の伝達関数は、負荷およびＰＣＳの双方の有効電力特性指数に基づく合成有効電力特性指数（α）と、電力系統の線路インピーダンス（Ｒ,Ｘ）と、ＰＣＳから出力された無効電力とによって、影響を受けた電力系統の周波数偏差（Δfc）を位相変化の時間微分で表現することを特徴とするＰＣＳの周波数フィードバック方式のハンチング現象の解析方法。
【数１】
0006466873.tif 000016
但しω0＝２πf0、R:線路インピーダンスの抵抗、X：線路インピーダンスのリアクタンスとする。
【請求項２】
前記現状解析用の開ループ伝達関数G1(s)と、時定数Thに基づいて一次遅れの特性を示す一次遅れ要素との積に基づく下記数式（２）で示すハンチング対策用の伝達関数G2(s)を、その変数ｓと一次遅れの時定数Thとに基づいて解析装置に演算させてグラフの中にナイキスト軌跡を描かせることを特徴とする請求項１に記載のＰＣＳの周波数フィードバック方式のハンチング現象の回避方法。
【数２】
0006466873.tif 000017
【請求項３】
太陽光発電用パワーコンディショナ（以下、「ＰＣＳ」という。）の単独運転検出に用いた周波数フィードバック方式によるハンチング現象を解析するために、周波数フィードバック方式の伝達関数と、ＰＣＳが接続される電力系統の伝達関数とに基づく下記の数式（３）の現状解析用の開ループ伝達関数G1(s)を簡略化した下記の数式（４）の安定限界関数AH（R,X）を、その変数Ｒ、Ｘに基づいて解析装置に演算させてグラフの中に安定限界曲線を描かせると共に、線路インピーダンスの抵抗ＲとリアクタンスＸとに基づく電力系統をグラフの中に描かせるものであって、
周波数フィードバック方式の伝達関数は、基準周波数の偏差（Δf0）に対する電力系統の周波数偏差（Δfc）の和である第１の周波数偏差の検出に要するむだ時間（τ）の特性を示すむだ時間要素と、電力系統の周波数を検出した２つの時期の時間差の特性をその時定数（Tr）に基づいて示すハイパスフィルタ要素と、２つの時期のうち最新時期における周波数の検出時間の特性をその時定数（Td）に基づいて示すローパスフィルタ要素と、２つの時期の周波数の差である第２の周波数偏差に応じてパワーコンディショナが出力する無効電力のゲイン（Ｇ）の特性を示すゲイン要素との積で表現し、
電力系統の伝達関数は、負荷およびＰＣＳの双方の有効電力特性指数に基づく合成有効電力特性指数（α）と、電力系統の線路インピーダンス（Ｒ,Ｘ）と、ＰＣＳから出力された無効電力とによって、影響を受けた電力系統の周波数偏差（Δfc）を位相変化の時間微分で表現し、
安定限界関数AH（R,X）は、現状解析用の開ループ伝達関数G1(s)のむだ時間要素、ハイパスフィルタ要素、変数ｓをそれぞれ１とした上で、その絶対値を１とする下記の数式（４）で表現することを特徴とするＰＣＳの周波数フィードバック方式のハンチング現象の解析方法。
【数３】
0006466873.tif 000018
【数４】
0006466873.tif 000019
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、太陽光発電用パワーコンディショナ（以下、「ＰＣＳ」と言う。）の単独運転検出方式に用いられる周波数フィードバック方式のハンチング現象の解析方法、および回避方法に関する。
【背景技術】
【０００２】
太陽光発電設備は、太陽光パネルと、太陽光パネルを電力会社の電力系統に接続するＰＣＳとを備えている。ＰＣＳは、太陽光パネルで発生した直流の電力を交流にして電力系統に供給する機能と、単独運転を検出する機能と、単独運転を検出した場合に電力の供給を停止する機能を備えている。単独運転とは、変電所から配電線への電力の送電が事故等で停止した場合に、当該配電線に太陽光発電設備から電力供給が行われることを言う。
近年、太陽光発電の普及により、配電線に複数台の太陽光発電設備が連系される事案が増大している。そして各々のＰＣＳの単独運転検出方式で発生させる信号が相互に干渉すると、ＰＣＳが単独運転を検出できなくなることから、単独運転検出方式を周波数フィードバック方式に統一することが日本電機工業会規格のＪＥＭ１４９８で定められている。
【０００３】
周波数フィードバック方式とは、１）電力系統の周波数を所定のサンプリング周期で取得し、２）図４に示すように２つの時期の周波数、つまり現在の周波数fnと過去の周波数fpを、それぞれ移動平均により算出し（この例では現在の周波数fnは現在から40ms前までの期間内の平均値、過去の周波数fpは現在よりも200ms前から80ms前までの期間内の平均値として算出し）、３）その算出した現在と過去の周波数の差、つまり周波数偏差Δf2を所定の周期ごとに検出して、周波数偏差Δf2の±方向および大きさに応じた無効電力を電力系統に供給することを繰り返すものである。そしてＰＣＳは、周波数フィードバック方式によって単独運転時には単独運転系統の周波数変動を生じさせ、単独運転を検出するものである。
【０００４】
また周波数フィードバック方式のＰＣＳの特性（周波数偏差Δf2に対する無効電力のゲイン特性）は図５に示すように、−0.01Hz≦Δf2≦＋0.01Hzの場合にはΔf2の±方向に応じて同じ方向の無効電力を少量供給し、それよりもΔf2が拡大して−0.01Hz≦Δf2≦＋0.01Hzの範囲外になった場合には−0.5Hz≦Δf2≦＋0.5Hzの範囲内で、その拡大に比例して上限値および下限値±0.25p.u.（ＰＣＳの定格容量を基準値とする。）とする大量の無効電力を供給し、周波数偏差Δf2がさらに拡大して−0.5Hz＞Δf2となった場合には、上限値の無効電力を供給し、＋0.5Hz＜Δf2となった場合には下限値の無効電力を供給するようになっている。つまりゲイン特性は、正確にはリミッタ付きゲイン特性と言える。
【先行技術文献】
【非特許文献】
【０００５】
【非特許文献１】日本電機工業会規格  ＪＥＭ１４９８  ステップ注入付周波数フィードバック方式（太陽光発電用パワーコンディショナの標準形能動的単独運転検出方式）  ２０１２年（平成２４年）８月２７日制定
【発明の概要】
【発明が解決しようとする課題】
【０００６】
このように周波数フィードバック方式などの能動的単独運転検出方式は電力系統を不安定化させる危険がある。またＰＣＳから供給された無効電力がハンチング現象を発生させたことが現実的に確認されているが、その仕組みは理論的に解明されていない。
【０００７】
本発明は上記実情を考慮して創作されたもので、その目的は、ＰＣＳがハンチング現象を起こす仕組みを理論的に解明し、その理論を適用した解析装置でハンチング現象を解析できるようにし、さらにはその解析装置でハンチング現象を回避できるようにすることを目的とする。
【課題を解決するための手段】
【０００８】
図１は、発電設備が連系された状態の電力系統を簡略化して示している。電力系統は、無限大母線１に線路２を介して、太陽光発電設備３と負荷４とが並列接続されたものである。また線路２のうち途中には開閉器５が介在しており、事故等の場合には開閉器５の箇所よりも下位系統を切り離すものである。そして単独系統になっていないにも係わらず、太陽光発電設備３内のＰＣＳから供給される無効電力の出力が周期的に変動するハンチング現象を生じさせることがある。
【０００９】
ハンチング現象は、無効電力の出力の±方向が非常に短い時間内に繰り返し反転する現象なので、負荷中のモータの回転数は変化せず、従って負荷は定インピーダンスに見える。また無効電力の初期値を０とすることにより、図１の電力系統では負荷は抵抗のみで模擬している。
【００１０】
図２は、ＰＣＳのハンチング現象を解析するための解析用の制御系を示すブロック線図である。解析用の制御系は、ＰＣＳの周波数フィードバック方式の制御系２０と、ＰＣＳが接続される電力系統の伝達系３０との、フィードバック接続によって表現される。
【００１１】
フィードバック接続は、加え合わせ点１０に対して電力系統の基準周波数の偏差Δf0（通常運転時は０）に電力系統の周波数偏差Δfcが加算されるように入力されるものである。これによってＰＣＳは、基準周波数の偏差Δf0に対する電力系統の周波数偏差Δfcの和である第１の周波数偏差Δf1を検出でき、電力系統が単独運転であるか否かが検出できる。
【００１２】
ＰＣＳの周波数フィードバック方式の制御系２０は、第１の周波数偏差Δf1の検出に要するむだ時間τの特性を示すむだ時間要素２１と、電力系統の周波数を検出した２つの時期の時間差の特性をその時定数Trに基づいて示すハイパスフィルタ要素（以下、「ＨＰＦ要素」とする。）２２と、２つの時期のうち最新時期における周波数の検出時間の特性をその時定数Tdに基づいて示すローパスフィルタ要素（以下、「ＬＰＦ要素」とする。）２３と、２つの時期の周波数の差である第２の周波数偏差Δf2に応じてＰＣＳが出力する無効電力のゲインＧの特性を示すゲイン要素（リミッタ付きゲイン要素）２４との積で表現される。したがってＰＣＳの周波数フィードバック方式の伝達関数ＧPCS(s)は、以下の数式（１）で表現される。
【数１】
0006466873.tif 000002
【００１３】
ＨＰＦ要素の２つの時期とは、周期毎に過去と最近の周波数を検出するにあたって用いた時期であって、過去の周波数の検出に用いた周波数の時系列データからなる所定期間における最後の時（図４では、80msの移動平均に用いた時系列データの最後の時）と、最新の周波数の検出に用いた周波数の時系列データからなる所定期間における最後の時（図４では、40msの移動平均に用いた時系列データの最後の時）のことである。
【００１４】
制御理論を適用するために以下のようにする。ＰＣＳのゲイン特性、つまり図５で示した周波数偏差Δf2に応じて出力させる無効電力との関係の特性を線形近似し、ゲイン要素のゲインＧは、Ｇ= -0.25/(0.5Hz/60Hz)= -30.0 とする。ＨＰＦ要素の時定数Tr=0.2s。ＬＰＦ要素の時定数Tdは、40msの移動平均の時定数であり、時定数Td=0.04s/2=0.02s。（分母の値２は経験値である。）また第１の周波数偏差Δf1の検出にはむだ時間e-τs（τ=0.5〜1サイクル程度）が伴うので、これも考慮する。
【００１５】
電力系統の伝達系３０は、以下のようにして求める。
周波数フィードバック方式により発生する無効電力ΔＱは、電力系統に供給され、電力系統の電圧の位相θ、および電圧の大きさＶを変化させる。位相変化をΔθ、電圧の大きさ変化をΔＶとする。また無効電力ΔＱと負荷および太陽光発電の電圧特性により、負荷の有効電力PL、および太陽光発電の有効電力PGが変化する。負荷の有効電力の変化をΔPL、太陽光発電の有効電力の変化をΔPGとする。そうすると、位相変化Δθ、電圧の大きさの変化ΔＶは、電力系統の有効電力をＰ、その変化をΔＰとすると、以下の数式（２），（３）で表される。
【数２】
0006466873.tif 000003
【数３】
0006466873.tif 000004
ここで太陽光発電および負荷の有効電力特性指数をαG,αLとし、合成有効電力特性指数をαとすると、ΔPG、ΔPL、ΔP、αは以下の数式（４）、（５）、（６）、（７）で表される。
【数４】
0006466873.tif 000005
PG0：太陽光発電の有効電力の初期値
【数５】
0006466873.tif 000006
PL0：負荷の有効電力の初期値
【数６】
0006466873.tif 000007
【数７】
0006466873.tif 000008
数式（２）、（３）、（６）、（７）から以下の数式（８）、（９）が得られる。
【数８】
0006466873.tif 000009
【数９】
0006466873.tif 000010
第１の周波数偏差Δｆ1は、位相変化の時間微分であるから以下の数式（１０）が得られる。
【数１０】
0006466873.tif 000011
ここでω0 ＝２πf0である。ω0：定格角周波数、f0：定格周波数。
【００１６】
上記した数式（１０）が、電力系統の伝達系３０での周波数偏差（Δfc）の元になるものを示すものであり、ブロック線図で表現すると、図２に示すようになる。
したがって電力系統の伝達関数Gx(s)は、以下の数式（１１）で表現される。
【数１１】
0006466873.tif 000012
電力系統の伝達系３０は、まずＰＣＳの制御系２０からの無効電力ΔＱが並列接続された線路インピーダンスの抵抗ＲとリアクタンスＸに別々に入力される。無効電力ΔＱとリアクタンスＸとの積である電圧ＸΔＱが加え合わせ点３１に対して加算されるように出力されると共に、抵抗の電圧−ＲΔＰが加え合わせ点３１に対して減算されるように出力され、加え合わせ点３１で電圧変化ΔＶが算出される。また電圧変化ΔＶと電力系統の有効電力特性指数αとの積αΔＶ＝−ΔＰが、引出し点３２を介して線路インピーダンスの抵抗Ｒに出力され、−ΔＰと抵抗Ｒとの積が電圧−ＲΔＰとなり、前記したように電圧−ＲΔＰが加え合わせ点３１に対して減算されるように出力される。いっぽう−ΔＰが引出し点３２を介して線路インピーダンスのリアクタンスＸに出力され、−ΔＰとリアクタンスＸとの積である位相−ＸΔＰが加え合わせ点３３に対して加算されるように出力され、ΔＱと線路インピーダンスの抵抗Ｒとの積である位相ＲΔＱが加え合わせ点３３に対して加算されるようにして出力される。加え合わせ点３３では位相変化Δθの負値が算出され、位相変化の負値−Δθと微分係数s/ω0との積が電力系統の周波数偏差Δfcとして算出され、フィードバック接続の加え合わせ点１０に出力される。
【００１７】
このように電力系統の伝達系３０が表現できることから、ＰＣＳの周波数フィードバック方式の制御系２０と電力系統の伝達系３０は、フィードバック制御系を形成する。本例では、ゲインＧ＝―３０という大きな負値をとっているので、フィードバック制御系が正帰還になり、不安定になりやすくなる。
【００１８】
フィードバック制御系は、加え合わせ点１０で切った場合における開ループ伝達関数（一巡伝達関数）Ｇ1(s)は、以下の数式（１２）で表される。
【数１２】
0006466873.tif 000013
【００１９】
フィードバック制御系が安定となるためには、高周波域（概ね７Hz以上）での（１１）式のゲインの絶対値AHが１未満でなければならない。またゲインの絶対値AHの数式（１１）を簡略化し、ゲインの絶対値AHが１となるようなＲ、Ｘを変数とするものが、安定限界関数AH（R,X）であり、下記の数式（１３）で示される。
【数１３】
0006466873.tif 000014
【００２０】
今、Ｇ＝-30、 Td＝0.02、αL＝２、αG＝１、PL＝PG＝1.0（p.u.:基準値は定格容量）の場合の安定限界関数AH（R,X）を安定限界曲線として示すＲ−Ｘ平面のグラフが、図６に示されている。グラフは、縦軸をＸ、横軸をＲとするもので、安定限界曲線は右肩下がりの曲線になる。安定限界曲線の内側にＲとＸが位置する電力系統の場合には安定であり、ハンチング現象が発生せず、安定限界曲線の外側にＲとＸが位置する電力系統の場合には不安定であり、ハンチング現象が発生する。
【００２１】
図２の解析用の制御系を電中研Ｙ法と称する電力系統安定度解析プログラム（以下、「Ｙ法」）を用いてシミュレーションした。なおＹ法ではむだ時間要素を用いることが現時点では困難であったので考慮しなかったが、原理的にはむだ時間要素を考慮しなくても安定限界は変わらないと考えられる。そこで電力系統の電線（電線サイズ80〜100sq）を想定し、R=X/2 としてＲ、Ｘすなわち配電線の亘長を変化させて、ハンチング現象を再現することとした。シミュレーション結果を図７に示す。太陽光発電の定格容量を基準としてR=0.17、 X=0.34 (p.u.)のとき、時刻t=3 秒以降でハンチング現象が発生した。またこの条件（R=0.17、 X=0.34)以外にも、ＲとＸの比を若干変化させて、Ｙ法でのハンチング現象の発生状況を確認した。
【００２２】
その確認結果を、図６に示すＲ−Ｘ平面のグラフに併記してある。ハンチング現象が発生した場合を黒三角で表記し、ハンチング現象が発生しなかった場合を白丸で表記してある。黒三角は、安定限界曲線の外側領域、つまり不安定領域に存在し、白丸の位置は、安定限界曲線の内側領域、つまり安定領域に存在している。したがってゲインの絶対値の数式（１２）が１未満のときに、フィードバック制御系が安定となるという説が図６において矛盾しないことが分かる。
【００２３】
ＰＣＳの制御系２０が不安定となる原因は、開ループ伝達関数Ｇ1(s)が高周波域で減衰しないためである。そこでハンチング対策として図３に示すように、ＰＣＳの制御系２０（現状解析用の制御系）に、時定数Thに基づいて一次遅れの特性を示す一次遅れ要素２５を付加したハンチング対策用の制御系２０Ａを作成し、高周波域でのゲインを落とすことを試みた。
ハンチング対策用の制御系２０Ａの伝達関数Ｇ2(s)は、現状解析用の開ループ伝達関数Ｇ1(s)と一次遅れ要素２５との積から、下記の数式（１４）で表現される。
【数１４】
0006466873.tif 000015
また開ループ伝達関数（一巡伝達関数）Ｇ1 (s)にｓ＝ｊω（ω：角周波数）を代入して周波数伝達関数Ｇ2(jω)を作成する。そして周波数伝達関数Ｇ2(jω)を用いてナイキスト線図のグラフを作成する。ナイキスト線図は、縦軸を虚数軸、横軸を実数軸とするものである。図８に示すナイキスト線図は、Th＝0、Th＝0.05の２通りの場合に、R=0.17、 X=0.34、Ｇ=-30.0、ＨＰＦ要素の時定数Tr=0.2s、ＬＰＦ要素の時定数Td=0.02s、むだ時間e-τs（τ=0.017）とし、角周波数ωを０〜無限大まで変化させて、周波数伝達関数Ｇ2(jω)の座標が動いた軌跡を描いたものである。Th＝０のときは現状解析用の制御系であり、このときのナイキスト軌跡はむだ時間要素により高周波域でも円を描き、不安定となっていることが分かる。いっぽうTh=0.05 の一次遅れを付加したときはハンチング対策用の制御系であり、このときのナイキスト軌跡は円が小さくなり、原点に収束する形になり、安定となっていることが分かる。
【００２４】
なおTh=0.05 の一次遅れを付加したハンチング対策用の制御系をＹ法に適用したシミュレーション結果を、図７に併記した。これによりハンチング対策用の制御系では、ハンチング現象が生じないことが確認された。したがってＰＣＳの制御回路に一次遅れ要素を付加することが、ハンチング対策には有効である。
【発明の効果】
【００２５】
本発明によれば、現状解析用の開ループ伝達関数G1(s)を、上記した数式で表現できることから、ＰＣＳがハンチング現象を起こす仕組みが理論的に解明され、解析装置でハンチング現象を解析できるようになった。また現状解析用の開ループ伝達関数G1(s)を利用したハンチング対策用の伝達関数G2(s)を解析装置で解析させることにより、ハンチング現象を回避できるようになった。
【図面の簡単な説明】
【００２６】
【図１】発電設備が連系された状態の電力系統を示す説明図である。
【図２】ハンチング現象が生じる現状の制御系を説明するブロック線図である。
【図３】ハンチング対策用の制御系を説明するブロック線図である。
【図４】周波数フィードバック方式によって無効電力を演算する基になる周波数偏差の検出方法を示す説明図である。
【図５】周波数フィードバック方式の無効電力のゲイン特性を示す説明図である。
【図６】Ｒ−Ｘ平面図のグラフである。
【図７】ハンチング現象のシミュレーション結果を示すグラフである。
【図８】ナイキスト線図のグラフである。
【発明を実施するための形態】
【００２７】
本発明のＰＣＳの周波数フィードバック方式のハンチング現象の解析方法、およびハンチング現象の回避方法は、コンピュータからなる解析装置で実行される。解析装置は、キーボードやマウス等からなる入力装置、ディスプレイ、プリンター等からなる出力装置、解析プログラムの命令を順番に実行するＣＰＵ、解析プログラムや解析プログラムの実行に必要なデータ及び計算結果等を保存する記憶装置を構成要素とする標準的なコンピュータである。解析プログラムをコンピュータに実行させることにより、解析プログラムとコンピュータとが協働した各種の手段として機能し、各種のステップが順番に行われる。解析プログラムには、グラフとしてのナイキスト線図を描かせるものと、グラフとしてのＲ−Ｘ平面図を描かせるものが存在する。
【００２８】
まずＲ−Ｘ平面図を描かせる解析プログラムを実行させた場合を説明する。
初めにデータ入力ステップが行われる。データ入力ステップでは、入力フォーム表示手段によって、入力フォームが記憶装置から読み込まれて、出力装置のディスプレイに表示される。入力フォームは、ゲインＧ、ＬＰＦ要素の時定数Td、合成有効電力特性指数α（αL、αG、PL、PG）、定格角周波数ω0、電力系統の線路インピーダンス（抵抗Ｒ、リアクタンスＸ）を入力するためのものである。そして入力フォームに対してこれら値を入力し、その入力操作を確定する操作をすると、これら入力値が記憶装置に保存されると共に、安定限界曲線を描画する安定限界曲線描画ステップが行われる。
【００２９】
安定限界曲線描画ステップでは、安定限界曲線描画手段が安定限界曲関数AH（R,X）を記憶装置から読み込んで、電力系統の線路インピ―ダンス以外の入力値を、安定限界関数AH（R,X）に代入し、その上で安定限界関数AH（R,X）中のＲとＸを変化させて、関数が１となるＲとＸを多数算出し、それらの値をグラフ中に安定限界曲線として出力装置のディスプレイに表示すると共に、記憶装置に保存する。
【００３０】
そして次に電力系統をグラフ中に描画する電力系統描画ステップが行われる。電力系統描画ステップでは、電力系統描画手段がその線路インピーダンスを記憶装置から読み込んで、グラフ中に電力系統を所定の印で描いて出力装置のディスプレイに表示する。以上でＲ−Ｘ平面図を描かせる解析プログラムの実行を終了する。
【００３１】
次にナイキスト線図を描かせる解析プログラムを実行させた場合を説明する。
初めにデータ入力ステップが行われる。データ入力ステップでは、入力フォーム表示手段によって、入力フォームが記憶装置から読み込まれて、出力装置のディスプレイに表示される。入力フォームは、むだ時間τ、ＨＰＦ要素の時定数Tr、ＬＰＦ要素の時定数Td、ＰＣＳ要素のゲインＧ、合成有効電力特性指数α（αL、αG、PL、PG）、電力系統の線路インピーダンス（抵抗Ｒ、リアクタンスＸ）、一次遅れの時定数Thを入力するためのものである。なお一次遅れの時定数Thには、複数の値が入力できるようになっており、たとえば０を入力したときにはハンチング対策用の伝達関数G2(s)が現状解析用の伝達関数G1(s)として機能し、Thに０よりも大きな値を入力したときにはハンチング対策用の伝達関数G2(s)がそのまま機能する。そして入力フォームに対してこれら値を入力し、その入力操作を確定する操作をすると、これら入力値が記憶装置に保存されると共に、ナイキスト軌跡を描画するナイキスト軌跡描画ステップが行われる。
【００３２】
ナイキスト軌跡描画ステップでは、ナイキスト軌跡描画手段がハンチング対策用の開ループ伝達関数G2(s)を読み込んで、変数ｓをｊωと置き換え、ω以外の値を開ループ伝達関数に代入し、その上でωを変数として０〜∞まで変化させたときの値を多数算出し、それらの算出した値（座標）を結んだナイキスト軌跡をグラフ中に描いて出力装置のディスプレイに表示すると共に、記憶装置に記憶する。以上でナイキスト線図を描かせる解析プログラムの実行を終了する。
【符号の説明】
【００３３】
１  無限大母線
２  線路
３  太陽光発電設備
４  負荷
５  開閉器
１０  加え合わせ点
２０  ＰＣＳの制御系
２０Ａ  ハンチング対策用の制御系
２１  むだ時間要素
２２  ＨＰＦ要素
２３  ＬＰＦ要素
２４  ゲイン要素
２５  一次遅れ要素
３０  電力系統の伝達系
３１  加え合わせ点
３２  引出し点
３３  加え合わせ点
【図１】
0006466873.tif 000020
【図２】
0006466873.tif 000021
【図３】
0006466873.tif 000022
【図４】
0006466873.tif 000023
【図５】
0006466873.tif 000024
【図６】
0006466873.tif 000025
【図７】
0006466873.tif 000026
【図８】
0006466873.tif 000027