28.02.2014 Views

The Development of Neural Network Based System Identification ...

The Development of Neural Network Based System Identification ...

The Development of Neural Network Based System Identification ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

164 CHAPTER 6 NEURAL NETWORK BASED PREDICTIVE CONTROL SYSTEM<br />

6.4 NON-MINIMAL STATE SPACE MODEL REALISATION<br />

<strong>The</strong>re are generally three types <strong>of</strong> models used in MPC design such as finite impulse<br />

response (FIR)/step response models, transfer function models and state space models<br />

[Rossiter, 2003]. <strong>The</strong> state space model is preferred over other model types in MPC<br />

controller design due to its simplicity and effective handling <strong>of</strong> multi-variable problems.<br />

<strong>The</strong> transfer function models obtained from instantaneous linearisation <strong>of</strong> the NN model<br />

can also be represented in terms <strong>of</strong> state space model realisation. This can be done<br />

by reformulating the state variables <strong>of</strong> the state space model to be identical to the<br />

feedback variables that have been used in the ARX model [Wang, 2009b, Ordys and<br />

Clarke, 1993]. Consider the general discrete time ARX model that describes the input<br />

and output relationship as:<br />

A(q −1 )(1 − q −1 )y(k) = B(q −1 )∆u(k) + q −1 ɛ(t) (6.16)<br />

where ɛ(t) is the input disturbance which is assumed to be a sequence <strong>of</strong> integrated<br />

white noise, A(q −1 ) and B(q −1 ) are polynomials in the time shift operator given by the<br />

following forms:<br />

A(q −1 ) = 1 + a 1 q −1 + a 2 q −2 · · · + a n q −n<br />

B(q −1 ) = b 1 q −1 + b 2 q −2 · · · + b n q −n (6.17)<br />

Let the polynomial A(q −1 )(1 − q −1 ) be referred as:<br />

A(q −1 )(1 − q −1 ) = 1 + ā 1 q −1 + ā 2 q −2 + · · · + ā n+1 q −(n+1) (6.18)<br />

Next, by choosing the state variable vector as:<br />

x(k) = [y(k) y(k − 1) · · · y(k − n − 1) ∆u(k − 1) ∆u(k − 2) · · · ∆u(k − n)] T (6.19)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!