A Hybrid Forecast of Exchange Rate based on Discrete Grey-Markov
and Grey Neural Network Model
Gol Kim a, Ri Suk Yun b
(a Center of Natural Science, University of Sciences, Pyongyang, DPR Korea,
E-mail: golkim124@yahoo.com
b

Foreign Economic General Bureau, Pyongyang, DPR Korea)

Abstract: We propose a hybrid forecast model based on discrete grey-fuzzy Markov and grey –neural network
model and show that our hybrid model can improve much more the performance of forecast than
traditional grey-Markov model and neural network models. Our simulation results are shown that our
hybrid forecast method with the combinational weight based on optimal grey relation degree method is
better than the hybrid model with combinational weight based minimization of error-squared criterion.
Keywords: Financial forecasting, Hybrid forecast, Discrete Grey-Fuzzy Markov，Grey –Neural Network.

1. Introduction
Investors have been trying to find a way to predict exchange rates and stock price accurately, but
haven’t been obtained the good results.Recently, artificial intelligence techniques like artificial
neural networks (ANN) and genetic algorithms (GA), and wavelet transform have been applied to
this area. However, the above-mentioned concern method showed that ANN had some limitations
in learning the patterns because foreign exchange rates data has tremendous noise and complex
dimensionality. Moreover, the sheer quantity of exchange rates data sometimes interferes with the
learning of patterns.
Hybrid forecast is a well-established and well-tested approach for improving the forecasting
accuracy. Therefore, the importance of hybrid forecast methods has steadily increased and it acts
still on time series forecasting. Hybrid forecast system using rough sets and artificial neural
networks has been proposed by M.Duraira, K. Meena (2011) and hybrid forecast method by
ARIMA and neural network model Has been researched by Reza Askari Moghadam, Mehrnaz
Keshmirpour (2011). Aobing Sun, Yubo Tan, Dexian Zhang (2008) have researched the hybrid
forecast model based on BP neural network.
Hybrid Forecast Model of ARIMA and Neural Network has been given by Šterba Ján, Hil’ovsk
Katarína (2010) and stacked heterogeneous neural networks for time series forecasting has been
suggested by Florin Leon, Mihai Horia Zaharia(2010).
Tarek Aboueldahab, Mahumod Fakhreldin (2011) have studied the forecast of stock market indices
using hybrid genetic algorithm/ particle swarm optimization with perturbation term.
Improvement method in forecasting accuracy using the hybrid model of ARFIMA and feed forward
neural network has been proposed by Cagdas Hakan Aladag, Erol Egrioglu, Cem Kadilar (2012).
Forecast of future stock close price using based on hybrid ANN Model of functional link fuzzy
logic neural model has been researched by Kumaran Kumar. J, Kailas (2012).
Forecasting of foreign exchange rate by grey- chaotic forecast hybrid has been researched by Kim Gol
(2008) and forecasting of stock price and index of oil by grey- fractal forecast hybrid model has
been studied by Kim Gol (2009).
Zhu Jiang-liang (2006) has studied the forecasting of the electric power load based on the grey
theory and BP neural network and Su Bo (2006) has considered the comparison and research of
gain production forecasting with methods of GM (1, N) grey system and BPNN. Method of the
wavelet neural network in the forecast of stock market has been suggested by,Yao Hong-Xing,
Shong Zhao-han, Chen Hong-xiang (2002) and Genetic algorithms approach to feature

1

discretization in artificial neural networks for the forecast of stock price has been researched by
Kyong-jac Kim, Ingoo Han (2000).
This paper deals with the development of an improved forecast model using grey-Markov chain
models and the grey-neural network and its application to financial time series forecasting. Our
method includes discrete grey- fuzzy weight Markov chain model, grey-neural network model, The
hybrid model with combinational weight based minimization of error-squared criterion, hybrid
model with combinational weight based on optimal grey relation degree method.
We demonstrate the advantage of our method compared with the traditional grey-Markov or neural
network methods through forecasting simulations.
2. Discrete grey- Fuzzy weight Markov model
2.1. Grey forecasting model
Grey forecasting model (GM) has three basic operation such as ;
(і) accumulated generation, (ⅱ) inverse accumulated generation, and (ⅲ) grey modeling.
The grey forecasting model uses the operations of accumulated generation to build differential
equations. Intrinsically speaking, it has the characteristics of requiring less data.
The GM (1, 1) grey model, i.e., a 1- variable first-order grey model, is summarized as follows
Step 1: Selection of the original sequence;
x (0) = ( x (0) (1), x (0) (2),⋯ x (0) (i ),⋯ , x (0) ( N ))
(5)
where x (0) (i ) is the time series data at time i .
Step 2: Accumulated generating;
Based on the time sequence x (0) , a new sequence x (1) is given by the accumulated generating
operation (AGO), where
(1)
x = ( x (1) (1), x (1) (2),⋯ x (1) (i ),⋯ , x (1) ( N )) and is x (1) (k ) derived as follows:
k

