Why LMI does not returne a positive matrix?

2 views (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)

Categories

Find more on Linear Matrix Inequalities in Help Center and File Exchange

Community Treasure Hunt

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

Start Hunting!