Ich vermute, dass die charakteristische Gleichung, von der Sie abweichen, sich von meiner unterscheidet. Lassen Sie mich in ein paar Schritten überprüfen, ob wir uns einig sind.
Betrachten Sie die Gleichung
λ2−ϕ1λ−ϕ2=0
Wenn eine Wurzel der "Standard" -Kennliniengleichung und , erhält die Anzeige durch Umschreiben der Standardgleichung folgende Ergebnisse:
Daher ist eine alternative Bedingung für die Stabilität eines dass sich alle Wurzeln der ersten Anzeige innerhalb des Einheitskreises befinden, .z1−ϕ1z−ϕ2z2=0z−1=λ1−ϕ1z−ϕ2z2⇒z−2−ϕ1z−1−ϕ2⇒λ2−ϕ1λ−ϕ2===000
AR(2)|z|>1⇔|λ|=|z−1|<1
Wir verwenden diese Darstellung, um das Stationaritätsdreieck eines -Prozesses abzuleiten, dh , dass ein stabil ist, wenn die folgenden drei Bedingungen erfüllt sind: AR(2)AR(2)
- ϕ2<1+ϕ1
- ϕ2<1−ϕ1
- ϕ2>−1
Denken Sie daran, dass Sie die Wurzeln der ersten Anzeige (falls vorhanden) als schreiben können, um die
zu finden erste zwei Bedingungen.λ1,2=ϕ1±ϕ21+4ϕ2−−−−−−−√2
Dann ist der stationär, wenn , also (wenn real ist):
Der größere der beiden ist begrenzt durch , oder:
Analog .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
Wenn komplex ist , dann und soDer Quadratmodul einer komplexen Zahl ist das Quadrat des Realteils plus das Quadrat des Imaginärteils. Daher ist
Dies ist stabil, wenn , also wenn oder , wie gezeigt werden sollte. (Die sich aus ergebende Einschränkung ist im Hinblick auf und redundant .)λ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
Das Zeichnen des Stationaritätsdreiecks, das auch die Linie angibt, die komplexe von realen Wurzeln trennt, erhalten wir
Produziert in R mit
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)