主要内容

。中文:中文:。

オンラインのARMAX多項式モデル推定

★★★★★★

スクリプトを開く

。1つの非線形化学反応プロセスについて、ARMAX モデルを 2 つ推定します。これらのモデルは、2 つの操作条件におけるプロセスの動作を再現します。このモデルの動作はオンラインで同定され、システムの操作時に適応 PI コントローラーのゲインを調整するために使用されます。

www.qqqq.com www.qqqq.com

。★★★★★★★★★

これは,Bequetteの著書“过程动力学:建模、分析和仿真”(普伦蒂斯·霍尔出版社,1998年)で詳細に記述されている被覆非断熱(非断熱的)タンク反応器です。容器は完全に混合されるものと仮定され,単一の1次不可逆発熱反応A - B >が発生します。我的意思是，我的意思是，我的意思是。攪拌後,最終製品が反応物一个と同じ流量で容器から流出します(反応タンクの容積は一定)。★★★★★★★★★★★★★★★★★★★★非断熱連続攪拌タンク反応器:仿金宝app真软件®でのシミュレーションによるMATLABファイルモデリング通吃通吃通吃通吃通吃

CSTR。

$$ \begin{array} {ll} u_1(t) = C_{Af}(t) \;& # 38;\textnormal{进料中A的浓度 流}[kgmol/m^3] \\ u_2(t) = T_f(t) \;& # 38;\textnormal{进料流温度}[K] \\ u_3(t) = T_j(t) \;& # 38;\textnormal{护套冷却液温度}[K] \\ \end{array} $$

