一、AP模型
AR(p)模型的均值、方差、自相关函数(ACF)和偏自相关函数(PACF)推导如下:
1. 均值(Mean)
AR(p)模型形式为:
Xt=c+ϕ1Xt−1+ϕ2Xt−2+⋯+ϕpXt−p+εt
其中 εt 是白噪声,满足 E[εt]=0,Var(εt)=σ2。
推导:
假设平稳性,均值为常数 E[Xt]=μ,则:
μ=c+ϕ1μ+ϕ2μ+⋯+ϕpμ
解得:
μ=c1−ϕ1−ϕ2−⋯−ϕp
均值:
E[Xt]=c1−ϕ1−ϕ2−⋯−ϕp
2. 方差(Variance)
推导:
对于平稳AR(p)模型,方差 γ0=Var(Xt) 满足:
γ0=ϕ1γ1+ϕ2γ2+⋯+ϕpγp+σ2
其中 γk 是滞后 k 的自协方差。通过Yule-Walker方程求解 γ0。
方差:
Var(Xt)=γ0
需通过Yule-Walker方程求解。
3. 自相关函数(ACF)
自协方差函数 γk 满足:
γk=ϕ1γk−1+ϕ2γk−2+⋯+ϕpγk−p,k>0
自相关系数 ρk=γkγ0 满足:
ρk=ϕ1ρk−1+ϕ2ρk−2+⋯+ϕpρk−p,k>0
ACF特性: ACF拖尾
指数或震荡衰减,无限延续。
4. 偏自相关函数(PACF)
PACF ϕkk 表示剔除中间 k−1 个变量影响后,Xt 与 Xt−k 的相关系数。
对于AR(p)模型,PACF在滞后 p 后截尾。
推导:
PACF ϕkk 可通过Yule-Walker方程或递归公式计算。对于AR(p)模型:
ϕkk=0,k>p
PACF特性:
- PACF截尾,当 k>p 时 ϕkk=0。
总结
统计量 |
AR(p) 模型特性 |
均值 |
E[Xt]=c1−ϕ1−ϕ2−⋯−ϕp |
方差 |
Var(Xt)=γ0,需通过Yule-Walker方程求解 |
ACF |
拖尾:指数或震荡衰减(无限延续) |
PACF |
截尾:当 k>p 时 ϕkk=0 |
示例(AR(1)): |
|
- 均值:μ=c1−ϕ1
- 方差:γ0=σ21−ϕ21
- ACF:ρk=ϕk1(指数衰减)
- PACF:ϕ11=ϕ1,ϕkk=0(k>1)
二、MA模型
MA(q)模型的均值、方差、自相关函数(ACF)和偏自相关函数(PACF)推导如下:
1. 均值(Mean)
MA(q)模型形式为:
Xt=μ+εt+θ1εt−1+θ2εt−2+⋯+θqεt−q
其中 εt 是白噪声,满足 E[εt]=0,Var(εt)=σ2。
推导:
E[Xt]=μ+E[εt]+θ1E[εt−1]+⋯+θqE[εt−q]=μ
均值:
E[Xt]=μ
2. 方差(Variance)
推导:
Var(Xt)=Var(εt+θ1εt−1+⋯+θqεt−q)
由于 εt 独立,方差为各分量方差之和:
Var(Xt)=σ2(1+θ21+θ22+⋯+θ2q)
方差:
Var(Xt)=σ2(1+q∑i=1θ2i)
3. 自相关函数(ACF)
自协方差函数为:
γk=Cov(Xt,Xt−k)=E[(Xt−μ)(Xt−k−μ)]
对于 k>q,Xt 和 Xt−k 无重叠的 ε 项,故 γk=0。
对于 1≤k≤q,协方差由共同的白噪声项贡献:
γk=σ2q−k∑j=0θjθj+k(其中 θ0=1)
自相关系数:
ρk=γkγ0=∑q−kj=0θjθj+k1+∑qi=1θ2i
ACF特性:
- 当 k>q 时,ρk=0(截尾现象)。
- 当 k≤q 时,ρk 由上述公式给出。
4. 偏自相关函数(PACF)
PACF ϕkk 表示剔除中间 k−1 个变量影响后,Xt 与 Xt−k 的相关系数。
对于MA(q)模型,PACF不截尾,而是拖尾(指数或震荡衰减),可通过以下方式理解:
- MA(q)可逆时,可表示为AR(∞),因此PACF无限延伸。
- 例(MA(1)): PACF满足
ϕkk=(−1)kθk(1−θ2)1−θ2(k+1)
(绝对值指数衰减)。
PACF特性:
- PACF拖尾,无闭合表达式,需通过Yule-Walker方程递归计算。
5. 总结
统计量 |
MA(q) 模型特性 |
均值 |
E[Xt]=μ |
方差 |
Var(Xt)=σ2(1+∑qi=1θ2i) |
ACF |
截尾:ρk=0(当 k>q) |
PACF |
拖尾:指数或震荡衰减(无限延续) |
示例(MA(1)): |
|
- ACF: ρ1=θ1+θ2,ρk=0(k>1)
- PACF: ϕkk∝(−θ)k(逐步衰减)
三、ARMA模型
ARMA(p, q)模型的均值、方差、自相关函数(ACF)和偏自相关函数(PACF)推导如下:
1. 均值(Mean)
ARMA(p, q)模型形式为:
Xt=c+ϕ1Xt−1+ϕ2Xt−2+⋯+ϕpXt−p+εt+θ1εt−1+θ2εt−2+⋯+θqεt−q
其中 εt 是白噪声,满足 E[εt]=0,Var(εt)=σ2。
推导:
假设平稳性,均值为常数 E[Xt]=μ,则:
μ=c+ϕ1μ+ϕ2μ+⋯+ϕpμ
解得:
μ=c1−ϕ1−ϕ2−⋯−ϕp
均值: E[Xt]=c1−ϕ1−ϕ2−⋯−ϕp
2. 方差(Variance)
推导:
对于平稳ARMA(p, q)模型,方差 γ0=Var(Xt) 满足:
γ0=p∑i=1ϕiγi+σ2(1+q∑j=1θ2j)
其中 γk 是滞后 k 的自协方差。通过Yule-Walker方程或递归公式求解 γ0。
方差: Var(Xt)=γ0,需通过递归或Yule-Walker方程求解。
3. 自相关函数(ACF)
自协方差函数 γk 满足:
γk=ϕ1γk−1+ϕ2γk−2+⋯+ϕpγk−p,k>q
自相关系数 ρk=γkγ0 满足:
ρk=ϕ1ρk−1+ϕ2ρk−2+⋯+ϕpρk−p,k>q
ACF特性:
4. 偏自相关函数(PACF)
PACF ϕkk 表示剔除中间 k−1 个变量影响后,Xt 与 Xt−k 的相关系数。
对于ARMA(p, q)模型,PACF在滞后 p 后拖尾。
推导:
PACF ϕkk 通过Yule-Walker方程或递归公式计算。对于ARMA(p, q)模型:
ϕkk 拖尾,无限延续
PACF特性:
- PACF拖尾,无闭合表达式,需通过递归计算。
总结
统计量 |
ARMA(p, q) 模型特性 |
均值 |
E[Xt]=c1−ϕ1−ϕ2−⋯−ϕp |
方差 |
Var(Xt)=γ0,需通过递归或Yule-Walker方程求解 |
ACF |
拖尾:指数或震荡衰减(无限延续) |
PACF |
拖尾:无限延续,无闭合表达式 |
补充说明
- ARMA(1,1) 示例:
- 均值:μ=c1−ϕ1
- 方差:γ0=1+θ21+2ϕ1θ11−ϕ21σ2
- ACF:ρ1=(1+ϕ1θ1)(ϕ1+θ1)1+θ21+2ϕ1θ1,ρk=ϕ1ρk−1(k>1)
- PACF:拖尾,无限延续。
- ARMA模型特性:
- ACF 和 PACF 均拖尾,无法通过简单截尾判断阶数,需结合其他方法(如AIC、BIC)确定 p 和 q。
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· TypeScript + Deepseek 打造卜卦网站:技术与玄学的结合
· 阿里巴巴 QwQ-32B真的超越了 DeepSeek R-1吗?
· 如何调用 DeepSeek 的自然语言处理 API 接口并集成到在线客服系统
· 【译】Visual Studio 中新的强大生产力特性
· 2025年我用 Compose 写了一个 Todo App