Mercurial > forge

changeset 9390:e090b67004a1 octave-forge

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
control: work on draft code
author paramaniac
date Tue, 31 Jan 2012 18:12:58 +0000
parents 97f77d4bd426
children ea88fce5f7ff
files main/control/devel/__frequency_vector_2__.m
diffstat 1 files changed, 4 insertions(+), 2 deletions(-) [+]
line wrap: on
line diff
--- a/main/control/devel/__frequency_vector_2__.m	Mon Jan 30 22:19:16 2012 +0000
+++ b/main/control/devel/__frequency_vector_2__.m	Tue Jan 31 18:12:58 2012 +0000
@@ -51,6 +51,10 @@
   dec_min = min (cell2mat (dec_min));
   dec_max = max (cell2mat (dec_max));
   zp = horzcat (zp{:});
+  
+  ## include zeros and poles for nice peaks in plots
+  idx = find (zp > 10^dec_min & zp < 10^dec_max);
+  zp = zp(idx);
 
   w = logspace (dec_min, dec_max, 500);
   w = unique ([w, zp]);                  # unique also sorts frequency vector
@@ -148,7 +152,5 @@
 
   ## include zeros and poles for nice peaks in plots
   zp = [abs(zer), abs(pol)];
-  idx = find (zp > 10^dec_min & zp < 10^dec_max);
-  zp = zp(idx);
 
 endfunction