Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- %%
- %preg 4b
- cl;
- A1=[0 5; -30 -1];
- A2=[0 5; -26 -1];
- A3=[-6 27; -150 -1];
- n=3;
- emptyCell = cell(1,10)
- emptyCell={A1,A2,A3}
- P=fun_4b(emptyCell)
- function [myP] = fun_4b(cel)
- nn = size(cel{1}, 1); % get matrix dimension
- setlmis([]); % reset lmi-system
- P = lmivar(1,[nn,1]); % create matrix-valued variable (symmetric block-diagonal, one block of size n)
- lmiterm([-1 1 1 P], 1,1);
- lmiterm([ 1 1 1 0],1);
- n=size(cel,2);
- for i=1:n
- lmiterm([ i+1 1 1 P], cel{i}',1,'s')
- cel{i}
- i
- end
- % define LMI system
- mylmi = getlmis;
- % solve LMI
- [tmin,xfeas]= feasp(mylmi);
- myP = dec2mat(mylmi, xfeas, P);
- %lmiterm([ 2 1 1 P], A',1,'s');
- %lmiterm([ 3 1 1 P], B',1,'s');
- end
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement