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

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Bézout<br />

equations<br />

Matrix<br />

polynomial<br />

equations<br />

To investigate the case of different minimal degree solutions, consider a right hand<br />

side of higher degree<br />

c = 15+15*s^4;<br />

As be<strong>for</strong>e, the call<br />

[x1,y1] = axbyc(a,b,c)<br />

x1 =<br />

y1 =<br />

8 - 13s + 15s^2<br />

7 + 12s + 2s^2<br />

results in the solution of minimal overall degree (in this case deg x1 deg y1<br />

2 ).<br />

A slightly different command<br />

[x2,y2] = axbyc(a,b,c,'minx')<br />

x2 =<br />

y2 =<br />

10.0000<br />

5 - 5s - 15s^2 - 15s^3<br />

returns another solution with the minimal degree of the first unknown. Finally,<br />

typing<br />

[x2,y2] = axbyc(a,b,c,'miny')<br />

x2 =<br />

y2 =<br />

10 - 15s + 15s^2<br />

5 + 10s<br />

produces the solution of minimal degree in the second unknown.<br />

Should the equation be unsolvable, the function returns NaNs..<br />

[x,y] = axbyc(s,s,1)<br />

x =<br />

y =<br />

NaN<br />

NaN<br />

A Diophantine equation with 1 on its right hand side is called a Bézout equation. It<br />

may look like<br />

a( s) x( s) b( s) y( s)<br />

1<br />

with a( s),<br />

b( s)<br />

given and x( s)<br />

, y( s)<br />

unknown.<br />

In the matrix case, the polynomial equation becomes a polynomial matrix equation.<br />

<strong>The</strong> basic matrix polynomial (or polynomial matrix) equations are<br />

and<br />

A( s) X ( s) B( s)<br />

X ( s) A( s) B( s)<br />

44

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

Saved successfully!

Ooh no, something went wrong!