1 ## Copyright (C) 2009 Lukas F. Reichlin
3 ## This file is part of LTI Syncope.
5 ## LTI Syncope is free software: you can redistribute it and/or modify
6 ## it under the terms of the GNU General Public License as published by
7 ## the Free Software Foundation, either version 3 of the License, or
8 ## (at your option) any later version.
10 ## LTI Syncope is distributed in the hope that it will be useful,
11 ## but WITHOUT ANY WARRANTY; without even the implied warranty of
12 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
13 ## GNU General Public License for more details.
15 ## You should have received a copy of the GNU General Public License
16 ## along with LTI Syncope. If not, see <http://www.gnu.org/licenses/>.
19 ## @deftypefn {Function File} pzmap (@var{sys})
20 ## @deftypefnx {Function File} {[@var{p}, @var{z}] =} pzmap (@var{sys})
21 ## Plot the poles and zeros of an LTI system in the complex plane.
22 ## If no output arguments are given, the result is plotted on the screen.
23 ## Otherwise, the poles and zeros are computed and returned.
34 ## Poles of @var{sys}.
36 ## Transmission zeros of @var{sys}.
40 ## Author: Lukas Reichlin <lukas.reichlin@gmail.com>
41 ## Created: November 2009
44 function [pol_r, zer_r] = pzmap (sys)
46 ## TODO: multiplot feature: pzmap (sys1, "b", sys2, "r", ...)
52 if (! isa (sys, "lti"))
53 error ("pzmap: argument must be an LTI system");
65 plot (pol_re, pol_im, "xb", zer_re, zer_im, "or")
67 title (["Pole-Zero Map of ", inputname(1)])
69 ylabel ("Imaginary Axis")