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
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