home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 10 / 10.iso / l / l455 / 7.ddi / ROBUST.DI$ / DINTEVA.M < prev    next >
Encoding:
Text File  |  1993-03-11  |  1.2 KB  |  49 lines
  1. %
  2. %  DINTEVA.M: Preparing singular value plots 
  3. %
  4.  
  5. % R. Y. Chiang & M. G. Safonov 
  6. % Copyright (c) 1988 by the MathWorks, Inc.
  7. % All Rights Reserved.
  8.  
  9. w = logspace(-3,3);
  10. %
  11. disp(' ')
  12. disp(' - - Computing the SV Bode plot of the plant open loop - -');
  13. svg = sigma(ag,bg,cg,dg,1,w);
  14. svg = 20*log10(svg);
  15. %
  16. disp(' ')
  17. disp(' - - Computing the SV Bode plot of Ty1u1 - -');
  18. svtt = sigma(acl,bcl,ccl,dcl,1,w); 
  19. svtt = 20*log10(svtt);
  20. %
  21. svw1i = bode(w1(2,:),w1(1,:),w); svw1i = 20*log10(svw1i);
  22. svw3i = bode(w3(2,:),w3(1,:),w); svw3i = 20*log10(svw3i);
  23. %
  24. disp(' ')
  25. disp(' - - Computing the SV Bode plot of Controller - -')
  26. svcp = sigma(acp,bcp,ccp,dcp,1,w); svcp = 20*log10(svcp);
  27. %
  28. [al,bl,cl,dl] = series(acp,bcp,ccp,dcp,ag,bg,cg,dg);
  29. %
  30. [as,bs,cs,ds] = feedbk(al,bl,cl,dl,1);
  31. disp(' ')
  32. disp(' - - Computing the SV Bode plot of S - - ')
  33. svs = sigma(as,bs,cs,ds,1,w);
  34. %[svs,temp] = sort(-svs); svs = -svs;     
  35. svs = 20*log10(svs);      
  36. %
  37. [at,bt,ct,dt] = feedbk(al,bl,cl,dl,2);
  38. disp(' ')
  39. disp(' - - Computing the SV Bode plot of I-S - - ')
  40. svt = sigma(at,bt,ct,dt,1,w);
  41. svt = 20*log10(svt);
  42. %
  43. disp(' ')
  44. disp(' - - Computing the time response - - ')
  45. t = 0:0.1:5;
  46. y = step(at,bt,ct,dt,1,t);
  47. %
  48. % -------- End of DINTEVA.M --- RYC/MGS %
  49.