Main Content

incidenceMatrix

方程式系の接続行列の検索

説明

m 個の方程式 eqsn 個の変数 vars に対する A = incidenceMatrix(eqs,vars) は、mn 列の行列 A を返します。ここで、eqs(i)vars(j) または vars(j) の導関数が含まれる場合は A(i,j) = 1 となります。A の他のすべての要素は 0 です。

接続行列

5 変数の 5 つの方程式から成る系の接続行列を求めます。

次の 5 つのシンボリックな微分方程式を含むシンボリック ベクトル eqs を作成します。

syms y1(t) y2(t) y3(t) y4(t) y5(t) c1 c3
eqs = [diff(y1(t),t) == y2(t),...
       diff(y2(t),t) == c1*y1(t) + c3*y3(t),...
       diff(y3(t),t) == y2(t) + y4(t),...
       diff(y4(t),t) == y3(t) + y5(t),...
       diff(y5(t),t) == y4(t)];

変数のベクトルを作成します。ここで、c1c3 はこの系のシンボリック パラメーターです (変数ではありません)。

vars = [y1(t), y2(t), y3(t), y4(t), y5(t)];

方程式 eqs および変数 vars についての接続行列 A を求めます。

A = incidenceMatrix(eqs, vars)
A =
     1     1     0     0     0
     1     1     1     0     0
     0     1     1     1     0
     0     0     1     1     1
     0     0     0     1     1

入力引数

すべて折りたたむ

方程式。シンボリックな方程式または式のベクトルとして指定します。

変数。x(t) など、シンボリック変数、シンボリック関数またはシンボリック関数呼び出しのベクトルとして指定します。

出力引数

すべて折りたたむ

接続行列。倍精度値の行列として返されます。

バージョン履歴

R2014b で導入