22.07.2013 Views

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

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!