Convert Riccati matrix inequalities into LMI

40 views
Skip to first unread message

stevens

unread,
Jun 28, 2018, 5:24:38 AM6/28/18
to YALMIP
Hey,

i am quite new in the topic LMIs... That's why I would like to ask, if someone can help me, solving these riccati matrix inequalities....
Can I form them into LMIs? 

Thank you very much!
Best regards,
Steven
Unbenannt.JPG

Johan Löfberg

unread,
Jun 28, 2018, 5:59:08 AM6/28/18
to YALMIP
yes, schur complements

stevens

unread,
Jun 28, 2018, 8:23:04 AM6/28/18
to YALMIP
Thank you very much  for that hint! 
With the schur complement I got these three LMIs. Is that correct?

so P1 and P2 should be symmetric?

Is it possible to define an equality (==) for the third equation in yalmip? 

Best regards,
Steven


Unbenannt.JPG

Johan Löfberg

unread,
Jun 28, 2018, 8:53:12 AM6/28/18
to YALMIP
You still have quadratic terms in P so you will have to apply a Schur complements once again

However, the equality is nonlinear and nonconvex, and cannot be linearized, so you have a pretty bad problem. Equalities are just == in YALMIP, but this is a nonlinear equality and hence an LMI solver will not be able to deal with this

stevens

unread,
Jun 29, 2018, 8:29:08 AM6/29/18
to YALMIP
Thank you very much Johan Löfberg!

stevens

unread,
Jun 29, 2018, 8:40:58 AM6/29/18
to YALMIP
Sorry, but I have a further question to a similar topic. I try to solve these three LMIs (22), (26a) and (26b) with the unkown matrices: X, Y, T, Wb, Wc, W_.
Unfortunately, I don't know how to handle the inequality W_ >= T*Y*T' (highlighted), because it is nonlinear.
Is there a possibility to implement this inequality?

Thank you for your support!
Best regards,
Steven


LMI_22.JPG
LMI_26.JPG

Johan Löfberg

unread,
Jun 29, 2018, 9:24:13 AM6/29/18
to YALMIP
W>=TYT is not possible to represent if you have all those variables as decision variables. Hence, I suspect you misunderstand some notation , and 22 is  unclear what they mean.  As the expression is non-symmetric, you would have to introduce a new variable Z, and add the constraints Z == [E 0;0 E]*[T I;I X], Z>=0
Reply all
Reply to author
Forward
0 new messages