home *** CD-ROM | disk | FTP | other *** search
/ OS/2 Shareware BBS: 10 Tools / 10-Tools.zip / octa21eb.zip / octave / SCRIPTS.ZIP / scripts.fat / control / zpout.m < prev   
Text File  |  1999-04-29  |  3KB  |  112 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 zpout(zer,pol,k,x)
  20. #
  21. # usage: zpout(zer,pol,k,[,x])
  22. #  
  23. # print formatted zero-pole form
  24. # to the screen
  25. # x defaults to the string "s"
  26. #
  27. #  SEE ALSO: polyval, polyvalm, poly, roots, conv, deconv, residue, 
  28. #    filter, polyderv, polyintg, polyout
  29.  
  30. # Written by A. Scottedward Hodel (scotte@eng.auburn.edu) June 1995)
  31.  
  32.   save_val = implicit_str_to_num_ok;
  33.   save_empty = empty_list_elements_ok;
  34.   empty_list_elements_ok = 1;
  35.  
  36.   implicit_str_to_num_ok = 1;
  37.  
  38.   if (nargin < 3 ) | (nargin > 4) | (nargout != 0 )
  39.     usage("zpout(zer,pol,k[,x])");
  40.   endif
  41.  
  42.   if( !(is_vec(zer) | isempty(zer)) | !(is_vec(pol) | isempty(pol)) )
  43.     error("zer, pol must be vectors or empty");
  44.   endif
  45.  
  46.   if(!is_scal(k))
  47.     error("zpout: argument k must be a scalar.")
  48.   endif
  49.  
  50.   if (nargin == 3)
  51.     x = 's';
  52.   elseif( ! isstr(x) )
  53.     error("zpout: third argument must be a string");
  54.   endif
  55.  
  56.   numstring = num2str(k);
  57.  
  58.   if(length(zer))
  59.     # find roots at z,s = 0
  60.     nzr = sum(zer == 0);
  61.     if(nzr)
  62.       if(nzr > 1)
  63.         numstring = [numstring,sprintf(" %s^%d",x,nzr)];
  64.       else
  65.         numstring = [numstring,sprintf(" %s",x)];
  66.       endif
  67.     endif
  68.     zer = sortcom(-zer);
  69.     for ii=1:length(zer)
  70.       if(zer(ii) != 0)
  71.         numstring = [numstring,sprintf(" (%s %s)",x,com2str(zer(ii),1) ) ];
  72.       endif
  73.     endfor
  74.   endif
  75.  
  76.   if(length(pol))
  77.     # find roots at z,s = 0
  78.     nzr = sum(pol == 0);
  79.     if(nzr)
  80.       if(nzr > 1)
  81.         denomstring = [sprintf("%s^%d",x,nzr)];
  82.       else
  83.         denomstring = [sprintf("%s",x)];
  84.       endif
  85.     else
  86.       denomstring = " ";
  87.     endif
  88.     pol = sortcom(-pol);
  89.     for ii=1:length(pol)
  90.       if(pol(ii) != 0)
  91.         denomstring = [denomstring,sprintf(" (%s %s)",x,com2str(pol(ii),1))];
  92.       endif
  93.     endfor
  94.   endif
  95.  
  96.   len = max(length(numstring),length(denomstring));
  97.   if(len > 0)
  98.     y = strrep(blanks(len)," ","-");
  99.     disp(numstring)
  100.     if(length(denomstring))
  101.       disp(y)
  102.       disp(denomstring)
  103.     endif
  104.   else
  105.     error('zpout: empty transfer function')
  106.   end
  107.  
  108.   implicit_str_to_num_ok = save_val;
  109.   empty_list_elements_ok = save_empty;
  110.  
  111. endfunction
  112.