Menu
▾
▴
[Octave-cvsupdate] SF.net SVN: octave:[10640] trunk/octave-forge/extra/control-devel/devel
|
From: <par...@us...> - 2012-06-19 11:41:11
|
Revision: 10640
http://octave.svn.sourceforge.net/octave/?rev=10640&view=rev
Author: paramaniac
Date: 2012-06-19 11:41:05 +0000 (Tue, 19 Jun 2012)
Log Message:
-----------
control-devel: compute initial states for examples
Modified Paths:
--------------
trunk/octave-forge/extra/control-devel/devel/GlassFurnaceARX.m
trunk/octave-forge/extra/control-devel/devel/pHarx.m
Modified: trunk/octave-forge/extra/control-devel/devel/GlassFurnaceARX.m
===================================================================
--- trunk/octave-forge/extra/control-devel/devel/GlassFurnaceARX.m 2012-06-19 10:43:36 UTC (rev 10639)
+++ trunk/octave-forge/extra/control-devel/devel/GlassFurnaceARX.m 2012-06-19 11:41:05 UTC (rev 10640)
@@ -49,12 +49,13 @@
dat = iddata (Y, U)
%[sys, x0] = ident (dat, 10, 5) % s=10, n=5
-sys = arx (dat, 5)
+%sys = arx (dat, 5)
+[sys, x0] = arx (dat, 5)
%[y, t] = lsim (sys, U, [], x0);
-[y, t] = lsim (sys(:, 1:3), U);
+%[y, t] = lsim (sys(:, 1:3), U);
+[y, t] = lsim (sys, U, [], x0);
-
err = norm (Y - y, 1) / norm (Y, 1)
figure (1)
Modified: trunk/octave-forge/extra/control-devel/devel/pHarx.m
===================================================================
--- trunk/octave-forge/extra/control-devel/devel/pHarx.m 2012-06-19 10:43:36 UTC (rev 10639)
+++ trunk/octave-forge/extra/control-devel/devel/pHarx.m 2012-06-19 11:41:05 UTC (rev 10640)
@@ -52,12 +52,14 @@
dat = iddata (Y, U)
% [sys, x0] = ident (dat, 15, 6) % s=15, n=6
-sys = arx (dat, 6) % normally na = nb
+% sys = arx (dat, 6) % normally na = nb
+[sys, x0] = arx (dat, 6) % normally na = nb
+
% [y, t] = lsim (sys, U, [], x0);
-[y, t] = lsim (sys(:, 1:2), U);
+% [y, t] = lsim (sys(:, 1:2), U);
+[y, t] = lsim (sys, U, [], x0);
-
err = norm (Y - y, 1) / norm (Y, 1)
figure (1)
This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site.
|
