Few modifications

Include the script for generating the feasible using brute force

Author: emammendes

Diff for: actual_build_lmi15_jianbo.m

@@ -8,12 +8,14 @@
 
 auxlmi15=[];
 
+auxm15=m15;
+
 for j1=1:3 % linhas 
     s=[];
     for j2=1:3 % colunas
         s=[s eval(m15{j1,j2})];
-    end;
+    end
     auxlmi15=[auxlmi15;s];
-end;
+end
 
 biglmi15= (auxlmi15>= 0);

Diff for: actual_build_lmi16_jianbo.m

@@ -31,10 +31,10 @@
             s='[';
             for j2=1:(M+1)
                 s=sprintf('%s m16{%d,%d}',s,j1,j2);
-            end;
+            end
             s=sprintf('%s];',s);
             auxlmi16=[auxlmi16;eval(s)];
-        end;
+        end
 
         biglmi16=[biglmi16;auxlmi16>=0];
     end

Diff for: actual_build_lmi20_jianbo.m

@@ -5,13 +5,6 @@
 
 % LMI (20) 
 
-% Build LMI 21
-
-% Rosileide Lopes - 03/13/2018
-% DELT - UFMG
-
-% LMI (21) 
-
 
 biglmi20=[];
 

Diff for: define_lmi_variables_jianbo.m

@@ -8,7 +8,7 @@
 
 gamma_1 = sdpvar(1);
 gamma_2 = sdpvar(1);
-uk      = sdpvar(1);assign(uk,0);   % Old - assign(uk,CIu);
+uk      = sdpvar(1);   % Old - assign(uk,CIu);
 
 % To build the other LMI variables we need to know the dimensions of the system and
 % control matrices.  This can be easily achieved by looking at the defintions

Diff for: script_generate_feasible_region.m

@@ -0,0 +1,98 @@
+%% Script - An attempt (with failure) to use yalmip plot to get the feasible region
+
+clear all;close all;clc;
+
+%% Call yalmpi and solvers - This is an example - I am running Matlab on a mac.
+
+addpath(genpath('~/Documents/MATLAB/yalmip'))
+addpath(genpath('~/Documents/MATLAB/cvx/sedumi'))
+addpath(genpath('~/Documents/MATLAB/cvx/sdpt3'))
+addpath(genpath('~/mosek/8/toolbox/r2014a'));
+
+%% Define Constants
+
+flagc=0;
+
+disp('***********************************************');
+
+if flagc == 1
+    disp('A Constrained problem will be solved');
+    disp('');
+    disp('LMIs 11, 15, 16, 20, 21, 23 and 25 will be used');
+elseif flagc == 2
+    disp('A Constrained problem will be solved');
+    disp('');
+    disp('LMIs 11, 15, 16, 20, 21, 23, 24, 25 and 26 will be used');
+else
+    disp('An unconstrained problem will be solved');
+    disp('');
+    disp('LMIs 11, 15 and 16 will be used');
+end
+
+opts = sdpsettings('verbose',0, 'warning',1,'solver','mosek');
+
+model_parameters_jianbo
+
+nq=nx;
+nr=nu;
+
+umax=1;
+xmax=1;
+
+rk=1;
+
+N=1;
+
+define_lmi_variables_jianbo
+
+xk=sdpvar(2,1);
+
+Obj=gamma_1+gamma_2;
+
+%% LMI 11
+
+pre_build_lmi11_jianbo
+actual_build_lmi11_jianbo
+
+%% LMI 15
+
+pre_build_lmi15_jianbo
+actual_build_lmi15_jianbo
+
+%% LMI 16
+
+pre_build_lmi16_jianbo
+actual_build_lmi16_jianbo
+
+%% Big LMI
+
+if flagc == 1 % Constrained on the input
+    pre_build_lmi20_jianbo
+    actual_build_lmi20_jianbo
+    pre_build_lmi21_jianbo
+    actual_build_lmi21_jianbo
+    pre_build_lmi23_jianbo
+    actual_build_lmi23_jianbo
+    pre_build_lmi25_jianbo
+    actual_build_lmi25_jianbo
+    LMIs_orig=[biglmi11;biglmi15;biglmi16;biglmi20;biglmi21;biglmi23;biglmi25,biglmi25_1];
+elseif flagc == 2 % Constrained on the input and states
+    pre_build_lmi20_jianbo
+    actual_build_lmi20_jianbo
+    pre_build_lmi21_jianbo
+    actual_build_lmi21_jianbo
+    pre_build_lmi23_jianbo
+    actual_build_lmi23_jianbo
+    pre_build_lmi24_jianbo
+    actual_build_lmi24_jianbo
+    pre_build_lmi25_jianbo
+    actual_build_lmi25_jianbo
+    pre_build_lmi26_jianbo
+    actual_build_lmi26_jianbo
+    LMIs_orig=[biglmi11;biglmi15;biglmi16;biglmi20;biglmi21;biglmi23;biglmi24;biglmi25,biglmi25_1;biglmi26,biglmi26_1];
+else % unconstrained
+    LMIs_orig=[biglmi11;biglmi15;biglmi16];
+end
+
+
+plot(LMIs_orig,xk,[],[],opts);