[Octave-cvsupdate] SF.net SVN: octave:[10352] trunk/octave-forge/main/control/inst/@ss/ __sys_inverse__.m

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

Revision: 10352
          http://octave.svn.sourceforge.net/octave/?rev=10352&view=rev
Author:   paramaniac
Date:     2012-05-01 20:43:31 +0000 (Tue, 01 May 2012)
Log Message:
-----------
control: don't compute inverse explicitly (where possible) for state-space inversion

Modified Paths:
--------------
    trunk/octave-forge/main/control/inst/@ss/__sys_inverse__.m

Modified: trunk/octave-forge/main/control/inst/@ss/__sys_inverse__.m
===================================================================
--- trunk/octave-forge/main/control/inst/@ss/__sys_inverse__.m	2012-05-01 20:31:07 UTC (rev 10351)
+++ trunk/octave-forge/main/control/inst/@ss/__sys_inverse__.m	2012-05-01 20:43:31 UTC (rev 10352)
@@ -1,4 +1,4 @@
-## Copyright (C) 2009, 2010   Lukas F. Reichlin
+## Copyright (C) 2009, 2010, 2012   Lukas F. Reichlin
 ##
 ## This file is part of LTI Syncope.
 ##
@@ -20,7 +20,7 @@
 
 ## Author: Lukas Reichlin <luk...@gm...>
 ## Created: October 2009
-## Version: 0.2
+## Version: 0.3
 
 function sys = __sys_inverse__ (sys)
 
@@ -49,12 +49,12 @@
 
   else                                   # proper ss
 
-    di = inv (d);
+    bid = b / d;
 
-    sys.a = a - b * di * c;
-    sys.b = -b * di;
-    sys.c = di * c;
-    sys.d = di;
+    sys.a = a - bid * c;
+    sys.b = -bid;
+    sys.c = d \ c;
+    sys.d = inv (d);
 
   endif
 

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