home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 10 / 10.iso / l / l455 / 7.ddi / ROBUST.DI$ / EVALSV.M < prev    next >
Encoding:
Text File  |  1993-03-11  |  1.8 KB  |  53 lines
  1. %
  2. %  EVALSV.M: Preparing singular value plots 
  3. %
  4. w = logspace(-3,3,100);
  5. %
  6. disp(' ')
  7. disp(' - - Computing the SV Bode plot of the plant open loop - -');
  8. svg = sigma(ag,bg,cg,dg,1,w);
  9. svg = 20*log10(svg);
  10. %
  11. disp(' ')
  12. disp(' - - Computing the SV Bode plot of Ty1u1 - -');
  13. svtt = sigma(acl,bcl,ccl,dcl,1,w); 
  14. %[svtt,temp] = sort(-svtt); svtt = -svtt;
  15. svtt = 20*log10(svtt);
  16. %
  17. svw1i = bode(1/Gam*w1(2,:),w1(1,:),w); svw1i = 20*log10(svw1i);
  18. svw3i = bode(w3(2,:),w3(1,:),w); svw3i = 20*log10(svw3i);
  19. %
  20. [al,bl,cl,dl] = series(acp,bcp,ccp,dcp,ag,bg,cg,dg);
  21. %
  22. [as,bs,cs,ds] = feedbk(al,bl,cl,dl,1);
  23. disp(' ')
  24. disp(' - - Computing the SV Bode plot of S - - ')
  25. svs = sigma(as,bs,cs,ds,1,w);
  26. %[svs,temp] = sort(-svs); svs = -svs;     
  27. svs = 20*log10(svs);      
  28. %
  29. [at,bt,ct,dt] = feedbk(al,bl,cl,dl,2);
  30. disp(' ')
  31. disp(' - - Computing the SV Bode plot of I-S - - ')
  32. svt = sigma(at,bt,ct,dt,1,w);
  33. %[svt,temp] = sort(-svt); svt = -svt;
  34. svt = 20*log10(svt);
  35. %
  36. [ali,bli,cli,dli] = series(ag,bg,cg,dg,acp,bcp,ccp,dcp);
  37. disp(' ');
  38. disp(' - - Computing the SV Bode plot of T at plant input - -');
  39. svti = sigma(ali,bli,cli,dli,4,w); 
  40. svti = -20*log10(svti);
  41. %
  42. %disp('                 << Plot variable names >>');
  43. %disp('---------------------------------------------------------------');
  44. %disp(' w ------ freqeuncy (rad/sec)');
  45. %disp(' svs ---- singular values of S(jw) (sensitivity)');
  46. %disp(' svt ---- singular values of I-S(jw) (complementary sensitivity)');
  47. %disp(' svtt --- singular values of the cost function (Ty1u1(jw))');
  48. %disp(' svw1i --- singular values of 1/W1(jw) weighting function');
  49. %disp(' svw3i --- singular values of 1/W3(jw) weighting function');
  50. %disp('---------------------------------------------------------------');
  51. %
  52. % -------- End of EVALSV.M --- RYC/MGS %
  53.