Mercurial > forge

changeset 10895:75b48e65ae10 octave-forge

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
control: use multiplot step command (4)
author paramaniac
date Sat, 22 Sep 2012 10:58:32 +0000
parents 02d87d375f99
children 1eb1dfc968ac
files main/control/inst/Anderson.m
diffstat 1 files changed, 12 insertions(+), 21 deletions(-) [+]
line wrap: on
line diff
--- a/main/control/inst/Anderson.m	Sat Sep 22 10:45:18 2012 +0000
+++ b/main/control/inst/Anderson.m	Sat Sep 22 10:58:32 2012 +0000
@@ -50,31 +50,22 @@
 L = lqe (G, W, V)
 
 % Coprime Factorization using Balanced Truncation Approximation
+Kr = arrayfun (@(k) cfconred (G, F, L, k), 8:-1:2, 'uniformoutput', false);  % 'method', 'bfsr-bta'
+T = cellfun (@(Kr) feedback (G*Kr), Kr, 'uniformoutput', false);
+
 figure (1)
-for k = 8:-1:2
-  Kr = cfconred (G, F, L, k);   % 'method', 'bfsr-bta'
-  T = feedback (G*Kr);
-  step (T, 200)
-  hold on
-endfor
-hold off
+step (T{:}, 200)
 
 % Coprime Factorization using Singular Perturbation Approximation
+Kr = arrayfun (@(k) cfconred (G, F, L, k, 'method', 'bfsr-spa'), 8:-1:2, 'uniformoutput', false);
+T = cellfun (@(Kr) feedback (G*Kr), Kr, 'uniformoutput', false);
+
 figure (2)
-for k = 8:-1:2
-  Kr = cfconred (G, F, L, k, 'method', 'bfsr-spa');
-  T = feedback (G*Kr);
-  step (T, 200)
-  hold on
-endfor
-hold off
+step (T{:}, 200)
 
 % Frequency-Weighted Coprime Factorization using BTA
+Kr = arrayfun (@(k) fwcfconred (G, F, L, k), 8:-1:2, 'uniformoutput', false);
+T = cellfun (@(Kr) feedback (G*Kr), Kr, 'uniformoutput', false);
+
 figure (3)
-for k = 8:-1:2
-  Kr = fwcfconred (G, F, L, k);
-  T = feedback (G*Kr);
-  step (T, 300)
-  hold on
-endfor
-hold off
+step (T{:}, 200)