Я предполагаю, что характерное уравнение, от которого вы уходите, отличается от моего. Позвольте мне сделать пару шагов, чтобы понять, согласны ли мы.
Рассмотрим уравнение
λ2−ϕ1λ−ϕ2=0
Если z является корнем «стандартного» характеристического уравнения 1−ϕ1z−ϕ2z2=0 и установки z−1=λ , отображение получается при переписывании стандартного уравнения следующим образом:
1−ϕ1z−ϕ2z2⇒z−2−ϕ1z−1−ϕ2⇒λ2−ϕ1λ−ϕ2===000
Следовательно, альтернативное условие устойчивостиAR(2)состоит в том, что все корни первого отображения находятсявнутриединичного круга,|z|>1⇔|λ|=|z−1|<1.
Мы используем это представление для получения треугольника стационарности процесса AR(2) , то есть, что AR(2) является стабильным, если выполняются следующие три условия:
- ϕ2<1+ϕ1
- ϕ2<1−ϕ1
- ϕ2>−1
Напомним, что вы можете записать корни первого отображения (если они действительны) как
λ1,2=ϕ1±ϕ21+4ϕ2−−−−−−−√2
чтобы найти первые два условия.
Тогда AR(2) является стационарным тогда и только тогда, когда |λ|<1 , следовательно (если λi действительны):
−1<ϕ1±ϕ21+4ϕ2−−−−−−−√2⇒−2<ϕ1±ϕ21+4ϕ2−−−−−−−√<<12
Больший из двухλiограниченϕ1+ϕ21+4ϕ2−−−−−−−√<2или:
ϕ1+ϕ21+4ϕ2−−−−−−−√⇒ϕ21+4ϕ2−−−−−−−√⇒ϕ21+4ϕ2⇒ϕ21+4ϕ2⇒ϕ2<<<<<22−ϕ1(2−ϕ1)24−4ϕ1+ϕ211−ϕ1
Аналогично находим, чтоϕ2<1+ϕ1.
Если λi сложен, то ϕ21<−4ϕ2 и так λ1,2=ϕ1/2±i−(ϕ21+4ϕ2)−−−−−−−−−−√/2.
Квадрат модуля комплексного числа представляет собой квадрат действительного плюс квадрат мнимой части. Следовательно,
λ2=(ϕ1/2)2+(−(ϕ21+4ϕ2)−−−−−−−−−−√/2)2=ϕ21/4−(ϕ21+4ϕ2)/4=−ϕ2.
Это стабильно, если|λ|<1, следовательно, если−ϕ2<1илиϕ2>−1, как должно было быть показано. (Ограничениеϕ2<1ϕ22<1ϕ2<1+ϕ1ϕ2<1−ϕ1
Построив треугольник стационарности, также указав линию, которая отделяет комплекс от реальных корней, мы получим
Произведено в R с использованием
phi1 <- seq(from = -2.5, to = 2.5, length = 51)
plot(phi1,1+phi1,lty="dashed",type="l",xlab="",ylab="",cex.axis=.8,ylim=c(-1.5,1.5))
abline(a = -1, b = 0, lty="dashed")
abline(a = 1, b = -1, lty="dashed")
title(ylab=expression(phi[2]),xlab=expression(phi[1]),cex.lab=.8)
polygon(x = phi1[6:46], y = 1-abs(phi1[6:46]), col="gray")
lines(phi1,-phi1^2/4)
text(0,-.5,expression(phi[2]<phi[1]^2/4),cex=.7)
text(1.2,.5,expression(phi[2]>1-phi[1]),cex=.7)
text(-1.75,.5,expression(phi[2]>1+phi[1]),cex=.7)