PID制御のプログラミング2

ヒータで層の温度を一定にする系で、実際にやったことがあるのでアドバイスします。

理論的なアドバイスでは有りませんが参考になればと思います。

2.について
参照先URLに紹介されている式では、I制御による変化量は、制御開始からずっと累積されてしまっているのがわかります。
系によると思いますが、おっしゃる通り、この式では積分時間やサンプリング時間を考慮していないので、この式をそのままプログラムに落とすのは、私の場合、あまり実用的では有りませんでした。

そこで、
ki=kp/Ti　，　kd=kp*Td　ですので、　　(Ti=積分時間、td=サンプル時間)
n=Ti/Td　　(n=サンプル回数)　として、
ei=　Σe／ｎ　（ei=積分時間あたりの偏差の平均）
ed=　en　-　en-1　(ed=前回偏差と今回偏差の差)
とし、
M= kp*( en + ei + ed )
と、操作量を決めてあげました。
このやり方は直感的にわかりやすく使いやすいです。

実際、操作量には上限がありますし、負の値を取れない系もありますので、
早速プログラムを組み始まってしまうのも一つの手だと思います。
ジーグラとニコルズの最適値は、オーバシュート量が25%に達しますので、
私の場合、カットアンドトライで最適値を決める必要がありました。
やってみて気付くことはとても多かったです。

1.についてですが、
LとTの単位はサンプリング時間と同じにするとやりやすいと思います。サンプリング時間の単位が秒であれば、時定数や無駄時間の単位も秒にします。
しかし、実際に計算してみると、これでも時間の単位の違いによって、パラメータは大きく異なってしまうのがわかりますが、これも実際やってみると、秒がいいのか分がいいのか、わかってくると思います。

#1に対して。

おっしゃる通り、Tdとtdは同じ物です。間違えました。

http://pelab.nagaokaut.ac.jp/kondolab/convenienc …

を参考にしました。ゲインの決め方ではなく、積分時間とサンプリング時間の決め方が書いてあります。参考になれば。

> 制御周期毎にPID演算を実行して辻褄が合うためには、

> 制御パラメータを制御周期で補正(?)しなければならないような気がして
その通りですね。
それが#2で述べた，積分器なら"*ΔT"，微分器なら"/ΔT"の項です。

> 制御パラーメータを求める際に使用する時間を制御周期を基準に
> (制御周期の何倍かを)代入すれば必要なパラメータが得られると
> 解釈しました。その理解で正しいでしょうか。
これも良いと思います。
ただ#2では基準時間は何でも良いと書きましたが，いや実際それで間違いはないのですが，
素直に1[s]としておいた方がいろいろと混乱を避けられるのも確かです。

あと#1さんへの補足にあった
> Ki＝0.6～0.7/KLやKd＝0.3～0.45T/Kを導き方
ですが，これらは経験則によって導かれたものですので，明快な導き方を示した資料はないのではないかと思います。

○積分器と微分器の実装方法が分かれば大方解決でしょうか。

いろいろ方法はありますが，
1. 積分器
・区分求積法使用
∫f(t)dt≒Σ{f(kΔT)*ΔT} とします(積分区間は省略)。
ここでΔTは制御周期(一般的にはサンプリング周期と同じ)，f(kΔT)はk回目のサンプリングにおける信号値です。
ΔTが十分に小さければこれで十分な精度が得られます。
さらに高精度化を図るには台形公式で近似します。
2. 微分器
・後退差分で近似
df(t)/dt≒(f(kΔT)-f((k-1)ΔT))/ΔT

これらの式で制御周期毎にPID演算を実行すればよいです。

.....とここまで書いて前の質問に記載されていたもう一つのURLを見てみたら同じ式が書かれてますね。
ただどちらのURLの内容も質問者様には帯に短しな内容みたいですので，もう少し資料を探された方が良いかと思います。
書籍でもWebサイトでも，詳しいのはいくらでもありますよ。

○単位について
どんな単位を使っても構わないのですが，混乱を避けるために同じ次元の数値には違う単位を使わいことをお勧めします。
例えば時間なら全て[秒]で統一し，[分]や[時]などは机上検討でも極力使わないようにするなどです。
[10msec]や[0.1msec]などでも構いません。
単位の変換に線型性が成立つ限り，実装時の単位は実装者の都合のいいように決定すればよいです。