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
Spectral factorization zpplot(1+2*i*s-s^2),grid One of quadratic equations in scalar and matrix polynomials frequently encountered is the polynomial spectral factorization A( s) X ( s) JX ( s) and the spectral co-factorization A( s) X ( s) JX ( s) . In either case, the given polynomial matrix A( s) satisfies A( s) A ( s) (we say it is para-Hermitian symmetric) and the unknown X ( s) is to be stable. The case of A( s) positive definite on the stability boundary results in J I . Spectral factorization with J = I is the main tool to design LQ and LQG controllers as well as Kalman filters. On the other hand, if A( s) is indefinite in matrix sense then J diag l+1, + 1, , + 1, 1, 1, 1q . This is the famous J -spectral factorization problem, which is an important tool for robust control and filter design based on H 50
norms. The Polynomial Toolbox provides two macros called spf and spcof to handle spectral factorization and co-factorization, respectively. By way of illustration consider the para-Hermitian matrix A = 34 - 56s^2 -13 - 22s + 60s^2 36 + 67s -13 + 22s + 60s^2 46 - 1e+002s^2 -42 - 26s + 38s^2 36 - 67s -42 + 26s + 38s^2 59 - 42s^2 Its spectral factorization follows by typing [X,J] = spf(A) X = J = 2.1 + 0.42s 5.2 + 0.39s -2 + 0.35s -5.5 + 4s 4.3 + 0.64s -7.4 - 5.5s 0.16 + 6.3s -0.31 - 10s -0.86 + 3.5s 1 0 0 0 1 0 0 0 1 while the spectral co-factorization of A is computed via [Xcof,J] = spcof(A) Xcof = J = 2.7 + 0.42s 4.8 + 4s 2 + 6.3s -1.6 + 0.39s 0.93 + 0.64s -6.5 - 10s 4.3 + 0.35s 2.7 - 5.5s 5.8 + 3.5s 1 0 0 0 1 0 0 0 1 The resulting J reveals that the given matrix A is positive-definite on the imaginary axis. On the other hand, the following matrix is indefinite B = 5 -6 - 18s -8 -6 + 18s -41 + 81s^2 -22 - 18s -8 -22 + 18s -13 Its spectral factorization follows as [Xf,J] = spf(B) Xf = J = 3.3 2.9 - 0.92s -1.6 1.8 6.6 + 0.44s 3.4 1.6 2.5 + 9s -2 51
- Page 2 and 3: PolyX, Ltd E-mail: info@polyx.com S
- Page 4 and 5: Rank ..............................
- Page 6 and 7: Resampling of polynomials in z ....
- Page 8 and 9: 1 Quick Start Initialization Help E
- Page 10 and 11: Typing the name of the matrix P now
- Page 12 and 13: Concatenation and working with subm
- Page 14 and 15: ans = 1 0 0 1 s^2 s^3 1 - s 1 The c
- Page 16 and 17: 2 Polynomial matrices Introduction
- Page 18 and 19: Changing the default indeterminate
- Page 20 and 21: Degrees and leading coefficients R
- Page 22 and 23: 3 + 8s integral(F) ans = 2s + 1.5s^
- Page 24 and 25: isunimod(U) ans = 1 Also the adjoin
- Page 26 and 27: Bases and null spaces ans = 1 confi
- Page 28 and 29: If P( s) is square then its roots a
- Page 30 and 31: and simply type or or hurwitz(p) an
- Page 32 and 33: Least common multiple If the only c
- Page 34 and 35: Greatest common left divisor 0 -3 +
- Page 36 and 37: Dual concepts M = lrm(A,B) M = 0.32
- Page 38 and 39: 0 1 s^2 -s^3 -s 0 Reduced and canon
- Page 40 and 41: Another triangular form Hermite for
- Page 42 and 43: Invariant polynomials The entries a
- Page 44 and 45: Bézout equations Matrix polynomial
- Page 46 and 47: Factorizations returns X = Y = 0.25
- Page 48 and 49: zpplot(1-s^2), grid zpplot(1+2*s^2+
- Page 52 and 53: Non-symmetric factorization or 1 0
- Page 54 and 55: Transformation to Kronecker canonic
- Page 56 and 57: common roots (including roots at in
- Page 58 and 59: Derivatives and integrals ans = ans
- Page 60 and 61: Conjugate transpose The operation o
- Page 62 and 63: Zeros of symmetric twosided polynom
- Page 64 and 65: For polynomial matrices with comple
- Page 66 and 67: Non-symmetric equation Discrete-tim
- Page 68 and 69: Resampling Sampling period Discrete
- Page 70 and 71: Resampling of two-sided polynomials
- Page 72 and 73: 4 The Polynomial Matrix Editor Intr
- Page 74 and 75: Matrix Pad window Matrix Pad button
- Page 76 and 77: 5 Polynomial matrix fractions Intro
- Page 78 and 79: The sdf command Comparison of fract
- Page 80 and 81: Reverse managed by switching the Po
- Page 82 and 83: The mdf command -1 + s 0.5s ------
- Page 84 and 85: Leftdenominatorfraction The ldf com
- Page 86 and 87: Gmdf==Grdf ans = Grdf==Gldf ans = 1
- Page 88 and 89: Entrywise and matrix division ans =
- Page 90 and 91: Coefficients and coefficient matric
- Page 92 and 93: Transfer function matrices or f2 =
- Page 94 and 95: 0 and then just types F=mdf(A,B,C,D
- Page 96 and 97: Conversion to the Control System To
- Page 98 and 99: y1 0 y2 0 Continuous-time model. Gt
Spectral<br />
factorization<br />
zpplot(1+2*i*s-s^2),grid<br />
One of quadratic equations in scalar and matrix polynomials frequently encountered<br />
is the polynomial spectral factorization<br />
A( s) X (<br />
s) JX ( s)<br />
and the spectral co-factorization<br />
A( s) X ( s) JX (<br />
s)<br />
.<br />
In either case, the given polynomial matrix A( s)<br />
satisfies A( s) A ( s)<br />
<br />
(we say it is<br />
para-Hermitian symmetric) and the unknown X ( s)<br />
is to be stable. <strong>The</strong> case of A( s)<br />
positive definite on the stability boundary results in J I .<br />
Spectral factorization with J = I is the main tool to design LQ and LQG controllers as<br />
well as Kalman filters. On the other hand, if A( s)<br />
is indefinite in matrix sense then<br />
J diag l+1, + 1, , + 1, 1, 1, 1q<br />
. This is the famous J -spectral factorization<br />
problem, which is an important tool <strong>for</strong> robust control and filter design based on H <br />
50