28.6 その他の関数
: y = poly (A)
: y = poly (x)
AがN行N列の正方行列である場合、はAの特性多項式 の 係数の行ベクトルです。 poly (A)det (z * eye (N) - A)
たとえば、次のコードは、の根であるAの固有値を見つけます。 poly (A)
roots (poly (eye (3)))
⇒ 1.00001 + 0.00001i
1.00001 - 0.00001i
0.99999 + 0.00000i
実際、3 つの固有値はすべて 1 であり、数値パフォーマンスのためにはeig関数を使用して固有値を計算する必要があることを強調しています。
xがベクトルの場合、 は、根がxの要素である多項式の係数のベクトルです。つまり、cが多項式の場合、 の要素は cに含まれます。ただし、ソートおよび数値エラーのため、ベクトルcとd は同一ではありません。 poly (x)d = roots (poly (c))
See also: roots, eig.
: polyout (c)
: polyout (c, x)
: str = polyout (…)
多項式cのフォーマットされたバージョンを表示します。
フォーマットされた多項式
c(x) = c(1) * x^n + ... + c(n) x + c(n+1)
がゼロの場合は文字列として返されるか、画面に書き込まれますnargout。
2 番目の引数x は、各用語に使用する変数名を指定し、デフォルトは文字列になります"s"。
See also: polyreduce.
: p = polyreduce (c)
先頭のゼロを削除して、多項式係数ベクトルを最小の項数に減らします。
See also: polyout.