白噪音及其检验

news/2024/12/29 18:30:52/

白噪音

  • 定义1:如果时间序列 { ε t , t = 1 , ⋯ T } \{\varepsilon_t, t = 1, \cdots T\} {εt,t=1,T}满足如下条件,则称该时间序列为白噪音序列
    • E ( ε t ) = 0 , V a r ( ε t ) = σ 2 E(\varepsilon_t) = 0, Var(\varepsilon_t) = \sigma^2 E(εt)=0,Var(εt)=σ2
    • E ( ε t × ε s ) = 0 , ∀ t ≠ s E(\varepsilon_t \times \varepsilon_s) = 0, \forall t \neq s E(εt×εs)=0,t=s。即当 s ≠ t s \neq t s=t时, ε t \varepsilon_t εt ε s \varepsilon_s εs不相关
  • 定义2:如果时间序列 { X t } \{X_t\} {Xt}满足如下性质,则称该序列为纯随机序列,也称为白噪音序列
    • E ( X t ) = μ , ∀ t ∈ T E(X_t) = \mu, \forall t \in T E(Xt)=μ,tT
    • r ( t , s ) = σ 2 , t = s , ∀ t , s ∈ T r(t, s) = \sigma^2, t = s, \forall t, s \in T r(t,s)=σ2,t=s,t,sT
    • r ( t , s ) = 0 , t ≠ s , ∀ t , s ∈ T r(t, s) = 0, t \neq s, \forall t, s \in T r(t,s)=0,t=s,t,sT
  • 性质
    • 纯随机性,无记忆: r ( k ) = 0 , ∀ k ≠ 0 r(k) = 0, \forall k \neq 0 r(k)=0,k=0
    • 方差齐性:序列中每个变量的方差都相等,即 D ( X t ) = r ( 0 ) = σ 2 D(X_t) = r(0) = \sigma^2 D(Xt)=r(0)=σ2。如果序列不满足方差齐性,则称序列具有异方差性质。
      • 根据马尔科夫定理,满足方差齐性时,用最小二乘得到的未知参数估计值是准确的、有效的。若不满足,最小二乘估计值不是方差最小线性无偏估计,拟合模型的精度会受到很大的影响
      • 模型拟合时,检验内容之一就是要检验拟合模型的残差是否满足方差齐性。若不满足,说明残差序列还不是白噪音序列,即拟合模型没有充分提取随机序列中的相关信息。此时通常需要使用适当的条件异方差模型来处理异方差信息。
      • 方差齐性检验

白噪音的检验

  • Barlett定理:如果一个时间序列是纯随机的,得到一个观察期数为 n n n的观察序列,那么该序列的延迟非零期的样本自相关系数将近似服从均为零,方差为序列观察期数倒数的正态分布: ρ k ∼ N ( 0 , 1 n ) , ∀ k ≠ 0 \rho_k \sim N(0, \frac{1}{n}), \forall k \neq 0 ρkN(0,n1)k=0
  • Q Q Q统计量/ Q B P Q_{BP} QBP统计量
    • 适用范围:大样本
    • 近似服从自由度为 m m m的卡方分布: Q = n × ∑ k = 1 m ρ k 2 ∼ χ 2 ( m ) Q = n \times \sum_{k=1}^m{\rho_k^2} \sim \chi^2(m) Q=n×k=1mρk2χ2(m)
      • 当统计量大于 χ 1 − α 2 ( m ) \chi_{1-\alpha}^2(m) χ1α2(m)分位点,或该统计量的 P P P值小于 α \alpha α时,则可以以 1 − α 1 - \alpha 1α的置信水平拒绝原假设,认为该序列为非白噪音序列。否则,为白噪音
  • L B LB LB统计量/ Q L B Q_{LB} QLB统计量
    • 适用范围:小样本
    • 近似服从自由度为 m m m的卡方分布: L B = n × ( n + 2 ) × ∑ k = 1 m ρ k 2 n − k LB = n \times (n + 2) \times \sum_{k=1}^m{\frac{\rho_k^2}{n - k}} LB=n×(n+2)×k=1mnkρk2
  • 对于一个长度为 T T T的白噪音序列,我们期望在95%的置信度下,它的自相关值(acf)属于区间 [ − 2 T , 2 T ] [\frac{-2}{\sqrt{T}}, \frac{2}{\sqrt{T}}] [T 2,T 2]

http://www.ppmy.cn/news/694631.html

相关文章

基于matlab根据特征进行全景图像拼接(附源码)

一、前言 此示例演示如何使用基于特征的图像配准技术自动创建全景图。 特征检测和匹配是许多计算机视觉应用(如图像配准、跟踪和对象检测)中使用的强大技术。在此示例中,基于特征的技术用于自动拼接一组图像。图像拼接过程是基于特征的图像…

1位十进制加法器设计

1位8421BCD码加法原理图如下: 修正逻辑:Adjust carry_out Z3∙Z2 Z3∙Z1

为什么1//0.1等于9.0,而1//-0.1=-10?

为什么1//0.1等于9.0,而1//-0.1-10? 在这之前我先纠正一下大家的一个普遍误区,讲解这个误区之前先来介绍一下Python中用于十进制数计算的模块decimal,首先引入这个模块(from decimal import Decimal)用它来测试一下0.1这个数究竟…

1到15的二进制、十进制、十六进制

十进制十六进制二进制11000122001033001144010055010166011077011188100099100110A101011B101112C110013D110114E111015F1111

夜深人静写算法(十四)- 0/1 背包

文章目录 一、前言二、0/1 背包问题1、状态设计2、状态转移方程3、初始状态4、非法状态5、状态初始化三、0/1 背包问题的实现1、背包物品结构设计2、状态数组3、状态转移4、状态初始化四、0/1 背包问题的扩展思考1、最大值问题2、最小值问题3、存在性问题4、方案数问题5、有顺序…

1093: 数1的个数

存限制 : 128 MB 题目描述 给定一个十进制正整数n(1≤n≤10000),写下从1到n的所有整数,然后数一下其中出现的数字“1”的个数。 例如当n2时,写下1,2。这样只出现了1个“1”;当n12时,写下1,2,…

如何用C语言实现十进制到二进制的转换并计算二进制中1的个数

<pre class"csharp" name"code">#define _CRT_SECURE_NO_WARNINGS 1 #include<stdio.h> #include<stdlib.h> int main() {int num, temp, i 0, k 0, arr[100];printf("请输入一个十进制的数:\n");scanf("%d", &am…

1-10的阶乘之和

对于十以内的阶乘之和代码如下&#xff1a; #include<stdio.h> int main() {int i 0;int sum_mul 1;int sum_add 0;for (i 1; i < 10; i){sum_mul sum_mul * i;sum_add sum_add sum_mul;}printf("%d\n", sum_add);return 0; }相较于1-10的相加&…