[Octave-cvsupdate] SF.net SVN: octave:[8862] trunk/octave-forge/extra/control-devel

[<par...@us...>]

Revision: 8862
          http://octave.svn.sourceforge.net/octave/?rev=8862&view=rev
Author:   paramaniac
Date:     2011-10-25 21:30:19 +0000 (Tue, 25 Oct 2011)
Log Message:
-----------
control-devel: compare different results

Modified Paths:
--------------
    trunk/octave-forge/extra/control-devel/inst/hnamodred.m

Added Paths:
-----------
    trunk/octave-forge/extra/control-devel/devel/compare_results_hnamodred.m

Added: trunk/octave-forge/extra/control-devel/devel/compare_results_hnamodred.m
===================================================================
--- trunk/octave-forge/extra/control-devel/devel/compare_results_hnamodred.m	                        (rev 0)
+++ trunk/octave-forge/extra/control-devel/devel/compare_results_hnamodred.m	2011-10-25 21:30:19 UTC (rev 8862)
@@ -0,0 +1,33 @@
+    Mo = [
+
+      -0.23915   0.30723   1.16297  -1.19671   1.04965
+      -2.97091  -0.23915   2.62702  -3.10273   3.70515
+       0.00000   0.00000  -0.51368   1.28421  -0.82227
+       0.00000   0.00000  -0.15189  -0.51368   0.74348
+       0.44660  -0.01427   0.47803  -0.20129   0.02190
+];
+
+    Me = [
+
+      -0.23910   0.30720   1.16300   1.19670  -1.04970
+      -2.97090  -0.23910   2.62700   3.10270  -3.70520
+       0.00000   0.00000  -0.51370  -1.28420   0.82230
+       0.00000   0.00000   0.15190  -0.51370   0.74350
+      -0.44660   0.01430  -0.47800  -0.20130   0.02190
+];
+
+syso = ss (Mo(1:4, 1:4), Mo(1:4, 5), Mo(5, 1:4), Mo(5, 5))
+
+syse = ss (Me(1:4, 1:4), Me(1:4, 5), Me(5, 1:4), Me(5, 5))
+
+figure (1)
+bode (syso)
+
+figure (2)
+bode (syse)
+
+figure (3)
+step (syso)
+
+figure (4)
+step (syse)
\ No newline at end of file

Modified: trunk/octave-forge/extra/control-devel/inst/hnamodred.m
===================================================================
--- trunk/octave-forge/extra/control-devel/inst/hnamodred.m	2011-10-25 15:20:09 UTC (rev 8861)
+++ trunk/octave-forge/extra/control-devel/inst/hnamodred.m	2011-10-25 21:30:19 UTC (rev 8862)
@@ -133,6 +133,8 @@
     endswitch
   endfor
   
+  ## TODO: handle jobv, jobw, (jobinv)
+  
   ## perform model order reduction
   [ar, br, cr, dr, nr] = slab09jd (a, b, c, d, dt, scaled, nr, ordsel, alpha, \
                                    jobv, av, bv, cv, dv, \

This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site.