【構造から設計まで】PID制御の基礎事項を詳細解説

この記事ではPID制御についてまとめます。以下の動画では、PID制御について約8分で解説しています。概要を述べただけの動画ですが、YouTubeでは1.6万回視聴されていて、制御工学チャンネル内でも最も視聴されている動画の一つとなっています。ここでは、ブロック線図の構造や、PID制御器の中身（連続・離散）、設計法、シミュレーションについて触れたいと思います。

youtu.be

PID制御とは
基本的なPID制御器設計法
制御シミュレーション
動画・関連記事
- PID制御の動画一覧（制御工学チャンネル）
- PID制御の解説記事

PID制御とは

PID制御は、PID制御は自動制御を行う際に最も多く利用されている制御手法の一つであり、参照入力（目標入力）と出力の偏差の各時刻での値に基づいて制御入力を決める手法です。Pは比例、Iは積分、Dは微分の頭文字であり、フィードバック制御方式において、偏差信号に比例した値、偏差を積分した値、偏差を微分した値の3つの和を取った値を制御に用います。

まず、PID制御の歴史について知りたい場合は荒木・須田の解説記事を読まれると良いと思います。

まずは、PID制御の基本ブロック線図を示します。

制御目的

制御対象に入力 $u$ を印加した結果出力 $y$ が得られますが、この $y$ の振る舞いは自由に決めることができません。 $y$ の目標としては目標信号を与え、例えばステップ上の目標値を与えてやりますと、ある一定の値に出力 $y$ を整定させたいというのが制御目的になります。その目的の達成度を評価するために各種の過渡応答に関する指標があります。

目標信号 $r$ と出力 $y$ の偏差が常にゼロとなればよいですが、必ずしもそうはならないため PID制御器をうまく調整することになります。制御原因を決定するときには3つの値を調整することになりますが、その際によく用いられる指標は以下の通りになります。

最大オーバーシュート
定常偏差
整定時間
立ち上がり時間

どのような応答波形なら実現できるかは、制御対象の特性に依存しますし、設計者のセンスにも依存します。様々な応答指標については以下の動画も併せてご覧頂ければと思います。

youtu.be

制御器の内部構造

その上で PID制御ではよく試行錯誤的にPIDのゲイン、すなわち係数を決定することになります。ここでは、PID制御器の構造について触れます。

比例動作

まず P 動作比例動作ですが、比例動作は出力と目標値との偏差を考えて、その偏差に比例係数をかけたものを制御入力として用います。

$\displaystyle u(t) = K_P e(t)$

ボード線図は次のようになります。ゲインが周波数によらず一定値となっていることがわかります。

微分動作

次にD動作です。 D動作では偏差の変化量に比例した入力を構成します。このためボード線図を見てやると低い周波数ではゲインが小さく、高い周波数ではゲインが大きくなります。

より変化に対して感度の高い入力を決定することになります。そしてこのD動作とP動作を組み合わせたPD制御がよく使われます。

$\displaystyle u(t) = K_P e(t) + K_D \dfrac{de(t)}{dt}$

PD制御では、ゲインはある周波数まではフラットとなって、ある周波数以降ではどんどん高い値になっていくボード線図の形状となります。このグラフの形状はPDゲインをどの値に設定するかによってシフトすることになります。

積分動作

最後にI動作です。 I動作では偏差信号を考え、その積分値を考えます。この積分値に比例した入力を決定します。このI動作に対してボード線図を見てやりますと、低い周波数では高いゲインとなり高い周波数になるに従って、そのゲインが下がっていくような構成になっています。

この積分器を構成する時にはその初期値を初期化することに注意する必要があります。この積分器は目標信号がステップの場合には、目標信号と出力との偏差をゼロにする働きがあります。

$\displaystyle u(t) = K_P e(t) + K_D \dfrac{de(t)}{dt} + K_I\int_0^t e(t) dt$

