|
3432
|
1 ## Copyright (C) 1996, 1998 Auburn University. All rights reserved. |
|
|
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, 59 Temple Place, Suite 330, Boston, MA 02111 USA. |
|
|
18 |
|
|
19 ## -*- texinfo -*- |
|
|
20 ## @deftypefn {Function File} {[@var{zer}, @var{pol}]=} pzmap (@var{sys}) |
|
|
21 ## Plots the zeros and poles of a system in the complex plane. |
|
|
22 ## @strong{Inputs} |
|
|
23 ## @var{sys} system data structure |
|
|
24 ## |
|
|
25 ## @strong{Outputs} |
|
|
26 ## if omitted, the poles and zeros are plotted on the screen. |
|
|
27 ## otherwise, pol, zer are returned as the system poles and zeros. |
|
|
28 ## (see sys2zp for a preferable function call) |
|
|
29 ## @end deftypefn |
|
|
30 |
|
|
31 function [zer, pol]=pzmap (sys) |
|
|
32 |
|
4460
|
33 save_warn_empty_list_elements = warn_empty_list_elements; |
|
|
34 unwind_protect |
|
|
35 warn_empty_list_elements = 0; |
|
3432
|
36 |
|
4460
|
37 if(nargin != 1) |
|
|
38 usage("pzmap(sys) or [zer,pol] = pzmap(sys)"); |
|
|
39 elseif (!isstruct(sys)); |
|
|
40 error("sys must be in system format"); |
|
|
41 endif |
|
|
42 |
|
|
43 [zer,pol] = sys2zp(sys); |
|
3432
|
44 |
|
4460
|
45 ## force to column vectors, split into real, imaginary parts |
|
|
46 zerdata = poldata = []; |
|
|
47 if(length(zer)) |
|
|
48 zer = reshape(zer,length(zer),1); |
|
|
49 zerdata = [real(zer(:,1)), imag(zer(:,1))]; |
|
|
50 endif |
|
|
51 if(length(pol)) |
|
|
52 pol = reshape(pol,length(pol),1); |
|
|
53 poldata = [real(pol(:,1)), imag(pol(:,1))]; |
|
|
54 endif |
|
3432
|
55 |
|
4460
|
56 ## determine continuous or discrete plane |
|
|
57 vars = "sz"; |
|
|
58 varstr = vars(is_digital(sys) + 1); |
|
3432
|
59 |
|
4460
|
60 ## Plot the data |
|
|
61 gset nologscale xy; |
|
|
62 if(is_siso(sys)) |
|
|
63 title(sprintf("Pole-zero map from %s to %s", ... |
|
|
64 sysgetsignals(sys,"in",1,1), sysgetsignals(sys,"out",1,1) )); |
|
|
65 endif |
|
|
66 xlabel(["Re(",varstr,")"]); |
|
|
67 ylabel(["Im(",varstr,")"]); |
|
|
68 grid; |
|
3432
|
69 |
|
4460
|
70 ## compute axis limits |
|
|
71 axis(axis2dlim([zerdata;poldata])); |
|
|
72 grid |
|
|
73 ## finally, plot the data |
|
|
74 if(length(zer) == 0) |
|
|
75 plot(poldata(:,1), poldata(:,2),"@12 ;poles (no zeros);"); |
|
|
76 elseif(length(pol) == 0) |
|
|
77 plot(zerdata(:,1), zerdata(:,2),"@31 ;zeros (no poles);"); |
|
|
78 else |
|
|
79 plot(zerdata(:,1), zerdata(:,2),"@31 ;zeros;", ... |
|
|
80 poldata(:,1), poldata(:,2),"@12 ;poles;"); |
|
|
81 endif |
|
|
82 replot |
|
3432
|
83 |
|
4460
|
84 unwind_protect_cleanup |
|
|
85 warn_empty_list_elements = save_warn_empty_list_elements; |
|
|
86 end_unwind_protect |
|
3432
|
87 |
|
|
88 endfunction |
