Advertisement
Guest User

Untitled

a guest
Jan 21st, 2017
102
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
MatLab 0.72 KB | None | 0 0
  1. %%
  2. %preg 4b
  3. cl;
  4. A1=[0 5; -30 -1];
  5. A2=[0 5; -26 -1];
  6. A3=[-6 27; -150 -1];
  7.  
  8. n=3;
  9. emptyCell = cell(1,10)
  10. emptyCell={A1,A2,A3}
  11.  
  12. P=fun_4b(emptyCell)
  13.  
  14.  
  15. function [myP] = fun_4b(cel)
  16.  
  17. nn = size(cel{1}, 1);       % get matrix dimension
  18. setlmis([]);          % reset lmi-system
  19.  
  20. P = lmivar(1,[nn,1]);  % create matrix-valued variable (symmetric block-diagonal, one block of size n)
  21.  
  22. lmiterm([-1 1 1 P], 1,1);
  23. lmiterm([ 1 1 1 0],1);
  24.  
  25. n=size(cel,2);
  26.  
  27. for i=1:n
  28.     lmiterm([ i+1 1 1 P], cel{i}',1,'s')
  29.     cel{i}
  30.     i
  31. end
  32.  
  33. % define LMI system
  34. mylmi = getlmis;
  35.  
  36. % solve LMI
  37. [tmin,xfeas]= feasp(mylmi);
  38.  
  39. myP = dec2mat(mylmi, xfeas, P);
  40. %lmiterm([ 2 1 1 P], A',1,'s');
  41. %lmiterm([ 3 1 1 P], B',1,'s');
  42.  
  43. end
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement