当前位置:   article > 正文

12灰色模型的检验_灰色预测模型事前检验

灰色预测模型事前检验

灰色模型的检验

   前两篇博客说了灰色预测的基本步骤,其中讲了两种检验模型的方法,本次再具体地说一下灰色模型的检验方法。

   设有 n n n 个原始数据 x ( 0 ) = ( x ( 0 ) ( 1 ) , x ( 0 ) ( 2 ) , . . . , , x ( 0 ) ( n ) ) x^{(0)}=(x^{(0)}(1), x^{(0)}(2),...,, x^{(0)}(n) ) x(0)=(x(0)(1),x(0)(2),...,,x(0)(n)),以及由灰色模型预测的值 x ^ ( 0 ) = ( x ^ ( 0 ) ( 1 ) , x ^ ( 0 ) ( 2 ) , . . . , , x ^ ( 0 ) ( n ) ) \hat x^{(0)}=(\hat x^{(0)}(1), \hat x^{(0)}(2),...,, \hat x^{(0)}(n) ) x^(0)=(x^(0)(1),x^(0)(2),...,,x^(0)(n))

1 残差检验
   令残差 ϵ ( 0 ) ( k ) \epsilon^{(0)}(k) ϵ(0)(k)为:

ϵ ( 0 ) ( k ) = x ( 0 ) ( k ) − x ^ ( 0 ) ( k ) , k = 1 , 2 , 3 , . . . , n (1) \epsilon^{(0)}(k)={x^{(0)}(k) - \hat x^{(0)}(k)},k=1,2,3,...,n \tag{1} ϵ(0)(k)=x(0)(k)x^(0)(k),k=1,2,3,...,n(1)

ϵ ( 0 ) = ( ϵ ( 0 ) ( 1 ) , ϵ ( 0 ) ( 2 ) , . . . , , ϵ ( 0 ) ( n ) ) (2) \epsilon^{(0)} = (\epsilon^{(0)}(1), \epsilon^{(0)}(2),...,, \epsilon^{(0)}(n) ) \tag{2} ϵ(0)=(ϵ(0)(1),ϵ(0)(2),...,,ϵ(0)(n))(2)

   令相对误差 Δ ( 0 ) ( k ) \Delta ^{(0)} (k) Δ(0)(k)为:

Δ ( 0 ) ( k ) = ∣ ϵ ( 0 ) ( k ) x ( 0 ) ( k ) ∣ = ∣ x ( 0 ) ( k ) − x ^ ( 0 ) ( k ) x ( 0 ) ( k ) ∣ , k = 1 , 2 , 3 , . . . , n (3) \Delta ^{(0)} (k)=\left | \frac {\epsilon^{(0)}(k)} {x^{(0)}(k)} \right|=\left | \frac {x^{(0)}(k) - \hat x^{(0)}(k)} {x^{(0)}(k)} \right|,k=1,2,3,...,n \tag{3} Δ(0)(k)=x(0)(k)ϵ(0)(k)=x(0)(k)x(0)(k)x^(0)(k),k=1,2,3,...,n(3)

Δ ( 0 ) = ( Δ ( 0 ) ( 1 ) , Δ ( 0 ) ( 2 ) , . . . , , Δ ( 0 ) ( n ) ) (4) \Delta ^{(0)} = (\Delta ^{(0)} (1), \Delta ^{(0)} (2),...,, \Delta ^{(0)} (n) ) \tag{4} Δ(0)=(Δ(0)(1),Δ(0)(2),...,,Δ(0)(n))(4)

   对于 k ≤ n k\leq n kn,定义 k k k点模拟相对误差为:

Δ k = ∣ ϵ ( 0 ) ( k ) x ( 0 ) ( k ) ∣ (5) \Delta _k=\left| \frac{\epsilon^{(0)}(k)} {x^{(0)}(k)} \right| \tag{5} Δk=x(0)(k)ϵ(0)(k)(5)

   定义平均相对误差:

Δ ˉ = 1 n ∑ k = 1 n Δ k (6) \bar \Delta= \frac{1}{n} \sum_{k=1}^n \Delta_k \tag{6} Δˉ=n1k=1nΔk(6)

   称 1 − Δ ˉ 1− \bar \Delta 1Δˉ 为平均相对精度,称 1 − Δ k 1− \Delta_k 1Δk k k k 点模拟精度;给定一个阈值 α \alpha α ,如果 Δ ˉ < α \bar \Delta < \alpha Δˉ<α Δ k < α \Delta_k < \alpha Δk<α成立时,称模型为残差合格模型。

