Why LMI does not returne a positive matrix?

3 Ansichten (letzte 30 Tage)
Amine
Amine am 30 Mai 2015
Kommentiert: Amine am 4 Mär. 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 Kommentare
Michael Hubatka
Michael Hubatka am 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 am 4 Mär. 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)

Melden Sie sich an, um zu kommentieren.

Antworten (0)

Kategorien

Mehr zu Linear Matrix Inequalities finden Sie in Help Center und File Exchange

Community Treasure Hunt

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

Start Hunting!

Translated by