The Polynomial Toolbox for MATLAB - DCE FEL ČVUT v Praze
The Polynomial Toolbox for MATLAB - DCE FEL ČVUT v Praze
The Polynomial Toolbox for MATLAB - DCE FEL ČVUT v Praze
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
Another<br />
triangular <strong>for</strong>m<br />
Hermite <strong>for</strong>m<br />
0 0.29 0.5<br />
0 -0.87 + 0.29s 0.5 + 0.5s<br />
-1 -0.29s^2 -0.5s^2<br />
If A( s)<br />
has not full row rank then T ( s)<br />
is in staircase <strong>for</strong>m. Similarly, an upper-right<br />
triangular (row staircase) <strong>for</strong>m is achieved by row (unimodular) operations. It results<br />
from the call tri(A,'row').<br />
If B( s)<br />
is a square polynomial matrix with nonsingular constant term then another<br />
upper-triangular <strong>for</strong>m may be obtained by the overloaded macro lu:<br />
B = [ 1 1 s; s+1 0 s; 1-s s 2+s]<br />
B =<br />
1 1 s<br />
1 + s 0 s<br />
1 - s s 2 + s<br />
[V,T] = lu(B)<br />
V =<br />
T =<br />
1 0.33s 0.11s<br />
1 + s 1 + 1.3s + 0.33s^2 0.44s + 0.11s^2<br />
1 - s 1 - 1.7s - 0.33s^2 1 - 0.56s - 0.11s^2<br />
1 1 + 0.33s 0.78s + 0.11s^3<br />
0 -1 -0.67s - s^2 + 0.33s^3<br />
0 0 2 + 2s + 2s^2 - s^3<br />
<strong>The</strong> triangular <strong>for</strong>ms described above are by no means unique. A canonical triangular<br />
<strong>for</strong>m is called the Hermite <strong>for</strong>m. An n m polynomial matrix A( s)<br />
of rank r is in<br />
column Hermite <strong>for</strong>m if it has the following properties:<br />
it is lower triangular<br />
the diagonal entries are all monic<br />
each diagonal entry has higher degree than any entry on its left<br />
in particular, if the diagonal element is constant then all off-diagonal elements in<br />
the same row are zero<br />
if n r then the last n r columns are zero<br />
<strong>The</strong> nomenclature in the literature is not consistent. Some authors (in particular<br />
Kailath, 1980) refer to this as the row Hermite <strong>for</strong>m. <strong>The</strong> polynomial matrix A is in<br />
row Hermite <strong>for</strong>m if it is the transpose of a matrix in column Hermite <strong>for</strong>m. <strong>The</strong><br />
command<br />
H = hermite(A)<br />
returns the column Hermite <strong>for</strong>m<br />
H =<br />
1 0 0<br />
1 + s s^2 0<br />
40