Why LMI does not returne a positive matrix?

1 view (last 30 days)
Amine
Amine on 30 May 2015
Commented: Amine on 4 Mar 2016
Hello! I have the following LMI problem: AX-BL+XA'-L'B'>0 And X>0 Where: P = inv(X) And: A, P and X are 4x4 matrices, B is 4x1 and L is 1x4 vectors. I wrote the following script. But P is not positive!? Is it a problem in my script or something else? Thanks in advance!
%%LMIs :
% LMI variables :
setlmis([]);
X = lmivar(1, [4 1]);
L = lmivar(2, [1 4]);
% First LMI :
lmiterm([-1 1 1 X],A,1,'s');
lmiterm([ 1 1 1 L],B,1,'s');
% Second LMI :
lmiterm([2 1 1 X],1,1);
% Get LMIs :
lmis = getlmis;
% P matrix:
[Tmin,Xfeas] = feasp(lmis)
X = dec2mat(lmis,Xfeas,X);
p = inv(X)
  4 Comments
Michael Hubatka
Michael Hubatka on 5 Feb 2016
The condition X > 0 is defined as
lmiterm([-2 1 1 X], 1, 1);
or
lmiterm([2 1 1 -X], 1, 1);
because lmiterm always uses the '<' condition.
Amine
Amine on 4 Mar 2016
First of all thanks for your interest in my question. Both code lines give the same result a negative P matrix (P > 0) where :
P = inv(X)

Sign in to comment.

Answers (0)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!