記述関数法シミュレーター — 非線形制御

リレー制御・オンオフ制御の設計：家庭用サーモスタット、冷蔵庫のコンプレッサ制御、電気ヒーターのオンオフ制御など、アクチュエータが「全開／全閉」しか取れない系は世の中に無数にあります。これらは本質的にリミットサイクルを持ち、室温やプロセス量が一定振幅で振れ続けます。記述関数法を使えば、ヒステリシス幅や線形部の時定数から「振れ幅」と「振動周期」を設計段階で見積もり、許容範囲に収めるための定数を選べます。

アクチュエータの飽和とリミットサイクル：サーボモータ、油圧バルブ、電力増幅器などはすべて出力に上限（飽和）があります。高ゲインのフィードバック系では、この飽和が原因で持続振動が起きることがあります。記述関数法で飽和要素を等価ゲイン N(A) に置き換えれば、ゲインをどこまで上げると振動に入るか、入った場合の振幅はどれくらいかを予測でき、ゲイン設計の指針になります。

機械系の摩擦・バックラッシュによる微小振動：歯車の遊び（バックラッシュ）や静摩擦・クーロン摩擦は不感帯やヒステリシスとしてモデル化できます。位置決めステージや工作機械の送り軸では、これらが原因で目標値付近の微小なリミットサイクル（リミットサイクル振動・スティックスリップ）が生じます。記述関数法は、こうした摩擦起因の振動が消えるか残るかを評価する古典的な手段です。

非線形制御の事前検討と教育：詳細な非線形シミュレーションや位相平面解析を行う前に、記述関数法で「リミットサイクルがありそうか」を素早く当たりづけします。計算が軽く、複素平面上の交点という直感的な描像が得られるため、制御工学の教育でも非線形系の入り口として広く使われます。本ツールのように L(jω) と −1/N(A) を同時に描くと、交点の意味が一目で理解できます。

まず最大の誤解が、「記述関数法は厳密解を与える」という思い込みです。記述関数法は出力の基本波成分だけを残し、第3・第5高調波などをすべて捨てる近似です。これが成り立つのは、線形部が高調波を十分に減衰させる「フィルタ仮説」が満たされる場合だけ。高調波の減衰が弱い系では、予測されたリミットサイクルの振幅・周波数が実際とずれます。本ツールの線形部は積分器＋二重ラグでローパス性が強いので近似は良好ですが、一般には記述関数法の結果は「目安」であり、最終確認は時間領域シミュレーションで行うのが鉄則です。

次に、「交点があれば必ず持続振動になる」と考えること。L(jω) と −1/N(A) の交点は「リミットサイクルの候補」であって、その振動が安定（外乱で振幅がずれても元に戻る）か不安定（ずれるとそのまま発散・消滅する）かは別問題です。安定性は、交点の近傍で振幅 A を少し増減させたとき、系がその振幅へ戻る向きに動くかどうかで判定します。飽和系では安定なリミットサイクル、不感帯系では不安定なリミットサイクルが現れやすいなど、非線形要素の性格によって挙動が変わります。本ツールは候補点を示すものであり、安定性判定は別途の検討が必要です。

最後に、「N(A) は実数とは限らない」という点。本ツールが扱うリレー（ヒステリシスなし）・飽和・不感帯はいずれも奇対称かつ単一値の静的非線形なので N(A) は実数になります。しかしヒステリシス付きリレーやバックラッシュのように「記憶」を持つ非線形要素では、N(A) は複素数になり、臨界曲線 −1/N(A) が実軸から離れて曲線を描きます。実在のリレーは必ず多少のヒステリシスを持つため、現実の設計では複素記述関数まで踏み込む必要があることを覚えておいてください。