incidenceMatrix

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

構文

A = incidenceMatrix(eqs,vars)

説明

m 個の方程式 eqs と n 個の変数 vars に対する A = incidenceMatrix(eqs,vars) は、m 行 n 列の行列 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)];

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

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