Main Content

hermiteForm

行列のエルミート型

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

構文

H = hermiteForm(A)
[U,H] = hermiteForm(A)
___ = hermiteForm(A,var)

説明

例

H = hermiteForm(A) は、行列 A のエルミート標準形を返します。A の要素は、整数または symvar(A,1) で決定される変数の多項式でなければなりません。エルミート型 H は上三角行列です。

例

[U,H] = hermiteForm(A) は、A のエルミート標準形と、H = U*A となるようなユニモジュラー変換行列 U を返します。

例

___ = hermiteForm(A,var) では、A の要素が、指定された変数 var の一変数多項式であると仮定します。A に他の変数が含まれる場合、hermiteForm ではそれらの変数をシンボリック パラメーターとして扱います。

前述のいずれの構文でも入力引数 var を使用できます。

A に var が含まれない場合、hermiteForm(A) と hermiteForm(A,var) は異なる結果を返します。

例

整数行列のエルミート型

逆ヒルベルト行列のエルミート型を求めます。

A = sym(invhilb(5))
H = hermiteForm(A)
A =
[    25,   -300,    1050,   -1400,    630]
[  -300,   4800,  -18900,   26880, -12600]
[  1050, -18900,   79380, -117600,  56700]
[ -1400,  26880, -117600,  179200, -88200]
[   630, -12600,   56700,  -88200,  44100]
 
H =
[ 5,  0, -210, -280,  630]
[ 0, 60,    0,    0,    0]
[ 0,  0,  420,    0,    0]
[ 0,  0,    0,  840,    0]
[ 0,  0,    0,    0, 2520]

一変数多項式の行列のエルミート型

2 行 2 列の行列を作成します。この行列の要素は、変数 x の多項式です。

syms x
A = [x^2 + 3, (2*x - 1)^2; (x + 2)^2, 3*x^2 + 5]
A =
[   x^2 + 3, (2*x - 1)^2]
[ (x + 2)^2,   3*x^2 + 5]

この行列のエルミート型を求めます。

H = hermiteForm(A)
H =
[ 1, (4*x^3)/49 + (47*x^2)/49 - (76*x)/49 + 20/49]
[ 0,            x^4 + 12*x^3 - 13*x^2 - 12*x - 11]

多変数多項式の行列のエルミート型

2 つの変数 x と y を含む 2 行 2 列の行列を作成します。

syms x y
A = [2/x + y, x^2 - y^2; 3*sin(x) + y, x]
A =
[      y + 2/x, x^2 - y^2]
[ y + 3*sin(x),         x]

この行列のエルミート型を求めます。多項式の変数を指定しない場合、hermiteForm では symvar(A,1) を使用して、多項式の変数が x であると決定します。3*sin(x) + y は x の多項式ではないため、hermiteForm はエラーをスローします。

H = hermiteForm(A)
Error using mupadengine/feval (line 163)
Cannot convert the matrix entries to integers or univariate polynomials.

A のすべての要素が変数 y の多項式であると指定して、A のエルミート型を求めます。

H = hermiteForm(A,y)
H =
[ 1, (x*y^2)/(3*x*sin(x) - 2) + (x*(x - x^2))/(3*x*sin(x) - 2)]
[ 0,     3*y^2*sin(x) - 3*x^2*sin(x) + y^3 + y*(- x^2 + x) + 2]

エルミート型と変換行列

逆ヒルベルト行列について、エルミート型とその対応する変換行列とを求めます。

A = sym(invhilb(3));
[U,H] = hermiteForm(A)
U =
[ 13,  9,  7]
[  6,  4,  3]
[ 20, 15, 12]
 
H =
[ 3,  0, 30]
[ 0, 12,  0]
[ 0,  0, 60]

H = U*A であることを検証します。

isAlways(H == U*A)
ans =
  3×3 logical array
     1     1     1
     1     1     1
     1     1     1

多項式の行列について、エルミート型とその対応する変換行列とを求めます。

syms x y
A = [2*(x - y), 3*(x^2 - y^2);
     4*(x^3 - y^3), 5*(x^4 - y^4)];
[U,H] = hermiteForm(A,x)
U =
[                   1/2,  0]
[ 2*x^2 + 2*x*y + 2*y^2, -1]
 
H =
[ x - y,         (3*x^2)/2 - (3*y^2)/2]
[     0, x^4 + 6*x^3*y - 6*x*y^3 - y^4]

H = U*A であることを検証します。

isAlways(H == U*A)
ans =
  2×2 logical array
     1     1
     1     1

整数行列の変数を指定する場合

行列に特定の変数が含まれておらず、その変数を 2 番目の引数として指定して hermiteForm を呼び出す場合、その結果はその変数を指定しない場合に得られる結果とは異なります。たとえば、変数を 1 つも含まない行列を作成します。

A = [9 -36 30; -36 192 -180; 30 -180 180]
A =
     9   -36    30
   -36   192  -180
    30  -180   180

変数 x を 2 番目の引数として指定して、hermiteForm を呼び出します。この場合、hermiteForm では A の要素が x の一変数多項式であると仮定します。

syms x
hermiteForm(A,x)
ans =
     1     0     0
     0     1     0
     0     0     1

変数を指定せずに hermiteForm を呼び出します。この場合、hermiteForm は A を整数行列として扱います。

hermiteForm(A)
ans =
     3     0    30
     0    12     0
     0     0    60

入力引数

すべて折りたたむ

入力行列。整数として指定します。この行列の要素は、整数または一変数多項式です。A の要素に複数の変数が含まれる場合、引数 var を使用して多項式変数を指定し、他のすべての変数をシンボリック パラメーターとして扱います。A が多変数で、var を指定していない場合、hermiteForm では symvar(A,1) を使用して多項式の変数を決定します。

多項式の変数。シンボリック変数として指定します。

出力引数

すべて折りたたむ

入力行列のエルミート標準形。シンボリック対角行列として返されます。行列のエルミート型は上三角行列です。

変換行列。ユニモジュラー シンボリック行列として返されます。A の要素が整数の場合、U の要素も整数です。また、det(U) = 1 であるか、det(U) = -1 です。A の要素が多項式の場合、U の要素は一変数多項式であり、det(U) は定数です。

詳細

すべて折りたたむ

エルミート標準形

整数の係数をもつ n 行 n 列の正方行列 A には、n 行 n 列の行列 H と、 A*U = H となるような n 行 n 列のユニモジュラー行列 U が存在します。ここで、H は A のエルミート標準形です。ユニモジュラー行列は、行列式が 1 または -1 に等しい実数正方行列です。A が多項式の行列の場合、U の行列式は定数です。

hermiteForm は、非特異整数正方行列 A のエルミート標準形を Hjj≥0 かつ −Hjj2<Hij≤Hjj2 となる上三角行列 H として返します。ここで、j>i です。A が正方行列でも特異行列でもない場合、行列 H は単純に上三角行列になります。

バージョン履歴

R2015b で導入

参考

|