Mercurial > forge
changeset 10454:52962c2ee902 octave-forge
control-devel: compute initial states for examples
author | paramaniac |
---|---|
date | Tue, 19 Jun 2012 11:41:05 +0000 |
parents | 05acd81505fd |
children | 34b4175ba644 |
files | extra/control-devel/devel/GlassFurnaceARX.m extra/control-devel/devel/pHarx.m |
diffstat | 2 files changed, 9 insertions(+), 6 deletions(-) [+] |
line wrap: on
line diff
--- a/extra/control-devel/devel/GlassFurnaceARX.m Tue Jun 19 10:43:36 2012 +0000 +++ b/extra/control-devel/devel/GlassFurnaceARX.m Tue Jun 19 11:41:05 2012 +0000 @@ -49,11 +49,12 @@ 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)
--- a/extra/control-devel/devel/pHarx.m Tue Jun 19 10:43:36 2012 +0000 +++ b/extra/control-devel/devel/pHarx.m Tue Jun 19 11:41:05 2012 +0000 @@ -52,11 +52,13 @@ 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)