25.01.2015 Views

Download Full Issue in PDF - Academy Publisher

Download Full Issue in PDF - Academy Publisher

Download Full Issue in PDF - Academy Publisher

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.

JOURNAL OF COMPUTERS, VOL. 8, NO. 6, JUNE 2013 1453<br />

z(t)<br />

30<br />

20<br />

10<br />

0<br />

-10<br />

-20<br />

-30<br />

20<br />

10<br />

0<br />

y(t)<br />

-10<br />

attractor of Lorenz<br />

-20 -20<br />

(a) three-dimensional map of Lorenz attractor<br />

y(t)<br />

20<br />

15<br />

10<br />

5<br />

0<br />

-5<br />

-10<br />

-15<br />

attractor of Lorenz<br />

-20<br />

-10 0 10 20 30 40 50<br />

x(t)<br />

(b) two-dimensional map of Lorenz attractor <strong>in</strong> the x-y-plan<br />

z(t)<br />

30<br />

20<br />

10<br />

0<br />

-10<br />

-20<br />

attractor of Lorenz<br />

-30<br />

-10 0 10 20 30 40 50<br />

x(t)<br />

(c) two-dimensional map of Lorenz attractor <strong>in</strong> the x-z-plan<br />

y(t)<br />

30<br />

20<br />

10<br />

0<br />

-10<br />

-20<br />

attractor of Lorenz<br />

-30<br />

-20 -15 -10 -5 0 5 10 15 20<br />

z(t)<br />

(d) two-dimensional map of Lorenz attractor <strong>in</strong> the z-y-plan<br />

0<br />

x(t)<br />

20<br />

40<br />

60<br />

Lorenz map:<br />

•<br />

⎧<br />

⎪<br />

x = σ ( y − x)<br />

⎪ •<br />

⎨ y = rx − y − xz<br />

(13)<br />

⎪ •<br />

⎪ z =− bz+<br />

xy<br />

⎩<br />

8<br />

Where σ = 10 , r = 28 , b = .The <strong>in</strong>itial value is<br />

3<br />

x (0) = 0 , y (0) = 5 , z (0) = − 5 ; and the fix<strong>in</strong>g step length<br />

of <strong>in</strong>itial value is 0.05s . Time series to the branch x with<br />

70s is produced by the Runge-Kutta algorithms and the<br />

total data is 1200. The embedded dimension of the<br />

sampl<strong>in</strong>g chaotic time series m is 8 by the G- P<br />

algorithms. The delay time is τ= 1 by the self-correlation<br />

function algorithms and the <strong>in</strong>put dimension of the<br />

adaptive RBF neural network filter<strong>in</strong>g is 8.The former<br />

1200 data is tra<strong>in</strong>ed and other 200 data is predicted by the<br />

adaptive RBF neural network filter<strong>in</strong>g predictive model.<br />

B. Evaluation of the Predictive Ability<br />

The model's predictive ability is generally measure the<br />

follow<strong>in</strong>g three <strong>in</strong>dicators: of MAPE (mean absolute<br />

percentage error), RMSE (root mean square error) and<br />

RMSPE (root mean square percentage error), they are<br />

calculated as follows:<br />

n<br />

1 yi<br />

− yi<br />

MAPE = 100<br />

n<br />

∑ × , (14)<br />

y<br />

i=<br />

1<br />

i<br />

n<br />

1 ⎛ yi<br />

− y ⎞<br />

i<br />

RMSPE = 100× ∑ ⎜ ⎟<br />

n I = 1 ⎝ y<br />

, (15)<br />

i ⎠<br />

n I = 1<br />

i i<br />

1 n<br />

RMSE = y − y<br />

∑( ) 2<br />

(16)<br />

where, y i<br />

is predictive value of the model; y i<br />

is the real<br />

value; n is prediction phases, and MAPE assess the<br />

predictive capability are as follows: less than or equal to<br />

10%, then predictive ability is excellent; 10% -20%, then<br />

the predictive ability is excellent; 20% -50%, more than<br />

50%, then the prediction is <strong>in</strong>accurate. For RMSPE, the<br />

prediction square vulnerable to the impact of outliers, for<br />

the larger error given greater weight, but still can be<br />

modeled on the MAPE to determ<strong>in</strong>e the model of the pros<br />

and cons. RMSPE values range from zero to <strong>in</strong>f<strong>in</strong>ity.<br />

MAPE and RMSPE are the relative <strong>in</strong>dicator, RMSE is<br />

the absolute <strong>in</strong>dicator. The RMSE is the smaller, the<br />

model predictive ability is the stronger.<br />

C. The Simulation Results<br />

That the experimental outcome of Lorenz chaotic<br />

sampl<strong>in</strong>g time series, the true value (real l<strong>in</strong>e) and the<br />

predictive value (star l<strong>in</strong>e) and the predictive error curve<br />

are showed <strong>in</strong> Figure 3., Figure 4. and Figure 5.<br />

2<br />

Figure 2. Lorenz attractor <strong>in</strong> the phase space reconstruction<br />

Consider<strong>in</strong>g Lorenz chaotic system<br />

© 2013 ACADEMY PUBLISHER

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

Saved successfully!

Ooh no, something went wrong!