Revision: 11081
http://octave.svn.sourceforge.net/octave/?rev=11081&view=rev
Author: paramaniac
Date: 2012-09-22 10:58:32 +0000 (Sat, 22 Sep 2012)
Log Message:
-----------
control: use multiplot step command (4)
Modified Paths:
--------------
trunk/octave-forge/main/control/inst/Anderson.m
Modified: trunk/octave-forge/main/control/inst/Anderson.m
===================================================================
--- trunk/octave-forge/main/control/inst/Anderson.m 2012-09-22 10:45:18 UTC (rev 11080)
+++ trunk/octave-forge/main/control/inst/Anderson.m 2012-09-22 10:58:32 UTC (rev 11081)
@@ -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)
This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site.
|