x (1) (k ) = ∑ x (0) (i )

(6)

i =1

Step 3: Making the differential equation;
The first-order differential equation holds true:
dx (1)
+ ax (1) = u
dt
Step 4: Obtaining the forecast model; We have by (7)
u
u
xˆ (1) (k + 1) = ( x (0) (1) − )e − ak +
a
a
(0)
(1)
(1)
xɵ (k + 1) = xɵ (k + 1) − xɵ (k )

(7)

(8)
(9)

where

[ a, u ]

T

= ( B T B ) −1 B T y N

 −0.5( x (1) (1) + x (1) (2))
1


−0.5( x (1) (2) + x (1) (3))
1

B=

⋮
⋮


(1)
(1)
x
n
x
N
−
0.5(
(
−
1)
+
(
))
1


(0)
(0)
(0)
yN = ( x (2), x (3),⋯ x (i ),⋯ , x (0) ( N ))T

(10)

(11)

(12)

xˆ (k + 1) is the predicted value of x (k + 1) at time k + 1 ,
The authors used the posterior test to evaluate the accuracy of the grey forecasting.
The forecasting error wear defined as q (0) ( k ) = x (0) ( k ) − xˆ (0) ( k ) , k = 1, 2,⋯ , N , the mean and the
(1)

(1)

2

standard deviation of the forecasting error is q and S2 .The mean and the standard deviation of
initial time series are
N

m=

∑x

(0)

(k )
, S1 =

k =1

N

∑ (x

(0)

( k ) − m) 2 / N − 1

N
k =1
The posterior ratio C is derived by dividing S2 by S1 .The lower C is, the better the model is. The
posterior ratio can indicate the change rate the forecasting error. Probability of small is defined as
P = prob.{ q (0) ( k ) − q ≤ 0.6745S1} , where k = 2, 3,⋯ , N . P This shows the probability that the
relative bias of the forecasting error is lower than 0.6745. P is commonly required to be larger than
0.95. The pairs of the forecasting indicators P and C can characterize four grades of forecasting
accuracy, as shown Table 1.
It is called that above-mentioned grey forecast model is traditional grey model.
Table 1. The Grades of Forecasting Accuracy
Forecasting index
Grade
Good
Qualified
Just
Unqualified

P

C

>0.95
>0.8
>0.7
≤ 0.7

<0.35
<0.5
<0.65

≥ 0.65

2.2. Non homogenous discrete grey model
Let x (0) = ( x (0) (1), x (0) (2),⋯ x (0) (i ),⋯ , x (0) ( N )) be an original sequence and
x (1) = ( x (1) (1), x (1) (2),⋯ x (1) (i ),⋯ , x (1) ( N )) be accumulated generating sequence by AGO.

Where x (1) (k ) is given such as:
k

x (1) (k ) = ∑ x (0) (i )
i =1

Then, we define non-homogenous discrete grey model (DGM) such as ([ 2]) ;
 xˆ (1) (k + 1) = β1 xˆ (1) (k ) + β 2 xˆ (1) (0) + β 3 k + β 4
 (1)
 xˆ (1) = ξ
where x̂ (1) (k) fitting value of original is sequence and β 1 , β 2 , β 3 , β 4 are parameters of system.
By least square method, the parameters of system β = [β1 , β2 , β3 , β4 ] are given by

β = [ β 1 , β 2 , β 3 , β 4 ] = ( B T B ) −1 B T Y
 x (1) (1)
 x (1) (2)
1 1
x ( 0) (1)
 (1)
 (1) 

2 1
x ( 2)
x ( 0 ) ( 2)
x (3) 

Here, B =
Y =

 ⋮ 
⋮
⋮
⋮
⋮
 (1)
 (1) 

(0)
 x (n − 1) x (n − 1) n − 1 1
 x (n)
The algorithm of non-homogenous discrete grey model (DGM) is given such as;
Step 1; Find the parameters of system β = [β1 , β2 , β3 , β4 ] .
Step 2; Put xˆ (1) (1) = ξ and find the simulation values x̂ ( 0) (k ), k = 2,3,⋯ , n .
xˆ ( 0 ) (k ) is the function of ξ .
n

