このページの内容は最新ではありません。最新版の英語を参照するには、ここをクリックします。

stepplot

ステップ応答をプロットし、プロット ハンドルを返す

ページ内をすべて折りたたむ

構文

h = stepplot(sys)
stepplot(sys,Tfinal)
stepplot(sys,t)
stepplot(sys1,sys2,...,sysN)
stepplot(sys1,sys2,...,sysN,Tfinal)
stepplot(sys1,sys2,...,sysN,t)
stepplot(AX,...)
stepplot(..., plotoptions)
stepplot(..., dataoptions)

説明

h = stepplot(sys) は、動的システム モデル sys のステップ応答をプロットします。また、プロット ハンドル h を返します。getoptions と setoptions コマンドを使ってプロットをカスタマイズするためにこのハンドルを使用することができます。次のように入力すると

help timeoptions

使用できるプロット オプションの一覧が得られます。

多入力モデルについては、個々の step コマンドが各入力チャネルに適用されます。時間範囲と計算点数は、自動的に選択されます。

stepplot(sys,Tfinal) は、t = 0 から最終時間 t = Tfinal までのステップ応答をシミュレーションします。システム時間単位で Tfinal を表し、sys の TimeUnit プロパティで指定されます。サンプル時間を指定しない離散時間システム (Ts = -1) では、stepplot は Tfinal を、シミュレートするサンプリング間隔の数として解釈します。

stepplot(sys,t) は、シミュレーションのためにユーザーが指定した時間ベクトル t を使用します。システム時間単位で t を表し、sys の TimeUnit プロパティで指定されます。離散時間モデルでは、t は Ti:Ts:Tf の形式です。ここで、Ts はサンプル時間です。連続時間モデルでは、t は Ti:dt:Tf の形式を取ります。ここで、dt は連続システムに対する離散近似のサンプル時間となります (step を参照してください)。stepplot コマンドは、Ti に関係なく、常に t=0 でステップ入力を適用します。

1 つのプロットで複数のモデル sys1,..., sys2 のステップ応答をプロットするには、次の構文を使用します。

stepplot(sys1,sys2,...,sysN)

stepplot(sys1,sys2,...,sysN,Tfinal)

stepplot(sys1,sys2,...,sysN,t)

システムごとに次のように、色、ラインスタイル、およびマーカーを指定することもできます。 

stepplot(sys1,'r',sys2,'y--',sys3,'gx')

stepplot(AX,...) は、ハンドル AX をもつ座標軸内にプロットします。

stepplot(..., plotoptions) は、オプション セットplotoptions を使用してプロットの外観をカスタマイズします。timeoptions を使用してオプション セットを作成します。

stepplot(..., dataoptions) は、オプション セット dataoptions を使用してステップ振幅や入力オフセットなどのオプションを指定します。stepDataOptions を使用してオプション セットを作成します。

例

すべて折りたたむ

この例では次を使用します。

ライブ スクリプトを開く

2 つの動的システムのステップ応答プロットを生成します。

sys1 = rss(3);
sys2 = rss(3);
h = stepplot(sys1,sys2);

Figure contains an axes. The axes contains 2 objects of type line. These objects represent sys1, sys2.

各ステップ応答は、異なる定常値で整定されます。プロット ハンドルを使用して、プロットされた応答を正規化します。

setoptions(h,'Normalize','on')

Figure contains an axes. The axes contains 2 objects of type line. These objects represent sys1, sys2.

ここで、応答は任意の単位で示された同じ値で整定されます。

この例では次を使用します。

ライブ スクリプトを開く

パラメトリックと同定されたモデルのステップ応答をノンパラメトリック (実測) モデルのステップ応答と比較し、3-σ の信頼領域も表示します (同定されたモデルには、System Identification Toolbox™ ソフトウェアが必要です)。

サンプル データからパラメトリックおよびノンパラメトリック モデルを同定します。

load iddata1 z1
sys1 = ssest(z1,4); 
sys2 = impulseest(z1);

両方の同定されたモデルのステップ応答をプロットします。プロット ハンドルを使用して、3-σ の信頼領域を表示します。

t = -1:0.1:5;
h = stepplot(sys1,sys2,t);
showConfidence(h,3)
legend('parametric','nonparametric')

Figure contains an axes. The axes with title From: u1 To: y1 contains 2 objects of type line. These objects represent parametric, nonparametric.

ノンパラメトリック モデル sys2 で不確かさが高まっているのがわかります。

ライブ スクリプトを開く

非線形の Hammerstein-Wiener モデルを推定するためのデータを読み込みます。

load(fullfile(matlabroot,'toolbox','ident','iddemos','data','twotankdata'));
z = iddata(y,u,0.2,'Name','Two tank system');

z は、入出力の推定データを格納する iddata オブジェクトです。

推定データを使用して次元が [1 5 3] の Hammerstein-Wiener モデルを推定します。入力の非線形性を区分的線形として指定し、出力の非線形性を 1 次元多項式として指定します。

sys = nlhw(z,[1 5 3],pwlinear,poly1d);

入力オフセットおよびステップ振幅レベルを指定するオプション セットを作成します。

opt = stepDataOptions('InputOffset',2,'StepAmplitude',0.5);

指定されたオプションを使用して 60 秒までのステップ応答をプロットします。

stepplot(sys,60,opt);

Figure contains an axes. The axes with title From: u1 To: y1 contains an object of type line. This object represents sys.

ヒント

単位などのプロットのプロパティを変更することができます。プロットのプロパティの変更方法の詳細については、を参照してください。

バージョン履歴

R2012a で導入

参考

| | |