home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 10 / 10.iso / l / l455 / 7.ddi / ROBUST.DI$ / MUDEMO1.M < prev    next >
Encoding:
Text File  |  1993-03-11  |  898 b   |  39 lines
  1. clc
  2. % MUSYNEX is a simple demo routine illustrating the
  3. % D-F iteration of mu-synthesis.  See also MUDEMO.
  4.  
  5. echo on
  6.  
  7. % Plant Data:
  8.  
  9. a=2;             b1=[.1 -1];          b2=-1;
  10. c1=[1 .01]';    d11=[.1 .2;.01 .01]; d12=[1 0]';
  11. c2=1;           d21=[0  1];          d22=3;
  12. tss_=mksys(a,b1,b2,c1,c2,d11,d12,d21,d22,'tss');
  13. w=logspace(-2,1);
  14.  
  15.  
  16. % H-infinity 
  17.  
  18. [gam0,ss_cp0,ss_cl0]=hinfopt(tss_);
  19. [mu0,logd0]=ssv(ss_cl0,w);
  20.  
  21. % Mu Synthesis -- Iteration No. 1 (CONSTANT D(s)):
  22.  
  23. [ss_d1,logd1]=fitd(logd0,w);  
  24. [gam1,ss_cp1,ss_cl1]=hinfopt(augd(tss_,ss_d1));
  25. [mu1,deltalogd]=ssv(ss_cl1,w);
  26.  
  27.  
  28.  
  29. % Mu Synthesis -- Iteration No. 2 (FIRST ORDER D(s)):
  30.  
  31. [ss_d2,logd2]=fitd(logd1+deltalogd,w,1);
  32. [gam2,ss_cp2,ss_cl2]=hinfopt(augd(tss_,ss_d2));
  33.  
  34. % Display Optimal hinf cost closed loop and SSV plot
  35.  
  36. clg
  37. loglog(w,max(sigma(ss_cl2,1,w))/gam2,w,ssv(ss_cl2,w)/gam2);
  38. echo off
  39. pause