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

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

98 CHAPTER 4 NEURAL NETWORK BASED SYSTEM IDENTIFICATION<br />

and outputs <strong>of</strong> the neural network are presented by m and n respectively. Similarly,<br />

the NNARX model can also be represented using HMLP network by introducing the<br />

lagged time input variables into the neural network model.<br />

<strong>The</strong> Elman network with modification in the network’s internal dynamics is found<br />

capable <strong>of</strong> identifying an n th order discrete dynamic system in Equation (4.23), through<br />

empirical findings in Pham and Liu [1993]. Furthermore, Pham and Liu [1996] suggested<br />

that in order to model the dynamic system in Equation (4.23) from experimental data<br />

using MLP network, a n y + n u input (regressor) nodes are needed to be included in<br />

the network. However, if an Elman network is used to model a NNARX model, only<br />

current measurement data are needed to be fed into the input nodes as shown in Figure<br />

4.9(c). <strong>The</strong>refore, the modified Elman network is significantly smaller in size compared<br />

with MLP or HMLP network when large time lags (model order) are used.<br />

4.3.2.1 Lag Space Selection for Feed-Forward MLP or HMLP <strong>Network</strong><br />

After deciding the input parameters to be used in the model, the number <strong>of</strong> past inputs<br />

and outputs fed into the MLP or HMLP neural networks were decided based on the<br />

calculation <strong>of</strong> the Lipschitz coefficient given in Norgaard [2000]. Using this coefficient,<br />

it is possible to determine the proper lag space via experimental data. <strong>The</strong> sizes <strong>of</strong><br />

output and input time regression vectors depend on the degree <strong>of</strong> non-linearity <strong>of</strong><br />

the Lipschitz coefficients where an insufficient number <strong>of</strong> regressors will result in high<br />

Lipschitz coefficients and small numerical values for extra regressors.<br />

<strong>The</strong> Lipschitz coefficient is calculated using the following formula for each input<br />

u i (t) and output y j (t) pairs:<br />

q ij =<br />

y (t i ) − y (t j )<br />

∣ϕ (t i ) − ϕ (t j ) ∣ i ≠ j (4.27)<br />

<strong>The</strong> outline <strong>of</strong> Lipschitz coefficient approach is given as follows:<br />

1. Determine the Lipschitz quotients for all combination <strong>of</strong> input-output pairs using<br />

(4.27) for a given choice <strong>of</strong> number <strong>of</strong> past outputs and inputs.

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

Saved successfully!

Ooh no, something went wrong!