didier henrion laas-cnrs toulouse, fr czech tech univ...
TRANSCRIPT
![Page 1: Didier HENRION LAAS-CNRS Toulouse, FR Czech Tech Univ ...cerfacs.fr/wp-content/uploads/2016/03/henrion.pdf · as a basis for polynomials of degree at most k ... 2.7603 5.1073 6.5115](https://reader033.vdocument.in/reader033/viewer/2022050410/5f86f5073f98be5f95111c28/html5/thumbnails/1.jpg)
Semidefinite programming
for polynomial optimization
and robust control
Didier HENRION
www.laas.fr/∼henrion
LAAS-CNRS Toulouse, FR
Czech Tech Univ Prague, CZ
CERFACS Sparse Days
Toulouse, October 2007
![Page 2: Didier HENRION LAAS-CNRS Toulouse, FR Czech Tech Univ ...cerfacs.fr/wp-content/uploads/2016/03/henrion.pdf · as a basis for polynomials of degree at most k ... 2.7603 5.1073 6.5115](https://reader033.vdocument.in/reader033/viewer/2022050410/5f86f5073f98be5f95111c28/html5/thumbnails/2.jpg)
Main points
Convex optimisation can be used to solve numerically non-convex
optimisation problems, specifically polynomial matrix inequalities
Numerical linear algebra is a key ingredient
Applications in many branches of engineering and applied maths,
in particular systems control
![Page 3: Didier HENRION LAAS-CNRS Toulouse, FR Czech Tech Univ ...cerfacs.fr/wp-content/uploads/2016/03/henrion.pdf · as a basis for polynomials of degree at most k ... 2.7603 5.1073 6.5115](https://reader033.vdocument.in/reader033/viewer/2022050410/5f86f5073f98be5f95111c28/html5/thumbnails/3.jpg)
PMIs and LMIs
Polynomial programming: objective function and constraints aremultivariate polynomials, or equivalently, we minimize a linearfunction on a (non-convex, non-connected, discrete..)semialgebraic set
Typically, the semialgebraic set is describedby polynomial matrix inequalities (PMIs)
Semidefinite programming (SDP): linear programming in the(convex but non-polyhedral) cone of positive semidefinitematrices, also called linear matrix inequality (LMI) optimisation
min cTxs.t. F (x) =
∑α Fαxα � 0
![Page 4: Didier HENRION LAAS-CNRS Toulouse, FR Czech Tech Univ ...cerfacs.fr/wp-content/uploads/2016/03/henrion.pdf · as a basis for polynomials of degree at most k ... 2.7603 5.1073 6.5115](https://reader033.vdocument.in/reader033/viewer/2022050410/5f86f5073f98be5f95111c28/html5/thumbnails/4.jpg)
An LMI set (and something more)
![Page 5: Didier HENRION LAAS-CNRS Toulouse, FR Czech Tech Univ ...cerfacs.fr/wp-content/uploads/2016/03/henrion.pdf · as a basis for polynomials of degree at most k ... 2.7603 5.1073 6.5115](https://reader033.vdocument.in/reader033/viewer/2022050410/5f86f5073f98be5f95111c28/html5/thumbnails/5.jpg)
Outline
1. Quasi-Hankel matrices
2. Generalized companion matrices
3. Systems control applications
![Page 6: Didier HENRION LAAS-CNRS Toulouse, FR Czech Tech Univ ...cerfacs.fr/wp-content/uploads/2016/03/henrion.pdf · as a basis for polynomials of degree at most k ... 2.7603 5.1073 6.5115](https://reader033.vdocument.in/reader033/viewer/2022050410/5f86f5073f98be5f95111c28/html5/thumbnails/6.jpg)
Outline
1. Quasi-Hankel matrices
2. Generalized companion matrices
3. Systems control applications
![Page 7: Didier HENRION LAAS-CNRS Toulouse, FR Czech Tech Univ ...cerfacs.fr/wp-content/uploads/2016/03/henrion.pdf · as a basis for polynomials of degree at most k ... 2.7603 5.1073 6.5115](https://reader033.vdocument.in/reader033/viewer/2022050410/5f86f5073f98be5f95111c28/html5/thumbnails/7.jpg)
Convex vs. non-convex optimisation
Build and solve a hierarchy of successive convex LMI relaxations
of increasing size, whose optima are guaranteed to converge
asymptotically to the global optimum
Relaxations are built thanks to results on the primal theory of
moments and the dual theory of sum-of-squares (SOS)
decompositions of multivariate polynomials
![Page 8: Didier HENRION LAAS-CNRS Toulouse, FR Czech Tech Univ ...cerfacs.fr/wp-content/uploads/2016/03/henrion.pdf · as a basis for polynomials of degree at most k ... 2.7603 5.1073 6.5115](https://reader033.vdocument.in/reader033/viewer/2022050410/5f86f5073f98be5f95111c28/html5/thumbnails/8.jpg)
LMI relaxations: an example
Quadratic problem
max x2s.t. 3− 2x2 − x1
2 − x22 ≥ 0
−x1 − x2 − x1x2 ≥ 01 + x1x2 ≥ 0
Non-convex feasible set delimited by circular and hyperbolic arcs
![Page 9: Didier HENRION LAAS-CNRS Toulouse, FR Czech Tech Univ ...cerfacs.fr/wp-content/uploads/2016/03/henrion.pdf · as a basis for polynomials of degree at most k ... 2.7603 5.1073 6.5115](https://reader033.vdocument.in/reader033/viewer/2022050410/5f86f5073f98be5f95111c28/html5/thumbnails/9.jpg)
Lifting
LMI relaxation built by replacing each monomial xi1x
j2 with
a lifting variable yij
For example, quadratic expression 3− 2x2 − x12 − x2
2 ≥ 0replaced with linear expression 3− 2y01 − y20 − y02 ≥ 0
Lifting variables yij satisfy non-convex relationssuch as y10y01 = y11 or y20 = y2
10
Relax these non-convex relations by enforcing LMI constraint
M11(y) =
1 y10 y01y10 y20 y11y01 y11 y02
� 0
Moment matrix of first order
![Page 10: Didier HENRION LAAS-CNRS Toulouse, FR Czech Tech Univ ...cerfacs.fr/wp-content/uploads/2016/03/henrion.pdf · as a basis for polynomials of degree at most k ... 2.7603 5.1073 6.5115](https://reader033.vdocument.in/reader033/viewer/2022050410/5f86f5073f98be5f95111c28/html5/thumbnails/10.jpg)
Projecting
First LMI relaxation
max y01
s.t.
1 y10 y01y10 y20 y11y01 y11 y02
� 0
3− 2y01 − y20 − y02 ≥ 0−y10 − y01 − y11 ≥ 01 + y11 ≥ 0
Projection onto the plane y10, y01 of first-order moments
LMI optimum = 2 = upper-bound on global optimum
![Page 11: Didier HENRION LAAS-CNRS Toulouse, FR Czech Tech Univ ...cerfacs.fr/wp-content/uploads/2016/03/henrion.pdf · as a basis for polynomials of degree at most k ... 2.7603 5.1073 6.5115](https://reader033.vdocument.in/reader033/viewer/2022050410/5f86f5073f98be5f95111c28/html5/thumbnails/11.jpg)
Higher order relaxation
To build second LMI relaxation, the moment matrix must capture
expressions of degrees up to 4
M22(y) =
1 y10 y01 y20 y11 y02y10 y20 y11 y30 y21 y12y01 y11 y02 y21 y12 y03y20 y30 y21 y40 y31 y22y11 y21 y12 y31 y22 y13y02 y12 y03 y22 y13 y04
Constraints are also lifted and relaxed with the help of
localizing matrices
Second LMI provides tighter relaxation
![Page 12: Didier HENRION LAAS-CNRS Toulouse, FR Czech Tech Univ ...cerfacs.fr/wp-content/uploads/2016/03/henrion.pdf · as a basis for polynomials of degree at most k ... 2.7603 5.1073 6.5115](https://reader033.vdocument.in/reader033/viewer/2022050410/5f86f5073f98be5f95111c28/html5/thumbnails/12.jpg)
max y01
s.t.
1 y10 y01 y20 y11 y02
y10 y20 y11 y30 y21 y12
y01 y11 y02 y21 y12 y03
y20 y30 y21 y40 y31 y22
y11 y21 y12 y31 y22 y13
y02 y12 y03 y22 y13 y04
� 0
3− 2y01 − y20 − y02 3y10 − 2y11 − y30 − y12 3y01 − 2y02 − y21 − y03
3y10 − 2y11 − y30 − y12 3y20 − 2y21 − y40 − y22 3y11 − 2y12 − y31 − y13
3y01 − 2y02 − y21 − y03 3y11 − 2y12 − y31 − y13 3y02 − 2y03 − y22 − y04
� 0
−y10 − y01 − y11 −y20 − y11 − y21 −y11 − y02 − y12
−y20 − y11 − y21 −y30 − y21 − y31 −y21 − y12 − y22
−y11 − y02 − y12 −y21 − y12 − y22 −y12 − y03 − y13
� 0
1 + y11 y10 + y21 y01 + y12
y10 + y21 y20 + y31 y11 + y22
y01 + y12 y11 + y22 y02 + y13
� 0
![Page 13: Didier HENRION LAAS-CNRS Toulouse, FR Czech Tech Univ ...cerfacs.fr/wp-content/uploads/2016/03/henrion.pdf · as a basis for polynomials of degree at most k ... 2.7603 5.1073 6.5115](https://reader033.vdocument.in/reader033/viewer/2022050410/5f86f5073f98be5f95111c28/html5/thumbnails/13.jpg)
Global optimality
Optimal value of 2nd LMI relaxation = 1.6180= global optimum within numerical accuracyNumerical certificate = moment matrix has rank one
First order moments
(y∗10, y∗01) = (−0.6180,1.6180)
provide optimal solution
of original problem
![Page 14: Didier HENRION LAAS-CNRS Toulouse, FR Czech Tech Univ ...cerfacs.fr/wp-content/uploads/2016/03/henrion.pdf · as a basis for polynomials of degree at most k ... 2.7603 5.1073 6.5115](https://reader033.vdocument.in/reader033/viewer/2022050410/5f86f5073f98be5f95111c28/html5/thumbnails/14.jpg)
Quasi-Hankel structure
Structure of moment matrices
univariate = Hankel
multivariate = sparse quasi-Hankel
1 x10 x01 x20 x11 x02x10 x20 x11 x30 x21 x12x01 x11 x02 x21 x12 x03x20 x30 x21 x40 x31 x22x11 x21 x12 x31 x22 x13x02 x12 x03 x22 x13 x04
Known sparsity pattern: LMI matrix coefficients
in pre-factored form (Cholesky)
![Page 15: Didier HENRION LAAS-CNRS Toulouse, FR Czech Tech Univ ...cerfacs.fr/wp-content/uploads/2016/03/henrion.pdf · as a basis for polynomials of degree at most k ... 2.7603 5.1073 6.5115](https://reader033.vdocument.in/reader033/viewer/2022050410/5f86f5073f98be5f95111c28/html5/thumbnails/15.jpg)
Quasi-Hankel structure
Efficient low-rank algebra to compute
gradient and Hessian in interior-point codes ?
Theoretical results by B. Mourrain and V. Pan
So far our implemention in GloptiPoly for Matlab
www.laas.fr/∼henrion/software/gloptipoly
uses a general purpose SDP solver (SeDuMi by default)
![Page 16: Didier HENRION LAAS-CNRS Toulouse, FR Czech Tech Univ ...cerfacs.fr/wp-content/uploads/2016/03/henrion.pdf · as a basis for polynomials of degree at most k ... 2.7603 5.1073 6.5115](https://reader033.vdocument.in/reader033/viewer/2022050410/5f86f5073f98be5f95111c28/html5/thumbnails/16.jpg)
Outline
1. Quasi-Hankel matrices
2. Generalized companion matrices
3. Systems control applications
![Page 17: Didier HENRION LAAS-CNRS Toulouse, FR Czech Tech Univ ...cerfacs.fr/wp-content/uploads/2016/03/henrion.pdf · as a basis for polynomials of degree at most k ... 2.7603 5.1073 6.5115](https://reader033.vdocument.in/reader033/viewer/2022050410/5f86f5073f98be5f95111c28/html5/thumbnails/17.jpg)
Detecting global optimalityGlobal optimisation problem
p? = minx g0(x)s.t. gi(x) ≥ 0, i = 1,2..
Let deg gi(x) = 2di − 1 or 2di and d = maxi di
LMI relaxation of order k
p∗k = miny∑
α(g0)αyα
s.t. Mk(y) � 0Mk−di
(giy) � 0, i = 1,2..
with solution y∗k
Hierarchy with convergence guarantee
p∗d ≤ p∗d+1 ≤ · · · ≤ p∗k∗ = p∗.
for (generally) small k∗
Sufficient condition for global optimalityrank of moment matrices
rank Mk(y∗k) = rank Mk−d(y∗k)
![Page 18: Didier HENRION LAAS-CNRS Toulouse, FR Czech Tech Univ ...cerfacs.fr/wp-content/uploads/2016/03/henrion.pdf · as a basis for polynomials of degree at most k ... 2.7603 5.1073 6.5115](https://reader033.vdocument.in/reader033/viewer/2022050410/5f86f5073f98be5f95111c28/html5/thumbnails/18.jpg)
Need for extraction procedures
Because rank condition is only sufficient
any other global optimality certificate is welcome
If there is a finite number of global optimisers,
they generate a zero-dimensional variety
Moment matrices are multiplication matrices in the quotient
ideal, a finite-dimension vector space
Key ideas:• Cholesky decomposition of moment matrix• column reduced echelon form• computation of common eigenvalues
![Page 19: Didier HENRION LAAS-CNRS Toulouse, FR Czech Tech Univ ...cerfacs.fr/wp-content/uploads/2016/03/henrion.pdf · as a basis for polynomials of degree at most k ... 2.7603 5.1073 6.5115](https://reader033.vdocument.in/reader033/viewer/2022050410/5f86f5073f98be5f95111c28/html5/thumbnails/19.jpg)
Cholesky factorization
Extract Cholesky factor V of moment matrix
Mk(y?k) = V V ′
Matrix V has r columns, corresponding to r globally optimal solutions x?j,
j = 1,2, . . . , r (provided global optimum was reached)
Choose e.g. v =[1 x1 x2 . . . xn x2
1 x1x2 . . . x1xn x22 x2x3 . . . x2
n . . . xkn
]′as a basis for polynomials of degree at most k
By definition of the moment matrix:
Mk(y?k) = V ?(V ?)′
where
V ? =[
v∗1 v∗2 · · · v∗r]
and v?j is polynomial basis v evaluated at solution x?
j
Extracting solutions amounts tofinding linear transformation
between Cholesky factors V and V ?
![Page 20: Didier HENRION LAAS-CNRS Toulouse, FR Czech Tech Univ ...cerfacs.fr/wp-content/uploads/2016/03/henrion.pdf · as a basis for polynomials of degree at most k ... 2.7603 5.1073 6.5115](https://reader033.vdocument.in/reader033/viewer/2022050410/5f86f5073f98be5f95111c28/html5/thumbnails/20.jpg)
Reduction to column echelon form
Next step is reduction of V into column echelon form
U =
1?0 10 0 1? ? ?
... . . .0 0 0 · · · 1? ? ? · · · ?
... ...? ? ? · · · ?
by Gaussian elimination with column pivotingEach row in U = monomial xα in basis vPivot entry in U = monomial xβj
in generating basis of the set of solutions
In other words, denoting w =[xβ1
xβ2. . . xβr
]′it holds
v = Uw
for all solutions x∗j, j = 1,2, . . . , r
![Page 21: Didier HENRION LAAS-CNRS Toulouse, FR Czech Tech Univ ...cerfacs.fr/wp-content/uploads/2016/03/henrion.pdf · as a basis for polynomials of degree at most k ... 2.7603 5.1073 6.5115](https://reader033.vdocument.in/reader033/viewer/2022050410/5f86f5073f98be5f95111c28/html5/thumbnails/21.jpg)
Multiplication matrices
For each first degree monomial xi extract from U the r-by-r
multiplication matrix Ni containing coefficients of product
monomials xixβjin generating basis w, i.e. such that
Niw = xiw i = 1,2, . . . , n
Given matrices Ni finding scalars xi is an eigenvalue problem
Extracting solutions amounts tosolving an eigenvalue problem
Eigenvectors w are shared by commuting matrices Ni so
it is a particular common eigenvalue problem
![Page 22: Didier HENRION LAAS-CNRS Toulouse, FR Czech Tech Univ ...cerfacs.fr/wp-content/uploads/2016/03/henrion.pdf · as a basis for polynomials of degree at most k ... 2.7603 5.1073 6.5115](https://reader033.vdocument.in/reader033/viewer/2022050410/5f86f5073f98be5f95111c28/html5/thumbnails/22.jpg)
Common eigenvalue problem
Build combination of multiplication matrices
N =n∑
i=1
λiNi
where λi are random positive numbers (summing up to one)
Compute ordered Schur decomposition
N = QTQ′
where Q =[
q1 q2 · · · qr
]is orthogonal and T upper triangular
Finally, due to orthogonality of vectors qi,ith entry in solution vector x?
j is given by
(x∗j)i = q′jNiqj, i = 1,2, . . . , n, j = 1,2, . . . , r
![Page 23: Didier HENRION LAAS-CNRS Toulouse, FR Czech Tech Univ ...cerfacs.fr/wp-content/uploads/2016/03/henrion.pdf · as a basis for polynomials of degree at most k ... 2.7603 5.1073 6.5115](https://reader033.vdocument.in/reader033/viewer/2022050410/5f86f5073f98be5f95111c28/html5/thumbnails/23.jpg)
Number of solutions
No easy way to control number of extracted solutions
in case of multiple global optima
Number of solutions = rank of moment matrix, but enforcing
rank in an LMI is a difficult non-convex problem
By default GloptiPoly minimizes the trace (sum of eigenvalues)
of the moment matrix, which may indirectly minimize the rank
(number of non-zero eigenvalues)
Practical experiments reveal that low rank moment matrices
ensure faster convergence of LMI relaxations to global optimum
![Page 24: Didier HENRION LAAS-CNRS Toulouse, FR Czech Tech Univ ...cerfacs.fr/wp-content/uploads/2016/03/henrion.pdf · as a basis for polynomials of degree at most k ... 2.7603 5.1073 6.5115](https://reader033.vdocument.in/reader033/viewer/2022050410/5f86f5073f98be5f95111c28/html5/thumbnails/24.jpg)
Example of extraction
Quadratic problem
p∗ = maxx (x1 − 1)2 + (x1 − x2)2 + (x2 − 3)2
s.t. (x1 − 1)2 ≤ 1(x1 − x2)
2 ≤ 1(x2 − 3)2 ≤ 1
First LMI relaxation yields p?1 = 3 and rank M1(y
∗) = 3,extraction algorithm fails due to incomplete monomial basis
Second LMI relaxation yields p?2 = 2 and
rank M1(y∗) = rank M2(y
∗) = 3
so rank condition ensures global optimality
![Page 25: Didier HENRION LAAS-CNRS Toulouse, FR Czech Tech Univ ...cerfacs.fr/wp-content/uploads/2016/03/henrion.pdf · as a basis for polynomials of degree at most k ... 2.7603 5.1073 6.5115](https://reader033.vdocument.in/reader033/viewer/2022050410/5f86f5073f98be5f95111c28/html5/thumbnails/25.jpg)
Cholesky factor
Moment matrix of order k = 2 reads
M2(y∗) =
1.0000 1.5868 2.2477 2.7603 3.6690 5.23871.5868 2.7603 3.6690 5.1073 6.5115 8.82452.2477 3.6690 5.2387 6.5115 8.8245 12.70722.7603 5.1073 6.5115 9.8013 12.1965 15.99603.6690 6.5115 8.8245 12.1965 15.9960 22.10845.2387 8.8245 12.7072 15.9960 22.1084 32.1036
Positive semidefinite with rank 3
Cholesky factor
V =
−0.9384 −0.0247 0.3447−1.6188 0.3036 0.2182−2.2486 −0.1822 0.3864−2.9796 0.9603 −0.0348−3.9813 0.3417 −0.1697−5.6128 −0.7627 −0.1365
satisfies
Mk(y?k) = V V ′
![Page 26: Didier HENRION LAAS-CNRS Toulouse, FR Czech Tech Univ ...cerfacs.fr/wp-content/uploads/2016/03/henrion.pdf · as a basis for polynomials of degree at most k ... 2.7603 5.1073 6.5115](https://reader033.vdocument.in/reader033/viewer/2022050410/5f86f5073f98be5f95111c28/html5/thumbnails/26.jpg)
Column echelon form
Gaussian elimination on V yields
U =
1 0 00 1 00 0 1
−2 3 0−4 2 2−6 0 5
1x1x2x21
x1x2x22
1 x1 x2
which means that solutions to be extracted satisfythe following system of polynomial equations
x21 = −2 + 3x1
x1x2 = −4 + 2x1 + 2x2x22 = −6 + 5x2
![Page 27: Didier HENRION LAAS-CNRS Toulouse, FR Czech Tech Univ ...cerfacs.fr/wp-content/uploads/2016/03/henrion.pdf · as a basis for polynomials of degree at most k ... 2.7603 5.1073 6.5115](https://reader033.vdocument.in/reader033/viewer/2022050410/5f86f5073f98be5f95111c28/html5/thumbnails/27.jpg)
Extraction
Multiplication matrices of monomials x1 and x2 in polynomial basis 1, x1, x2
are extracted from U :
N1 =
0 1 0−2 3 0−4 2 2
, N2 =
0 0 1−4 2 2−6 0 5
Random linear combination
N = 0.6909N1 + 0.3091N2
Schur decomposition of N = QTQ′ yields
Q =
0.4082 0.1826 −0.89440.4082 −0.9129 −0.00000.8165 0.3651 0.4472
Projections of orthogonal columns of Q onto N yield the 3 expectedglobally optimal solutions
x∗1 =
[12
]x∗2 =
[22
]x∗3 =
[23
]
![Page 28: Didier HENRION LAAS-CNRS Toulouse, FR Czech Tech Univ ...cerfacs.fr/wp-content/uploads/2016/03/henrion.pdf · as a basis for polynomials of degree at most k ... 2.7603 5.1073 6.5115](https://reader033.vdocument.in/reader033/viewer/2022050410/5f86f5073f98be5f95111c28/html5/thumbnails/28.jpg)
Outline
1. Quasi-Hankel matrices
2. Generalized companion matrices
3. Systems control applications
![Page 29: Didier HENRION LAAS-CNRS Toulouse, FR Czech Tech Univ ...cerfacs.fr/wp-content/uploads/2016/03/henrion.pdf · as a basis for polynomials of degree at most k ... 2.7603 5.1073 6.5115](https://reader033.vdocument.in/reader033/viewer/2022050410/5f86f5073f98be5f95111c28/html5/thumbnails/29.jpg)
Robust stability
Robust stability of linear system
x(t) = A(q)x(t)
where uncertain parameter q belongs to a compact set Q
Characteristic polynomial
p(s) = det(sI −A(q)) =∑k
skpk(q)
has roots in the open left half-plane iff
H(p) =∑i
∑j
pi(q)pj(q)Hij � 0
with Hermite matrix H(p) (Bezoutian)
![Page 30: Didier HENRION LAAS-CNRS Toulouse, FR Czech Tech Univ ...cerfacs.fr/wp-content/uploads/2016/03/henrion.pdf · as a basis for polynomials of degree at most k ... 2.7603 5.1073 6.5115](https://reader033.vdocument.in/reader033/viewer/2022050410/5f86f5073f98be5f95111c28/html5/thumbnails/30.jpg)
Robust stability
Assuming nominal system stable (q = 0), uncertain system
remains stable if and only if the multivariate polynomial
detH(p) =∑α
qαhα
remains positive for all q ∈ Q
General framework of µ-analysis (Doyle, Safonov)
Assessing robust stability of linear systems
amounts to polynomial minimisation
over compact sets (typically the unit box)
![Page 31: Didier HENRION LAAS-CNRS Toulouse, FR Czech Tech Univ ...cerfacs.fr/wp-content/uploads/2016/03/henrion.pdf · as a basis for polynomials of degree at most k ... 2.7603 5.1073 6.5115](https://reader033.vdocument.in/reader033/viewer/2022050410/5f86f5073f98be5f95111c28/html5/thumbnails/31.jpg)
Gain scheduling
Linear system
z(s) =b(s, q)
a(s, q)u(s)
whose transfer function depends on
an exogeneous parameter q ∈ Q
Find a linear controller
u(s) =y(s, q)
x(s, q)z(s)
also depending on q (scheduled in q) and ensuring
closed-loop stability and performance (H2 or H∞ norm)
![Page 32: Didier HENRION LAAS-CNRS Toulouse, FR Czech Tech Univ ...cerfacs.fr/wp-content/uploads/2016/03/henrion.pdf · as a basis for polynomials of degree at most k ... 2.7603 5.1073 6.5115](https://reader033.vdocument.in/reader033/viewer/2022050410/5f86f5073f98be5f95111c28/html5/thumbnails/32.jpg)
Gain scheduling
Controller design can be reformuled as
a parametrized PMI problem
F (k, q) � 0
where k is a variable to be found (parametrizing the controller)
and q can be anywhere in Q
Difficult semi-infinite non-convex optimisation problem for which
a hierarchy of finite-dimensional LMI problems can be devised
![Page 33: Didier HENRION LAAS-CNRS Toulouse, FR Czech Tech Univ ...cerfacs.fr/wp-content/uploads/2016/03/henrion.pdf · as a basis for polynomials of degree at most k ... 2.7603 5.1073 6.5115](https://reader033.vdocument.in/reader033/viewer/2022050410/5f86f5073f98be5f95111c28/html5/thumbnails/33.jpg)
Aircraft turbofan engines
Research contract with Snecma since 2000
Inputs: fuel flow WF32, nozzle area A8Outputs: fan speed XN2, compressor speed XN25Scheduling parameter q: pressure in combustion chamber
![Page 34: Didier HENRION LAAS-CNRS Toulouse, FR Czech Tech Univ ...cerfacs.fr/wp-content/uploads/2016/03/henrion.pdf · as a basis for polynomials of degree at most k ... 2.7603 5.1073 6.5115](https://reader033.vdocument.in/reader033/viewer/2022050410/5f86f5073f98be5f95111c28/html5/thumbnails/34.jpg)
Main points
Convex optimisation can be used to solve numerically non-convex
optimisation problems, specifically polynomial matrix inequalities
Numerical linear algebra is a key ingredient
Applications in many branches of engineering and applied maths,
in particular systems control
![Page 35: Didier HENRION LAAS-CNRS Toulouse, FR Czech Tech Univ ...cerfacs.fr/wp-content/uploads/2016/03/henrion.pdf · as a basis for polynomials of degree at most k ... 2.7603 5.1073 6.5115](https://reader033.vdocument.in/reader033/viewer/2022050410/5f86f5073f98be5f95111c28/html5/thumbnails/35.jpg)
Current limitations
In general, no bounds available on the size of the LMI problem
to be solved - we can construct worst-case exponential instances
No working implementation of SDP solver exploting the sparse
quasi-Hankel structure
Numerical stability of solution extraction mechanism remains
unclear = unknown sensitivity w.r.t. problem formulation