home *** CD-ROM | disk | FTP | other *** search
/ OS/2 Shareware BBS: 10 Tools / 10-Tools.zip / octa21eb.zip / octave / SCRIPTS.ZIP / scripts.fat / control / rldemo.m < prev    next >
Text File  |  1999-04-29  |  10KB  |  297 lines
  1. # Copyright (C) 1996 A. Scottedward Hodel 
  2. #
  3. # This file is part of Octave. 
  4. #
  5. # Octave is free software; you can redistribute it and/or modify it 
  6. # under the terms of the GNU General Public License as published by the 
  7. # Free Software Foundation; either version 2, or (at your option) any 
  8. # later version. 
  9. # 
  10. # Octave is distributed in the hope that it will be useful, but WITHOUT 
  11. # ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or 
  12. # FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License 
  13. # for more details.
  14. # 
  15. # You should have received a copy of the GNU General Public License 
  16. # along with Octave; see the file COPYING.  If not, write to the Free 
  17. # Software Foundation, 675 Mass Ave, Cambridge, MA 02139, USA. 
  18.  
  19. function rldemo()
  20. # Octave Controls toolbox demo: Root Locus demo
  21. # Written by David Clem August 15, 1994
  22. # Updated by John Ingram December 1996
  23.  
  24.   while (1)
  25.     clc
  26.     k = menu("Octave Root Locus Demo", ...
  27.         "Display continuous system's open loop poles and zeros (pzmap)", ...
  28.         "Display discrete system's open loop poles and zeros (pzmap)", ...
  29.         "Display root locus diagram of SISO continuous system (rlocus)", ...
  30.         "Display root locus diagram of SISO discrete system (rlocus)", ...    
  31.         "Return to main demo menu");
  32.     gset autoscale
  33.     if (k == 1)
  34.       clc
  35.       help pzmap
  36.       prompt
  37.  
  38.       clc
  39.       disp("Display continuous system's open loop poles and zeros (pzmap)\n");
  40.       disp("Example #1, Consider the following continuous transfer function:");
  41.       cmd = "sys1 = tf2sys([1.5 18.5 6],[1 4 155 302 5050]);";
  42.       disp(cmd);
  43.       eval(cmd);
  44.       cmd ="sysout(sys1);";
  45.       disp(cmd);
  46.       eval(cmd);
  47.       disp("\nPole-zero form can be obtained as follows:");
  48.       cmd = "sysout(sys1,""zp"");";
  49.       disp(cmd);
  50.       eval(cmd);     
  51.       disp("View the system's open loop poles and zeros with the command:")
  52.       cmd = "pzmap(sys1);";
  53.       run_cmd
  54.       prompt     
  55.  
  56.       clc
  57.       disp("Example #2, Consider the following set of poles and zeros:");
  58.       cmd = "sys2 = zp2sys([-1 5 -23],[-1 -10 -7+5i -7-5i],5);";
  59.       disp(cmd);
  60.       eval(cmd);
  61.       cmd = "sysout(sys2);";
  62.       disp(cmd);
  63.       eval(cmd);
  64.       disp("\nThe pzmap command for the zp form is the same as the tf form:")
  65.       cmd = "pzmap(sys2);";
  66.       run_cmd;
  67.       disp("\nThe internal representation of the system is not important;");
  68.       disp("pzmap automatically sorts it out internally.");
  69.       prompt;
  70.  
  71.       clc
  72.       disp("Example #3, Consider the following state space system:\n");
  73.       cmd = "sys3=ss2sys([0 1; -10 -11], [0;1], [0 -2], 1);";
  74.       disp(cmd); 
  75.       eval(cmd);
  76.       cmd = "sysout(sys3);";
  77.       disp(cmd); 
  78.       eval(cmd);
  79.       disp("\nPole-zero form can be obtained as follows:");
  80.       cmd = "sysout(sys3,""zp"");";
  81.       disp(cmd);
  82.       eval(cmd); 
  83.       disp("\nOnce again, the pzmap command is the same:");
  84.       cmd = "pzmap(sys3);";
  85.       run_cmd;
  86.       prompt;
  87.  
  88.       closeplot
  89.       clc
  90.     
  91.     elseif (k == 2)
  92.       clc
  93.       help pzmap
  94.       prompt
  95.  
  96.       clc
  97.       disp("\nDisplay discrete system's open loop poles and zeros (pzmap)\n");
  98.       disp("First we must define a sampling time, as follows:\n");
  99.       cmd = "Tsam = 1;";
  100.       run_cmd;
  101.       disp("Example #1, Consider the following discrete transfer function:");
  102.       cmd = "sys1 = tf2sys([1.05 -0.09048],[1 -2 1],Tsam);";
  103.       disp(cmd);
  104.       eval(cmd);
  105.       cmd ="sysout(sys1);";
  106.       disp(cmd);
  107.       eval(cmd);
  108.       disp("\nPole-zero form can be obtained as follows:");
  109.       cmd = "sysout(sys1,""zp"");";
  110.       disp(cmd);
  111.       eval(cmd);     
  112.       disp("View the system's open loop poles and zeros with the command:")
  113.       cmd = "pzmap(sys1);";
  114.       run_cmd
  115.       prompt     
  116.  
  117.       clc
  118.       disp("Example #2, Consider the following set of discrete poles and zeros:");
  119.       cmd = "sys2 = zp2sys([-0.717],[1 -0.368],3.68,Tsam);";
  120.       disp(cmd);
  121.       eval(cmd);
  122.       cmd = "sysout(sys2);";
  123.       disp(cmd);
  124.       eval(cmd);
  125.       disp("\nThe pzmap command for the zp form is the same as the tf form:")
  126.       cmd = "pzmap(sys2);";
  127.       run_cmd;
  128.       disp("\nThe internal representation of the system is not important;");
  129.       disp("pzmap automatically sorts it out internally.");
  130.       prompt;
  131.  
  132.       clc
  133.       disp("Example #3, Consider the following discrete state space system:\n");
  134.       cmd = "sys3=ss2sys([1 0.0952;0 0.905], [0.00484; 0.0952], [1 0], 0, Tsam);";
  135.       disp(cmd); 
  136.       eval(cmd);
  137.       cmd = "sysout(sys3);";
  138.       disp(cmd); 
  139.       eval(cmd);
  140.       disp("\nPole-zero form can be obtained as follows:");
  141.       cmd = "sysout(sys3,""zp"");";
  142.       disp(cmd);
  143.       eval(cmd); 
  144.       disp("\nOnce again, the pzmap command is the same:");
  145.       cmd = "pzmap(sys3);";
  146.       run_cmd;
  147.       prompt;
  148.  
  149.       closeplot
  150.       clc
  151.  
  152.     elseif (k == 3)
  153.       clc
  154.       help rlocus
  155.       prompt;
  156.  
  157.       clc
  158.       disp("Display root locus of a continuous SISO system (rlocus)\n")
  159.       disp("Example #1, Consider the following continuous transfer function:");
  160.       cmd = "sys1 = tf2sys([1.5 18.5 6],[1 4 155 302 5050]);";
  161.       disp(cmd);
  162.       eval(cmd);
  163.       cmd ="sysout(sys1);";
  164.       disp(cmd);
  165.       eval(cmd);
  166.       disp("\nPole-zero form can be obtained as follows:");
  167.       cmd = "sysout(sys1,""zp"");";
  168.       disp(cmd);
  169.       eval(cmd); 
  170.       disp("\nWhen using rlocus, inital system poles are displayed as X's.")
  171.       disp("Moving poles are displayed as diamonds.  Zeros are displayed as")
  172.       disp("boxes.  The simplest form of the rlocus command is as follows:")
  173.       cmd = "rlocus(sys1);";
  174.       run_cmd
  175.       disp("\nrlocus automatically selects the minimum and maximum gains based")
  176.       disp("on the real-axis locus breakpoints.  The plot limits are chosen")
  177.       disp("to be no more than 10 times the maximum magnitude of the open")
  178.       disp("loop poles/zeros.");
  179.       prompt
  180.  
  181.       clc
  182.       disp("Example #2, Consider the following set of poles and zeros:");
  183.       cmd = "sys2 = zp2sys([],[0 -20 -2 -0.1],5);";
  184.       disp(cmd);
  185.       eval(cmd);
  186.       cmd = "sysout(sys2);";
  187.       disp(cmd);
  188.       eval(cmd);
  189.       disp("\nThe rlocus command for the zp form is the same as the tf form:")
  190.       cmd = "rlocus(sys2);";
  191.       run_cmd;
  192.       disp("\nThe internal representation of the system is not important;");
  193.       disp("rlocus automatically sorts it out internally.");
  194.       prompt;
  195.  
  196.       clc
  197.       disp("Example #3, Consider the following state space system:\n");
  198.       cmd = "sys3=ss2sys([0 1; -10 -11], [0;1], [0 -2], 0);";
  199.       disp(cmd); 
  200.       eval(cmd);
  201.       cmd = "sysout(sys3);";
  202.       disp(cmd); 
  203.       eval(cmd);
  204.       disp("\nPole-zero form can be obtained as follows:");
  205.       cmd = "sysout(sys3,""zp"");";
  206.       disp(cmd);
  207.       eval(cmd); 
  208.       disp("\nOnce again, the rlocus command is the same:");
  209.       cmd = "rlocus(sys3);";
  210.       run_cmd;
  211.  
  212.       disp("\nNo matter what form the system is in, the rlocus command works the");
  213.       disp("the same.");
  214.       prompt;
  215.  
  216.       closeplot
  217.       clc
  218.  
  219.     elseif (k == 4)
  220.       clc
  221.       help rlocus
  222.       prompt
  223.  
  224.       clc
  225.       disp("Display root locus of a discrete SISO system (rlocus)\n")
  226.       disp("First we must define a sampling time, as follows:\n");
  227.       cmd = "Tsam = 1;";
  228.       run_cmd;
  229.       disp("Example #1, Consider the following discrete transfer function:");
  230.       cmd = "sys1 = tf2sys([1.05 -0.09048],[1 -2 1],Tsam);";
  231.       disp(cmd);
  232.       eval(cmd);
  233.       cmd ="sysout(sys1);";
  234.       disp(cmd);
  235.       eval(cmd);
  236.       disp("\nPole-zero form can be obtained as follows:");
  237.       cmd = "sysout(sys1,""zp"");";
  238.       disp(cmd);
  239.       eval(cmd);     
  240.       disp("\nWhen using rlocus, inital system poles are displayed as X's.")
  241.       disp("Moving poles are displayed as diamonds.  Zeros are displayed as")
  242.       disp("boxes.  The simplest form of the rlocus command is as follows:")
  243.       cmd = "rlocus(sys1);";
  244.       run_cmd
  245.       disp("\nrlocus automatically selects the minimum and maximum gains based")
  246.       disp("on the real-axis locus breakpoints.  The plot limits are chosen")
  247.       disp("to be no more than 10 times the maximum magnitude of the open")
  248.       disp("loop poles/zeros.");
  249.       prompt
  250.  
  251.       clc
  252.       disp("Example #2, Consider the following set of discrete poles and zeros:");
  253.       cmd = "sys2 = zp2sys([-0.717],[1 -0.368],3.68,Tsam);";
  254.       disp(cmd);
  255.       eval(cmd);
  256.       cmd = "sysout(sys2);";
  257.       disp(cmd);
  258.       eval(cmd);
  259.       disp("\nThe rlocus command for the zp form is the same as the tf form:")
  260.       cmd = "rlocus(sys2);";
  261.       run_cmd;
  262.       disp("\nThe internal representation of the system is not important;");
  263.       disp("rlocus automatically sorts it out internally.  Also, it does not");
  264.       disp("matter if the system is continuous or discrete.  rlocus also sorts");
  265.       disp("this out automatically");
  266.       prompt;
  267.  
  268.       clc
  269.       disp("Example #3, Consider the following discrete state space system:\n");
  270.       cmd = "sys3=ss2sys([1 0.0952;0 0.905], [0.00484; 0.0952], [1 0], 0, Tsam);";
  271.       disp(cmd); 
  272.       eval(cmd);
  273.       cmd = "sysout(sys3);";
  274.       disp(cmd); 
  275.       eval(cmd);
  276.       disp("\nPole-zero form can be obtained as follows:");
  277.       cmd = "sysout(sys3,""zp"");";
  278.       disp(cmd);
  279.       eval(cmd); 
  280.       disp("\nOnce again, the rlocus command is the same:");
  281.       cmd = "rlocus(sys3);";
  282.       run_cmd;
  283.  
  284.       disp("\nNo matter what form the system is in, the rlocus command works the");
  285.       disp("the same.");
  286.  
  287.       prompt;
  288.  
  289.       closeplot
  290.       clc
  291.  
  292.     elseif (k == 5)
  293.       return
  294.     endif
  295.   endwhile  
  296. endfunction
  297.