2 级比偏差值检验
   计算级比 λ ( 0 ) ( k ) \lambda^{(0)}(k) λ(0)(k)(可参考之前博客),已经在模型建立阶段计算了发展系数 a a a,那么计算相应的级比偏差:

ρ ( k ) = 1 − ( 1 − 0.5 a 1 + 0.5 a ) λ ( 0 ) ( k ) (7) \rho(k)=1- \left ( \frac{1-0.5a}{1+0.5a} \right) \lambda^{(0)}(k) \tag{7} ρ(k)=1(1+0.5a10.5a)λ(0)(k)(7)

   如果所有的 ρ ( k ) \rho(k) ρ(k)均小于0.2,则认为模型可行,当然 ρ ( k ) \rho(k) ρ(k)越小越好。如果所有的 ρ ( k ) \rho(k) ρ(k)均小于0.1,则认为模型达到较高的要求。

3 关联度合格检验

   之前的博客介绍了怎么计算两个数列的关联度,如果两个数列差别越小,那么关联度越高。现在我们有原始数列以及由灰色模型预测的数列,我们的目标是让预测的数列和原始数列尽可能一样,那也就是说两个数列的关联度要越大越好。设 g ( 0 ) g^{(0)} g(0)为原始数列 x ( 0 ) x^{(0)} x(0)和预测数列 x ^ ( 0 ) \hat x^{(0)} x^(0)的绝对关联度,若 g ( 0 ) g^{(0)} g(0)大于给定的 g 0 ( 0 ) g^{(0)}_0 g0(0),则称模型为关联度合格模型。

4 均方差比合格模型

   定义 x ( 0 ) x^{(0)} x(0)的均值和方差分别为:

x ˉ = 1 n ∑ k = 1 n x ( 0 ) ( k ) (8) \bar x=\frac{1}{n}\sum_{k=1}^nx^{(0)}(k) \tag{8} xˉ=n1k=1nx(0)(k)(8)

S 1 2 = 1 n ∑ k = 1 n ( x ( 0 ) ( k ) − x ˉ ) 2 (9) S_1^2=\frac{1}{n}\sum_{k=1}^n(x^{(0)}(k)-\bar x)^2 \tag{9} S12=n1k=1n(x(0)(k)xˉ)2(9)

   定义 ϵ ( 0 ) \epsilon ^{(0)} ϵ(0)的均值和方差分别为:

ϵ ˉ = 1 n ∑ k = 1 n ϵ ( 0 ) ( k ) (10) \bar \epsilon=\frac{1}{n}\sum_{k=1}^n\epsilon^{(0)}(k) \tag{10} ϵˉ=n1k=1nϵ(0)(k)(10)

S 2 2 = 1 n ∑ k = 1 n ( ϵ ( 0 ) ( k ) − ϵ ˉ ) 2 (11) S_2^2=\frac{1}{n}\sum_{k=1}^n(\epsilon^{(0)}(k)-\bar \epsilon)^2 \tag{11} S22=n1k=1n(ϵ(0)(k)ϵˉ)2(11)

   定义均方差比值为:

C = S 2 S 1 (12) C=\frac{S_2}{S_1} \tag{12} C=S1S2(12)

   对于给定的 C 0 > 0 C_0 > 0 C0>0 ,当 C < C 0 C < C_0 C<C0 时,称模型为均方差比合格模型。

5 小误差概率合格模型

   定义小概率为:
p = P { ∣ ϵ ( 0 ) ( k ) − ϵ ˉ ∣ < 0.6745 S 1 } (13) p=P\{ \left | \epsilon^{(0)}(k)-\bar \epsilon \right|<0.6745S_1 \} \tag{13} p=P{ϵ(0)(k)ϵˉ<0.6745S1}(13)

   对于给定的 p 0 > 0 p_0 > 0 p0>0 ,当 p < p 0 p < p_0 p<p0 时,称模型为小误差概率合格模型。

   下标列出了灰色模型精度检验等级。

等级相对误差 Δ \Delta Δ绝对关联度 g ( 0 ) g^{(0)} g(0)均方值比值 C C C小误差概率 p p p
0.010.900.350.95
0.050.800.500.80
0.100.700.650.10
0.200.600.800.60

   等级越高(等级一最高),说明模型越好。

声明:本文内容由网友自发贡献,不代表【wpsshop博客】立场,版权归原作者所有,本站不承担相应法律责任。如您发现有侵权的内容,请联系我们。转载请注明出处:https://www.wpsshop.cn/w/weixin_40725706/article/detail/265949?site
推荐阅读
相关标签
  

闽ICP备14008679号