Step 3; Calculate Q = ∑ [ xˆ ( 0) (k ) − x ( 0 ) (k )]2 . Q is the function of ξ .
k =1

3

dQ
= 0 and calculate the value ξ which Q have to take minimum.
dξ
Step 5; Calculate the value xˆ (1) (k ) corresponding to ξ obtained from step 4.
Step 6; calculate the forecast value xˆ ( 0) (k ), k = 2,3,⋯ , n by inverse accumulated generating
operator (IAGO). That is, xˆ ( 0) (k ) = xˆ (1) (k + 1) − xˆ (1) (k ) .
Step 4; Put

2.3. Fuzzy weight Markov model
We assume that X̂ t is the fitting curve which is obtained through forecasting for the time series Yt
by non homogenous discrete grey model.
By considering the actual meaning of the original sequence, we have generated a new time series
Z t = (Yt − Xˆ t ) / Yt −1 (t = 2,3, ⋯ , N ) which the fitting curve X̂ t takes as reference
Here Yt −1 is one step of time delay value.
In accordance with the distribution of the random sequence Zt , it is divided by k pieces of state and
is taken its partition value. That is,
e1 = [m0 , m1 ] , e2 = [m1 , m2 ] , , ⋯ , ek −1 = [mk − 2 , mk −1 ] , ek = [mk −1 , mk ]
Here mi ≥ mi −1 .
A one-step transition probability P is associated with each possible transition from stat ei to stat e j ,
and P can be estimated using Pij = M ij / M i (i, j = 1,2,⋯ , m) . M i means the numbers of divided
pieces whose residuals are stat ei , and M ij is number of transition from ei to e j that have occurred
by passed one step. These Pij values can be presented as a transition matrix R .
Then, we accept the statistical quantity such as;
k
k
k
k
Pij
P0 j = ∑ M ij / ∑ M i , χ 2 = 2∑∑ M ij log
P0 j
i =1 j =1
i =1
i =1

Then, when N is comparatively large, the statistical quantity χ 2 is according to χ 2 -distribution
with the degree of freedom (k − 1) 2 .
When giving the confidence degree α , if χ 2 > χα2 (m − 1) 2 , then we confirm that the random
sequence Zt have Markov’s property, unless, the sequence haven’t Markov’s property.
Suppose U is the range which random variable of Markov chain takes the value.
We construct the fuzzy state set S1 , S 2 ,⋯, S l .If for arbitrary u ∈ U the condition
l

∑µ
m =1

Sm

(u ) = 1

is satisfied, and then µ Sm (u ) is called the membership degree of the fuzzy state S M for numerical
valueU .
Suppose that µ Si ( Z t ) ⋅ µ S j ( Z t +1 ) is the fuzzy state transition coefficient from the state Si to S j when
time is turned from t to t + 1 .
Then,
N −1

aij = ∑ µ Si ( Z t ) ⋅ µ S j ( Z t +1 )
t =1

is called fuzzy transition frequency number from the state Si to the state S j .
When the state Zt belongs to Si with degree of member µ Si ( Z t ) and belongs to S j with
membership degree µ S j ( Z t ) , the transition order is only expressed by the product between
membership degrees µ Si ( Z t ) ⋅ µ S j ( Z t +1 ) .

4

Owing to the fuzzy transition probability from the state Si to the state S j is denoted by

Pij =

aij

