概述
z3 v: I8 J$ ?8 N众多外汇和加密货币交易者都希望从短期走势当中获益,但缺乏可以帮助他们奋斗的基本面信息,因而备受困扰。 而这恰是标准时间序列技术可以提供帮助的地方。 George Box 和 Gwilym Jenkins 开发出可以说是最受尊崇的时间序列预测方法。 尽管已经取得了许多大幅进步,改进了原始方法,但其基本原则在今天依旧适用。" u* O W/ A8 E4 E( }) \2 A- F
他们方法的衍生物之一是自回归综合移动平均线(ARIMA),其已成为时间序列预测的流行方法。 它是一类模型,用于捕获数据序列中的时间依赖关系,并提供了针对非稳态时间序列进行建模的框架。 在本文中,我们将采用函数最小化 Powells 方法作为基础创建 ARIMA 训练算法,并利用 mql5 编程语言实现。赫兹量化交易软件
5 l( o. O" R6 O9 }& AARIMA 概览) h4 u& l0 b3 R+ f! x$ K: c% o6 A6 e/ K
Box 和 Jenkins 表示,大多数时间序列可以由两个框架中之一个或两个来建模。 其一是自回归(AR),这意味着序列的值可以由与其先前的值相关性来解释,而当中的恒定偏移量和微小的差异,通常称为创新或噪声。 请注意,在本文中,我们将噪声或误差部分称为创新。 创新解释了无法解释的随机变化。赫兹量化交易软件
( H3 d, a: A5 d2 U, {: T4 IARIMA 模型的第二个框架是移动平均线(MA)。该模型指出,序列的值是特定数量的前期创新项、当前创新、和恒定偏移量的比例和。 还有许多其它统计条件可以定义这些模型,但我们不会深入研究细节。 网上有许多资源可以提供更多信息。我们对它们的应用更感兴趣。$ j8 G5 ^8 [$ ^2 X
我们不会受限于纯 MA 和 AR 模型,我们可以将它们组合在一起,生成称为自回归移动平均模型(ARMA)的混合模型。 在 ARMA 模型当中,除了恒定偏移量和当前创新项外,我们还指定了有限数量的滞后序列和噪声项。赫兹量化交易软件 W. J* I; c, g: L% j" ?
影响所有这些框架应用的基本要求之一是,被建模的序列必须是稳态的。 取决于您对稳态定义的严格程度,到目前为止所描述的模型在技术上都不适合应用于金融时间序列。 这就是 ARIMA 的用武之地。 数学积分是微分的逆向。 当非稳态时间序列相差一次或更多次时,成果序列通常具有更好的稳态。 首先对序列进行微分,可以将这些模型应用于成果序列。ARIMA 中的 “I” 是指逆向(积分)所应用的差分的要求,以便将建模的序列返回到其原始域。赫兹量化交易软件. S/ G- o8 P4 x6 [" z
自回归模型标注法
& G' x m& @! Y有一个标准标注法来管控模型的描述。 AR 项(不包括常数项)的个数通常称为 p。 MA 项表示为 q,且 d 表示原始序列差异的次数。 使用这些术语,ARIMA 模型被指定为 ARIMA(p,d,q)。 纯过程可以描述为 MA(q) 和 AR(p)。 没有差分的混合模型写为 ARMA(p,q)。 此标注法假定这些术语是连续的。 例如,ARMA(4,2) 表示序列可以用 4 个连续的 AR 项,和两个连续的创新项来描述。 采用 ARIMA,我们可以通过将 p、q 或 d 指定为零值来描绘纯过程。 例如,ARIMA(1,0,0) 简化为纯 AR(1) 模型。赫兹量化交易软件
( c* f; q1 Y' \大多数自回归模型指定各自的项是连续的,分别是 AR 项从滞后 1 到滞后 p,及 MA 项滞后 q。 将要演示的算法将允许指定 MA 和/或 AR 项的非连续滞后。 该算法将引入的另一个灵活性是能够指定模型是否带有恒定偏移量。赫兹量化交易软件: w7 z: D5 h3 L
例如,可以构建由以下函数定义的模型:- ]5 v$ [& P% @+ V' ~
y(t) = AR1* y(t-4) + AR2*y(t-6) + E(t) (1)9 a+ Y! q0 D$ x' m( w0 C
上面的函数描述了一个纯 AR(2) 过程,该过程没有恒偏移量,当前值由前期第 4 和第 6 时隙的序列值定义。 标准标注法没有提供指定这种模型的方法,但我们不必被这些限制所束缚。
: ]/ I( K2 U! p- [; S6 ^- E计算模型系数和恒定偏移量5 f4 r- _5 G8 ? p
模型可以具有必须经计算的 p+q 系数。 为此,我们使用模型的规范对已知序列值进行预测,然后将预测值与已知值进行比较,并计算误差的平方和。 最优系数将是产生最小误差平方和的系数。
2 e1 W- T/ x4 ^. e: n在进行预测时必须小心,无限延伸的数据并不可用,因此会施加限制。 如果模型规范有任何 AR 项,我们只能跳过与所有 AR 项的最大滞后值对应的数值个数之后,再开始进行预测。5 e* Z& X& L: F( J& Z" m
使用上面(1)指定的示例,我们只能从时隙 7 开始进行预测。 因为任何超前的预测都会在序列开始之前引用未知值。
- ^; G( d# c1 T" o( t& K6 C应该注意的是,如果(1)有任何 MA 项,此时该模型将被视为纯自回归,因为我们还没有创新序列。 随着预测的推进,一系列的创新数值就会不断积累。 回到这个示例,第 7 处时隙的第一个预测将利用人工初始 AR 系数进行计算。
, ~. k- k/ O& z计算出的预测值与第 7 处时隙的已知值之间的差值就是该时隙的创新值。 如果指定了任何 MA 项,则当相应的创新滞后值已知时,在计算预测时将包含它们在内。 否则,MA 项将清零。 在纯 MA 模型的情况下,遵循类似的过程,除了这次,如果应包含恒定偏移量,则将其初始化为序列的均值。" v$ v. X8 ~' u( a( J
刚才讲述的方法只有一个明显的局限性。 已知序列需要包含与其所应用模型按顺序等量的数值。 更多的项和/或滞后越大,我们需要的值就越多,才能有效地拟合模型。 然后,通过应用相应的全局最小化算法来优化系数,从而完善训练过程。 我们将采用的最小化预测误差的算法是 Powells 方法 。 该实现于此应用的详情记录在文章利用指数平滑进行时间序列预测之中。7 c6 Q ^/ ?% u# s7 n: t
CArima 类6 H- ]1 ?: |, x- P, i7 u) a% Z, p
ARIMA 训练算法将被包含在 Arima.mqh 文件定义的 CArima 类当中。 该类有两个构造函数,每个构造函数初始化一个自回归模型。 默认构造函数创建一个具有恒定偏移量的纯 AR(1) 模型。
! l) e# s6 V B4 VCArima::CArima(void)- v. ~+ a2 X& W# s2 h
{7 T/ O. }- u1 a
m_ar_order=1;" M% f" l( e! q+ W/ E1 c
//---
/ k) d5 i! R9 a% rArrayResize(m_arlags,m_ar_order);1 v0 I! o, ?. |5 H8 [2 _
for(uint i=0; i<m_ar_order; i++)( o/ Z. k9 q5 a2 `$ y k/ l9 c
m_arlags=i+1;
. a, p9 T$ K0 h! ~* m//---
1 u7 _7 v4 z6 X8 i. k* zm_ma_order=m_diff_order=0;
* a, T6 P8 X! a0 `+ t* N4 bm_istrained=false;
2 k, c3 y. s) H& d$ {! p* wm_const=true;
& R s D" e) A; s6 N- A) EArrayResize(m_model,m_ar_order+m_ma_order+m_const);
, {! ]0 a9 k, }, q4 F+ vArrayInitialize(m_model,0);
2 Q" T$ z O I, ]& W. e}
; |! _+ u; ?) Q1 r4 h5 }参数化构造函数允许在指定模型时进行更多控制。 它需要下面列出的四个参数:- t& ]3 j( A9 ?4 Y1 I: ]
参数) m5 s; q1 p7 D
参数类型- W$ G) t, e+ X
参数说明$ K9 d0 C2 P u
p+ E$ J6 `+ H& }" V
无符号整数1 W( T% C( ?5 w2 |3 C! q
指定了模型的 AR 项数1 C( B3 L+ |/ T, ?: i+ K/ P8 P
d
w- t) {; k, B( F. n; M# o无符号整数/ D) h6 b/ z9 r& q+ x
指定要应用于正在建模的序列的差异程度
: I2 s/ y. |9 k+ ]* W# Pq- V2 o1 [# [; c$ s# U) U4 k
无符号整数
; g0 |# C2 n* @" a- u指示模型应包含的 MA 项数
F/ s; f: U3 N# J; K- ]0 uuse_const_term
0 p% @0 _$ Y- Q, W! p布尔值4 C. H$ L; _. ~
设置启用模型中恒定偏移量
7 M; f+ R0 a' P* v( I3 NCArima::CArima(const uint p,const uint d,const uint q,bool use_const_term=true)7 Z" r6 G8 F7 i5 c! J) ]! L* ]
{( X* p& A( Q& o: r: k! o! n
m_ar_order=m_ma_order=m_diff_order=0;. K, w, C* f8 I$ T4 Y) o8 g
if(d). D- F a0 s3 K ^: {6 Q# h/ Y% \
m_diff_order=d;5 H4 o! ?* ?+ i5 x; L' z/ {
if(p)
7 y t6 g7 [$ y1 E0 e. r{3 t+ Q4 `: E$ y) z
m_ar_order=p;
; R0 C* ~7 s( n# G& s" E3 jArrayResize(m_arlags,p);
J/ p' E. n, w- s. r" J% ^for(uint i=0; i<m_ar_order; i++)
8 c' `% f7 R) H! j; D$ nm_arlags=i+1;
& m! @6 M8 L/ {5 r}
* m& K3 |+ k# X/ O; y$ V0 \: W$ L+ iif(q): _( M: K, C! Z; O4 T/ B
{
7 h4 l# z4 F2 Fm_ma_order=q;# G; l* s3 |1 r3 H% ?$ ^
ArrayResize(m_malags,q);2 K' a) U( S+ v( e) z/ s; y
for(uint i=0; i<m_ma_order; i++). h* m8 O( c- o+ L7 Y' O; N' A
m_malags=i+1;
& O! [8 V0 b, H% Z}
! P, I, y! b" Rm_istrained=false;
3 _6 G2 Y' c1 a* nm_const=use_const_term;
' ~* l8 _! }, k+ HArrayResize(m_model,m_ar_order+m_ma_order+m_const);, I% O& L0 p( J/ Y
ArrayInitialize(m_model,0);7 g# _, Z# I* B9 ~/ f9 B, l
}4 s, Y# E$ I4 o5 `2 U( ~
除了提供的两个构造函数外,还可以使用重载的 Fit() 方法之一来指定模型。 这两种方法都把建模的数据序列作为其第一个参数。 第一个 Fit() 方法只有一个参数,而第二个方法需要更多的四个参数,所有这些参数都与上表中所记录的参数相同。
1 X* n$ h$ s- L% p+ Jbool CArima::Fit(double &input_series[])
6 I! W) E$ h0 Q" K6 e{
( V/ b# L1 L: b" Q+ C- ]3 _1 ouint input_size=ArraySize(input_series);# {$ C8 V6 U7 |* P6 H
uint in = m_ar_order+ (m_ma_order*2);8 a) p, x5 v, {1 z
if(input_size<=0 || input_size<in)
1 M7 M7 R- w! S+ J3 Zreturn false;/ C T( L) t/ t0 x0 E
if(m_diff_order)
( E& V* \: A0 W9 L/ Odifference(m_diff_order,input_series,m_differenced,m_leads);! V( w$ L0 _+ g3 ^4 S
else
2 q& V: G6 {# e% x7 w" p, T2 {ArrayCopy(m_differenced,input_series);2 \- r5 P4 A( q/ [ \
ArrayResize(m_innovation,ArraySize(m_differenced));
. m& S4 X: S6 ^- J& l1 w. e3 Bdouble parameters[];
$ G0 ~+ ]8 e+ Y! {+ W' ]6 ?& qArrayResize(parameters,(m_const)?m_ar_order+m_ma_order+1:m_ar_order+m_ma_order);
) R- e( \# C' _5 L0 LArrayInitialize(parameters,0.0);
( Z& U( V/ H; ^& r% rint iterations = Optimize(parameters);
2 C# B# A3 |* x, ~if(iterations>0)7 u U, }- g! T- }4 G
m_istrained=true;! [6 x3 }! r- X3 t9 P" `" y- `
else# f! w! E* \0 \8 |+ E3 A. a( e& L+ e
return false;
R. O9 ?, t: u% v- ~m_sse=PowellsMethod::GetFret();( D$ t L8 N! u( i5 R
ArrayCopy(m_model,parameters);/ G1 h" K" g, ~8 c, L. m
return true;
$ P6 h1 Z2 V3 W5 E9 j$ D}9 @6 }+ Q3 s1 z) g& a
调用更多参数的方法会覆盖以前指定的任何模型,并且还假定每一项的滞后是彼此相邻的。 两者都假定第一个参数指定的数据序列并无差别。 因此,如果由模型参数指定,则会应用差别。 这两种方法都返回一个布尔值,表示模型训练过程的成功或失败。
0 Q2 F( i) @5 _& l6 O( Nbool CArima::Fit(double&input_series[],const uint p,const uint d,const uint q,bool use_const_term=true)
5 m; j* b2 G2 F H4 W{5 S. ]& [& H0 C8 R* W
m_ar_order=m_ma_order=m_diff_order=0;
9 i M( k, g4 uif(d)) M) T% f$ u! i( t' j
m_diff_order=d;& ^, E* X. Y8 L$ B ^
if(p)3 O0 i( ?' r" e6 U9 ^2 H
{
: d+ O1 O: x7 ^4 B& I- {; rm_ar_order=p;
+ W3 H$ _% I" m+ P' e2 L$ gArrayResize(m_arlags,p);
6 Z; o; `4 Q) P! Ifor(uint i=0; i<m_ar_order; i++)- y4 @1 X) y6 r( ?8 I7 C7 A
m_arlags=i+1;& T2 {/ ~2 w3 k2 }- s+ a- w
}
) U, `% r' q, k& f: k; B' cif(q)
0 y l% d6 F$ y+ x9 }2 ]: Y' q$ a{; F7 r2 b' k+ Z; e6 v9 N2 G* s
m_ma_order=q;
8 b; C# w! g3 @# ~2 eArrayResize(m_malags,q);2 w5 \$ K! k' Y; ?. ]2 W' A: P
for(uint i=0; i<m_ma_order; i++)
1 `$ [8 a: H" Vm_malags=i+1;; T3 T( x6 B6 @
}6 G1 V- X+ p" I" F. y( e; w
m_istrained=false;
6 e( `( _+ N. z7 ]2 p3 W8 Q! J4 H& K: fm_const=use_const_term;
7 D& R6 s# g% U/ S% ^/ o1 \ZeroMemory(m_innovation);1 w# Z: Q7 Q( m" u3 f
ZeroMemory(m_model);
+ B+ U3 ]% e6 n3 U& }8 hZeroMemory(m_differenced);& j* m5 x$ [ E( q% a& p% r1 F$ \
ArrayResize(m_model,m_ar_order+m_ma_order+m_const);
) K1 C7 t" S' q: IArrayInitialize(m_model,0);1 |) q' P/ a# N
return Fit(input_series);- j" g L7 {3 X$ r# z
} |