Mercurial > forge

changeset 8676:6277883c4428 octave-forge

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
control-devel: compare different results
author paramaniac
date Tue, 25 Oct 2011 21:30:19 +0000
parents 8be778dfb306
children f70999737e92
files extra/control-devel/devel/compare_results_hnamodred.m extra/control-devel/inst/hnamodred.m
diffstat 2 files changed, 35 insertions(+), 0 deletions(-) [+]
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/extra/control-devel/devel/compare_results_hnamodred.m	Tue Oct 25 21:30:19 2011 +0000
@@ -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
--- a/extra/control-devel/inst/hnamodred.m	Tue Oct 25 15:20:09 2011 +0000
+++ b/extra/control-devel/inst/hnamodred.m	Tue Oct 25 21:30:19 2011 +0000
@@ -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, \