home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 10 / 10.iso / l / l455 / 7.ddi / ROBUST.DI$ / MRDEMO.M < prev    next >
Encoding:
Text File  |  1993-03-11  |  4.1 KB  |  117 lines
  1. %
  2. % MRDEMO Robust model reduction demonstration
  3. %
  4.  
  5. % R. Y. Chiang & M. G. Safonov 3/88
  6. % Copyright (c) 1988 by the MathWorks, Inc.
  7. % All Rights Reserved.
  8. % ----------------------------------------------------------------
  9. clc
  10. disp(' ')
  11. disp(' ')
  12. disp('          <<<  Robust Basis Free Model Reduction Techniques Demo  >>>') 
  13. disp('    ')
  14. disp('   ')
  15. disp('                 1). Balanced model reduction techniques')
  16. disp('                     (Schur method vs. Square-root method)')
  17. disp(' ')
  18. disp('                 2). Optimal Descriptor Hankel model reduction & ')
  19. disp('                     anticausal projection w/o balancing')
  20. disp(' ') 
  21. disp('                 3). Schur Balanced Stochastic Trunction /')
  22. disp('                     Relative Error Method')
  23. disp(' ')
  24. disp('                 4). Comparison of 1, 2, and 3.')
  25. disp('     ')
  26. disp('                 5). Go to the main menu')
  27. disp(' ')
  28. disp('                 0). Quit ..')
  29. disp('  ')
  30. Demono = input('      Select one of the above options: ');
  31. % ----------------------------------------------------------------------
  32. if Demono == 0
  33.    clc
  34.    disp(' ')
  35.    disp(' ')
  36.    disp('  ')
  37.    disp('  ')
  38.    disp('                   * * * * * * * * * * * * * * * * *')
  39.    disp('                   *  End of MRDEMO  ......        *')
  40.    disp('                   *                               *')
  41.    disp('                   * Good luck with your design !! *')
  42.    disp('                   * * * * * * * * * * * * * * * * *')
  43.    break
  44. end
  45. if Demono == 5
  46.    rctdemo1
  47. end
  48. % ----------------------------------------------------------------------
  49. clc
  50. disp(' ' )
  51. disp('                  << Motivation of Robust Model Reduction >>')
  52. disp('    ')
  53. disp('    Robust model reduction techniques have become more and more important.')
  54. disp('  ')
  55. disp('    This is due to the following:')
  56. disp('     ')
  57. disp('          1). Complexity of the design problem (e.g. large size plant)')
  58. disp('  ')
  59. disp('          2). Complexity of the design algorithm')
  60. disp('  ')
  61. disp('          3). Large size of the controller generated by modern control')
  62. disp('              techniques -- e.g. LQG, LQG/LTR, H2, H-inf ..etc..')
  63. disp('    ')
  64. disp('    Model reduction algorithms must be both numerically robust and')
  65. disp(' ')
  66. disp('    be able to address closed-loop "robustness" issues.')
  67. disp('   ')
  68. disp('  ')
  69. disp('                                 (strike a key to continue ...)')
  70. pause
  71. clc
  72. disp('    The model reduction algorithms included here have the following features:')
  73. disp(' ')
  74. disp('       1). They all by-pass the ill-conditioned "balanced transformation"')
  75. disp(' ')
  76. disp('       2). They employ Schur decomposition as an intermediate step')
  77. disp('    ')
  78. disp('           for robust computation of orthonormal bases for eigenspaces.')
  79. disp('    ')
  80. disp('       3). They all have infinity-norm error bounds:')
  81. disp('             ')
  82. disp('                Methods                     Infinity-Norm Error Bounds')
  83. disp('     ----------------------------         -------------------------------')
  84. disp('      Schur/Square-root Balanced           twice of the tails of Hankel')
  85. disp('           (schur.m & balmr.m)             singular values')
  86. disp('  ')
  87. disp('      Optimal Descriptor Hankel                  same as above')
  88. disp('      without balancing (ohklmr.m)')
  89. disp('                                            -1   ^            n    hsv')
  90. disp('      Balanced Stochastic Truncation/     |G  (G-G) |   <= 2 Sum -------')
  91. disp('      Relative Error Method (reschmr.m)              inf     k+1  1-hsv')
  92. disp('  ')
  93. disp('                                 (strike a key to continue ...)')
  94. pause
  95. % ---------------------------------------------------------------------
  96. if Demono == 1
  97.    baldemo
  98.    mrdemo
  99. end
  100. % ----------------------------------------------------------------------
  101. if Demono == 2
  102.    ohkdemo
  103.    mrdemo
  104. end
  105. % ----------------------------------------------------------------------
  106. if Demono == 3
  107.    remdemo
  108.    mrdemo
  109. end
  110. % ----------------------------------------------------------------------
  111. if Demono == 4
  112.    bhrdemo
  113.    mrdemo
  114. end
  115. %
  116. % --------- End of MRDEMO.M --- RYC/MGS %
  117.