home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 10 / 10.iso / l / l455 / 7.ddi / ROBUST.DI$ / BILEXP.M < prev    next >
Encoding:
Text File  |  1993-03-11  |  970 b   |  31 lines
  1. %
  2. k = 1; m1 = 1; m2 = 1;
  3. agg = [0      0     1     0
  4.       0      0     0     1
  5.    -k/m1   k/m1    0     0
  6.     k/m2  -k/m2    0     0];
  7. bgg = [0 0 1/m1 0]'; cgg = [0 1 0 0]; dgg = 0;
  8. lamg = eig(agg);
  9. [ta1,ta2,ta3,ta4] = bilin(agg,bgg,cgg,dgg,1,'Sft_jw',[-12 -2]);
  10. lamtg = eig(ta1);
  11. xaxis = 2:0.2:12;
  12. [tmp1,tmp2]=size(xaxis);
  13. y = sqrt(25-(xaxis-7*ones(tmp1,tmp2)).^2);
  14. y = [y;-y];
  15. subplot(111);
  16. plot(xaxis,y,'--',-xaxis,y,'--',real(lamg),imag(lamg),'x',...
  17. real(lamtg),imag(lamtg),'+',-real(lamtg),imag(lamtg),'*',...
  18. -2*real(lamtg),imag(lamtg),'o')
  19. text(0.2,0.35,'1) x -- original poles','sc')
  20. text(0.2,0.3,'2) + -- shifted poles','sc')
  21. text(0.2,0.25,'3) * -- H-Inf closed-loop poles in s~-plane','sc')
  22. text(0.2,0.20,'4) o -- H-Inf closed-loop poles in s-plane','sc')
  23. text(-3,0,'p1')
  24. text(3,0,'-p1')
  25. text(-13,0,'p2')
  26. text(13,0,'-p2')
  27. xlabel('(p1 = -2, p2 = -12)')
  28. grid
  29. title('Example of Bilinear Mapping: s~ = (-s + p1) / (s/p2 -1)')
  30. pause
  31.