(。

数组$ $ \开始{}{你}& # xA; y_1 (t) = x_1 (t) = C_{一}\ (t);& # 38;\textnormal{反应器槽内A浓度}[kgmol/m^3] \\ y_2(t) = x_2(t) = t (t) \;& # 38;\textnormal{反应器温度}[K] \\ \end{array} $$

【中文翻译美元$ C_{一}(t)を,時間とともに変化する目的のレベル美元$ C_ {Aref} (t)草皮，草皮，草皮你知道吗美元美元T_f (t)我的意思是，我的意思是，我的意思是，我的意思是美元美元T_j (t)。PI美元C_ {Aref} (t) -C_A (t)美元啊哈!你说得对美元$ C_ {Af} (t)。金宝appSimulinkiddemo_cstrは,装运箱プラントをブロック装运箱> > > > > >

open_system (“iddemo_cstr”）;

你知道吗美元美元T_f (t)。0.0075 [k ^2]。★★★★★★★★美元$ C_ {Aref} (t)【中文翻译】

この例の信号美元$ C_ {Aref} (t)は,t = 400美元1.5 [kgmol/m^3] 2 [kgmol/m^3]★★★★★★★★★★★★★★★★★★★★★★★★★★★★美元$ C_ {Aref} (t)。0.015。。閉ループ同定で参照信号が十分に励起されない場合,一意のモデルを同定するために十分な情報が得られないことがあります。美元$ C_ {Aref} (t)中文:iddemo_cstr / CA参考。

適応制御のオンライン推定

【中文翻译美元$ C_{一}(t)は美元$ C_{一}(t)【翻译美元美元T_j (t)> > > > > > > > > > > > > > > > > > > >この感度の変化を検出するために,递归的多项式模型估计ブロックを使用します。【翻译】美元$ C_{一}(t)。これは,不安定の要因となり得る高ゲイン制御ループを使用しないためです。

递归多项式模型估计器美元美元T_j (t)から美元$ C_{一}(t)。。追従誤差美元C_ {Aref} (t) -C_A (t)美元。この正則化はπコントローラーの設計用途である,初期操作点における追従誤差でゲイン1を得ることにより行われます。# # # # # # # # # # # # # # # # # # # # # # # # # # #。この適応コントローラーはiddemo_cstr/自适应PI控制器。

递归多项式模型估计器

“递归多项式模型估计器”，Simulink金宝app系统识别工具箱/估计器。。阿玛斯，阿玛斯，阿玛斯，阿玛斯，阿玛斯。

数组$ $ & # xA; \开始{}{1}{y} (q) \酒吧(t) = B (q) \酒吧{你}(t) + C (q) \酒吧}{e (t) \ \ [0.1] & # xA; (q) = 1 + a_1z ^ {1} + a_2z ^ {2} + \ cdots +现代z ^ {na} {na} \ \ & # xA; B (q) = (b_01 + b_1z ^ {1} + b_2z ^ {2} + \ cdots +现代z ^ {nb-1} {nb + 1}) z ^ {nk} \ \ & # xA; C (q) = 1 + c_1z ^ {1} + c_2z ^ {2} + \ cdots + c_ z ^{数控}{数控}\ \ & # xA;结束\{数组}& # xA; $ $

  • 再帰多項式モデル推定器ブロックの输入?输出你说得对，你说得对酒吧\{你}(t)美元と酒吧\ {y} (t)美元? ? ? ?装运箱モデルの$ {y} $ \酒吧と{你}$ $ \酒吧【中文翻译】:中文翻译:;$ {y} \酒吧= C_A (t) - C{} \酒吧(t)美元、${你}\酒吧= T_j酒吧(t) - \ {t} _j (t)美元……ピーク振幅が1になるよう{你}$ $ \酒吧と$ {y} $ \酒吧。平衡化操作点$ {C} \酒吧(t)美元?$ {T} \酒吧_j (T)美元。★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★iddemo_cstr /预处理Tjブロックとiddemo_cstr / CA进行预处理。

open_system (“iddemo_cstr /预处理Tj”）;

  • 递归多项式模型估计器启用。。启用。【翻译】t \[200400)美元とt \[600800)美元★★★★★★★★★★★★★★★★★★★美元美元T_j (t)。

递归多项式模型估计器

。[]。

  • 模型结构:ARMAX。我的意思是，我的意思是，我的意思是美元美元T_f (t)【中文译文美元美元C_A (t)> > > > > > > > > > > > > > > > > > >

  • 最初的估计:一个也没有。。

  • A(q) (na)：2。非線形モデルは2つの状態をとります。

  • B(q) (nb)：2。

  • C(q) (nc)：2。Na, nb, nc, ni, ni, ni, ni, ni, ni, ni, ni, ni, ni, ni, ni, ni, ni, ni, ni, ni, ni, ni, ni, ni, ni, ni, ni。

  • (英文)：1。多くの物理システムと同様に,装运箱システムには直達がありません。“i / o”，“i / o”，“i / o”，“i / o”。

  • パラメーター共分散行列：1 e4。中文意思是:“我的意思是，我的意思是，我的意思是，我的意思是，我的意思是。”

  • 【翻译：0.1。装运箱モデルは約0.25赫兹の帯域幅をもつことがわかっています。T_s = 0.1美元を選択するのは,1/0.1がこの帯域幅の20倍(5赫兹)より大きくなるからです。

[算法和块选项]。

  • 推定法：遗忘因子。

  • 忘却係数：1-5e-3。推定パラメーターが操作点とともに変化することが予測されるため,忘却係数を1未満の値に設定します。$\lambda = 1-5e-3$? ? ? ?★★★★記憶の時定数$ T_0 = \压裂{Ts}{1 - \λ}= 100美元> > > > > >記憶の時間を100秒に設定することにより,各操作点で200秒の同定期間から大量のデータが取得され,同定に使用されます。

  • [输出估计误差]。。

  • [添加启用端口]。。。

  • 这是真的吗：没有一个。

递归多项式模型估计器

【中文译文】:【中文译文一个(q)美元、B (q)美元、C (q)美元★★★★$ {e} $ \酒吧? ? ? ?この多項式モデル推定器ブロックの错误【翻译酒吧\ {e} (t)美元。ブロックの参数。【中文】na = 2美元、nb = 2美元、数控= 2美元、nk = 1美元) ? ? ?参数。

$$ \begin{array} {l} A(t) = [1 \;a_1 (t) \;B(t) = [0 \;b_0 (t) \;C(t) = [1 \;c₁(t) \;\\ \end{array} $$

多項式(q), B (q), C (q)の推定パラメーターは,シミュレーション中に次のように変化します。

sim卡(“iddemo_cstr”）;close_system (“iddemo_cstr /预处理Tj”）;open_system (“iddemo_cstr / ABC”）;

初期のパラメーター共分散行列に高い値が選択されているため,パラメーターの推定値は初期値0から即座に変化します。多項式一个(q)美元?B (q)美元のパラメーターは,t = 200美元★★★★★★★【翻译C (q)美元。(1)、(1)、(1)、(1)、(1)美元美元C_A (t)这是一个很好的例子美元美元T_f (t)。中文:我的意思是，我的意思是，我的意思是，我的意思是，我的意思是C (q)美元★★★★★★★★★★★★美元美元T_j (t)から美元美元C_A (t)への関係が伝達関数$ \压裂{B (q)} {(q)} $。★★★★★★★★★C (q)美元。

t \[200400)美元() () () (启用。パラメーター推定は,装运箱タンクが新規の操作点に切り替わるt = 400美元★★★★★★一个(q)美元?B (q)美元のパラメーターは,t = 600美元★★★★★★★★★启用■■■■■■■■■■■■■■■■この操作点での一个(q)美元とB (q)美元。この遅い収束はt = 400におけるパラメーター共分散行列の固有値がt = 0で設定された初期の1 e4値と比較して小さいことが原因です。。t = 0における推定の信頼度が低かったため,アルゴリズムはパラメーター推定値をすばやく変更しました。改善されたパラメーター推定値によりシステムの動作が的確に把握され,その結果,パラメーター共分散行列の1ステップ先の予測誤差と固定値が小さくなります(信頼性が向上します)。i =400}}}}}}}。递归多项式模型估计ブロックの[外部リセット] オプションを使用して、t=400 におけるパラメーター共分散の新しい値を指定できます。パラメーター共分散の値を表示するには、Recursive Polynomial Model Estimator ブロックの[输出参数协方差矩阵]。

英文翻译

递归多项式模型估计ブロックの错误1。

close_system (“iddemo_cstr / ABC”）;open_system ('iddemo_cstr/提前1步预测错误'）;

【翻译美元美元T_j (t)チャネルに追加の摂動を与えない場合,1ステップ先の誤差はより大きい値になります。誤差が大きくなる理由として,推定器ブロックが依存する入力チャネル美元美元T_j (t)★★★★★★★★★★★★★★★() () () (美元美元C_A (t)。。

推定モデルのより厳密な検査では,推定モデルのシミュレーションを実行し,実際のモデル出力と比較します。iddemo_cstr /时变ARMAXブロックは,在线多项式模型估计ブロックで推定した時変ARMAXモデルを実装しています。中文翻译:中文翻译:中文翻译:中文翻译:

close_system ('iddemo_cstr/提前1步预测错误'）;open_system (“iddemo_cstr /模拟误差”）;

★★★★★★★★美元美元C_A (t)。。

中文:matlab。操作点C_A = 1.5美元kgmol / m ^ 3美元?C_A = 2美元kgmol / m ^ 3美元★★★★★★★★$t = 200$?$t = 600$。。

Ts = 0.1;= find(t>=200,1);P200 = idpoly (AHat (:,:, tidx) BHat (:,:, tidx),聊天(:,:,tidx), 1, 1, [], Ts);= find(t bb0 =600,1);P600 = idpoly (AHat (:,:, tidx) BHat (:,:, tidx),聊天(:,:,tidx), 1, 1, [], Ts);bodemag (P200“b”P600,“r——”{10 ^ 1 20});传奇(在C_A=1.5 [kgmol/m^3]时的估计模型，…估算模型C_A=2.0 [kgmol/m^3]，…“位置”，“最佳”）;

。★★★★非線形 CSTR プラントに関する事前知識と一致しています。美元C_A (t) = 2 kgmol / m ^ 3美元における伝達関数は,低周波数においてゲインが6 db美元。

∕∕∕

2つのARMAXモデルを推定して,2つの操作条件における非線形装运箱プランドの動作を再現しました。。推定結果を検証するために,2つの信号を確認しました。それらは1ステップ先の予測誤差および装运箱プラント出力と推定モデルのシミュレーション間の誤差です。これらの誤差信号は,装运箱プラント出力と比較すると,一定の範囲内にあり,小さい値でした。。

bdclose (“iddemo_cstr”）;

系统识别工具箱

超级超级超级超级

MATLAB,金宝app Simulink，

無料評価版のダウンロード