システム制御情報学会 学会誌 「システム/制御/情報」 56 巻,6 号,pp.???-???,2012 年

『初学者のための図解でわかる制御工学 II』 特集号

市原:むだ時間系の状態予測制御

 M ファイル
spc.zip
 spc.zip は,以下の M ファイルを含みます.
 ・dd_proper.m
 ・dd_strictproper.m
 ・delay.m
 ・loop_gain_proper.m
 ・loop_gain_strictproper.m
 ・rotatepend.m
 ・spc_sample.m
これら M ファイルはすべて同じフォルダに置いてください.
ユーザが書き変えて使うのは,おもに spc_sample.mrotatepend.m です.
実行に際して,Control System Toolbox が必要になります.

 spc_sample.m では,
  - 状態予測制御による時間応答(解説の第 6 図)および周波数応答(解説の第 7 図)
が計算できます.分布遅れ制御則に関しては,spc_sample.m において
  1. プロパな近似
     % 制御器 (どちらかを選択)
      % プロパ近似
      [input,c_mem] = dd_proper(A,B,Ft,h,N,x0,c_mem);
      % 厳密プロパ近似
      % [input,c_mem] = dd_strictproper(A,B,Ft,h,N,x0,c_mem,del,tau);
    
    % 制御器 (どちらかを選択)
     % プロパ近似
     loop_gain_proper;
     % 厳密プロパ近似
     % loop_gain_strictproper;
    
  2. 厳密にプロパな近似
     % 制御器 (どちらかを選択)
      % プロパ近似
      % [input,c_mem] = dd_proper(A,B,Ft,h,N,x0,c_mem);
      % 厳密プロパ近似
      [input,c_mem] = dd_strictproper(A,B,Ft,h,N,x0,c_mem,del,tau);
    
    % 制御器 (どちらかを選択)
     % プロパ近似
     % loop_gain_proper;
     % 厳密プロパ近似
     loop_gain_strictproper;
    

のいずれかの方法を選択して計算してください.また,厳密にプロパな近似に用いる調整パラメータ(tau:)は rotatepend.m の 22 行目
tau = 0.03;
で定義されており,ユーザは書き換えることができます.