From: <par...@us...> - 2012-09-22 10:58:38
|
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. |