, ( (i, j = 1, 2, ⋯ , M )

∑a

ij

j =1

, the time series time series which we are going to establish is given such as;
k

Yˆt = Xˆ t +

∑
i =1

k

µ Si (Z t −1 )

1

∑ 2 (m

i −1

+ mi ) Pij Yt −1 , t = 1,2,⋯ , N

j =1

Here, X̂ t is the predicted value obtained by using the discrete grey model.
The forecasting model obtained by above methods is called the discrete grey -fuzzy weight Markov
model. We denote this model by DGM-FMarkov.

3. The forecast of JPY/USD exchange rate based on the discrete grey -fuzzy
weight Markov model.
We take JPY/USD exchange rate of total 278 barter period from September 1994 to July 1997. The
data of this paper is taken from G. Shinarist, P. Waimiso (1997). The concrete data is omitted.
We assume that ŷ k is the forecasted value of the original data obtained by DGM model in k
time. The curve yˆ k = xˆ k( 0+)1 deflects the total change tendency of the original data.
Based on the real phase of the sample data and considering the real meaning, let divide the state
6’ th item, that is,
⊗ = [⊗ , ⊗ ] = [ 
y + a x (0) , 
y + b x (0) ] (i = 1, 2,⋯ , 6)
i

1i

2i

k

i k

k

i k

where ai = (a1 , a2 ,⋯ , am ) = (−0.09, −0.025, −0.01, 0, 0.01, 0.025) ,

bi = (b1 , b2 ,⋯ , bm ) = (−0.025, −0.01, 0, 0.01, 0.025, 0.09)
⊗ :[ 
y − 0.09 x (0) , 
y − 0.025 x (0) ] , ⊗ :[ 
y − 0.025 x (0) , 
y − 0.01x (0) ]
1

k

k

k

k

k

2

k

k

k

⊗3 :[ 
yk − 0.01xk(0) , 
yk ] , ⊗4 :[ 
yk , 
yk + 0.01xk(0) ]
⊗ :[ 
y + 0.01x (0) , 
y + 0.025 x (0) ] , ⊗ :[ 
y + 0.025 x (0) , 
y + 0.09 x (0) ]
5

k

k

k

k

6

k

k

k

k

where ŷ k is the predict value of the foreign exchange rate obtained by DGM model in k time.
The state shows the degree of deviation of the initial data deviated from the curve ŷ k . Where ⊗1 denote that the disadvantage is putted within 2~9% range and ⊗2 disadvantage is putted within
1~2.5% range. And ⊗3 and ⊗4 denote that the profit and disadvantage are putted within 1% range
and ⊗5 denote that the profit is putted within 1~2.5% range and ⊗6 denote that the profit is putted
within 2.5% range.
The states are ⊗1 , ⊗2 , ⊗3 , ⊗4 , ⊗5 , ⊗6 . The numbers of sample point lay at each state are
respectively n1 = 73, n2 = 42, n3 = 29, n4 = 19, n5 = 44, n6 = 71 . The 1-step state matrix composed by
each state is given by
53 14 3 0 2 1 
13 18 9 2 0 0 


 4 9 12 3 1 0 
(nij )6×6 = 

0 2 6 8 3 0
 1 0 0 9 13 10 


 0 0 0 0 23 48

5

The last state transition of the initial sample is uncertain. Therefore, when calculating the state
transition matrix, last data have to exclude. From the formula (15), we can obtain
M 0 j = [0.2626, 0.1511, 0.1043, 0.0683, 0.1583, 0.2554]T
From the formula (16), we can obtain the state transition matrix by
0
 0.7260 0.1918 0.9411
 0.3095 0.4286 0.2143 0.0476

 0.1319 0.1303 0.4138 0.1034
( M ij )6×6 = 
0.1053 0.3158 0.4211
 0
 0.0227
0
0
0.2645

0
0
0
 0
Calculating the statistical quantity, we have
m m
M
χ 2 = 2∑∑ nij log ij ≈ 280.89
M0 j
i =1 j =1

0.0274 0.0137 
0
0 
0.0345
0 

0.1579
0 
0.2955 0.5000 

0.3239 0.7324 

Taking the degree of confidence α = 0.01 and inspecting the tables, then we have
2
2
χ 0.01
(25) = 44.3, χ 2 ≥ χ 0.01
(25) .Therefore, the 1-step transition matrix ( M ij )6×6 then we confirm
that the sequence is coincidence to Markov’s property, unless, the sequence is not coincidence to
Markov’s property. Then, based on this transition matrix, we can predict the foreign exchange rates
which would be putted on the next 1-period or k -period. Therefore, the forecast value of the
exchange rates of the 279’th period is given by
(0)
(0)
(0)
(0)
(0)
1
1
xɵ 279 = [ ɵy 278 + (0.01xɵ 278 + 0.025 xɵ 278 ) p65 + (0.025 xɵ 278 + 0.09 xɵ 278 ) p66 ] = 117.49
2
2
.
From the step-by-step formula the forecast value of the stock price of the period which would be
putted on afterward the state 279’th period is given by
(0)
1
xɵ k +1 = ɵy k + xˆk(0) [0.115 p66 (k + 1) + 0.025 p65 (k + 1)
2
−0.01 p63 (k + 1) − 0.035 p62 (k + 1) − 0.115 p61 (k + 1)] , k ≥ 279

Accuracy was the most important criterion, followed by the cost savings generated from improved
decisions. In particular, execution issues such as ease of interpretation and ease of use were also
highly rated. In this study, there are three criteria used to evaluate forecasting models.
The first evaluate criteria is mean square error, MSE:
N

MSE = ∑ ( xˆ (0) (k ) − x (0) (k ))2 / N

(17)

k =1

where xˆ (0) (k ) is the predicted value at time k , x (0) (k ) is the actual value at time k , and N is the
number of forecasts.
The second criterion is mean absolute error, MAE:
N

N

k =1

k =1

MAE = ∑ e(0) (k ) / N ≡ ∑ xˆ (0) (k ) − x (0) (k ) / N

(18)

The third criterion is mean absolute percent error, MAPE:

1 N (0)
1
q (k ) × 100(%) ≡
∑
N i =1
N
The fourth criterion is Theil Coefficient, µ :
MAPE =

6

N

∑
i =1

xˆ (0) (k ) − x (0) (k )
× 100(%)
x (0) (k )

(19)

µ2 =

[∑ ( x (0) (k ) − xˆ (0) (k ))2 ] / N
[∑ ( xˆ (0) (k )) 2 ] / N

(20)

If x (0) (k ) = xˆ (0) (k ) (1 ≤ k ≤ N ) , i.e., forecast value is coinciding with actual value perfectly, we
have µ 2 = 0 . Therefore, if µ 2 is approached to 0 , then forecast method become finer.

4. Grey Neural Network
The Grey GM (1, 1) model was established forecast model by using new data accumulated
generator. The accumulated generator data can make weak the randomness a certain degree and
can easily find the data change rule. The grey GM (1, 1) model has the advantage which demands
the small sample data.
The neural network has the capacity of self-learning, nonlinear mapping, and parallel distribution
processing, etc.
These tow forecast method have already used in foreign exchange rate predicting.
The system of foreign exchange rate forecast is a complicated system which have the great
randomness and there are many index influenced in that system.
Combining grey system idea and neural network, we can construct the grey neural network (GNN)
and sufficiently demonstrate the advantage of the model by using the modeling method of grey
system with the small data and the characters of nonlinear mapping of neural network.
Thus, we can more enhance forecast accuracy.
Here we discuss the method of combing the GM model and neural network model. There are three
kinds of forecasting model structure. That is, Parallel grey neural network (PGNN), series grey
neural network (SGNN), and inlaid grey neural network (IGNN).

4.1. Parallel grey neural network (PGNN)
PGNN uses grey model and neural network to predict separately, then uses neural network to
combine the predicting results. Fig. 1 shows the principle scheme of PGNN
GM(1,1)

model

Hybrid
Forecast
model

Initial data
Neural network model

Forecast results

Fig.1 Parallel grey neural network model

PGNN model is essential is the hybrid model. The aim of making the PGNN model is to decrees the
randomness of data and avoids data lack used by single model, and enhances the forecast accuracy
degree by using data offered from various kinds methods synthetically.
On the ground of hybrid forecast principle, there are various hybrid methods that is, arithmetic
mean hybrid method, geometric mean hybrid method, harmonic mean hybrid method, etc. Those
hybrid formulas are given by
xɵ (k ) = w xɵ 1 (k ) + w xɵ 2 (k ) , k = 1, 2,⋯ , N
1

2

xɵ (k ) = xɵ 1 (k ) w1 + xɵ 2 (k ) w2 , k = 1, 2,⋯ , N
1
xɵ (k ) =
, k = 1, 2,⋯ , N
w / xɵ 1 (k ) + w / xɵ 2 (k )
1

2

7

Where N is total forecast data number, xɵ 1 (k ) and xɵ 2 (k ) are the forecast value by GM(1,1) model
and neural network model, respectively, and w1 and w2 are weights of two kind of forecast method.
xɵ (k ) is actual value.
Now, we consider the conception of effective degree. The conception of effective degree has some
rationality because it reflects the efficient of the methods.
Its idea is as follow;
We put
x(k ) − xˆ (k )
x(k ) − ( w1 xɵ 1 (k ) + w2 xɵ 2 (k ))
= 1−
A(k ) = 1 −
x(k )
x(k )
Then, A(k ) is the accuracy series of hybrid forecast. We denote the mean value E and mean square
deviation of A(k ) respectively by
1 N
1 N
E = ∑ A(k ) , σ = (∑ ( A(k ) − E ) 2 )1/ 2 .
N k =1
N k =1
We define effective degree of hybrid forecast method by
S = E (1 − σ )
If S is more great, then forecast accuracy is more raised and forecast error is more safety. It says
that the model is more effective.
To find the weights w1 , w2 , we can use optimization method by
inf S ( w1 , w2 ) ≡ S = E (1 − σ )
{ w1 , w2 }

But, this method is very complicated.
Now, we would use a simple method by using of physical meaning.
Let A1 (k ) and A2 (k ) are the series of forecast accuracy by grey GM (1, 1) model method and neural
network model method respectively, that is,
x(k ) − xˆi (k )
A1 (k ) = 1 −
, i = 1, 2; k = 1, 2,⋯ , N
x (k )
Let Si are the effective degree by GM (1, 1) Model and neural network model respectively. Then,
we can defined weights w1 and w2 by
S
wi = 2 i (i = 1, 2) .
∑Sj
j =1

4.2. Series grey neural network (SGNN)
Series grey neural network employ grey model to predict, then use neural network to combine the
predicting results.
Established forecast model GM(1,1) with each other data for already given the same series, then the
gained forecast results are different each other.
To obtain forecast results approached to actual value, for various grey forecast results we can
combine its results by using neural network. This is immediately SGNN model.
Fig. 2 shows the principle scheme of SGNN.
The input nerve element number of neural network is the number of each other GM model and the
number of output nerve element number of neural network is single element. The hidden nerve
element number is confirmed by test methods. The learning of neural network is progressed by error
back-propagation method.
PGNN and SGNN are all essential hybrid forecast models.
Theoretically, we can prove that the forecast results of hybrid model are advantage than single GM
model or neural network model.
8

In SGNN model, we find the combined weights of several kinds grey model by using the nonlinear
matching ability of neural network. The combination of SGNN is nonlinear, the other side, the
combination of PGNN is linear.
GM1 model

Neural Network
Model

GM 2 model

Initial data

Forecast results

⋮
GM m model
Fig.2 Series grey neural network model

4.3. Inlaid grey neural network (IGNN)
IGNN model is built by adding a grey layer before neural input layer and a white layer after
neural output layer. Fig. 3 shows the principle scheme of IGNN

Initial data

Grey layer

Neural network

white layer

Forecast results

Fig.3 Inlaid grey neural network model

The action of grey layer is to weaken the randomness of initial data. Generally, the role of grey
layer is to make training sample of neural network with new data gained from 1-AGO for initial
data. In that case, the approximation in nonlinear excitation function of neural network is easy and
studying time of network is very short, and convergence process is quick and forecast accuracy is
raised because accumulated operator data have monotone increasing tendency.

5. The hybrid forecast method
5.1 The determination of combinational weight based minimization of error-squared
criterion
In the actual problem, we can predict for forecasting problem by using the various forecast
methods. Generally speaking, the forecast result forecasted by the distinctive forecast method is
different each other. To obtain the more improved forecast result, it is necessary to apply the
various kinds of hybrid forecast methods.
Let x (0) (k ) (k = 1, 2,⋯ , N ) are an actual observation value and xˆi (k ) (k = 1, 2,⋯ , p ) are the
predict value by the i ' th method. Let wi is the weight of the i ' th method and ei (k ) = x (0) (k ) − xˆi (k )
Is the forecast error. From this notation, we have
p

x (0) (k ) = ∑ ( wi xˆi (k ) + ei (k )) (k = 1, 2,⋯ , N ) .
i =1

Now, let wˆ i is the calculated value of the weight and xˆ (0) (k ) is the hybrid forecast value. Then we
can obtain
p

(0)
 i xɵ i (k )
xɵ (k ) = ∑ w
i =1

*
i

To find the optimal weight w , we can use various optimization methods. For example,
find optimal weight wi* for optimization problem;

9

N

p

k =1

i =1

J [ wi* ] = min J [ wi ] = ∑ ( x (()) (k ) − ∑ wi xˆi (k )) 2
wi

(21)

p

∑w

i

= 1 , wi ≥ 0 , i = 1, 2,⋯ , p

i =1

To obtain the solution of optimization problem (21), we can apply various optimization methods
, for example, quadratic programming, genetic algorithm, neural network method, etc.
(0)
(0)
For instance, xˆ1 (k ) and xˆ2 (k ) are the predict value by two methods for the same phenomena.
Then,
(0)
(0)
(0)
(0)
(0)
xɵ (k ) = w xɵ 1 (k ) + w xɵ 2 (k ) = ρ xɵ 1 (k ) + (1 − ρ ) xɵ 2 (k )
1

2

(0)
(0)
is called the hybrid forecast model for two models xˆ1 (k ) and xˆ2 (k ) .Where 0 ≤ ρ ≤ 1 .

Now, we put e1 = e1(0) (k ) , e2 = e2(0) (k ) , from the hybrid model, we have e = ρ e1 + (1 − ρ )e2 .
The square deviation of e is given by
D (e) = ρ 2 D (e1 ) + (1 − ρ 2 ) D (e2 ) + 2 ρ (1 − ρ ) cov(e1 , e2 ) .
dD (e)
= 0.
Where cov(e1 , e2 ) is covariance of e1 and e2 . To obtain least value of D (e) , we put
dρ
Then, we can obtain the least value of D (e) by
D (e1 ) − cov(e1 , e2 )
ρ* =
.
D (e) = D (e1 ) + D (e2 ) − 2 cov(e1 , e2 )
That is, when ρ = ρ * , D (e) is achieved to the least value. In the actual calculation, we can
approvable that e1 and e2 are the mutual independence random quantity. Therefore, we
have cov(e1 , e2 ) = 0 . Accordingly, we have
D (e1 )
ρ* =
(22)
D (e) = D (e1 ) + D (e2 )
The final forecast model is given by
(0)
(0)
(0)
xɵ (k ) = ρ * xɵ 1 (k ) + (1 − ρ * ) xɵ 2 (k )
(23)

5.2. The determination of combinational weight based on optimal grey relation degree
method
Let { y (k ) k = 1,2,⋯, N } be original value of time series, yˆ j (k ) (i = 1,2,⋯, m, k = 1,2,⋯, N ) be predicted
value of the j ’th forecast method in k time ( j = 1,2, ⋯, m, t = 1,2,⋯, N ) .
We put
γ0j

1
=
N

N

∑
t =1

min min e j (t ) + ρ max max e j (t )

1≤ j ≤ m 1≤t ≤ N

1≤ j ≤ m 1≤t ≤ N

e j (t ) + ρ max max e j (t )
1≤ j ≤ m 1≤t ≤ N

Then, γ 0 j is called the grey relation degree of the predicted value sequence
Here ρ ∈ (0,1) is called the identification coefficient. Generally we take ρ =0.5.
e j (t ) = y (t ) − yˆ j (t ) is the predicted error of k time for i ’th forecast method.
Let ŷ( t ) =

m

∑w

j

yˆ j (t ) be the predicted value of y (t ) by the hybrid forecast method.

j =1

Here w1 , w２, ⋯，wm is the weight coefficient of m kind of forecast method, which satisfies
ｍ

∑ w ＝１, w j ≥ 0
j

( j = 1,2,⋯ , m) .

j＝１

The combinational weight based on the optimal grey relation degree is determined such as;

10

1
max γ (W ) =
N

N

∑
t =1

min min e j (t ) + ρ max max e j (t )

1≤ j ≤ m 1≤ t ≤ N

1≤ j ≤ m 1≤ t ≤ N

m

∑w e
j

j

(t ) + ρ max max e j (t )
1≤ j ≤ m 1≤ t ≤ N

j =1

ｍ
 ∑ w j＝１
s.t.  j＝１
 w ≥ 0, j = 1,2, ⋯ , m
 j

By solving this optimal problem, we can determine the combinational weight
W = {w 1 , w 2 , ⋯ , w m }T based on the optimal grey relation degree.

6. The hybrid forecast of exchange rates by DGM-FMarkov model and Grey Neural Network
(0)
(0)
Now, let xˆ1 (k ) is the predict value obtained by DGM-FMarkov model and xˆ2 (k ) is the
forecasted value obtained by inlaid grey neural network model.
We have to obtain the calculating results by DGM-FMarkov model and grey neural network model
for exchange rates.
The structure of neural network is taken 4 × 4 × 1 BP network. The transfer function of input layer
and hidden layer are sigmoid type and transfer function of output is linear function.
In GNN forecast, we have applied IGNN mode. We have progressed 1-AGO smoothing handling
by adding a grey layer before input layer of BP neural network and have progressed 1-IAGO
reduction handling for white layer after neural output layer.
We have calculated according to formula (22), (23) by unbiased grey-Markov- grey neural network
hybrid forecast value for foreign exchange rates.
The mean absolute percent error MAPE of DGM-FMarkov model is given by
1 N xˆ ( 0 ) (k ) − x ( 0) (k )
MAPE = ∑
× 100% ≈ 3.65(%)
N i =1
x ( 0) (k )
Theil coefficient of DGM-FMarkov model is given by µ = 0.0053
The mean absolute percent error MAPE by hybrid model of DGM-FMarkov model - grey neural
network on the basis of the determination of combinational weight based minimization of errorsquared criterion is given by
1 N xˆ ( 0 ) (k ) − x ( 0) (k )
MAPE = ∑
× 100% ≈ 2.17(%)
N i =1
x ( 0 ) (k )

Then, Theil coefficient is given by µ = 0.0032
The mean absolute percent error MAPE by hybrid model of DGM-FMarkov model - grey neural
network on the basis of the determination of combinational weight based on optimal grey relation
degree method is given by
1 N xˆ ( 0 ) (k ) − x ( 0) (k )
MAPE = ∑
× 100% ≈ 1.85(%)
N i =1
x ( 0) (k )
Then, Theil coefficient is given by µ = 0.0021

11

7. Conclusion
As known from the hybrid forecast model, its model can use the information of the original data
sufficiently and absorb the characters and the advantages of two models, and avoid the limitation of
single model.
The forecast accuracy of hybrid model by DGM-FMarkov model - grey neural network is more
higher than the DGM-FMarkov model. Especially, The forecast accuracy of hybrid model by
DGM-FMarkov - grey neural network model on the basis of the determination of combinational
weight based on optimal grey relation degree method is more higher than the DGM-FMarkov-grey
neural network model on the basis of the determination of combinational weight based
minimization of error-squared criterion
It is shown that this model has high forecast accuracy for forecast of exchange rates with tendency
and fluctuation.
In this paper we have researched an hybrid forecast model by using DGM-FMarkov model and
grey neural network forecast model and have exhibited that this method is efficient to the exchange
rates forecast
In the hybrid forecast model, we can use some forecast models, fore example; not only forecast not
only DGM-FMarkov forecast, grey-regression forecast, but also grey- chaotic forecast and greyfractal forecast ([13, 14]), ect.

References
[1] M. Duraira, K. Meena (2011), A Hybrid Forecast System Using Rough Sets and Artificial
Neural Networks, International Journal of Innovative Technology & Creative Engineering,1(7)
16-23
[2] Reza Askari Moghadam, Mehrnaz Keshmirpour (2011), Hybrid ARIMA and Neural Network
Model for Measurement Estimation in Energy-Efficient Wireless Sensor Networks,
Communications in Computer and Information Science, 253, 35-48
[3] Aobing Sun, Yubo Tan, Dexian Zhang(2008) , Hybrid Forecast Model Based on BP Neural
Network for Lung Cancer, Proceedings of 2008 IEEE International Symposium on IT in
Medicine and Education, 532-535
[4] Cagdas Hakan Aladag, Erol Egrioglu, Cem Kadilar(2012), Improvement in Forecasting
Accuracy Using the Hybrid Model of ARFIMA and Feed Forward Neural Network, American
Journal of Intelligent Systems, 2(2): 12-17
[5] Kumaran Kumar. J, Kailas (2012). A Forecast of Future Stock Close Price using Proposed
Hybrid ANN Model of Functional Link Fuzzy Logic Neural Model, IAES International Journal
of Artificial Intelligence,1(1), 25~30
[6] Tarek Aboueldahab, Mahumod Fakhreldin(2011), Forecast of Stock Market Indices using
Hybrid Genetic Algorithm/ Particle Swarm Optimization with Perturbation Term, ICSI 2011:
International Conference on swarm intelligence Cergy, France, 1-10
[7] Florin Leon, Mihai Horia Zaharia (2010), Stacked Heterogeneous Neural Networks for
Time Series Forecasting, Mathematical Problems in Engineering, 20, 1-19
[8] Šterba Ján, Hil’ovsk Katarína (2010), The Implementation of Hybrid ARIMA-Neural Network
Forecast Model for Aggregate Water Consumption Forecast, Aplimat – Journal of Applied
Mathematics, 3(3), 123-130
[9] Zhu Jiang-liang (2006), The forecasting of the electric power load based on the grey theory and
BP neural network, Electric Machines and Control, 10(4), 440-442
[10] Su Bo (2006), Comparison and research of gain production forecasting with methods of
GM (1, N) grey system and BPNN, Journal of China Agricultural University, 11(4), 99-104,
[11] Yao Hong-Xing, Shong Zhao-han, Chen Hong-xiang(2002), Method of the wavelet neural
network in the forecast of stock market, Systems Engineering-Theory & Practice, 6, 35-38,

12

[12] Kyong-jac Kim, Ingoo Han(2000), Genetic algorithms approach to feature discretization in
artificial neural networks for the forecast of stock price, Expert Systems with Applications, 19,
125-132,
[13] Kim Gol (2008), Forecasting of foreign exchange rate by grey- chaotic forecast hybrid
model, Journal of University of Science, 4, 33-40

[14] Kim Gol (2009), Research of forecasting of stock price and index of oil by grey- fractal
forecast hybrid model, Information Science, 2, 28-35
[15] G. Shinarist, P. Waimiso (1997), Excepting effect of sample in admission fluctuation foreign
exchange rate model, International finance papers, 301, 16-21

13

