时间序列:白噪声

白噪声

一、白噪声定义及性质

在时间序列中,最简单的平稳过程(纯随机过程)就是白噪声过程(White Noise),具体如下:
{ Z t {Z_{t}} Zt} 是白噪声过程,如果满足:
E ( Z t ) = 0 E(Z_{t})=0 E(Zt)=0 V a r ( Z t ) = σ z 2 Var(Z_{t})=\sigma_{z}^{2} Var(Zt)=σz2 C o v ( Z t , Z s ) = 0 , t ≠ s Cov(Z_{t},Z_{s})=0 ,t\neq s Cov(Zt,Zs)=0,t=s也就是均值为0,方差为 σ z 2 \sigma_{z}^{2} σz2 ,协方差为0 (无自相关性) 的序列,简单记为 Z t ∼ W N ( 0 , σ z 2 ) Z_{t}\sim WN(0,\sigma_{z}^{2}) ZtWN(0,σz2)
从白噪声序列的协方差为0可以得到,其ACF除在0处之外均为0,即 ρ k = C o v ( Z t , Z t + k ) V a r ( Z t ) V a r ( Z t + k ) = 0 , k ≠ 0 \rho_{k}=\frac{Cov(Z_{t},Z_{t+k})}{\sqrt{Var(Z_{t})Var(Z_{t+k})}}=0 ,k\neq 0 ρk=Var(Zt)Var(Zt+k) Cov(Zt,Zt+k)=0,k=0
只有当序列为白噪声序列才有上述的关系,容易出错的是,很多人往往计算时会下意识默认序列为平稳序列,于是直接将分母处的 σ Z t 2 \sigma_{Z_{t}}^{2} σZt2 σ Z t + k 2 \sigma_{Z_{t+k}}^{2} σZt+k2 默认为相等,如果序列是平稳的,那没有问题;否则,请按照原先公式进行计算。

二、白噪声检验

现在我们需要解决的问题是得到的平稳序列是否是一个完全随机的序列。如果序列是完全随机的,那无论怎么利用这个序列也无法进行预测;反之,我们可以找出内在联系,为预测的准确性打下基础。

其中最简单的方法就是观察序列的ACF图,完全随机的序列应该是没有明显的自相关性。

而较为严格的一种方法就是:混成检验:
H 0 : ρ 1 = ρ 2 ⋯ = ρ k = 0 V S H 1 : ∃ ρ i ≠ 0 , 1 ≤ i ≤ k H_{0}: \rho_{1}=\rho_{2} \cdots =\rho_{k}=0 \quad VS \quad H_{1}: \exists \rho_{i}\neq 0, 1\leq i \leq k H0:ρ1=ρ2=ρk=0VSH1:ρi=0,1ik

对于一列完全随机的序列,有近似分布 ρ i ^ ∼ N ( 0 , 1 T ) \hat{\rho_{i}} \sim N(0,\frac{1}{T}) ρi^N(0,T1) ,对于任意的 k k k,下面定义检验统计量为: Q ( k ) = T ∑ i = 1 k ρ k 2 ^ Q(k)= T \sum_{i=1}^{k}\hat{\rho^{2}_{k}} Q(k)=Ti=1kρk2^
在完全随机的假设下有:
Q ( k ) ≈ χ k 2 Q(k)\approx \chi_{k}^{2} Q(k)χk2
Q ( k ) ≈ χ k 2 > χ k 2 ( 1 − α ) Q(k)\approx \chi_{k}^{2} > \chi_{k}^{2}(1-\alpha) Q(k)χk2>χk2(1α) 时,拒绝原假设,即不能认为序列是纯随机的。而这种检验方式被称为 B o x − P i e r c e ′ s Box-Pierce's BoxPierces 混成检验。

下面换了一个检验统计量—— L j u n g − B o x ′ s Ljung-Box's LjungBoxs 检验统计量:
L B ( k ) = T ( T + 2 ) ∑ i = 1 k ρ i 2 ^ T − k LB(k)=T(T+2)\sum_{i=1}^{k}\frac{\hat{\rho_{i}^{2}}}{T-k} LB(k)=T(T+2)i=1kTkρi2^
对于修正的 L j u n g − B o x ′ s Ljung-Box's LjungBoxs 检验统计量的分布在样本量较小时比上面的 B o x − P i e r c e ′ s Box-Pierce's BoxPierces 混成检验统计量更接近 χ k 2 \chi_{k}^{2} χk2 分布,而在样本量较大时两者差别却并不大。所以在实际中 L j u n g − B o x ′ s Ljung-Box's LjungBoxs 检验统计量更常用。

在R语言中使用 B o x . t e s t Box.test Box.test 函数可以实现上述两个检验。

  • 2
    点赞
  • 22
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 1
    评论
在Matlab中进行时间序列预测时,可以使用噪声检验来检验时间序列数据是否具有噪声的特性。噪声是指一个时间序列的观测值之间没有任何相关性。在Matlab中,可以使用函数`archtest`进行噪声检验。这个函数可以用来对给定的时间序列进行ARCH效应检验,ARCH效应是指时间序列的波动是否存在异方差性。如果时间序列具有噪声特性,那么满足条件ARCH效应为0。因此,可以使用`archtest`函数检验时间序列数据的ARCH效应,并根据检验结果判断时间序列数据是否具有噪声特性。 以下是一个使用`archtest`函数进行噪声检验的示例代码: ```matlab % 假设你已经有了一个名为data的时间序列数据 % 这里的data是一个列向量,其中存储了观测值 % 进行噪声检验 [pValue, h, stat = archtest(data); % 根据检验结果判断时间序列数据是否具有噪声特性 if h == 0 disp('时间序列数据通过噪声检验,具有噪声特性。'); else disp('时间序列数据未通过噪声检验,不具有噪声特性。'); end ``` 在这个示例代码中,`archtest`函数返回的`pValue`是检验结果的p值,`h`是假设检验的结果,如果h为0表示通过检验,如果h为1表示未通过检验,`stat`是检验统计量。根据判断结果,可以确定时间序列数据是否具有噪声特性。<span class="em">1</span><span class="em">2</span><span class="em">3</span> #### 引用[.reference_title] - *1* *2* *3* [【LSTM时序预测】基于matlab LSTM时间序列神经网络预测【含Matlab源码 2267期】](https://blog.csdn.net/TIQCmatlab/article/details/128243832)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_2"}}] [.reference_item style="max-width: 100%"] [ .reference_list ]

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

梦什

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值