这是 AR(1) 与 AR(2) 案例的部分答案。
在 AR(1) 的情况下,方差是 (设置 )
在 AR(2) 的情况下,可以证明
这个表达式可能有助于激发我为什么不能上来沿着这条路线给出一般答案,因为在一般 AR(p) 情况下,方差的表达式只会变得更加复杂。σ2=1
γ0=11−ϕ2
γ0=(1−ϕ2)1−ϕ2−ϕ21−ϕ21ϕ2−ϕ22(1−ϕ2)
您的限制意味着(仅限于静止情况),,,因此也是。0<ϕ<1ϕ1,ϕ2>0ϕ1+ϕ2=ϕϕ>ϕ1
将其映射到 R 进行数值评估给了我
phiAR1 <- seq(0.01,0.99,by=0.01)
phi_1 <- seq(0.01,0.99,by=0.01)
gamma0AR1 <- function(phiAR1) 1/(1-phiAR1^2)
gamma0AR2 <- function(phiAR1,phi_1){
phi_2 <- phiAR1 - phi_1
return(ifelse(phi_2>0,(1-phi_2)/(1-phi_2-phi_1^2 -phi_1^2*phi_2 - phi_2^2*(1-phi_2)),NA))
}
Vardiff <- function(phiAR1,phi_1) gamma0AR1(phiAR1)-gamma0AR2(phiAR1,phi_1)
Vardiffs <- outer(phiAR1,phi_1,Vardiff)
persp(phiAR1,phi_1,Vardiffs)
min(Vardiffs[!is.na(Vardiffs)])
这似乎是一个非负函数:

我还试图让 MAPLE 通过分析验证这一点
restart;
assume(0<phi[AR1],phi[AR1]<1,phi[1]>0,phi[AR1]>phi[1]);
is(1/(1-phi[AR1]^2) - (1-phi[2])/(1-phi[2]-(phi[1])^2 -(phi[1])^2*phi[2] - (phi[2])^2*(1-phi[2]))>0);
但这似乎不是一个正确的方法。