これについては内部モデル原理の動画をご覧下さい。そしてPID制御器を構成すると次のようなボード線図の形状となります。

積分器の初期値は通常0としますが、その設定を考えずに実装すると想定外の応答になることがありますので注意してください。

コントローラーを設計する際には、まず出力 $y$ がどのような応答になればよいかということを考える必要があります。その上でコントローラーの係数をどのような組み合わせにするかを考えることで、制御性能を良くしていくことになります。出力の結果がどのような応答波形となるかは、制御対象の特性によっても依存しますしゲインの係数によっても依存して変わってくることになります。以下で設計法について触れます。

離散時間系のPID制御

連続時間と同様に離散時間の場合もPID制御系を構成することができます。実装は、和分と差分を組み合わせて行います。

$\displaystyle u(k) = K_P e(k) + K_D (e(k)-e(k-1)) + K_I\sum_{i=0}^k e(i)$

基本的なPID制御器設計法

制御を実装するためには、 $K_P, K_D, K_I$ をそれぞれ定める必要があります。パラメータ設計法は様々であり、モデルベース手法やモデルレス手法があります。ここではいくつか紹介します。

限界感度法

初期の設計方法として限界感度法が知られています。限界感度法は、比例動作のみでフィードバックを施し、ゲインを徐々に大きくした上で安定限界となった際のゲインに基づいてPID制御パラメータを決定する手法です。モデルレスの設計法ですが、高次の複雑な制御対象や不安定化すると危険な制御対象にはお勧めできません。どのような対象に有用であるかは、教科書等を含めた各種文献をご覧ください。

ステップ応答法

ステップ入力を印加したときの波形に基づいてゲインを設定する手法がステップ応答法です。制御対象をむだ時間プラス一次遅れ系として表現し、その係数に基づいてPIDゲインを定めます。どのような対象に有用であるかは、教科書等を含めた各種文献をご覧ください。

北森法

北森法は、部分モデルマッチング法とも呼ばれ、PIDゲインを規範モデルに基づいて決定する手法です。モデルとなる入出力伝達関数（規範モデルの伝達関数）と、制御対象、PIDゲインの組において、PIDゲインを設定した制御系が低周波数から規範モデルに合わせこむ手法となっており、規範モデルさえ与えられていれば機械的にゲインを定めることができます。

《PID制御》第4回: PID制御方式の原理・実用化と設計法

周波数領域での設計

推定された周波数応答を使用した PID コントローラーの設計 - MATLAB & Simulink - MathWorks 日本

Measurement of the Loop Gain Frequency Response of Digitally Controlled Power Converters | IEEE Journals & Magazine | IEEE Xplore

周波数特性分析器技術解説集 - エヌエフ回路設計ブロック (nfcorp.co.jp)

VRFT, FRIT

制御データを用いて直接的に制御器設計を行う手法として、IFTやVRFT、FRITなどが知られています。PID制御器設計においても、これらの制御データに基づく手法が多く研究されています。

閉ループデータに基づく直接的PID調整とその不安定プロセス制御への適用 (jst.go.jp)

Pythonによるデータ駆動制御VRFTによるPID制御器オートチューニング #データサイエンス - Qiita

制御性能評価を用いたオンライン型FRITに基づくフィードフォワードコントローラの調整法 (jst.go.jp)

非反証制御による設計

PID制御器のパラメータ空間設計法(<特集>PID制御のパラメータチューニング)

過渡応答データに基づくPID制御器の直接設計 (jst.go.jp)

GKYP

周波数応答に基づいて設計を行う場合、一つのアプローチとしてGKYPが挙げられます。GKYPは一般化KYP補題であり、周波数帯域ごとに制約を与えて設計することが可能です。PID制御は、係数を決める問題ですので、GKYPとの親和性が高いです。GKYPに基づいた設計法については、以下の解説記事を参照ください。

一般化KYP補題に基づくPID制御器の設計(<特集>PID制御のパラメータチューニング)