赞
踩
前两篇博客说了灰色预测的基本步骤,其中讲了两种检验模型的方法,本次再具体地说一下灰色模型的检验方法。
设有 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 k≤n,定义 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=1∑nΔ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.5a1−0.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=1∑nx(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=1∑n(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=1∑nϵ(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=1∑n(ϵ(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.01 | 0.90 | 0.35 | 0.95 |
二 | 0.05 | 0.80 | 0.50 | 0.80 |
三 | 0.10 | 0.70 | 0.65 | 0.10 |
四 | 0.20 | 0.60 | 0.80 | 0.60 |
等级越高(等级一最高),说明模型越好。
Copyright © 2003-2013 www.wpsshop.cn 版权所有,并保留所有权利。